SolutionTransfer derived class which is specifically for transferring solutions back and forth between a VolumeMesh and its associated BoundaryMesh. More...

#include <boundary_volume_solution_transfer.h>

Inheritance diagram for libMesh::BoundaryVolumeSolutionTransfer:

[legend]

Public Member Functions
	BoundaryVolumeSolutionTransfer (const Parallel::Communicator &comm_in)

virtual	~BoundaryVolumeSolutionTransfer ()=default

virtual void	transfer (const Variable &from_var, const Variable &to_var) override
	Transfer values from a Variable in a System associated with a volume mesh to a Variable in a System associated with the corresponding BoundaryMesh, or vice-versa. More...

const Parallel::Communicator &	comm () const

processor_id_type	n_processors () const

processor_id_type	processor_id () const

Protected Attributes
const Parallel::Communicator &	_communicator

Private Member Functions
void	transfer_volume_boundary (const Variable &from_var, const Variable &to_var)
	Transfer values from volume mesh to boundary mesh. More...

void	transfer_boundary_volume (const Variable &from_var, const Variable &to_var)
	Transfer values from boundary mesh to volume mesh. More...

Detailed Description

SolutionTransfer derived class which is specifically for transferring solutions back and forth between a VolumeMesh and its associated BoundaryMesh.

That is, when calling transfer() below, we must have either:

.) Transfer from the boundary of a VolumeMesh to its associated BoundaryMesh, or .) Transfer from a BoundaryMesh to its corresponding VolumeMesh.

Because of this assumption, we can use the interior_parent() to directly define the mapping between solutions. The transfer() function internally determines which direction the transfer is going, and calls the appropriate algorithm internally.

Author: Xikai Jiang; John W. Peterson

Date: 2016

Definition at line 47 of file boundary_volume_solution_transfer.h.

Constructor & Destructor Documentation

◆ BoundaryVolumeSolutionTransfer()

libMesh::BoundaryVolumeSolutionTransfer::BoundaryVolumeSolutionTransfer ( const Parallel::Communicator & comm_in )

inline

Definition at line 50 of file boundary_volume_solution_transfer.h.

                                                                         :
     SolutionTransfer(comm_in)
   {}

◆ ~BoundaryVolumeSolutionTransfer()

virtual libMesh::BoundaryVolumeSolutionTransfer::~BoundaryVolumeSolutionTransfer ( )

virtualdefault

Member Function Documentation

◆ comm()

const Parallel::Communicator& libMesh::ParallelObject::comm ( ) const

inlineinherited

Returns: A reference to the Parallel::Communicator object used by this mesh.

Definition at line 97 of file parallel_object.h.

References libMesh::ParallelObject::_communicator.

98 { return _communicator; }

libMesh::ParallelObject::_communicator

const Parallel::Communicator & _communicator

Definition: parallel_object.h:120

◆ n_processors()

processor_id_type libMesh::ParallelObject::n_processors ( ) const

inlineinherited

Returns: The number of processors in the group.

Definition at line 103 of file parallel_object.h.

References libMesh::ParallelObject::_communicator, libMesh::libmesh_assert(), and TIMPI::Communicator::size().

   {
     processor_id_type returnval =
       cast_int<processor_id_type>(_communicator.size());
     libmesh_assert(returnval); // We never have an empty comm
     return returnval;
   }

◆ processor_id()

processor_id_type libMesh::ParallelObject::processor_id ( ) const

inlineinherited

Returns: The rank of this processor in the group.

Definition at line 114 of file parallel_object.h.

References libMesh::ParallelObject::_communicator, and TIMPI::Communicator::rank().

115 { return cast_int<processor_id_type>(_communicator.rank()); }

TIMPI::Communicator::rank

processor_id_type rank() const

libMesh::ParallelObject::_communicator

const Parallel::Communicator & _communicator

Definition: parallel_object.h:120

◆ transfer()

void libMesh::BoundaryVolumeSolutionTransfer::transfer	(	const Variable &	from_var,
		const Variable &	to_var
	)

overridevirtual

Transfer values from a Variable in a System associated with a volume mesh to a Variable in a System associated with the corresponding BoundaryMesh, or vice-versa.

Implements libMesh::SolutionTransfer.

Definition at line 25 of file boundary_volume_solution_transfer.C.

