zeromodel.hh

#ifndef __DUNE_ACFEM_MODELS_MODULES_ZEROMODEL_HH__
#define __DUNE_ACFEM_MODELS_MODULES_ZEROMODEL_HH__
 
#include "../../expressions/terminal.hh"
 
#include "../modelbase.hh"
#include "../modeltraits.hh"
 
namespace Dune {
 
  namespace ACFem::PDEModel {
 
    template<class DomainFunctionSpace, class RangeFunctionSpace = DomainFunctionSpace>
    class ZeroModel
      : public ModelBase<DomainFunctionSpace, RangeFunctionSpace>
      , public Expressions::SelfExpression<ZeroModel<DomainFunctionSpace, RangeFunctionSpace> >
      , public ZeroExpression
    {
     public:
      ZeroModel(const std::string& name = "")
        : name_(name == "" ? "0" : name)
      {}
 
      std::string name() const
      {
        return name_;
      }
 
     protected:
      const std::string name_;
    };
 
    template<class T, class F = Expressions::Closure,
             std::enable_if_t<IsPDEModel<T>::value, int> = 0>
    auto zeroModel(const T& t, const std::string& name, F closure = F{})
    {
      return closure(ZeroModel<typename T::DomainFunctionSpaceType, typename T::RangeFunctionSpaceType>(name));
    }
 
    template<class T, class F = Expressions::Closure, std::enable_if_t<IsPDEModel<T>::value, int> = 0>
    auto zeroModel(const T&, F closure = F{})
    {
      return closure(ZeroModel<typename T::DomainFunctionSpaceType, typename T::RangeFunctionSpaceType>());
    }
 
    template<class T, class F = Expressions::Closure, std::enable_if_t<IsPDEModel<T>::value, int> = 0>
    auto zeroModel(F closure = F{})
    {
      return closure(ZeroModel<typename T::DomainFunctionSpaceType, typename T::RangeFunctionSpaceType>());
    }
 
 
 
 
  } // namespace ACFem::PDEModel
 
  namespace ACFem {
 
    using PDEModel::zeroModel;
 
  }
 
}  //Namespace Dune
 
 
#endif // __DUNE_ACFEM_MODELS_ZEROMODEL_HH__