gsMaterialMatrixEval.h

Go to the documentation of this file.

 
#pragma once
 
#include <gsCore/gsFunction.h>
#include <gsKLShell/src/gsMaterialMatrixContainer.h>
 
namespace gismo
{
 
template <class T, enum MaterialOutput out> class gsMaterialMatrixEvalSingle;
 
template <class T, enum MaterialOutput out>
class gsMaterialMatrixEval : public gsFunction<T>
{
public:
    gsMaterialMatrixEval(   const gsMaterialMatrixContainer<T> & materialMatrices,
                            const gsFunctionSet<T> * deformed,
                            const gsMatrix<T> z);
 
    gsMaterialMatrixEval(   gsMaterialMatrixBase<T> * materialMatrix,
                            const gsFunctionSet<T> * deformed,
                            const gsMatrix<T> z);
 
    gsMaterialMatrixEval(   typename gsMaterialMatrixBase<T>::uPtr & materialMatrix,
                            const gsFunctionSet<T> * deformed,
                            const gsMatrix<T> z);
 
    gsMaterialMatrixEval(  const gsMaterialMatrixContainer<T> & materialMatrices,
                            const gsFunctionSet<T> * undeformed,
                            const gsFunctionSet<T> * deformed,
                            const gsMatrix<T> z);
 
    gsMaterialMatrixEval(  gsMaterialMatrixBase<T> * materialMatrix,
                            const gsFunctionSet<T> * undeformed,
                            const gsFunctionSet<T> * deformed,
                            const gsMatrix<T> z);
 
    gsMaterialMatrixEval(   typename gsMaterialMatrixBase<T>::uPtr & materialMatrix,
                            const gsFunctionSet<T> * undeformed,
                            const gsFunctionSet<T> * deformed,
                            const gsMatrix<T> z);
 
    ~gsMaterialMatrixEval();
 
    short_t domainDim() const {return 2;}
 
    short_t targetDim() const { return this->piece(0).targetDim(); }
 
    const gsFunction<T> & piece(const index_t p) const { return *m_pieces[p]; }
 
    void eval_into(const gsMatrix<T>& /*u*/, gsMatrix<T>& /*result*/) const
    { GISMO_NO_IMPLEMENTATION; }
 
protected:
    void _makePieces();
 
    void _makePieces(const gsFunctionSet<T> * undeformed);
 
protected:
    gsMaterialMatrixContainer<T> m_materialMatrices;
    const gsFunctionSet<T> * m_deformed;
    gsMatrix<T> m_z;
    mutable std::vector<gsMaterialMatrixEvalSingle<T,out> *> m_pieces;
};
 
template <class T, enum MaterialOutput out>
class gsMaterialMatrixEvalSingle : public gsFunction<T>
{
public:
 
    gsMaterialMatrixEvalSingle( index_t patch,
                                gsMaterialMatrixBase<T> * materialMatrix,
                                const gsFunctionSet<T> * deformed,
                                const gsMatrix<T> z);
 
    gsMaterialMatrixEvalSingle( index_t patch,
                                gsMaterialMatrixBase<T> * materialMatrix,
                                const gsFunctionSet<T> * undeformed,
                                const gsFunctionSet<T> * deformed,
                                const gsMatrix<T> z);
 
    gsMaterialMatrixEvalSingle(){};
 
    short_t domainDim() const {return 2;}
 
