Detailed Description

template<int DisplacementDim>
struct MaterialLib::Solids::MFront::detail::SetGradient< DisplacementDim >

A helper struct filling MFront's gradient data (or thermodynamic forces) with values taken from a MaterialPropertyLib::VariableArray.

Definition at line 181 of file MFrontGeneric.h.

#include <MFrontGeneric.h>

Collaboration diagram for MaterialLib::Solids::MFront::detail::SetGradient< DisplacementDim >:

[legend]

Public Member Functions
template<typename Grad>
void	operator() (Grad)

Public Attributes
MaterialPropertyLib::VariableArray const &	variable_array
std::optional< MathLib::KelvinVector::KelvinMatrixType< DisplacementDim > > const &	Q
double *	target

Member Function Documentation

◆ operator()()

template<int DisplacementDim>

template<typename Grad>

void MaterialLib::Solids::MFront::detail::SetGradient< DisplacementDim >::operator() ( Grad )

inline

Definition at line 189 of file MFrontGeneric.h.

    {
        auto constexpr num_comp = Grad::template size<DisplacementDim>();
 
        if constexpr (Grad::type == mgis::behaviour::Variable::Type::SCALAR)
        {
            *target = variable_array.*Grad::mpl_var;
        }
        else if constexpr (Grad::type ==
                           mgis::behaviour::Variable::Type::STENSOR)
        {
            using MPLType = MapToMPLType_t<DisplacementDim, Grad::type>;
            auto const& grad_ogs =
                std::get<MPLType>(variable_array.*Grad::mpl_var);
 
            auto const grad_mfront =
                Q ? eigenSwap45View(Q->transpose() * grad_ogs).eval()
                  : eigenSwap45View(grad_ogs).eval();
            std::copy_n(grad_mfront.data(), num_comp, target);
        }
        else if constexpr (Grad::type ==
                           mgis::behaviour::Variable::Type::TENSOR)
        {
            using MPLType = MapToMPLType_t<DisplacementDim, Grad::type>;
            auto const& grad_ogs =
                std::get<MPLType>(variable_array.*Grad::mpl_var);
 
            if (Q.has_value())
            {
                OGS_FATAL("Rotations of tensors are not implemented.");
            }
 
            Eigen::Map<Eigen::Vector<double, num_comp>>{target} =
                ogsTensorToMFrontTensor<DisplacementDim>(grad_ogs);
        }
        else
        {
            OGS_FATAL("Unsupported gradient type {}.",
                      varTypeToString(Grad::type));
        }
 
        target += num_comp;
    }