Detailed Description

Definition at line 37 of file DeactivatedSubdomainDirichlet.h.

#include <DeactivatedSubdomainDirichlet.h>

Inheritance diagram for ProcessLib::DeactivatedSubdomainDirichlet:

Collaboration diagram for ProcessLib::DeactivatedSubdomainDirichlet:

Public Member Functions
	DeactivatedSubdomainDirichlet (std::vector< std::size_t > const *active_element_ids, MathLib::PiecewiseLinearInterpolation time_interval, ParameterLib::Parameter< double > const &parameter, bool const set_outer_nodes_dirichlet_values, DeactivatedSubdomainMesh const &subdomain, NumLib::LocalToGlobalIndexMap const &dof_table_bulk, int const variable_id, int const component_id)

void	getEssentialBCValues (const double t, GlobalVector const &x, NumLib::IndexValueVector< GlobalIndexType > &bc_values) const override
	Writes the values of essential BCs to `bc_values`. More...

Public Member Functions inherited from ProcessLib::BoundaryCondition
virtual void	applyNaturalBC (const double, std::vector< GlobalVector * > const &, int const, GlobalMatrix &, GlobalVector &, GlobalMatrix *)

virtual void	preTimestep (const double, std::vector< GlobalVector * > const &, int const)

virtual void	postTimestep (const double, std::vector< GlobalVector * > const &, int const)

virtual	~BoundaryCondition ()=default

Private Member Functions
void	config (NumLib::LocalToGlobalIndexMap const &dof_table_bulk)

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

DeactivatedSubdomainMesh const &	_subdomain

std::unique_ptr< NumLib::LocalToGlobalIndexMap const >	_dof_table_boundary

int const	_variable_id

int const	_component_id

MathLib::PiecewiseLinearInterpolation const	_time_interval

std::vector< std::size_t > const *	_active_element_ids = nullptr

bool const	_set_outer_nodes_dirichlet_values

Constructor & Destructor Documentation

◆ DeactivatedSubdomainDirichlet()

ProcessLib::DeactivatedSubdomainDirichlet::DeactivatedSubdomainDirichlet	(	std::vector< std::size_t > const *	active_element_ids,
		MathLib::PiecewiseLinearInterpolation	time_interval,
		ParameterLib::Parameter< double > const &	parameter,
		bool const	set_outer_nodes_dirichlet_values,
		DeactivatedSubdomainMesh const &	subdomain,
		NumLib::LocalToGlobalIndexMap const &	dof_table_bulk,
		int const	variable_id,
		int const	component_id
	)

Definition at line 20 of file DeactivatedSubdomainDirichlet.cpp.

     : _parameter(parameter),
       _subdomain(subdomain),
       _variable_id(variable_id),
       _component_id(component_id),
       _time_interval(std::move(time_interval)),
       _active_element_ids(active_element_ids),
       _set_outer_nodes_dirichlet_values(set_outer_nodes_dirichlet_values)
 {
     config(dof_table_bulk);
 }

References config().

Member Function Documentation

◆ config()

void ProcessLib::DeactivatedSubdomainDirichlet::config ( NumLib::LocalToGlobalIndexMap const & dof_table_bulk )

private

Definition at line 39 of file DeactivatedSubdomainDirichlet.cpp.

 {
     checkParametersOfDirichletBoundaryCondition(
         *_subdomain.mesh, dof_table_bulk, _variable_id, _component_id);
  
     std::vector<MeshLib::Node*> const& bc_nodes = _subdomain.mesh->getNodes();
     MeshLib::MeshSubset subdomain_mesh_subset(*_subdomain.mesh, bc_nodes);
  
     // Create local DOF table from the BC mesh subset for the given variable
     // and component id.
     _dof_table_boundary.reset(dof_table_bulk.deriveBoundaryConstrainedMap(
         _variable_id, {_component_id}, std::move(subdomain_mesh_subset)));
 }

References _component_id, _dof_table_boundary, _subdomain, _variable_id, ProcessLib::checkParametersOfDirichletBoundaryCondition(), NumLib::LocalToGlobalIndexMap::deriveBoundaryConstrainedMap(), and ProcessLib::DeactivatedSubdomainMesh::mesh.

Referenced by DeactivatedSubdomainDirichlet().

◆ getEssentialBCValues()

void ProcessLib::DeactivatedSubdomainDirichlet::getEssentialBCValues	(	const double	,
		GlobalVector const &	,
		NumLib::IndexValueVector< GlobalIndexType > &
	)		const

overridevirtual

Writes the values of essential BCs to bc_values.

Reimplemented from ProcessLib::BoundaryCondition.

Definition at line 54 of file DeactivatedSubdomainDirichlet.cpp.

 {
     auto const& bulk_element_ids =
         *_subdomain.mesh->getProperties()
              .template getPropertyVector<std::size_t>(
                  "bulk_element_ids", MeshLib::MeshItemType::Cell, 1);
  
     auto is_inactive = [&](MeshLib::Element const* const e)
     {
         return none_of(begin(*_active_element_ids), end(*_active_element_ids),
                        [&](auto const id)
                        { return id == bulk_element_ids[e->getID()]; });
     };
  
     std::vector<MeshLib::Node*> inactive_nodes_in_bc_mesh;
     std::copy_if(begin(_subdomain.inner_nodes), end(_subdomain.inner_nodes),
                  back_inserter(inactive_nodes_in_bc_mesh),
                  [&](MeshLib::Node* const n)
                  {
                      const auto& connected_elements =
                          _subdomain.mesh->getElementsConnectedToNode(*n);
  
                      return std::all_of(begin(connected_elements),
                                         end(connected_elements), is_inactive);
                  });
  
     if (_set_outer_nodes_dirichlet_values)
     {
         std::copy_if(begin(_subdomain.outer_nodes), end(_subdomain.outer_nodes),
                      back_inserter(inactive_nodes_in_bc_mesh),
                      [&](MeshLib::Node* const n)
                      {
                          const auto& connected_elements =
                              _subdomain.mesh->getElementsConnectedToNode(*n);
  
                          return std::all_of(begin(connected_elements),
                                             end(connected_elements),
                                             is_inactive);
                      });
     }
  
     auto time_interval_contains = [&](double const t)
     {
         return _time_interval.getSupportMin() <= t &&
                t <= _time_interval.getSupportMax();
     };
     if (time_interval_contains(t))
     {
         getEssentialBCValuesLocal(
             _parameter, *_subdomain.mesh, inactive_nodes_in_bc_mesh,
             *_dof_table_boundary, _variable_id, _component_id, t, x, bc_values);
         return;
     }
  
     bc_values.ids.clear();
     bc_values.values.clear();
 }