References libMesh::System::get_mesh(), libMesh::MeshBase::mesh_dimension(), libMesh::Variable::system(), transfer_boundary_volume(), and transfer_volume_boundary().

 {
   // Determine which direction the transfer is in
   System * from_sys = from_var.system();
   System * to_sys = to_var.system();
 
   unsigned int
     from_dimension = from_sys->get_mesh().mesh_dimension(),
     to_dimension = to_sys->get_mesh().mesh_dimension();
 
   // Sanity check
   libmesh_error_msg_if(from_dimension == to_dimension,
                        "Error: Transfer must be from volume mesh to its boundary or vice-versa!");
 
   if (from_dimension > to_dimension)
     this->transfer_volume_boundary(from_var, to_var);
   else
     this->transfer_boundary_volume(from_var, to_var);
 }

◆ transfer_boundary_volume()

void libMesh::BoundaryVolumeSolutionTransfer::transfer_boundary_volume	(	const Variable &	from_var,
		const Variable &	to_var
	)

private

Transfer values from boundary mesh to volume mesh.

Called when needed from transfer().

Definition at line 144 of file boundary_volume_solution_transfer.C.

References libMesh::System::current_local_solution, libMesh::DofMap::dof_indices(), libMesh::DofObject::dof_number(), libMesh::System::get_dof_map(), libMesh::System::get_mesh(), libMesh::Elem::interior_parent(), libMesh::Variable::n_components(), libMesh::Elem::node_ref_range(), libMesh::Variable::number(), libMesh::System::number(), libMesh::System::solution, libMesh::Variable::system(), libMesh::TOLERANCE, and value.

Referenced by transfer().

 {
   // Get references to the Systems from the Variables
   System * from_sys = from_var.system(); // boundary system
   System * to_sys = to_var.system();     // volume system
 
   // Get reference to BoundaryMesh.
   const MeshBase & from_mesh = from_sys->get_mesh();
 
   // DofMap for BoundaryMesh
   const DofMap & dof_map = from_sys->get_dof_map();
 
   // Get system number and variable numbers
   const unsigned int to_sys_number = to_sys->number();
   const unsigned int to_var_number = to_var.number();
   const unsigned int from_var_number = from_var.number();
   const unsigned int to_n_comp = to_var.n_components(to_sys->get_mesh());
 
   // In order to get solution vectors from BoundaryMesh
   std::vector<dof_id_type> from_dof_indices;
   std::vector<Number> value;
 
   // Loop through all boundary elements.
   for (const auto & from_elem : from_mesh.active_local_element_ptr_range())
     {
       const Elem * to_elem = from_elem->interior_parent();
 
       libmesh_error_msg_if(!to_elem, "Error, transfer must be between a Mesh and its associated BoundaryMesh.");
 
       // Get dof indices for phi2 for all nodes on this boundary element
       dof_map.dof_indices(from_elem, from_dof_indices, from_var_number);
 
       // Get values from BoundaryMesh for this element
       from_sys->current_local_solution->get(from_dof_indices, value);
 
       // loop through all nodes in each boundary element.
       for (auto node : from_elem->node_index_range())
         {
           // Node in boundary element.
           const Node * from_node = from_elem->node_ptr(node);
 
           // Nodes in interior_parent element.
           for (const Node & to_node : to_elem->node_ref_range())
             {
               // Match BoundaryNode & VolumeNode.
               if (to_node.absolute_fuzzy_equals(*from_node, TOLERANCE))
                 {
                   // Global dof_index for node in VolumeMesh.
                   const dof_id_type to_dof = to_node.dof_number(to_sys_number,
                                                                 to_var_number,
                                                                 to_n_comp - 1);
 
                   // Assign values to boundary in VolumeMesh.
                   to_sys->solution->set(to_dof, value[node]);
                 }
             }
         }
     }
 }

◆ transfer_volume_boundary()

void libMesh::BoundaryVolumeSolutionTransfer::transfer_volume_boundary	(	const Variable &	from_var,
		const Variable &	to_var
	)

private

Transfer values from volume mesh to boundary mesh.

Called when needed from transfer().

Definition at line 49 of file boundary_volume_solution_transfer.C.

References libMesh::System::get_mesh(), libMesh::MeshTools::Generation::Private::idx(), libMesh::Elem::interior_parent(), libMesh::Variable::n_components(), libMesh::Elem::node_ref_range(), libMesh::Variable::number(), libMesh::System::number(), libMesh::System::solution, libMesh::Variable::system(), and libMesh::TOLERANCE.

