fn_intpoints_gauss_tri7.hpp

Go to the documentation of this file.

/*
 * BELFEM -- The Berkeley Lab Finite Element Framework
 * Copyright (c) 2026, The Regents of the University of California,
 * through Lawrence Berkeley National Laboratory (subject to receipt of any required
 * approvals from the U.S. Dept. of Energy).  All rights reserved.
 *
 * Developers: Christian Messe, Gregory Giard
 *
 * See the top-level LICENSE file for the complete license and disclaimer.
 */
 
#ifndef BELFEM_FN_INTPOINTS_GAUSS_TRI7_HPP
#define BELFEM_FN_INTPOINTS_GAUSS_TRI7_HPP
 
#include "typedefs.hpp"
#include "cl_Vector.hpp"
#include "cl_Matrix.hpp"
#include "fn_sum.hpp"
 
namespace belfem
{
    namespace integration
    {
// ----------------------------------------------------------------------------

        inline void
        gauss_tri7(
                Vector< real > & aWeights,
                Matrix< real > & aPoints )
        {
            aPoints.set_size( 3, 7 );
 
            aPoints( 0, 0 ) = 1./3.;
            aPoints( 1, 0 ) = 1./3.;
 
            aPoints( 0, 1 ) = 0.10128650732345634;
            aPoints( 1, 1 ) = 0.10128650732345634;
 
            aPoints( 0, 2 ) = 0.4701420641051151;
            aPoints( 1, 2 ) = 0.4701420641051151;
 
            aPoints( 0, 3 ) = 0.10128650732345634;
            aPoints( 1, 3 ) = 0.7974269853530873;
 
            aPoints( 0, 4 ) = 0.4701420641051151;
            aPoints( 1, 4 ) = 0.05971587178976982;
 
            aPoints( 0, 5 ) = 0.7974269853530873;
            aPoints( 1, 5 ) = 0.10128650732345634;
 
            aPoints( 0, 6 ) = 0.05971587178976982;
            aPoints( 1, 6 ) = 0.4701420641051151;
 
            for( uint k=0; k<7; ++k )
            {
                aPoints( 2, k ) =  1.0
                        - aPoints( 0, k )
                        - aPoints( 1, k );
            }
 
            aWeights.set_size( 7 );
 
            aWeights( 0 ) = 0.1125;
            aWeights( 0 ) = 0.0 ;
            aWeights( 1 ) = 0.06296959027241358;
            aWeights( 2 ) = 0.0661970763942531;
            aWeights( 3 ) = 0.06296959027241358;
            aWeights( 4 ) = 0.0661970763942531;
            aWeights( 5 ) = 0.06296959027241358;
            aWeights( 6 ) = 0.0661970763942531;
 
            aWeights( 0 ) = 0.0;
            aWeights( 0 ) = 0.5 - sum( aWeights );
        }
 
// ----------------------------------------------------------------------------
    } /* namespace integration */
} /* end namespace belfem */
 
#endif  // BELFEM_FN_INTPOINTS_GAUSS_TRI7_HPP