FEMfaceEvaluation.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_FEMfaceEvaluation_H_
#define model_FEMfaceEvaluation_H_
 
#include <Eigen/Dense>
#include <FEMbaseEvaluation.h>
namespace model
{
    
    /**************************************************************************/
    /**************************************************************************/
    template <typename ElementType,int rows,int cols>
    struct FEMfaceEvaluation : public FEMbaseEvaluation<ElementType,rows,cols>
    {
        constexpr static int dim=ElementType::dim;
        const ElementType& ele;
        const Eigen::Matrix<double,dim+1,1> domainBary;
        const int boundaryFace;
        const double& weight;
//        const Eigen::Matrix<double,dim,1> P;
        
        /**********************************************************************/
        FEMfaceEvaluation(const ElementType& _ele,
                                 const Eigen::Matrix<double,dim+1,1>& _domainBary,
                                 const int& _boundaryFace,
                                 const double& _weight) :
        /* init */ FEMbaseEvaluation<ElementType,rows,cols>(_ele.position(_domainBary))
//        /* init */ Eigen::Matrix<double,rows,cols>(Eigen::Matrix<double,rows,cols>::Zero()),
        /* init */,ele(_ele)
        /* init */,domainBary(_domainBary)
        /* init */,boundaryFace(_boundaryFace)
        /* init */,weight(_weight)
//        /* init */ P(ele.position(domainBary))
        {
 
        }
        
        /**********************************************************************/
        const FEMfaceEvaluation& operator=(const Eigen::Matrix<double,rows,cols>& val)
        {
            static_cast<Eigen::Matrix<double,rows,cols>*>(this)->operator=(val);
            return *this;
        }
        
    };
    
}       // close namespace
#endif