IntegrationTypeTensorProduct.h

Go to the documentation of this file.

#pragma once

#include <vector>
#include "nuto/mechanics/integrationtypes/IntegrationTypeBase.h"
#include "nuto/math/shapes/Line.h"
#include "nuto/math/shapes/Quadrilateral.h"
#include "nuto/math/shapes/Hexahedron.h"

namespace NuTo
{

enum class eIntegrationMethod
{
    GAUSS,
    LOBATTO
};


template <int TDim>
class IntegrationTypeTensorProduct : public IntegrationTypeBase
{
public:
    IntegrationTypeTensorProduct(size_t numIps, NuTo::eIntegrationMethod method);

    Eigen::VectorXd GetLocalIntegrationPointCoordinates(int rIpNum) const override;

    int GetNumIntegrationPoints() const override;

    double GetIntegrationPointWeight(int rIpNum) const override;

    const Shape& GetShape() const override
    {
        return mShape;
    }

private:
    std::vector<double> mIPts1D;
    std::vector<Eigen::Matrix<double, TDim, 1>> mIPts;
    std::vector<double> mWeights;

    using ShapeType =
            typename std::conditional<TDim == 1, Line,
                                      typename std::conditional<TDim == 2, Quadrilateral, Hexahedron>::type>::type;

    ShapeType mShape;
};
} // namespace