codac2_AnalyticExpr.h

Go to the documentation of this file.

 
#pragma once
 
#include <map>
#include <memory>
#include <utility>
#include "codac2_ExprBase.h"
#include "codac2_Domain.h"
#include "codac2_FunctionArgsList.h"
#include "codac2_AnalyticType.h"
 
namespace codac2
{
  using ValuesMap = std::map<ExprID,std::shared_ptr<AnalyticTypeBase>>;
 
  template<typename T>
  class AnalyticExpr : public ExprBase
  {
    public:
 
      AnalyticExpr<T>& operator=(const AnalyticExpr<T>& x) = delete;
 
      virtual T fwd_eval(ValuesMap& v, Index total_input_size, bool natural_eval) const = 0;
      virtual void bwd_eval(ValuesMap& v) const = 0;
 
      T init_value(ValuesMap& v, const T& x) const
      {
        auto it = v.find(unique_id());
 
        if(it == v.end())
        {
          v[unique_id()] = std::make_shared<T>(x);
          return x;
        }
 
        *std::dynamic_pointer_cast<T>(it->second) = x;
        return *std::dynamic_pointer_cast<T>(it->second);
      }
 
      T& value(ValuesMap& v) const
      {
        assert(v.find(unique_id()) != v.end() && "argument cannot be found");
        return *std::dynamic_pointer_cast<T>(v[unique_id()]);
      }
 
      virtual bool belongs_to_args_list(const FunctionArgsList& args) const = 0;
  };
 
  template<typename C,typename Y,typename... X>
  class AnalyticOperationExpr : public AnalyticExpr<Y>, public OperationExprBase<AnalyticExpr<X>...>
  {
    public:
 
      AnalyticOperationExpr(std::shared_ptr<AnalyticExpr<X>>... x)
        : OperationExprBase<AnalyticExpr<X>...>(x...)
      { }
 
      AnalyticOperationExpr(const AnalyticOperationExpr<C,Y,X...>& e)
        : OperationExprBase<AnalyticExpr<X>...>(e)
      { }
 
      std::shared_ptr<ExprBase> copy() const
      {
        return std::make_shared<AnalyticOperationExpr<C,Y,X...>>(*this);
      }
 
      void replace_expr(const ExprID& old_expr_id, const std::shared_ptr<ExprBase>& new_expr)
      {
        return OperationExprBase<AnalyticExpr<X>...>::replace_expr(old_expr_id, new_expr);
      }
 
      Y fwd_eval(ValuesMap& v, Index total_input_size, bool natural_eval) const
      {
        return std::apply(
          [this,&v,total_input_size,natural_eval](auto &&... x)
          {
            if(natural_eval)
              return AnalyticExpr<Y>::init_value(v, C::fwd_natural(x->fwd_eval(v, total_input_size, natural_eval)...));
 
            else
              return AnalyticExpr<Y>::init_value(v, C::fwd_centered(x->fwd_eval(v, total_input_size, natural_eval)...));
          },
        this->_x);
      }
 
      void bwd_eval(ValuesMap& v) const
      {
        auto y = AnalyticExpr<Y>::value(v);
 
        std::apply([&v,y](auto &&... x)
        {
          C::bwd(y.a, x->value(v).a...);
        }, this->_x);
 
        std::apply([&v](auto &&... x)
        {
          (x->bwd_eval(v), ...);
        }, this->_x);
      }
 
      virtual bool belongs_to_args_list(const FunctionArgsList& args) const
      {
        bool b = true;
 
        std::apply([&b,args](auto &&... x)
        {
          ((b &= x->belongs_to_args_list(args)), ...);
        }, this->_x);
 
        return b;
      }
  };
}