Referenced by transfer().

 {
   // Get references to the Systems from the Variables
   System * from_sys = from_var.system(); // volume system
   System * to_sys = to_var.system();     // boundary system
 
   // Get reference to the BoundaryMesh.  Note: we always loop over the
   // BoundaryMesh since, by definition, it has fewer dofs than the
   // volume mesh.
   const MeshBase & to_mesh = to_sys->get_mesh();
 
   // Get system number and variable numbers
   const unsigned int from_sys_number = from_sys->number();
   const unsigned int from_var_number = from_var.number();
 
   const unsigned int to_sys_number = to_sys->number();
   const unsigned int to_var_number = to_var.number();
 
   // Get a constant reference to variables, get their number of components
   const unsigned int from_n_comp = from_var.n_components(from_sys->get_mesh());
   const unsigned int to_n_comp = to_var.n_components(to_mesh);
 
   // Sanity check that the variables have the same number of components
   libmesh_assert_equal_to(from_n_comp, to_n_comp);
 
   // Construct map from "from" dofs to "to" dofs.
   typedef std::map<numeric_index_type, numeric_index_type> DofMapping;
   DofMapping dof_mapping;
 
   // Loop through all boundary elements.
   for (const auto & to_elem : to_mesh.active_local_element_ptr_range())
     {
       const Elem * from_elem = to_elem->interior_parent();
 
       libmesh_error_msg_if(!from_elem, "Error, transfer must be between a Mesh and its associated BoundaryMesh.");
 
       // loop through all nodes in each boundary element.
       for (const Node & to_node : to_elem->node_ref_range())
         {
           // Nodes in interior_parent element.
           for (const Node & from_node : from_elem->node_ref_range())
             {
               const dof_id_type from_dof = from_node.dof_number(from_sys_number,
                                                                 from_var_number,
                                                                 from_n_comp - 1);
 
               // If we've already mapped this dof, we don't need to map
               // it again or do floating point comparisons.
               if (!dof_mapping.count(from_dof) &&
                   from_node.absolute_fuzzy_equals(to_node, TOLERANCE))
                 {
                   // Global dof_index for node in BoundaryMesh.
                   const dof_id_type to_dof = to_node.dof_number(to_sys_number,
                                                                 to_var_number,
                                                                 to_n_comp - 1);
 
                   // Keep track of the volume system dof index which is needed.
                   dof_mapping[from_dof] = to_dof;
                 }
             }
         }
     }
 
   // Construct a vector of the indices needed from the Volume system's
   // global solution vector on this processor.
   std::vector<numeric_index_type> needed_indices;
   needed_indices.reserve(dof_mapping.size());
   {
     DofMapping::iterator
       it = dof_mapping.begin(),
       end = dof_mapping.end();
 
     for (; it!=end; ++it)
       needed_indices.push_back(it->first);
   }
 
   // Communicate just the required values without making a copy of the
   // global solution vector.
   std::vector<Number> needed_values;
   from_sys->solution->localize(needed_values, needed_indices);
 
   // Loop over DofMapping again, assigning values in the
   // Boundary System solution vector.
   {
     DofMapping::iterator
       it = dof_mapping.begin(),
       end = dof_mapping.end();
 
     for (unsigned idx=0; it!=end; ++it, ++idx)
       to_sys->solution->set(it->second, needed_values[idx]);
   }
 }

Member Data Documentation

◆ _communicator

const Parallel::Communicator& libMesh::ParallelObject::_communicator

protectedinherited

Definition at line 120 of file parallel_object.h.

Referenced by libMesh::EquationSystems::build_parallel_elemental_solution_vector(), libMesh::EquationSystems::build_parallel_solution_vector(), libMesh::StaticCondensation::close(), libMesh::ParallelObject::comm(), libMesh::CondensedEigenSystem::initialize_condensed_matrices(), libMesh::ParallelObject::n_processors(), libMesh::ParallelObject::operator=(), libMesh::ParallelObject::processor_id(), libMesh::BoundaryInfo::regenerate_id_sets(), and libMesh::StaticCondensationDofMap::reinit().

The documentation for this class was generated from the following files:

include/solution_transfer/boundary_volume_solution_transfer.h
src/solution_transfer/boundary_volume_solution_transfer.C

Public Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ BoundaryVolumeSolutionTransfer()

◆ ~BoundaryVolumeSolutionTransfer()

Member Function Documentation

◆ comm()

◆ n_processors()

◆ processor_id()

◆ transfer()

◆ transfer_boundary_volume()

◆ transfer_volume_boundary()

Member Data Documentation

◆ _communicator