Enumerations
enum	{ New, Exists, Not_found }
	Enumerators for element status in location procedure. More...

Functions
void	send_and_receive_missing_zetas (Problem *problem_pt)
	Helper function to send any "missing" zeta coordinates to the next process and receive any coordinates from previous process. More...

void	send_and_receive_located_info (int &iproc, Mesh const &external_mesh_pt, Problem problem_pt)
	Helper function to send back any located information. More...

void	add_external_haloed_node_to_storage (int &iproc, Node nod_pt, Problem problem_pt, Mesh *const &external_mesh_pt, int &n_cont_inter_values)
	Helper function to add external haloed nodes, including any masters. More...

void	recursively_add_masters_of_external_haloed_node (int &iproc, Node nod_pt, Problem problem_pt, Mesh *const &external_mesh_pt, int &n_cont_inter_values)
	Recursively add any master nodes (and their master nodes etc) of external haloed nodes. More...

void	add_external_haloed_node_helper (int &iproc, Node nod_pt, Problem problem_pt, Mesh *const &external_mesh_pt, int &n_cont_inter_values)
	Helper to add external haloed node that is not a master. More...

void	add_external_haloed_master_node_helper (int &iproc, Node master_nod_pt, Problem problem_pt, Mesh *const &external_mesh_pt, int &n_cont_inter_values)
	Helper function to add external haloed node that is a master. More...

void	get_required_nodal_information_helper (int &iproc, Node nod_pt, Problem problem_pt, Mesh *const &external_mesh_pt, int &n_cont_inter_values)
	Helper function to get the required nodal information from an external haloed node so that a fully-functional external halo node (and therefore element) can be created on the receiving process. More...

void	get_required_master_nodal_information_helper (int &iproc, Node master_nod_pt, Problem problem_pt, Mesh *const &external_mesh_pt, int &n_cont_inter_values)
	Helper function to get the required master nodal information from an external haloed master node so that a fully-functional external halo master node (and possible element) can be created on the receiving proc. More...

void	add_external_halo_node_helper (Node &new_nod_pt, Mesh const &external_mesh_pt, unsigned &loc_p, unsigned &node_index, FiniteElement const &new_el_pt, int &n_cont_inter_values, Problem problem_pt)
	Helper functiono to add external halo node that is not a master. More...

void	construct_new_external_halo_node_helper (Node &new_nod_pt, unsigned &loc_p, unsigned &node_index, FiniteElement const &new_el_pt, Mesh const &external_mesh_pt, Problem problem_pt)
	Helper function which constructs a new external halo node (on an element) with the information sent from the haloed process. More...

void	locate_zeta_for_missing_coordinates (int &iproc, Mesh const &external_mesh_pt, Problem problem_pt, Vector< MeshAsGeomObject *> &mesh_geom_obj_pt)
	Locate zeta for current set of missing coordinates; vector-based version. More...

void	locate_zeta_for_local_coordinates (const Vector< Mesh * > &mesh_pt, Mesh const &external_mesh_pt, Vector< MeshAsGeomObject > &mesh_geom_obj_pt, const unsigned &interaction_index)
	Helper function to locate "local" zeta coordinates This is the vector-based version which operates simultaenously on the meshes contained in the Vectors. More...

void	get_dim_helper (Problem problem_pt, Mesh const &mesh_pt, Mesh *const &external_mesh_pt)
	Helper function that computes the dimension of the elements within each of the specified meshes (and checks they are the same) Stores result in Dim. More...

void	clean_up ()
	Helper function that clears all the intermediate information used during the external storage creation at the end of the procedure. More...

bool	first_closer_than_second (const std::pair< FiniteElement , Vector< double > > &p1, const std::pair< FiniteElement , Vector< double > > &p2)

template<class BULK_ELEMENT , unsigned DIM>
void	setup_bulk_elements_adjacent_to_face_mesh (Problem problem_pt, Vector< unsigned > &boundary_in_bulk_mesh, Mesh const &bulk_mesh_pt, Vector< Mesh *> &face_mesh_pt, const unsigned &interaction=0)
	Identify the `FaceElements` (stored in the mesh pointed to by `face_mesh_pt`) that are adjacent to the bulk elements next to the `boundary_in_bulk_mesh` -th boundary of the mesh pointed to by `bulk_mesh_pt`. The `FaceElements` must be derived from the `ElementWithExternalElement` base class and the adjacent bulk elements are stored as their external elements. More...

template<class BULK_ELEMENT , unsigned DIM>
void	setup_bulk_elements_adjacent_to_face_mesh (Problem problem_pt, const unsigned &boundary_in_bulk_mesh, Mesh const &bulk_mesh_pt, Mesh *const &face_mesh_pt, const unsigned &interaction=0)
	Identify the `FaceElements` (stored in the mesh pointed to by `face_mesh_pt`) that are adjacent to the bulk elements next to the `boundary_in_bulk_mesh` -th boundary of the mesh pointed to by `bulk_mesh_pt`. The `FaceElements` must be derived from the `ElementWithExternalElement` base class and the adjacent bulk elements are stored as their external elements. More...

template<class ELEMENT_0 , class ELEMENT_1 >
void	setup_multi_domain_interactions (Problem problem_pt, Mesh const &first_mesh_pt, Mesh *const &second_mesh_pt, const unsigned &first_interaction=0, const unsigned &second_interaction=0)
	Set up the two-way multi-domain interactions for the problem pointed to by `problem_pt`. Use this for cases where first_mesh_pt and second_mesh_pt occupy the same physical space and are populated by ELEMENT_0 and ELEMENT_1 respectively, and are combined to solve a single problem. The elements in two meshes interact both ways the elements in each mesh act as "external elements" for the elements in the "other" mesh. The interaction indices allow the specification of which interaction we're setting up in the two meshes. They default to zero, which is appropriate if there's only a single interaction. More...

template<class EXT_ELEMENT >
void	setup_multi_domain_interaction (Problem problem_pt, Mesh const &mesh_pt, Mesh *const &external_mesh_pt, const unsigned &interaction_index=0)
	Function to set up the one-way multi-domain interaction for problems where the meshes pointed to by `mesh_pt` and `external_mesh_pt` occupy the same physical space, and the elements in `external_mesh_pt` act as "external elements" for the `ElementWithExternalElements` in `mesh_pt` (but not vice versa): More...

template<class EXT_ELEMENT , class FACE_ELEMENT_GEOM_OBJECT >
void	setup_multi_domain_interaction (Problem problem_pt, Mesh const &mesh_pt, Mesh const &external_mesh_pt, Mesh const &external_face_mesh_pt, const unsigned &interaction_index=0)
	Function to set up the one-way multi-domain interaction for FSI-like problems. More...

template<class EXT_ELEMENT , class FACE_ELEMENT_GEOM_OBJECT >
void	setup_multi_domain_interaction (Problem problem_pt, const Vector< Mesh > &mesh_pt, Mesh const &external_mesh_pt, const Vector< Mesh > &external_face_mesh_pt, const unsigned &interaction_index=0)
	Function to set up the one-way multi-domain interaction for FSI-like problems. More...

template<class EXT_ELEMENT , class GEOM_OBJECT >
void	aux_setup_multi_domain_interaction (Problem problem_pt, Mesh const &mesh_pt, Mesh const &external_mesh_pt, const unsigned &interaction_index, Mesh const &external_face_mesh_pt=0)
	Auxiliary helper function. More...

template<class EXT_ELEMENT , class GEOM_OBJECT >
void	aux_setup_multi_domain_interaction (Problem problem_pt, const Vector< Mesh > &mesh_pt, Mesh const &external_mesh_pt, const unsigned &interaction_index, const Vector< Mesh > &external_face_mesh_pt)
	Auxiliary helper function. More...

template<class EXT_ELEMENT >
void	create_external_halo_elements (int &iproc, const Vector< Mesh > &mesh_pt, Mesh const &external_mesh_pt, Problem *problem_pt, const unsigned &interaction_index)
	Create external (halo) elements on the loop process based on the information received from each locate_zeta call on other processes This is the vector-based function which operates simultaneously on the meshes contained in the vectors. More...

