Detailed Description

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

Definition at line 28 of file RobinBoundaryConditionLocalAssembler.h.

#include <RobinBoundaryConditionLocalAssembler.h>

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

[legend]

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

[legend]

Public Member Functions
	RobinBoundaryConditionLocalAssembler (MeshLib::Element const &e, std::size_t const local_matrix_size, NumLib::GenericIntegrationMethod const &integration_method, bool is_axially_symmetric, RobinBoundaryConditionData const &data)
void	assemble (std::size_t const id, NumLib::LocalToGlobalIndexMap const &dof_table_boundary, double const t, std::vector< GlobalVector * > const &xs, int const process_id, GlobalMatrix K, GlobalVector &b, GlobalMatrix Jac) override
Public Member Functions inherited from ProcessLib::GenericNaturalBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >
	GenericNaturalBoundaryConditionLocalAssembler (MeshLib::Element const &e, bool is_axially_symmetric, NumLib::GenericIntegrationMethod const &integration_method)
Public Member Functions inherited from ProcessLib::GenericNaturalBoundaryConditionLocalAssemblerInterface
virtual	~GenericNaturalBoundaryConditionLocalAssemblerInterface ()=default

Public Attributes
	EIGEN_MAKE_ALIGNED_OPERATOR_NEW

Private Types
using	Base
using	ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, GlobalDim>

Private Attributes
RobinBoundaryConditionData const &	_data
Base::NodalMatrixType	_local_K
Base::NodalVectorType	_local_rhs

Additional Inherited Members
Protected Types inherited from ProcessLib::GenericNaturalBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >
using	ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, GlobalDim>
using	NodalMatrixType = typename ShapeMatricesType::NodalMatrixType
using	NodalVectorType = typename ShapeMatricesType::NodalVectorType
Protected Attributes inherited from ProcessLib::GenericNaturalBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >
NumLib::GenericIntegrationMethod const &	_integration_method
std::vector< NAndWeight, Eigen::aligned_allocator< NAndWeight > > const	_ns_and_weights
MeshLib::Element const &	_element

Member Typedef Documentation

◆ Base

template<typename ShapeFunction, int GlobalDim>

using ProcessLib::RobinBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >::Base

private

Initial value:

GenericNaturalBoundaryConditionLocalAssembler<ShapeFunction, GlobalDim>

ProcessLib::GenericNaturalBoundaryConditionLocalAssembler::GenericNaturalBoundaryConditionLocalAssembler

GenericNaturalBoundaryConditionLocalAssembler(MeshLib::Element const &e, bool is_axially_symmetric, NumLib::GenericIntegrationMethod const &integration_method)

Definition GenericNaturalBoundaryConditionLocalAssembler.h:79

Definition at line 32 of file RobinBoundaryConditionLocalAssembler.h.

◆ ShapeMatricesType

template<typename ShapeFunction, int GlobalDim>

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

private

Definition at line 34 of file RobinBoundaryConditionLocalAssembler.h.

Constructor & Destructor Documentation

◆ RobinBoundaryConditionLocalAssembler()

template<typename ShapeFunction, int GlobalDim>

ProcessLib::RobinBoundaryConditionLocalAssembler< ShapeFunction, GlobalDim >::RobinBoundaryConditionLocalAssembler	(	MeshLib::Element const &	e,
		std::size_t const	local_matrix_size,
		NumLib::GenericIntegrationMethod const &	integration_method,
		bool	is_axially_symmetric,
		RobinBoundaryConditionData const &	data )

inline

Definition at line 37 of file RobinBoundaryConditionLocalAssembler.h.

        : Base(e, is_axially_symmetric, integration_method),
          _data(data),
          _local_K(local_matrix_size, local_matrix_size),
          _local_rhs(local_matrix_size)
    {
    }

References _data, _local_K, and _local_rhs.

Member Function Documentation

◆ assemble()

template<typename ShapeFunction, int GlobalDim>

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

inlineoverridevirtual

Implements ProcessLib::GenericNaturalBoundaryConditionLocalAssemblerInterface.

Definition at line 50 of file RobinBoundaryConditionLocalAssembler.h.

    {
        _local_K.setZero();
        _local_rhs.setZero();
 
        auto& x = *xs[process_id];
 
        auto const indices = NumLib::getIndices(id, dof_table_boundary);
        auto const local_x = x.get(indices);
        auto const u =
            MathLib::toVector<Eigen::Matrix<double, ShapeFunction::NPOINTS, 1>>(
                local_x, ShapeFunction::NPOINTS);
 
        unsigned const n_integration_points =
            Base::_integration_method.getNumberOfPoints();
 
        typename Base::NodalVectorType const alpha =
            _data.alpha.getNodalValuesOnElement(Base::_element, t)
                .template topRows<ShapeFunction::MeshElement::n_all_nodes>();
        typename Base::NodalVectorType const u_0 =
            _data.u_0.getNodalValuesOnElement(Base::_element, t)
                .template topRows<ShapeFunction::MeshElement::n_all_nodes>();
 
        for (unsigned ip = 0; ip < n_integration_points; ++ip)
        {
            auto const& ip_data = Base::_ns_and_weights[ip];
            auto const& N = ip_data.N;
            auto const& w = ip_data.weight;
 
            ParameterLib::SpatialPosition const position{
                std::nullopt, Base::_element.getID(),
                MathLib::Point3d(
                    NumLib::interpolateCoordinates<ShapeFunction,
                                                   ShapeMatricesType>(
                        Base::_element, N))};
 
            double integral_measure = 1.0;
            if (_data.integral_measure)
            {
                integral_measure = (*_data.integral_measure)(t, position)[0];
            }
 
            double const a = alpha.dot(N) * w * integral_measure;
 
            // The local K matrix is used for both, the Newton and Picard
            // methods.
            _local_K.noalias() += N.transpose() * N * a;
 
            if (Jac != nullptr)
            {
                _local_rhs.noalias() -= N.transpose() * (u - u_0).dot(N) * a;
            }
            else
            {
                _local_rhs.noalias() += N.transpose() * (u_0.dot(N) * a);
            }
        }
 
        b.add(indices, _local_rhs);
        if (Jac != nullptr)
        {
            Jac->add(NumLib::LocalToGlobalIndexMap::RowColumnIndices(indices,
                                                                     indices),
                     _local_K);
        }
        else if (K != nullptr)
        {
            K->add(NumLib::LocalToGlobalIndexMap::RowColumnIndices(indices,
                                                                   indices),
                   _local_K);
        }
        else
        {
            OGS_FATAL(
                "In Robin boundary condition assembler, both, the Jacobian and "
                "the K-matrices are null, but one matrix must be provided.");
        }
    }