Constant.h

Go to the documentation of this file.

/* This file is part of MoDELib, the Mechanics Of Defects Evolution Library.
 *
 *
 * MoDELib is distributed without any warranty under the
 * GNU General Public License (GPL) v2 <http://www.gnu.org/licenses/>.
 */
 
#ifndef model_Constatnt_H_
#define model_Constatnt_H_
 
#include <memory>
#include <utility> // for std::move
 
#include <Eigen/Dense>
//#include <EvalExpression.h>
#include <EvalFunction.h>
 
 
namespace model
{
    template <typename T, int _rows, int _cols>
        struct Constant : public EvalFunction<Constant<T,_rows,_cols> >
    {
        constexpr static int rows=_rows;
        constexpr static int cols=_cols;
   
        const T c;
        Constant(const T& c_in) : c(c_in)
        {
        }
        
        template<typename ElementType, typename BaryType>
        const T& operator() (const ElementType&, const BaryType&) const
        {
            return c;
        }
    };
    
    // Operators
    template <int rows, int cols>
    auto make_constant(const Eigen::Matrix<double,rows,cols>& c)
    {
        return std::make_shared<Constant<Eigen::Matrix<double,rows,cols>,rows,cols>>(c);
    }
 
    static inline auto make_constant(const double& c)
    {
        return std::make_shared<Constant<double,1,1>>(c);
    }
     
}       // close namespace
#endif