template<class EXT_ELEMENT >
void	add_external_halo_node_to_storage (Node &new_nod_pt, Mesh const &external_mesh_pt, unsigned &loc_p, unsigned &node_index, FiniteElement const &new_el_pt, int &n_cont_inter_values, Problem problem_pt)
	Helper function to add external halo nodes, including any masters, based on information received from the haloed process. More...

template<class EXT_ELEMENT >
void	recursively_add_masters_of_external_halo_node_to_storage (Node &new_nod_pt, Mesh const &external_mesh_pt, unsigned &loc_p, unsigned &node_index, FiniteElement const &new_el_pt, int &n_cont_inter_values, Problem problem_pt)
	Recursively add masters of external halo nodes (and their masters, etc) based on information received from the haloed process. More...

template<class EXT_ELEMENT >
void	add_external_halo_master_node_helper (Node &new_master_nod_pt, Node &new_nod_pt, Mesh const &external_mesh_pt, unsigned &loc_p, int &n_cont_inter_values, Problem problem_pt)
	Helper function to add external halo node that is a master. More...

template<class EXT_ELEMENT >
void	construct_new_external_halo_master_node_helper (Node &new_master_nod_pt, Node &nod_pt, unsigned &loc_p, Mesh const &external_mesh_pt, Problem problem_pt)
	Helper function which constructs a new external halo master node with the information sent from the haloed process. More...

Variables
std::ofstream	Doc_boundary_coordinate_file
	Output file to document the boundary coordinate along the mesh boundary of the bulk mesh during call to setup_bulk_elements_adjacent_to_face_mesh(...) More...

bool	Accept_failed_locate_zeta_in_setup_multi_domain_interaction =false
	Boolean to indicate that failure in setup multi domain functions is acceptable; defaults to false. If set to true external element pointers are set to null for those elements for which external elements couldn't be located. More...

unsigned	Dim
	Dimension of zeta tuples (set by get_dim_helper) – needed because we store the scalar coordinates in flat-packed form. More...

Vector< Vector< unsigned > >	External_element_located
	Lookup scheme for whether a local element's integration point has had an external element assigned to it – essentially boolean. External_element_located[e][ipt] = {0,1} if external element for ipt-th integration in local element e {has not, has} been found. Used locally to ensure that we're not searching for the same elements over and over again when we go around the spirals. More...

Vector< double >	Flat_packed_zetas_not_found_locally
	Vector of flat-packed zeta coordinates for which the external element could not be found during current local search. These will be sent to the next processor in the ring-like parallel search. The zeta coordinates come in groups of Dim (scalar) coordinates. More...

Vector< double >	Received_flat_packed_zetas_to_be_found
	Vector of flat-packed zeta coordinates for which the external element could not be found on another processor and for which we're currently searching here. Whatever can't be found here, gets written into Flat_packed_zetas_not_found_locally and then passed on to the next processor during the ring-like parallel search. The zeta coordinates come in groups of Dim (scalar) coordinates. More...

Vector< int >	Proc_id_plus_one_of_external_element
	Proc_id_plus_one_of_external_element[i] contains the processor id (plus one) of the processor on which the i-th zeta coordinate tuple received from elsewhere (in the order in which these are stored in Received_flat_packed_zetas_to_be_found) was located; it's zero if it wasn't found during the current stage of the ring-like parallel search. More...

Vector< unsigned >	Located_element_status
	Vector to indicate (to another processor) whether a located element (that will have to represented as an external halo element on that processor) should be newly created on that processor (2), already exists on that processor (1), or is not on the current processor either (0). More...

Vector< double >	Flat_packed_located_coordinates
	Vector of flat-packed local coordinates for zeta tuples that have been located. More...

Vector< double >	Flat_packed_doubles
	Vector of flat-packed doubles to be communicated with other processors. More...

unsigned	Counter_for_flat_packed_doubles
	Counter used when processing vector of flat-packed doubles – this is really "private" data, declared here to avoid having to pass it (and the associated array) between the various helper functions. More...

Vector< unsigned >	Flat_packed_unsigneds
	Vector of flat-packed unsigneds to be communicated with other processors – this is really "private" data, declared here to avoid having to pass the array between the various helper functions. More...

Vector< std::string >	Flat_packed_unsigneds_string

unsigned	Counter_for_flat_packed_unsigneds
	Counter used when processing vector of flat-packed unsigneds – this is really "private" data, declared here to avoid having to pass it (and the associated array) between the various helper functions. More...

bool	Use_bulk_element_as_external =false
	Boolean to indicate when to use the bulk element as the external element. Defaults to false, you must have set up FaceElements properly first in order for it to work. More...

bool	Allow_use_of_halo_elements_as_external_elements =true
	Boolean to indicate if we're allowed to use halo elements as external elements. Can drastically reduce the number of external halo elements – currently not aware of any problems therefore set to true by default but retention of this flag allows easy return to previous implementation. More...

bool	Allow_use_of_halo_elements_as_external_elements_for_projection =true
	Indicate whether we are allowed to use halo elements as external elements for projection, possibly only required in parallel unstructured mesh generation during the projection stage. Default set to true. More...

bool	Doc_timings =false
	Boolean to indicate whether to doc timings or not. More...

bool	Doc_stats =false
	Boolean to indicate whether to output basic info during setup_multi_domain_interaction() routines. More...

bool	Doc_full_stats =false
	Boolean to indicate whether to output further info during setup_multi_domain_interaction() routines. More...

Vector< double >	Zeta_coords_for_further_away_comparison

Enumeration Type Documentation

◆ anonymous enum

anonymous enum

Enumerators for element status in location procedure.

Enumerator
New
Exists
Not_found

Definition at line 153 of file multi_domain.h.

Function Documentation

◆ add_external_halo_master_node_helper()

template<class EXT_ELEMENT >

void oomph::Multi_domain_functions::add_external_halo_master_node_helper	(	Node *&	new_master_nod_pt,
		Node *&	new_nod_pt,
		Mesh *const &	external_mesh_pt,
		unsigned &	loc_p,
		int &	n_cont_inter_values,
		Problem *	problem_pt
	)

Helper function to add external halo node that is a master.

Definition at line 2011 of file multi_domain.template.cc.

References construct_new_external_halo_master_node_helper(), Counter_for_flat_packed_unsigneds, oomph::Mesh::external_halo_node_pt(), Flat_packed_unsigneds, and oomph::oomph_info.

Referenced by recursively_add_masters_of_external_halo_node_to_storage().

◆ add_external_halo_node_helper()

void oomph::Multi_domain_functions::add_external_halo_node_helper	(	Node *&	new_nod_pt,
		Mesh *const &	external_mesh_pt,
		unsigned &	loc_p,
		unsigned &	node_index,
		FiniteElement *const &	new_el_pt,
		int &	n_cont_inter_values,
		Problem *	problem_pt
	)

Helper functiono to add external halo node that is not a master.

Helper function to add external halo node that is not a master.

Definition at line 1286 of file multi_domain.cc.

References construct_new_external_halo_node_helper(), Counter_for_flat_packed_unsigneds, oomph::Mesh::external_halo_node_pt(), oomph::FiniteElement::node_pt(), and oomph::oomph_info.

Referenced by get_required_master_nodal_information_helper().

◆ add_external_halo_node_to_storage()

template<class EXT_ELEMENT >

void oomph::Multi_domain_functions::add_external_halo_node_to_storage	(	Node *&	new_nod_pt,
		Mesh *const &	external_mesh_pt,
		unsigned &	loc_p,
		unsigned &	node_index,
		FiniteElement *const &	new_el_pt,
		int &	n_cont_inter_values,
		Problem *	problem_pt
	)

Helper function to add external halo nodes, including any masters, based on information received from the haloed process.

Helper function to add external halo nodes, including any masters, based on information received from the haloed process

Definition at line 1925 of file multi_domain.template.cc.

References oomph::Missing_masters_functions::add_external_halo_node_helper(), and recursively_add_masters_of_external_halo_node_to_storage().

Referenced by create_external_halo_elements().

◆ add_external_haloed_master_node_helper()

void oomph::Multi_domain_functions::add_external_haloed_master_node_helper	(	int &	iproc,
		Node *	master_nod_pt,
		Problem *	problem_pt,
		Mesh *const &	external_mesh_pt,
		int &	n_cont_inter_values
	)

