IntegrationGaussLegendreTri.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 "BaseLib/Error.h"
#include "MathLib/Integration/GaussLegendreTri.h"
#include "MathLib/WeightedPoint.h"
 
namespace NumLib
{
class IntegrationGaussLegendreTri
{
public:
    explicit IntegrationGaussLegendreTri(unsigned order = 2) : _order(order)
    {
        this->setIntegrationOrder(order);
    }

    void setIntegrationOrder(unsigned order)
    {
        _order = order;
        _n_sampl_pt = getNumberOfPoints(order);
    }

    unsigned getIntegrationOrder() const { return _order; }

    unsigned getNumberOfPoints() const { return _n_sampl_pt; }

    MathLib::WeightedPoint getWeightedPoint(unsigned const igp) const
    {
        return getWeightedPoint(getIntegrationOrder(), igp);
    }

    static MathLib::WeightedPoint getWeightedPoint(unsigned const order, unsigned const igp)
    {
        switch (order)
        {
            case 1:
                return getWeightedPoint<MathLib::GaussLegendreTri<1>>(igp);
            case 2:
                return getWeightedPoint<MathLib::GaussLegendreTri<2>>(igp);
            case 3:
                return getWeightedPoint<MathLib::GaussLegendreTri<3>>(igp);
            case 4:
                return getWeightedPoint<MathLib::GaussLegendreTri<4>>(igp);
        }
        OGS_FATAL("Integration order {:d} not implemented for triangles.",
                  order);
    }
 
    template <typename Method>
    static MathLib::WeightedPoint getWeightedPoint(unsigned const igp)
    {
        return MathLib::WeightedPoint(Method::X[igp], 0.5 * Method::W[igp]);
    }

    static unsigned getNumberOfPoints(unsigned order)
    {
        switch (order)
        {
            case 1:
                return MathLib::GaussLegendreTri<1>::NPoints;
            case 2:
                return MathLib::GaussLegendreTri<2>::NPoints;
            case 3:
                return MathLib::GaussLegendreTri<3>::NPoints;
            case 4:
                return MathLib::GaussLegendreTri<4>::NPoints;
        }
        OGS_FATAL("Integration order {:d} not implemented for triangles.",
                  order);
    }
 
private:
    unsigned _order;
    unsigned _n_sampl_pt{0};
};
 
}  // namespace NumLib