ShapeMatrices.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include <Eigen/Core>
#include <iosfwd>
 
namespace NumLib
{
enum class ShapeMatrixType
{
    N,       
    DNDR,    
    N_J,     
    DNDR_J,  
    DNDX,    
    ALL      
};

template <class T_N, class T_DNDR, class T_J, class T_DNDX>
struct ShapeMatrices
{
    using ShapeType = T_N;
    using DrShapeType = T_DNDR;
    using JacobianType = T_J;
    using DxShapeType = T_DNDX;
 
    ShapeType N;        
    DrShapeType dNdr;   
    JacobianType J;     
    double detJ;        
    JacobianType invJ;  
    DxShapeType dNdx;   
    double integralMeasure;

    ShapeMatrices() = delete;

    ShapeMatrices(std::size_t local_dim, std::size_t global_dim,
                  std::size_t n_nodes)
        : N(n_nodes),
          dNdr(local_dim, n_nodes),
          J(local_dim, local_dim),
          detJ(.0),
          invJ(local_dim, local_dim),
          dNdx(global_dim, n_nodes),
          integralMeasure(0.0)
    {
        setZero();
    }

    void setZero();

    template <ShapeMatrixType T_SHAPE_MATRIX_TYPE>
    void setZero();

    void write(std::ostream& out) const;
 
    EIGEN_MAKE_ALIGNED_OPERATOR_NEW;
};  // ShapeMatrices
 
}  // namespace NumLib
 
#include "ShapeMatrices-impl.h"