GMatrixPolicy.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 "MathLib/VectorizedTensor.h"
#include "NumLib/Fem/ShapeMatrixPolicy.h"
 
namespace ProcessLib
{
template <typename ShapeFunction, int DisplacementDim>
class GMatrixPolicyType
{
private:
    template <int N>
    using VectorTypeFixedSize = typename EigenFixedShapeMatrixPolicy<
        ShapeFunction, DisplacementDim>::template VectorType<N>;

    template <int N, int M>
    using MatrixType = typename ShapeMatrixPolicyType<
        ShapeFunction, DisplacementDim>::template MatrixType<N, M>;
 
    // Dimensions of specific b-matrix for n-points and displacement dimension.
    static int const _number_of_dof = ShapeFunction::NPOINTS * DisplacementDim;
 
public:
    // For the 2D case the 33-component is needed (and the four entries
    // of the non-symmetric matrix); In 3d there are nine entries.
    using GradientMatrixType =
        MatrixType<MathLib::VectorizedTensor::size(DisplacementDim),
                   _number_of_dof>;
    using GradientVectorType =
        VectorTypeFixedSize<MathLib::VectorizedTensor::size(DisplacementDim)>;
};
}  // namespace ProcessLib