Helper function to add external haloed node that is a master.

Definition at line 590 of file multi_domain.cc.

References oomph::Mesh::add_external_haloed_node_pt(), get_required_master_nodal_information_helper(), and oomph::Mesh::nexternal_haloed_node().

Referenced by recursively_add_masters_of_external_haloed_node().

◆ add_external_haloed_node_helper()

void oomph::Multi_domain_functions::add_external_haloed_node_helper	(	int &	iproc,
		Node *	nod_pt,
		Problem *	problem_pt,
		Mesh *const &	external_mesh_pt,
		int &	n_cont_inter_values
	)

Helper to add external haloed node that is not a master.

Helper function to add external haloed node that is not a master.

Definition at line 532 of file multi_domain.cc.

References oomph::Mesh::add_external_haloed_node_pt(), get_required_nodal_information_helper(), and oomph::Mesh::nexternal_haloed_node().

Referenced by add_external_haloed_node_to_storage().

◆ add_external_haloed_node_to_storage()

void oomph::Multi_domain_functions::add_external_haloed_node_to_storage	(	int &	iproc,
		Node *	nod_pt,
		Problem *	problem_pt,
		Mesh *const &	external_mesh_pt,
		int &	n_cont_inter_values
	)

Helper function to add external haloed nodes, including any masters.

Helper function to add external haloed nodes, inc. masters of external haloed nodes.

Definition at line 446 of file multi_domain.cc.

References add_external_haloed_node_helper(), and recursively_add_masters_of_external_haloed_node().

Referenced by get_required_master_nodal_information_helper(), and locate_zeta_for_missing_coordinates().

◆ aux_setup_multi_domain_interaction() [1/2]

template<class EXT_ELEMENT , class GEOM_OBJECT >

void oomph::Multi_domain_functions::aux_setup_multi_domain_interaction	(	Problem *	problem_pt,
		Mesh *const &	mesh_pt,
		Mesh *const &	external_mesh_pt,
		const unsigned &	interaction_index,
		Mesh *const &	external_face_mesh_pt = `0`
	)

Auxiliary helper function.