    short_t targetDim() const { return targetDim_impl<out>(); }
 
 
private:
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Density   , short_t>::type targetDim_impl() const { return 1; };
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::VectorN ||
                            _out==MaterialOutput::CauchyVectorN ||
                            _out==MaterialOutput::VectorM ||
                            _out==MaterialOutput::CauchyVectorM ||
                            _out==MaterialOutput::Generic  , short_t>::type targetDim_impl() const { return 3; };
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::MatrixA ||
                            _out==MaterialOutput::MatrixB ||
                            _out==MaterialOutput::MatrixC ||
                            _out==MaterialOutput::MatrixD   , short_t>::type targetDim_impl() const { return 9; };
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Stretch   , short_t>::type targetDim_impl() const { return 3; };
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::PStress  ||
                            _out==MaterialOutput::PStressN ||
                            _out==MaterialOutput::PCauchyStressN ||
                            _out==MaterialOutput::PStressM ||
                            _out==MaterialOutput::PCauchyStressM ||
                            _out==MaterialOutput::PStrainN ||
                            _out==MaterialOutput::PStrainM   , short_t>::type targetDim_impl() const { return 2; };
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::StretchDir, short_t>::type targetDim_impl() const { return 9; };
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::PStressDir, short_t>::type targetDim_impl() const { return 9; };
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Spec2CovTransform ||
                            _out==MaterialOutput::Spec2ConTransform ||
                            _out==MaterialOutput::Cov2CartTransform ||
                            _out==MaterialOutput::Con2CartTransform, short_t>::type targetDim_impl() const { return 9; };
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::TensionField, short_t>::type targetDim_impl() const { return 1; };
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::StretchTransform ||
                            _out==MaterialOutput::PStressTransform, short_t>::type targetDim_impl() const { return 9; };
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Theta ||
                            _out==MaterialOutput::Gamma   , short_t>::type targetDim_impl() const { return 1; };
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Strain  ||
                            _out==MaterialOutput::StrainN ||
                            _out==MaterialOutput::StrainM   , short_t>::type targetDim_impl() const { return 3; };
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Stress  ||
                            _out==MaterialOutput::StressN ||
                            _out==MaterialOutput::StressM ||
                            _out==MaterialOutput::CauchyStressN ||
                            _out==MaterialOutput::CauchyStressM   , short_t>::type targetDim_impl() const { return 3; };
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Thickness, short_t>::type targetDim_impl() const { return 1; };
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Parameters, short_t>::type targetDim_impl() const { return m_materialMat->numParameters(); };
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Deformation, short_t>::type targetDim_impl() const { return 9; };
 
protected:
    void setPatch(index_t p) {m_pIndex = p; }
 
public:
 
    void eval_into(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
private:
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Density   , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::VectorN ||
                            _out==MaterialOutput::VectorM ||
                            _out==MaterialOutput::Generic , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::CauchyVectorN ||
                            _out==MaterialOutput::CauchyVectorM , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::MatrixA ||
                            _out==MaterialOutput::MatrixB ||
                            _out==MaterialOutput::MatrixC ||
                            _out==MaterialOutput::MatrixD   , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::PStress  ||
                            _out==MaterialOutput::PStressN ||
                            _out==MaterialOutput::PStressM  , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::PCauchyStressN ||
                            _out==MaterialOutput::PCauchyStressM  , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::PStrainN ||
                            _out==MaterialOutput::PStrainM  , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Stretch   , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::StretchDir, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::PStressDir, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::StretchTransform, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::PStressTransform, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Spec2CovTransform, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Spec2ConTransform, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Cov2CartTransform, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Con2CartTransform, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::TensionField, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Theta, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Gamma, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Strain  ||
                            _out==MaterialOutput::StrainN ||
                            _out==MaterialOutput::StrainM   , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Stress  ||
                            _out==MaterialOutput::StressN ||
                            _out==MaterialOutput::StressM   , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::CauchyStress  ||
                            _out==MaterialOutput::CauchyStressN ||
                            _out==MaterialOutput::CauchyStressM   , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Thickness, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Parameters, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
    template<enum MaterialOutput _out>
    typename std::enable_if<_out==MaterialOutput::Deformation, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
 
protected:
    index_t m_pIndex;
    gsMaterialMatrixBase<T> * m_materialMat;
    gsMatrix<T> m_z;
};
 
} // namespace
 
 
#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsMaterialMatrixEval.hpp)
#endif