TriQuadrature.cpp

Go to the documentation of this file.

#include "TriQuadrature.hpp"
#include "LineQuadrature.hpp"
 
#include <vector>
#include <cassert>
#include <cmath>
#include <stdexcept>
 
namespace polyfem
{
    namespace quadrature
    {
        namespace
        {
            void get_weight_and_points(const int order, const bool use_corner_quadrature, Eigen::MatrixXd &points, Eigen::VectorXd &weights)
            {
                if (use_corner_quadrature)
                {
                    switch (order)
                    {
#include <polyfem/autogen/auto_triangle_corner.ipp>
 
                    default:
                        throw std::runtime_error("TriQuadrature: unsupported corner quadrature order " + std::to_string(order));
                    };
                }
                else
                {
                    switch (order)
                    {
#include <polyfem/autogen/auto_triangle.ipp>
 
                    default:
                        throw std::runtime_error("TriQuadrature: unsupported quadrature order " + std::to_string(order));
                    };
                }
            }
        } // namespace
 
        TriQuadrature::TriQuadrature(bool use_corner_quadrature) : use_corner_quadrature_(use_corner_quadrature)
        {
        }
 
        void TriQuadrature::get_quadrature(const int order, Quadrature &quad)
        {
            Quadrature tmp;
 
            get_weight_and_points(order, use_corner_quadrature_, quad.points, quad.weights);
 
            assert(use_corner_quadrature_ || fabs(quad.weights.sum() - 1) < 1e-14);
            assert(quad.points.minCoeff() >= 0 && quad.points.maxCoeff() <= 1);
 
            assert(quad.points.rows() == quad.weights.size());
 
            quad.weights /= 2;
        }
    } // namespace quadrature
} // namespace polyfem