This routine calls the locate_zeta routine (simultaneously on each processor for each individual processor's element set if necessary) and sets up the external (halo) element and node storage as necessary. The locate_zeta procedure here works for all multi-domain problems where either two meshes occupy the same physical space but have differing element types (e.g. a Boussinesq convection problem where AdvectionDiffusion elements interact with Navier-Stokes type elements) or two meshes interact along some boundary of the external mesh, represented by a "face mesh", such as an FSI problem.

Definition at line 496 of file multi_domain.template.cc.

Referenced by setup_multi_domain_interaction().

◆ aux_setup_multi_domain_interaction() [2/2]

template<class EXT_ELEMENT , class GEOM_OBJECT >

void oomph::Multi_domain_functions::aux_setup_multi_domain_interaction	(	Problem *	problem_pt,
		const Vector< Mesh *> &	mesh_pt,
		Mesh *const &	external_mesh_pt,
		const unsigned &	interaction_index,
		const Vector< Mesh *> &	external_face_mesh_pt
	)

Auxiliary helper function.

This routine calls the locate_zeta routine (simultaneously on each processor for each individual processor's element set if necessary) and sets up the external (halo) element and node storage as necessary. The locate_zeta procedure here works for all multi-domain problems where either two meshes occupy the same physical space but have differing element types (e.g. a Boussinesq convection problem where AdvectionDiffusion elements interact with Navier-Stokes type elements) or two meshes interact along some boundary of the external mesh, represented by a "face mesh", such as an FSI problem.

Vector-based version operates simultaneously on the meshes contained in the vectors.

If it's is now zero then break out of the spirals loop

Definition at line 533 of file multi_domain.template.cc.

References Accept_failed_locate_zeta_in_setup_multi_domain_interaction, clean_up(), oomph::Problem::communicator_pt(), create_external_halo_elements(), NonRefineableBinArray::current_max_spiral_level(), NonRefineableBinArray::current_min_spiral_level(), oomph::SamplePointContainerParameters::Default_nsample_points_generated_per_element, Dim, oomph::FiniteElement::dim(), oomph::Missing_masters_functions::Doc_full_stats, oomph::Missing_masters_functions::Doc_stats, oomph::Missing_masters_functions::Doc_timings, e, oomph::Mesh::element_pt(), External_element_located, oomph::ElementWithExternalElement::external_element_pt(), oomph::Mesh::finite_element_pt(), Flat_packed_zetas_not_found_locally, i, oomph::ElementWithExternalElement::initialise_external_element_storage(), oomph::FiniteElement::integral_pt(), oomph::FiniteElement::interpolated_x(), oomph::GeneralisedElement::is_halo(), oomph::Integral::knot(), locate_zeta_for_local_coordinates(), locate_zeta_for_missing_coordinates(), BinArray::max_bin_dimension(), NonRefineableBinArray::n_spiral_chunk(), oomph::Node::ndim(), oomph::Mesh::nelement(), oomph::Mesh::nexternal_halo_element(), oomph::Mesh::nexternal_halo_node(), oomph::Mesh::nexternal_haloed_element(), oomph::Mesh::nexternal_haloed_node(), oomph::Mesh::nhalo_node(), oomph::Mesh::nhaloed_node(), oomph::Mesh::nnode(), oomph::FiniteElement::node_pt(), oomph::Mesh::nroot_halo_element(), oomph::Mesh::nroot_haloed_element(), oomph::Integral::nweight(), oomph::oomph_info, oomph::FiniteElement::output(), BinArray::output_bins(), oomph::Problem::problem_has_been_distributed(), s, send_and_receive_located_info(), send_and_receive_missing_zetas(), t, oomph::TimingHelpers::timer(), Use_bulk_element_as_external, oomph::UseCGALSamplePointContainer, oomph::UseNonRefineableBinArray, and oomph::UseRefineableBinArray.

◆ clean_up()

void oomph::Multi_domain_functions::clean_up ( )

Helper function that clears all the intermediate information used during the external storage creation at the end of the procedure.

Helper function that clears all the information used during the external storage creation

Definition at line 2276 of file multi_domain.cc.

Referenced by aux_setup_multi_domain_interaction().

◆ construct_new_external_halo_master_node_helper()

template<class EXT_ELEMENT >

void oomph::Multi_domain_functions::construct_new_external_halo_master_node_helper	(	Node *&	new_master_nod_pt,
		Node *&	nod_pt,
		unsigned &	loc_p,
		Mesh *const &	external_mesh_pt,
		Problem *	problem_pt
	)

Helper function which constructs a new external halo master node with the information sent from the haloed process.

Helper function which constructs a new external halo master node with the required information sent from the haloed process

The first entry of All_unsigned_values is the default node update id

... so for the specified update_id, call add_node_update_info

Now call update_node_update

Definition at line 2050 of file multi_domain.template.cc.

References oomph::Mesh::add_boundary_node(), oomph::Mesh::add_external_halo_element_pt(), oomph::Mesh::add_external_halo_node_pt(), oomph::AlgebraicNode::add_node_update_info(), Counter_for_flat_packed_doubles, Counter_for_flat_packed_unsigneds, oomph::FiniteElement::dim(), oomph::Mesh::external_halo_element_pt(), Flat_packed_doubles, Flat_packed_unsigneds, oomph::AlgebraicMesh::geom_object_list_pt(), oomph::MacroElementNodeUpdateMesh::geom_object_vector_pt(), i, oomph::BoundaryNodeBase::index_of_first_value_assigned_by_face_element_pt(), oomph::FiniteElement::local_coordinate_of_node(), oomph::MacroElementNodeUpdateMesh::macro_domain_pt(), oomph::Domain::macro_element_pt(), oomph::Node::ndim(), oomph::FiniteElement::nnode(), oomph::FiniteElement::node_pt(), oomph::MacroElementNodeUpdateNode::node_update_element_pt(), oomph::TimeStepper::ntstorage(), oomph::Data::nvalue(), oomph::oomph_info, oomph::QElementBase::s_macro_ll(), oomph::QElementBase::s_macro_ur(), oomph::FiniteElement::set_macro_elem_pt(), oomph::MacroElementNodeUpdateNode::set_node_update_info(), oomph::Data::set_value(), t, oomph::Problem::time_stepper_pt(), oomph::AlgebraicMesh::update_node_update(), oomph::SolidNode::variable_position_pt(), and oomph::Node::x().

Referenced by add_external_halo_master_node_helper().

◆ construct_new_external_halo_node_helper()

void oomph::Multi_domain_functions::construct_new_external_halo_node_helper	(	Node *&	new_nod_pt,
		unsigned &	loc_p,
		unsigned &	node_index,
		FiniteElement *const &	new_el_pt,
		Mesh *const &	external_mesh_pt,
		Problem *	problem_pt
	)

Helper function which constructs a new external halo node (on an element) with the information sent from the haloed process.

Helper function which constructs a new external halo node (on new element) with the required information sent from the haloed process

The first entry of All_alg_nodal_info contains the default node update id e.g. for the quarter circle there are "Upper_left_box", "Lower right box" etc...

again we need the size of this vector as it varies between meshes; we also need some indication as to which geometric object should be used...

For the received update_id, ref_value, geom_object call add_node_update_info

Now call update_node_update

Definition at line 1333 of file multi_domain.cc.

References oomph::Mesh::add_boundary_node(), oomph::Mesh::add_external_halo_node_pt(), oomph::AlgebraicNode::add_node_update_info(), oomph::FiniteElement::construct_boundary_node(), oomph::FiniteElement::construct_node(), Counter_for_flat_packed_unsigneds, oomph::MacroElementNodeUpdateMesh::geom_object_vector_pt(), i, oomph::BoundaryNodeBase::index_of_first_value_assigned_by_face_element_pt(), oomph::FiniteElement::local_coordinate_of_node(), locate_zeta_for_missing_coordinates(), oomph::Node::ndim(), oomph::TimeStepper::ntstorage(), oomph::Data::nvalue(), oomph::oomph_info, oomph::Node::resize(), oomph::MacroElementNodeUpdateNode::set_node_update_info(), oomph::Data::set_value(), t, oomph::Problem::time_stepper_pt(), oomph::AlgebraicMesh::update_node_update(), oomph::SolidNode::variable_position_pt(), and oomph::Node::x().

Referenced by add_external_halo_node_helper().

◆ create_external_halo_elements()

template<class EXT_ELEMENT >

void oomph::Multi_domain_functions::create_external_halo_elements	(	int &	iproc,
		const Vector< Mesh *> &	mesh_pt,
		Mesh *const &	external_mesh_pt,
		Problem *	problem_pt,
		const unsigned &	interaction_index
	)

Create external (halo) elements on the loop process based on the information received from each locate_zeta call on other processes This is the vector-based function which operates simultaneously on the meshes contained in the vectors.

Creates external (halo) elements on the loop process based on the information received from each locate_zeta call on other processes. vector based version

Definition at line 1683 of file multi_domain.template.cc.

References add_external_halo_node_to_storage(), oomph::Problem::communicator_pt(), Counter_for_flat_packed_doubles, Counter_for_flat_packed_unsigneds, oomph::FiniteElement::dim(), e, External_element_located, oomph::ElementWithExternalElement::external_element_pt(), Flat_packed_doubles, Flat_packed_located_coordinates, Flat_packed_unsigneds, i, oomph::FiniteElement::integral_pt(), oomph::GeneralisedElement::is_halo(), Located_element_status, oomph::MacroElementNodeUpdateMesh::macro_domain_pt(), New, oomph::FiniteElement::nnode(), oomph::Integral::nweight(), oomph::oomph_info, Proc_id_plus_one_of_external_element, oomph::QElementBase::s_macro_ll(), oomph::QElementBase::s_macro_ur(), and oomph::FiniteElement::set_macro_elem_pt().

Referenced by aux_setup_multi_domain_interaction().

◆ first_closer_than_second()

bool oomph::Multi_domain_functions::first_closer_than_second	(	const std::pair< FiniteElement *, Vector< double > > &	p1,
		const std::pair< FiniteElement *, Vector< double > > &	p2
	)

Comparison function for sorting entries in bin: Returns true if point identified by p1 (comprising pointer to finite element and vector of local coordinates within that element) is closer to Zeta_coords_for_further_away_comparison than p2

Definition at line 2296 of file multi_domain.cc.

References Dim, i, oomph::FiniteElement::interpolated_zeta(), and s.

◆ get_dim_helper()

void oomph::Multi_domain_functions::get_dim_helper	(	Problem *	problem_pt,
		Mesh *const &	mesh_pt,
		Mesh *const &	external_mesh_pt
	)

Helper function that computes the dimension of the elements within each of the specified meshes (and checks they are the same) Stores result in Dim.

Helper function that computes the dimension of the elements within each of the specified meshes (and checks they are the same). Stores result in Dim.

Definition at line 2212 of file multi_domain.cc.

References oomph::Problem::communicator_pt(), oomph::Mesh::element_pt(), and oomph::Mesh::nelement().

Referenced by setup_multi_domain_interaction().

◆ get_required_master_nodal_information_helper()

void oomph::Multi_domain_functions::get_required_master_nodal_information_helper	(	int &	iproc,
		Node *	master_nod_pt,
		Problem *	problem_pt,
		Mesh *const &	external_mesh_pt,
		int &	n_cont_inter_values
	)

Helper function to get the required master nodal information from an external haloed master node so that a fully-functional external halo master node (and possible element) can be created on the receiving proc.

Helper function to get the required master nodal information from an external haloed master node so that a fully-functional external halo master node (and possible element) can be created on the receiving process

Definition at line 899 of file multi_domain.cc.

References add_external_halo_node_helper(), add_external_haloed_node_to_storage(), oomph::FiniteElement::dim(), oomph::AlgebraicMesh::geom_object_list_pt(), oomph::AlgebraicNode::geom_object_pt(), i, oomph::BoundaryNodeBase::index_of_first_value_assigned_by_face_element_pt(), oomph::Node::is_on_boundary(), oomph::FiniteElement::macro_elem_pt(), oomph::MacroElement::macro_element_number(), oomph::Mesh::nboundary(), oomph::Node::ndim(), oomph::AlgebraicNode::ngeom_object(), oomph::AlgebraicMesh::ngeom_object_list_pt(), oomph::SolidNode::nlagrangian(), oomph::SolidNode::nlagrangian_type(), oomph::FiniteElement::nnode(), oomph::FiniteElement::node_pt(), oomph::MacroElementNodeUpdateNode::node_update_element_pt(), oomph::AlgebraicNode::node_update_fct_id(), oomph::Node::nposition_type(), oomph::AlgebraicNode::nref_value(), oomph::Problem::ntime_stepper(), oomph::TimeStepper::ntstorage(), oomph::Data::nvalue(), oomph::AlgebraicNode::ref_value(), oomph::QElementBase::s_macro_ll(), oomph::QElementBase::s_macro_ur(), t, oomph::Data::time_stepper_pt(), oomph::Problem::time_stepper_pt(), oomph::Node::value(), oomph::SolidNode::variable_position_pt(), and oomph::Node::x().

Referenced by add_external_haloed_master_node_helper(), and get_required_nodal_information_helper().

◆ get_required_nodal_information_helper()

void oomph::Multi_domain_functions::get_required_nodal_information_helper	(	int &	iproc,
		Node *	nod_pt,
		Problem *	problem_pt,
		Mesh *const &	external_mesh_pt,
		int &	n_cont_inter_values
	)

Helper function to get the required nodal information from an external haloed node so that a fully-functional external halo node (and therefore element) can be created on the receiving process.

Helper function to get the required nodal information from an external haloed node so that a fully-functional external halo node (and therefore element) can be created on the receiving process

Definition at line 643 of file multi_domain.cc.

References oomph::AlgebraicMesh::geom_object_list_pt(), oomph::AlgebraicNode::geom_object_pt(), get_required_master_nodal_information_helper(), i, oomph::BoundaryNodeBase::index_of_first_value_assigned_by_face_element_pt(), oomph::Node::is_on_boundary(), oomph::Mesh::nboundary(), oomph::Node::ndim(), oomph::AlgebraicNode::ngeom_object(), oomph::AlgebraicMesh::ngeom_object_list_pt(), oomph::AlgebraicNode::node_update_fct_id(), oomph::AlgebraicNode::nref_value(), oomph::Problem::ntime_stepper(), oomph::TimeStepper::ntstorage(), oomph::Data::nvalue(), oomph::AlgebraicNode::ref_value(), t, oomph::Data::time_stepper_pt(), oomph::Problem::time_stepper_pt(), oomph::Node::value(), oomph::SolidNode::variable_position_pt(), and oomph::Node::x().

Referenced by add_external_haloed_node_helper().

◆ locate_zeta_for_local_coordinates()

void oomph::Multi_domain_functions::locate_zeta_for_local_coordinates	(	const Vector< Mesh *> &	mesh_pt,
		Mesh *const &	external_mesh_pt,
		Vector< MeshAsGeomObject *> &	mesh_geom_obj_pt,
		const unsigned &	interaction_index
	)

Helper function to locate "local" zeta coordinates This is the vector-based version which operates simultaenously on the meshes contained in the Vectors.

locate zeta for current set of "local" coordinates vector-based version

Definition at line 2022 of file multi_domain.cc.

References oomph::FaceElement::bulk_element_pt(), Dim, oomph::FiniteElement::dim(), e, oomph::ElementWithExternalElement::external_element_local_coord(), oomph::ElementWithExternalElement::external_element_pt(), oomph::FaceElement::get_local_coordinate_in_bulk(), i, oomph::FiniteElement::integral_pt(), oomph::FiniteElement::interpolated_zeta(), oomph::GeneralisedElement::is_halo(), oomph::Integral::knot(), and oomph::Integral::nweight().

Referenced by aux_setup_multi_domain_interaction(), and locate_zeta_for_missing_coordinates().

◆ locate_zeta_for_missing_coordinates()

void oomph::Multi_domain_functions::locate_zeta_for_missing_coordinates	(	int &	iproc,
		Mesh *const &	external_mesh_pt,
		Problem *	problem_pt,
		Vector< MeshAsGeomObject *> &	mesh_geom_obj_pt
	)

Locate zeta for current set of missing coordinates; vector-based version.

Helper function to locate these "missing" zeta coordinates. This is the vector-based function which operates simultaneously on the meshes contained in the vectors.

Definition at line 1677 of file multi_domain.cc.

References oomph::Mesh::add_external_haloed_element_pt(), add_external_haloed_node_to_storage(), oomph::Problem::communicator_pt(), Dim, oomph::FiniteElement::dim(), Exists, oomph::FaceElement::get_local_coordinate_in_bulk(), i, oomph::GeneralisedElement::is_halo(), locate_zeta_for_local_coordinates(), oomph::FiniteElement::macro_elem_pt(), oomph::MacroElement::macro_element_number(), New, oomph::FiniteElement::nnode(), oomph::FiniteElement::node_pt(), Not_found, oomph::QElementBase::s_macro_ll(), and oomph::QElementBase::s_macro_ur().

Referenced by aux_setup_multi_domain_interaction(), and construct_new_external_halo_node_helper().

◆ recursively_add_masters_of_external_halo_node_to_storage()

template<class EXT_ELEMENT >

void oomph::Multi_domain_functions::recursively_add_masters_of_external_halo_node_to_storage	(	Node *&	new_nod_pt,
		Mesh *const &	external_mesh_pt,
		unsigned &	loc_p,
		unsigned &	node_index,
		FiniteElement *const &	new_el_pt,
		int &	n_cont_inter_values,
		Problem *	problem_pt
	)

Recursively add masters of external halo nodes (and their masters, etc) based on information received from the haloed process.

Recursively add masters of external halo nodes (and their masters, etc) based on information received from the haloed process

Definition at line 1951 of file multi_domain.template.cc.

References add_external_halo_master_node_helper(), Counter_for_flat_packed_doubles, Counter_for_flat_packed_unsigneds, Flat_packed_doubles, Flat_packed_unsigneds, oomph::oomph_info, oomph::Node::set_hanging_pt(), and oomph::HangInfo::set_master_node_pt().

Referenced by add_external_halo_node_to_storage().

◆ recursively_add_masters_of_external_haloed_node()

void oomph::Multi_domain_functions::recursively_add_masters_of_external_haloed_node	(	int &	iproc,
		Node *	nod_pt,
		Problem *	problem_pt,
		Mesh *const &	external_mesh_pt,
		int &	n_cont_inter_values
	)

Recursively add any master nodes (and their master nodes etc) of external haloed nodes.

Recursively add any master nodes (and their master nodes etc) of external nodes

Definition at line 468 of file multi_domain.cc.

References add_external_haloed_master_node_helper(), oomph::Node::hanging_pt(), oomph::Node::is_hanging(), oomph::HangInfo::master_node_pt(), oomph::HangInfo::master_weight(), and oomph::HangInfo::nmaster().

Referenced by add_external_haloed_node_to_storage().

◆ send_and_receive_located_info()

void oomph::Multi_domain_functions::send_and_receive_located_info	(	int &	iproc,
		Mesh *const &	external_mesh_pt,
		Problem *	problem_pt
	)

Helper function to send back any located information.

Send location information from current process; Received location information from (current process + iproc) modulo (nproc)

Definition at line 248 of file multi_domain.cc.

References oomph::Problem::communicator_pt(), Dim, i, and oomph::oomph_info.

Referenced by aux_setup_multi_domain_interaction(), and send_and_receive_missing_zetas().

◆ send_and_receive_missing_zetas()

void oomph::Multi_domain_functions::send_and_receive_missing_zetas ( Problem * problem_pt )

Helper function to send any "missing" zeta coordinates to the next process and receive any coordinates from previous process.

Send the zeta coordinates from the current process to the next process; receive from the previous process

Definition at line 184 of file multi_domain.cc.

References oomph::Problem::communicator_pt(), and send_and_receive_located_info().

Referenced by aux_setup_multi_domain_interaction().

◆ setup_bulk_elements_adjacent_to_face_mesh() [1/2]

template<class BULK_ELEMENT , unsigned DIM>

void oomph::Multi_domain_functions::setup_bulk_elements_adjacent_to_face_mesh	(	Problem *	problem_pt,
		Vector< unsigned > &	boundary_in_bulk_mesh,
		Mesh *const &	bulk_mesh_pt,
		Vector< Mesh *> &	face_mesh_pt,
		const unsigned &	interaction = `0`
	)

Identify the FaceElements (stored in the mesh pointed to by face_mesh_pt) that are adjacent to the bulk elements next to the boundary_in_bulk_mesh -th boundary of the mesh pointed to by bulk_mesh_pt. The FaceElements must be derived from the ElementWithExternalElement base class and the adjacent bulk elements are stored as their external elements.

This is the vector-based version which deals with multiple bulk mesh boundaries at the same time.

Identify the FaceElements (stored in the mesh pointed to by face_mesh_pt) that are adjacent to the bulk elements next to the boundary_in_bulk_mesh -th boundary of the mesh pointed to by bulk_mesh_pt. The FaceElements must be derived from the ElementWithExternalElement base class and the adjacent bulk elements are stored as their external elements.

This is the vector-based version which deals with multiple bulk mesh boundaries at the same time.

Definition at line 78 of file multi_domain.template.cc.

References Doc_boundary_coordinate_file, e, oomph::FiniteElement::get_s_plot(), i, oomph::FaceElement::interpolated_x(), oomph::FiniteElement::interpolated_zeta(), oomph::FiniteElement::nplot_points(), s, oomph::FaceElement::set_boundary_number_in_bulk_mesh(), setup_multi_domain_interaction(), oomph::FiniteElement::tecplot_zone_string(), and oomph::FiniteElement::write_tecplot_zone_footer().

Referenced by oomph::FSI_functions::setup_fluid_load_info_for_solid_elements(), and oomph::FSI_functions::setup_solid_elements_for_displacement_bc().

◆ setup_bulk_elements_adjacent_to_face_mesh() [2/2]

template<class BULK_ELEMENT , unsigned DIM>

void oomph::Multi_domain_functions::setup_bulk_elements_adjacent_to_face_mesh	(	Problem *	problem_pt,
		const unsigned &	boundary_in_bulk_mesh,
		Mesh *const &	bulk_mesh_pt,
		Mesh *const &	face_mesh_pt,
		const unsigned &	interaction = `0`
	)

Identify the FaceElements (stored in the mesh pointed to by face_mesh_pt) that are adjacent to the bulk elements next to the boundary_in_bulk_mesh -th boundary of the mesh pointed to by bulk_mesh_pt. The FaceElements must be derived from the ElementWithExternalElement base class and the adjacent bulk elements are stored as their external elements.

Definition at line 216 of file multi_domain.template.cc.

References setup_multi_domain_interactions().

◆ setup_multi_domain_interaction() [1/3]

template<class EXT_ELEMENT >

void oomph::Multi_domain_functions::setup_multi_domain_interaction	(	Problem *	problem_pt,
		Mesh *const &	mesh_pt,
		Mesh *const &	external_mesh_pt,
		const unsigned &	interaction_index = `0`
	)

Function to set up the one-way multi-domain interaction for problems where the meshes pointed to by mesh_pt and external_mesh_pt occupy the same physical space, and the elements in external_mesh_pt act as "external elements" for the ElementWithExternalElements in mesh_pt (but not vice versa):

mesh_pt points to the mesh of ElemenWithExternalElements for which the interaction is set up.
external_mesh_pt points to the mesh that contains the elements of type EXT_ELEMENT that act as "external elements" for the ElementWithExternalElements in \ mesh_pt.
The interaction_index parameter defaults to zero and must be otherwise set by the user if there is more than one mesh that provides sources for the Mesh pointed to by mesh_pt.

Function to set up the one-way multi-domain interaction for problems where the meshes pointed to by mesh_pt and external_mesh_pt occupy the same physical space, and the elements in external_mesh_pt act as "external elements" for the ElementWithExternalElements in mesh_pt (but not vice versa):
mesh_pt points to the mesh of ElemenWithExternalElements for which the interaction is set up.
external_mesh_pt points to the mesh that contains the elements of type EXT_ELEMENT that act as "external elements" for the ElementWithExternalElements in mesh_pt.
The interaction_index parameter defaults to zero and must be otherwise set by the user if there is more than one mesh that provides sources for the Mesh pointed to by mesh_pt.

Definition at line 286 of file multi_domain.template.cc.

References Dim, get_dim_helper(), and Use_bulk_element_as_external.

Referenced by oomph::ProjectionProblem< PROJECTABLE_ELEMENT >::project(), setup_bulk_elements_adjacent_to_face_mesh(), setup_multi_domain_interaction(), and setup_multi_domain_interactions().

◆ setup_multi_domain_interaction() [2/3]

template<class EXT_ELEMENT , class FACE_ELEMENT_GEOM_OBJECT >

void oomph::Multi_domain_functions::setup_multi_domain_interaction	(	Problem *	problem_pt,
		Mesh *const &	mesh_pt,
		Mesh *const &	external_mesh_pt,
		Mesh *const &	external_face_mesh_pt,
		const unsigned &	interaction_index = `0`
	)

Function to set up the one-way multi-domain interaction for FSI-like problems.

mesh_pt points to the mesh of ElemenWithExternalElements for which the interaction is set up. In an FSI example, this mesh would contain the FSIWallElements (either beam/shell elements or the FSISolidTractionElements that apply the traction to a "bulk" solid mesh that is loaded by the fluid.)
external_mesh_pt points to the mesh that contains the elements of type EXT_ELEMENT that provide the "source" for the ElementWithExternalElements. In an FSI example, this mesh would contain the "bulk" fluid elements.
external_face_mesh_pt points to the mesh of FaceElements attached to the external_mesh_pt. The mesh pointed to by external_face_mesh_pt has the same dimension as mesh_pt. The elements contained in external_face_mesh_pt are of type FACE_ELEMENT_GEOM_OBJECT. In an FSI example, these elements are usually the FaceElementAsGeomObjects (templated by the type of the "bulk" fluid elements to which they are attached) that define the FSI boundary of the fluid domain.
The interaction_index parameter defaults to zero and must otherwise be set by the user if there is more than one mesh that provides "external elements" for the Mesh pointed to by mesh_pt (e.g. in the case when a beam or shell structure is loaded by fluid from both sides.)

Function to set up the one-way multi-domain interaction for FSI-like problems.

mesh_pt points to the mesh of ElemenWithExternalElements for which the interaction is set up. In an FSI example, this mesh would contain the FSIWallElements (either beam/shell elements or the FSISolidTractionElements that apply the traction to a "bulk" solid mesh that is loaded by the fluid.)
external_mesh_pt points to the mesh that contains the elements of type EXT_ELEMENT that provide the "source" for the ElementWithExternalElements. In an FSI example, this mesh would contain the "bulk" fluid elements.
external_face_mesh_pt points to the mesh of FaceElements attached to the external_mesh_pt. The mesh pointed to by external_face_mesh_pt has the same dimension as mesh_pt. The elements contained in external_face_mesh_pt are of type FACE_ELEMENT_GEOM_OBJECT. In an FSI example, these elements are usually the FaceElementAsGeomObjects (templated by the type of the "bulk" fluid elements to which they are attached) that define the FSI boundary of the fluid domain.
The interaction_index parameter defaults to zero and must otherwise be set by the user if there is more than one mesh that provides "external elements" for the Mesh pointed to by mesh_pt (e.g. in the case when a beam or shell structure is loaded by fluid from both sides.)

Definition at line 343 of file multi_domain.template.cc.

References aux_setup_multi_domain_interaction(), Dim, get_dim_helper(), setup_multi_domain_interaction(), and Use_bulk_element_as_external.

◆ setup_multi_domain_interaction() [3/3]

template<class EXT_ELEMENT , class FACE_ELEMENT_GEOM_OBJECT >

void oomph::Multi_domain_functions::setup_multi_domain_interaction	(	Problem *	problem_pt,
		const Vector< Mesh *> &	mesh_pt,
		Mesh *const &	external_mesh_pt,
		const Vector< Mesh *> &	external_face_mesh_pt,
		const unsigned &	interaction_index = `0`
	)

Function to set up the one-way multi-domain interaction for FSI-like problems.

mesh_pt points to the mesh of ElemenWithExternalElements for which the interaction is set up. In an FSI example, this mesh would contain the FSIWallElements (either beam/shell elements or the FSISolidTractionElements that apply the traction to a "bulk" solid mesh that is loaded by the fluid.)
external_mesh_pt points to the mesh that contains the elements of type EXT_ELEMENT that provide the "source" for the ElementWithExternalElements. In an FSI example, this mesh would contain the "bulk" fluid elements.
external_face_mesh_pt points to the mesh of FaceElements attached to the external_mesh_pt. The mesh pointed to by external_face_mesh_pt has the same dimension as mesh_pt. The elements contained in external_face_mesh_pt are of type FACE_ELEMENT_GEOM_OBJECT. In an FSI example, these elements are usually the FaceElementAsGeomObjects (templated by the type of the "bulk" fluid elements to which they are attached) that define the FSI boundary of the fluid domain.
The interaction_index parameter defaults to zero and must otherwise be set by the user if there is more than one mesh that provides "external elements" for the Mesh pointed to by mesh_pt (e.g. in the case when a beam or shell structure is loaded by fluid from both sides.)

This is the vector-based version which operates simultaneously on the meshes contained in the Vector arguments.

Function to set up the one-way multi-domain interaction for FSI-like problems.

mesh_pt points to the mesh of ElemenWithExternalElements for which the interaction is set up. In an FSI example, this mesh would contain the FSIWallElements (either beam/shell elements or the FSISolidTractionElements that apply the traction to a "bulk" solid mesh that is loaded by the fluid.)
external_mesh_pt points to the mesh that contains the elements of type EXT_ELEMENT that provide the "source" for the ElementWithExternalElements. In an FSI example, this mesh would contain the "bulk" fluid elements.
external_face_mesh_pt points to the mesh of FaceElements attached to the external_mesh_pt. The mesh pointed to by external_face_mesh_pt has the same dimension as mesh_pt. The elements contained in external_face_mesh_pt are of type FACE_ELEMENT_GEOM_OBJECT. In an FSI example, these elements are usually the FaceElementAsGeomObjects (templated by the type of the "bulk" fluid elements to which they are attached) that define the FSI boundary of the fluid domain.
The interaction_index parameter defaults to zero and must otherwise be set by the user if there is more than one mesh that provides "external elements" for the Mesh pointed to by mesh_pt (e.g. in the case when a beam or shell structure is loaded by fluid from both sides.)

Vector-based version operates simultaneously on the meshes contained in the vectors.

Definition at line 404 of file multi_domain.template.cc.

References aux_setup_multi_domain_interaction(), Dim, get_dim_helper(), i, and Use_bulk_element_as_external.

◆ setup_multi_domain_interactions()

template<class ELEMENT_0 , class ELEMENT_1 >

void oomph::Multi_domain_functions::setup_multi_domain_interactions	(	Problem *	problem_pt,
		Mesh *const &	first_mesh_pt,
		Mesh *const &	second_mesh_pt,
		const unsigned &	first_interaction = `0`,
		const unsigned &	second_interaction = `0`
	)

Set up the two-way multi-domain interactions for the problem pointed to by problem_pt. Use this for cases where first_mesh_pt and second_mesh_pt occupy the same physical space and are populated by ELEMENT_0 and ELEMENT_1 respectively, and are combined to solve a single problem. The elements in two meshes interact both ways the elements in each mesh act as "external elements" for the elements in the "other" mesh. The interaction indices allow the specification of which interaction we're setting up in the two meshes. They default to zero, which is appropriate if there's only a single interaction.

Definition at line 252 of file multi_domain.template.cc.

References oomph::Mesh::delete_all_external_storage(), and setup_multi_domain_interaction().

Referenced by setup_bulk_elements_adjacent_to_face_mesh().

Variable Documentation

◆ Accept_failed_locate_zeta_in_setup_multi_domain_interaction

bool oomph::Multi_domain_functions::Accept_failed_locate_zeta_in_setup_multi_domain_interaction =false

Boolean to indicate that failure in setup multi domain functions is acceptable; defaults to false. If set to true external element pointers are set to null for those elements for which external elements couldn't be located.

Definition at line 62 of file multi_domain.cc.

Referenced by aux_setup_multi_domain_interaction().

◆ Allow_use_of_halo_elements_as_external_elements

bool oomph::Multi_domain_functions::Allow_use_of_halo_elements_as_external_elements =true

Boolean to indicate if we're allowed to use halo elements as external elements. Can drastically reduce the number of external halo elements – currently not aware of any problems therefore set to true by default but retention of this flag allows easy return to previous implementation.

Definition at line 155 of file multi_domain.cc.

◆ Allow_use_of_halo_elements_as_external_elements_for_projection

bool oomph::Multi_domain_functions::Allow_use_of_halo_elements_as_external_elements_for_projection =true

Indicate whether we are allowed to use halo elements as external elements for projection, possibly only required in parallel unstructured mesh generation during the projection stage. Default set to true.

Definition at line 161 of file multi_domain.cc.

◆ Counter_for_flat_packed_doubles

unsigned oomph::Multi_domain_functions::Counter_for_flat_packed_doubles

Counter used when processing vector of flat-packed doubles – this is really "private" data, declared here to avoid having to pass it (and the associated array) between the various helper functions.

Definition at line 118 of file multi_domain.cc.

Referenced by construct_new_external_halo_master_node_helper(), oomph::RefineableTriangleMesh< ELEMENT >::construct_new_halo_node_helper(), oomph::RefineableTriangleMesh< ELEMENT >::construct_new_node_load_balance_helper(), create_external_halo_elements(), oomph::RefineableTriangleMesh< ELEMENT >::load_balance(), recursively_add_masters_of_external_halo_node_to_storage(), oomph::RefineableTriangleMesh< ELEMENT >::reset_halo_haloed_scheme(), and oomph::RefineableTriangleMesh< ELEMENT >::reset_halo_haloed_scheme_helper().

◆ Counter_for_flat_packed_unsigneds

unsigned oomph::Multi_domain_functions::Counter_for_flat_packed_unsigneds

Counter used when processing vector of flat-packed unsigneds – this is really "private" data, declared here to avoid having to pass it (and the associated array) between the various helper functions.

Definition at line 140 of file multi_domain.cc.

◆ Dim

unsigned oomph::Multi_domain_functions::Dim

Dimension of zeta tuples (set by get_dim_helper) – needed because we store the scalar coordinates in flat-packed form.

Definition at line 66 of file multi_domain.cc.

Referenced by aux_setup_multi_domain_interaction(), oomph::HelmholtzBCElementBase< ELEMENT >::compute_contribution_to_fourier_components(), oomph::HelmholtzDtNBoundaryElement< ELEMENT >::compute_gamma_contribution(), oomph::LinearisedAxisymPoroelasticBJS_FSIElement< FLUID_BULK_ELEMENT, POROELASTICITY_BULK_ELEMENT >::contribution_to_total_porous_flux(), oomph::HelmholtzBCElementBase< ELEMENT >::d_shape_and_test_local(), oomph::NetFluxControlElement::dim(), oomph::QSUPGAdvectionDiffusionElement< DIM, NNODE_1D >::dshape_and_dtest_eulerian_at_knot_adv_diff(), oomph::NavierStokesFluxControlElement< ELEMENT >::fill_in_generic_residual_contribution_fluid_traction(), oomph::HelmholtzAbsorbingBCElement< ELEMENT >::fill_in_generic_residual_contribution_helmholtz_abc(), oomph::HelmholtzDtNBoundaryElement< ELEMENT >::fill_in_generic_residual_contribution_helmholtz_DtN_bc(), first_closer_than_second(), oomph::PseudoElasticPreconditioner::lagrange_multiplier_preconditioner_solve(), oomph::PseudoElasticPreconditionerOld::lagrange_multiplier_preconditioner_solve(), locate_zeta_for_local_coordinates(), locate_zeta_for_missing_coordinates(), oomph::NavierStokesFluxControlElement< ELEMENT >::NavierStokesFluxControlElement(), oomph::LinearisedAxisymPoroelasticBJS_FSIElement< FLUID_BULK_ELEMENT, POROELASTICITY_BULK_ELEMENT >::output(), oomph::HelmholtzFluxFromNormalDisplacementBCElement< HELMHOLTZ_BULK_ELEMENT, ELASTICITY_BULK_ELEMENT >::output(), oomph::FourierDecomposedHelmholtzFluxFromNormalDisplacementBCElement< HELMHOLTZ_BULK_ELEMENT, ELASTICITY_BULK_ELEMENT >::output(), oomph::LinearisedFSIAxisymmetricNStNoSlipBCElementElement< FLUID_BULK_ELEMENT, SOLID_BULK_ELEMENT >::output(), oomph::PMLHelmholtzFluxFromNormalDisplacementBCElement< HELMHOLTZ_BULK_ELEMENT, ELASTICITY_BULK_ELEMENT >::output(), oomph::StorableShapeElementBase::pre_compute_d2shape_local_at_knots(), oomph::StorableShapeElementBase::pre_compute_dshape_local_at_knots(), oomph::RefineableNavierStokesFluxControlElement< ELEMENT >::refineable_fill_in_generic_residual_contribution_fluid_traction(), send_and_receive_located_info(), oomph::PseudoElasticPreconditioner::setup(), oomph::PseudoElasticPreconditionerOld::setup(), setup_multi_domain_interaction(), and oomph::PVDEquationsWithPressure< DIM >::solid_pshape_at_knot().

◆ Doc_boundary_coordinate_file

std::ofstream oomph::Multi_domain_functions::Doc_boundary_coordinate_file

Output file to document the boundary coordinate along the mesh boundary of the bulk mesh during call to setup_bulk_elements_adjacent_to_face_mesh(...)

Definition at line 53 of file multi_domain.cc.

Referenced by setup_bulk_elements_adjacent_to_face_mesh().

◆ Doc_full_stats

bool oomph::Multi_domain_functions::Doc_full_stats =false

Boolean to indicate whether to output further info during setup_multi_domain_interaction() routines.

Boolean to indicate whether to document further info (to screen) during setup_multi_domain_interaction() routines.

Definition at line 172 of file multi_domain.cc.

Referenced by oomph::RefineableTriangleMesh< ELEMENT >::adapt().

◆ Doc_stats

bool oomph::Multi_domain_functions::Doc_stats =false

Boolean to indicate whether to output basic info during setup_multi_domain_interaction() routines.

Boolean to indicate whether to document basic info (to screen) during setup_multi_domain_interaction() routines.

Definition at line 168 of file multi_domain.cc.

Referenced by oomph::RefineableTriangleMesh< ELEMENT >::adapt().

◆ Doc_timings

bool oomph::Multi_domain_functions::Doc_timings =false

Boolean to indicate whether to doc timings or not.

Definition at line 164 of file multi_domain.cc.

Referenced by oomph::RefineableTriangleMesh< ELEMENT >::adapt().

◆ External_element_located

Vector< Vector< unsigned > > oomph::Multi_domain_functions::External_element_located

Lookup scheme for whether a local element's integration point has had an external element assigned to it – essentially boolean. External_element_located[e][ipt] = {0,1} if external element for ipt-th integration in local element e {has not, has} been found. Used locally to ensure that we're not searching for the same elements over and over again when we go around the spirals.

Definition at line 74 of file multi_domain.cc.

Referenced by aux_setup_multi_domain_interaction(), and create_external_halo_elements().

◆ Flat_packed_doubles

Vector< double > oomph::Multi_domain_functions::Flat_packed_doubles

Vector of flat-packed doubles to be communicated with other processors.

Definition at line 112 of file multi_domain.cc.

Referenced by construct_new_external_halo_master_node_helper(), oomph::RefineableTriangleMesh< ELEMENT >::construct_new_halo_node_helper(), oomph::RefineableTriangleMesh< ELEMENT >::construct_new_node_load_balance_helper(), create_external_halo_elements(), oomph::RefineableTriangleMesh< ELEMENT >::get_required_nodal_information_helper(), oomph::RefineableTriangleMesh< ELEMENT >::get_required_nodal_information_load_balance_helper(), oomph::RefineableTriangleMesh< ELEMENT >::load_balance(), recursively_add_masters_of_external_halo_node_to_storage(), oomph::RefineableTriangleMesh< ELEMENT >::reset_halo_haloed_scheme(), oomph::RefineableTriangleMesh< ELEMENT >::reset_halo_haloed_scheme_helper(), and oomph::RefineableTriangleMesh< ELEMENT >::send_and_receive_elements_nodes_info().

◆ Flat_packed_located_coordinates

Vector< double > oomph::Multi_domain_functions::Flat_packed_located_coordinates

Vector of flat-packed local coordinates for zeta tuples that have been located.

Definition at line 108 of file multi_domain.cc.

Referenced by create_external_halo_elements().

◆ Flat_packed_unsigneds

Vector< unsigned > oomph::Multi_domain_functions::Flat_packed_unsigneds

Vector of flat-packed unsigneds to be communicated with other processors – this is really "private" data, declared here to avoid having to pass the array between the various helper functions.

Definition at line 124 of file multi_domain.cc.

◆ Flat_packed_unsigneds_string

Vector< std::string > oomph::Multi_domain_functions::Flat_packed_unsigneds_string

Definition at line 131 of file multi_domain.cc.

◆ Flat_packed_zetas_not_found_locally

Vector< double > oomph::Multi_domain_functions::Flat_packed_zetas_not_found_locally

Vector of flat-packed zeta coordinates for which the external element could not be found during current local search. These will be sent to the next processor in the ring-like parallel search. The zeta coordinates come in groups of Dim (scalar) coordinates.

Definition at line 80 of file multi_domain.cc.

Referenced by aux_setup_multi_domain_interaction().

◆ Located_element_status

Vector< unsigned > oomph::Multi_domain_functions::Located_element_status

Vector to indicate (to another processor) whether a located element (that will have to represented as an external halo element on that processor) should be newly created on that processor (2), already exists on that processor (1), or is not on the current processor either (0).

Definition at line 104 of file multi_domain.cc.

Referenced by create_external_halo_elements().

◆ Proc_id_plus_one_of_external_element

Vector< int > oomph::Multi_domain_functions::Proc_id_plus_one_of_external_element

Proc_id_plus_one_of_external_element[i] contains the processor id (plus one) of the processor on which the i-th zeta coordinate tuple received from elsewhere (in the order in which these are stored in Received_flat_packed_zetas_to_be_found) was located; it's zero if it wasn't found during the current stage of the ring-like parallel search.

Definition at line 97 of file multi_domain.cc.

Referenced by create_external_halo_elements().

◆ Received_flat_packed_zetas_to_be_found

Vector< double > oomph::Multi_domain_functions::Received_flat_packed_zetas_to_be_found

Vector of flat-packed zeta coordinates for which the external element could not be found on another processor and for which we're currently searching here. Whatever can't be found here, gets written into Flat_packed_zetas_not_found_locally and then passed on to the next processor during the ring-like parallel search. The zeta coordinates come in groups of Dim (scalar) coordinates.

Definition at line 88 of file multi_domain.cc.

◆ Use_bulk_element_as_external

bool oomph::Multi_domain_functions::Use_bulk_element_as_external =false

Boolean to indicate when to use the bulk element as the external element. Defaults to false, you must have set up FaceElements properly first in order for it to work.

Definition at line 148 of file multi_domain.cc.

Referenced by aux_setup_multi_domain_interaction(), and setup_multi_domain_interaction().

◆ Zeta_coords_for_further_away_comparison

Vector<double> oomph::Multi_domain_functions::Zeta_coords_for_further_away_comparison

Vector of zeta coordinates that we're currently trying to locate; used in sorting of bin entries in further_away() comparison function

Definition at line 2290 of file multi_domain.cc.

Enumerations

Functions

Variables

Enumeration Type Documentation

◆ anonymous enum

Function Documentation

◆ add_external_halo_master_node_helper()

◆ add_external_halo_node_helper()

◆ add_external_halo_node_to_storage()

◆ add_external_haloed_master_node_helper()

◆ add_external_haloed_node_helper()

◆ add_external_haloed_node_to_storage()

◆ aux_setup_multi_domain_interaction() [1/2]

◆ aux_setup_multi_domain_interaction() [2/2]

◆ clean_up()

◆ construct_new_external_halo_master_node_helper()

◆ construct_new_external_halo_node_helper()

◆ create_external_halo_elements()

◆ first_closer_than_second()

◆ get_dim_helper()

◆ get_required_master_nodal_information_helper()

◆ get_required_nodal_information_helper()

◆ locate_zeta_for_local_coordinates()

◆ locate_zeta_for_missing_coordinates()

◆ recursively_add_masters_of_external_halo_node_to_storage()

◆ recursively_add_masters_of_external_haloed_node()

◆ send_and_receive_located_info()

◆ send_and_receive_missing_zetas()

◆ setup_bulk_elements_adjacent_to_face_mesh() [1/2]

◆ setup_bulk_elements_adjacent_to_face_mesh() [2/2]

◆ setup_multi_domain_interaction() [1/3]

◆ setup_multi_domain_interaction() [2/3]

◆ setup_multi_domain_interaction() [3/3]

◆ setup_multi_domain_interactions()

Variable Documentation

◆ Accept_failed_locate_zeta_in_setup_multi_domain_interaction

◆ Allow_use_of_halo_elements_as_external_elements

◆ Allow_use_of_halo_elements_as_external_elements_for_projection

◆ Counter_for_flat_packed_doubles

◆ Counter_for_flat_packed_unsigneds

◆ Dim

◆ Doc_boundary_coordinate_file

◆ Doc_full_stats

◆ Doc_stats

◆ Doc_timings

◆ External_element_located

◆ Flat_packed_doubles

◆ Flat_packed_located_coordinates

◆ Flat_packed_unsigneds

◆ Flat_packed_unsigneds_string

◆ Flat_packed_zetas_not_found_locally

◆ Located_element_status

◆ Proc_id_plus_one_of_external_element

◆ Received_flat_packed_zetas_to_be_found

◆ Use_bulk_element_as_external

◆ Zeta_coords_for_further_away_comparison