Detailed Description

template<typename ShapeFunction, int GlobalDim>
class ProcessLib::VolumetricSourceTermLocalAssembler< ShapeFunction, GlobalDim >

Definition at line 37 of file VolumetricSourceTermFEM.h.

#include <VolumetricSourceTermFEM.h>

Inheritance diagram for ProcessLib::VolumetricSourceTermLocalAssembler< ShapeFunction, GlobalDim >:

[legend]

Collaboration diagram for ProcessLib::VolumetricSourceTermLocalAssembler< ShapeFunction, GlobalDim >:

[legend]

Public Member Functions
	VolumetricSourceTermLocalAssembler (MeshLib::Element const &element, std::size_t const local_matrix_size, NumLib::GenericIntegrationMethod const &integration_method, bool const is_axially_symmetric, ParameterLib::Parameter< double > const &volumetric_source_term)

void	integrate (std::size_t const id, NumLib::LocalToGlobalIndexMap const &source_term_dof_table, double const t, GlobalVector &b) override

Public Member Functions inherited from ProcessLib::VolumetricSourceTermLocalAssemblerInterface
virtual	~VolumetricSourceTermLocalAssemblerInterface ()=default

Private Types
using	ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, GlobalDim>

using	LocalAssemblerTraits

using	NodalVectorType = typename LocalAssemblerTraits::LocalVector

using	NodalRowVectorType = typename ShapeMatricesType::NodalRowVectorType

Private Attributes
ParameterLib::Parameter< double > const &	_volumetric_source_term

NumLib::GenericIntegrationMethod const &	_integration_method

std::vector< SourceTermIntegrationPointData< NodalRowVectorType >, Eigen::aligned_allocator< SourceTermIntegrationPointData< NodalRowVectorType > > >	_ip_data

MeshLib::Element const &	_element

NodalVectorType	_local_rhs

Static Private Attributes
static const unsigned	NUM_NODAL_DOF = 1

Member Typedef Documentation

◆ LocalAssemblerTraits

template<typename ShapeFunction , int GlobalDim>

using ProcessLib::VolumetricSourceTermLocalAssembler< ShapeFunction, GlobalDim >::LocalAssemblerTraits

private

Initial value:

ProcessLib::LocalAssemblerTraits<

ShapeMatricesType, ShapeFunction::NPOINTS, NUM_NODAL_DOF, GlobalDim>

ProcessLib::VolumetricSourceTermLocalAssembler::ShapeMatricesType

ShapeMatrixPolicyType< ShapeFunction, GlobalDim > ShapeMatricesType

Definition VolumetricSourceTermFEM.h:42

ProcessLib::VolumetricSourceTermLocalAssembler::NUM_NODAL_DOF

static const unsigned NUM_NODAL_DOF

Definition VolumetricSourceTermFEM.h:40

ProcessLib::detail::LocalAssemblerTraitsFixed

Definition LocalAssemblerTraits.h:37

Definition at line 44 of file VolumetricSourceTermFEM.h.

◆ NodalRowVectorType

template<typename ShapeFunction , int GlobalDim>

using ProcessLib::VolumetricSourceTermLocalAssembler< ShapeFunction, GlobalDim >::NodalRowVectorType = typename ShapeMatricesType::NodalRowVectorType

private

Definition at line 48 of file VolumetricSourceTermFEM.h.

◆ NodalVectorType

template<typename ShapeFunction , int GlobalDim>

using ProcessLib::VolumetricSourceTermLocalAssembler< ShapeFunction, GlobalDim >::NodalVectorType = typename LocalAssemblerTraits::LocalVector

private

Definition at line 47 of file VolumetricSourceTermFEM.h.

◆ ShapeMatricesType

template<typename ShapeFunction , int GlobalDim>

using ProcessLib::VolumetricSourceTermLocalAssembler< ShapeFunction, GlobalDim >::ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, GlobalDim>

private

Definition at line 42 of file VolumetricSourceTermFEM.h.

Constructor & Destructor Documentation

◆ VolumetricSourceTermLocalAssembler()

template<typename ShapeFunction , int GlobalDim>

ProcessLib::VolumetricSourceTermLocalAssembler< ShapeFunction, GlobalDim >::VolumetricSourceTermLocalAssembler	(	MeshLib::Element const &	element,
		std::size_t const	local_matrix_size,
		NumLib::GenericIntegrationMethod const &	integration_method,
		bool const	is_axially_symmetric,
		ParameterLib::Parameter< double > const &	volumetric_source_term )

inline

Definition at line 51 of file VolumetricSourceTermFEM.h.

        : _volumetric_source_term(volumetric_source_term),
          _integration_method(integration_method),
          _element(element),
          _local_rhs(local_matrix_size)
    {
        unsigned const n_integration_points =
            _integration_method.getNumberOfPoints();
 
        auto const shape_matrices =
            NumLib::initShapeMatrices<ShapeFunction, ShapeMatricesType,
                                      GlobalDim>(_element, is_axially_symmetric,
                                                 _integration_method);
 
        for (unsigned ip = 0; ip < n_integration_points; ip++)
        {
            _ip_data.emplace_back(
                shape_matrices[ip].N,
                _integration_method.getWeightedPoint(ip).getWeight() *
                    shape_matrices[ip].integralMeasure *
                    shape_matrices[ip].detJ);
        }
    }

References ProcessLib::VolumetricSourceTermLocalAssembler< ShapeFunction, GlobalDim >::_element, ProcessLib::VolumetricSourceTermLocalAssembler< ShapeFunction, GlobalDim >::_integration_method, ProcessLib::VolumetricSourceTermLocalAssembler< ShapeFunction, GlobalDim >::_ip_data, NumLib::GenericIntegrationMethod::getNumberOfPoints(), MathLib::WeightedPoint::getWeight(), NumLib::GenericIntegrationMethod::getWeightedPoint(), and NumLib::initShapeMatrices().

Member Function Documentation

◆ integrate()

template<typename ShapeFunction , int GlobalDim>

void ProcessLib::VolumetricSourceTermLocalAssembler< ShapeFunction, GlobalDim >::integrate	(	std::size_t const	id,
		NumLib::LocalToGlobalIndexMap const &	source_term_dof_table,
		double const	t,
		GlobalVector &	b )

inlineoverridevirtual

Implements ProcessLib::VolumetricSourceTermLocalAssemblerInterface.

Definition at line 80 of file VolumetricSourceTermFEM.h.

    {
        _local_rhs.setZero();
 
        unsigned const n_integration_points =
            _integration_method.getNumberOfPoints();
 
        for (unsigned ip = 0; ip < n_integration_points; ip++)
        {
            auto const& N = _ip_data[ip].N;
            auto const& w = _ip_data[ip].integration_weight;
 
            ParameterLib::SpatialPosition const pos{
                std::nullopt, _element.getID(), ip,
                MathLib::Point3d(
                    NumLib::interpolateCoordinates<ShapeFunction,
                                                   ShapeMatricesType>(_element,
                                                                      N))};
            auto const st_val = _volumetric_source_term(t, pos)[0];
 
            _local_rhs.noalias() += N.transpose() * st_val * w;
        }
        auto const indices = NumLib::getIndices(id, source_term_dof_table);
        b.add(indices, _local_rhs);
    }