Detailed Description

template<typename ShapeFunction, typename LowerOrderShapeFunction, int GlobalDim>
class ProcessLib::PythonBoundaryConditionLocalAssembler< ShapeFunction, LowerOrderShapeFunction, GlobalDim >

Definition at line 21 of file PythonBoundaryConditionLocalAssembler.h.

#include <PythonBoundaryConditionLocalAssembler.h>

Inheritance diagram for ProcessLib::PythonBoundaryConditionLocalAssembler< ShapeFunction, LowerOrderShapeFunction, GlobalDim >:

[legend]

Collaboration diagram for ProcessLib::PythonBoundaryConditionLocalAssembler< ShapeFunction, LowerOrderShapeFunction, GlobalDim >:

[legend]

Public Member Functions
	PythonBoundaryConditionLocalAssembler (MeshLib::Element const &e, std::size_t const, NumLib::GenericIntegrationMethod const &integration_method, bool is_axially_symmetric, PythonBcData const &data)
void	assemble (std::size_t const boundary_element_id, NumLib::LocalToGlobalIndexMap const &dof_table_boundary, double const t, std::vector< GlobalVector * > const &xs, int const process_id, GlobalMatrix , GlobalVector &b, GlobalMatrix const Jac) override
double	interpolate (unsigned const local_node_id, NumLib::LocalToGlobalIndexMap const &dof_table_boundary, GlobalVector const &x, int const var, int const comp) const override
Public Member Functions inherited from ProcessLib::GenericNaturalBoundaryConditionLocalAssemblerInterface
virtual	~GenericNaturalBoundaryConditionLocalAssemblerInterface ()=default

Private Types
using	LocAsmImpl
using	Traits = typename LocAsmImpl::Traits

Private Member Functions
Traits::LowerOrderShapeMatrix	computeLowerOrderShapeMatrix (unsigned const local_node_id) const

Private Attributes
LocAsmImpl const	impl_

Member Typedef Documentation

◆ LocAsmImpl

template<typename ShapeFunction, typename LowerOrderShapeFunction, int GlobalDim>

using ProcessLib::PythonBoundaryConditionLocalAssembler< ShapeFunction, LowerOrderShapeFunction, GlobalDim >::LocAsmImpl

private

Initial value:

 ProcessLib::BoundaryConditionAndSourceTerm::Python::
        BcAndStLocalAssemblerImpl<PythonBcData, ShapeFunction,
                                  LowerOrderShapeFunction, GlobalDim>

Definition at line 24 of file PythonBoundaryConditionLocalAssembler.h.

◆ Traits

template<typename ShapeFunction, typename LowerOrderShapeFunction, int GlobalDim>

using ProcessLib::PythonBoundaryConditionLocalAssembler< ShapeFunction, LowerOrderShapeFunction, GlobalDim >::Traits = typename LocAsmImpl::Traits

private

Definition at line 27 of file PythonBoundaryConditionLocalAssembler.h.

Constructor & Destructor Documentation

◆ PythonBoundaryConditionLocalAssembler()

template<typename ShapeFunction, typename LowerOrderShapeFunction, int GlobalDim>

ProcessLib::PythonBoundaryConditionLocalAssembler< ShapeFunction, LowerOrderShapeFunction, GlobalDim >::PythonBoundaryConditionLocalAssembler	(	MeshLib::Element const &	e,
		std::size_t const	,
		NumLib::GenericIntegrationMethod const &	integration_method,
		bool	is_axially_symmetric,
		PythonBcData const &	data )

inline

Definition at line 30 of file PythonBoundaryConditionLocalAssembler.h.

        : impl_{e, integration_method, is_axially_symmetric, data}
    {
    }

References impl_.

Member Function Documentation

◆ assemble()

template<typename ShapeFunction, typename LowerOrderShapeFunction, int GlobalDim>

void ProcessLib::PythonBoundaryConditionLocalAssembler< ShapeFunction, LowerOrderShapeFunction, GlobalDim >::assemble	(	std::size_t const	boundary_element_id,
		NumLib::LocalToGlobalIndexMap const &	dof_table_boundary,
		double const	t,
		std::vector< GlobalVector * > const &	xs,
		int const	process_id,
		GlobalMatrix *	,
		GlobalVector &	b,
		GlobalMatrix *const	Jac )

inlineoverridevirtual

Implements ProcessLib::GenericNaturalBoundaryConditionLocalAssemblerInterface.

Definition at line 40 of file PythonBoundaryConditionLocalAssembler.h.

    {
        impl_.assemble(boundary_element_id, dof_table_boundary, t,
                       *xs[process_id], b, Jac);
    }

References impl_.

◆ computeLowerOrderShapeMatrix()

template<typename ShapeFunction, typename LowerOrderShapeFunction, int GlobalDim>

Traits::LowerOrderShapeMatrix ProcessLib::PythonBoundaryConditionLocalAssembler< ShapeFunction, LowerOrderShapeFunction, GlobalDim >::computeLowerOrderShapeMatrix ( unsigned const local_node_id ) const

inlineprivate

Definition at line 76 of file PythonBoundaryConditionLocalAssembler.h.

    {
        using HigherOrderMeshElement = typename ShapeFunction::MeshElement;
 
        assert(local_node_id < impl_.element.getNumberOfNodes());
 
        std::array natural_coordss{MathLib::Point3d{NumLib::NaturalCoordinates<
            HigherOrderMeshElement>::coordinates[local_node_id]}};
 
        bool const is_axially_symmetric = false;  // does not matter for N
        auto const shape_matrices = NumLib::computeShapeMatrices<
            LowerOrderShapeFunction,
            typename Traits::LowerOrderShapeMatrixPolicy, GlobalDim,
            NumLib::ShapeMatrixType::N>(impl_.element, is_axially_symmetric,
                                        natural_coordss);
        return shape_matrices.front().N;
    }

References NumLib::computeShapeMatrices(), impl_, and NumLib::N.

Referenced by interpolate().

◆ interpolate()

template<typename ShapeFunction, typename LowerOrderShapeFunction, int GlobalDim>

double ProcessLib::PythonBoundaryConditionLocalAssembler< ShapeFunction, LowerOrderShapeFunction, GlobalDim >::interpolate	(	unsigned const	local_node_id,
		NumLib::LocalToGlobalIndexMap const &	dof_table_boundary,
		GlobalVector const &	x,
		int const	var,
		int const	comp ) const

inlineoverridevirtual

Interpolates the given component of the given variable to the given local_node_id.

The local_node_id is the number of the node within the current boundary element.

Implements ProcessLib::PythonBoundaryConditionLocalAssemblerInterface.

Definition at line 50 of file PythonBoundaryConditionLocalAssembler.h.

    {
        if constexpr (ShapeFunction::ORDER < 2 ||
                      LowerOrderShapeFunction::ORDER > 1)
        {
            return std::numeric_limits<double>::quiet_NaN();
        }
        else
        {
            auto const N = computeLowerOrderShapeMatrix(local_node_id);
 
            auto const nodal_values_base_node =
                ProcessLib::BoundaryConditionAndSourceTerm::Python::
                    collectDofsToMatrixOnBaseNodesSingleComponent(
                        impl_.element,
                        impl_.bc_or_st_data.bc_or_st_mesh.getID(),
                        dof_table_boundary, x, var, comp);
 
            return N * nodal_values_base_node;
        }
    }