#include <pseudo_elastic_fsi_preconditioner.h>

Inheritance diagram for oomph::PseudoElasticFSIPreconditioner:

Public Member Functions
	PseudoElasticFSIPreconditioner (const unsigned &dim, Problem *problem_pt)

virtual	~PseudoElasticFSIPreconditioner ()

	PseudoElasticFSIPreconditioner (const PseudoElasticFSIPreconditioner &)=delete
	Broken copy constructor. More...

void	clean_up_memory ()
	Broken assignment operator. More...

void	setup ()
	Setup the precoonditioner. More...

void	preconditioner_solve (const DoubleVector &r, DoubleVector &z)
	Apply the preconditioner. More...

void	set_fluid_and_pseudo_elastic_mesh_pt (Mesh *mesh_pt)
	specify the mesh containing the combined fluid/pseudo solid elements More...

void	set_solid_mesh_pt (Mesh *mesh_pt)
	specify the mesh containing the solid elements More...

void	set_lagrange_multiplier_mesh_pt (Mesh *mesh_pt)
	specify the mesh containing the lagrange multiplier elements More...

void	set_solid_preconditioner (Preconditioner *prec_pt)

PseudoElasticPreconditioner *const	pseudo_elastic_preconditioner_pt ()
	Access function to the pseudo elastic subsidiary preconditioner. More...

NavierStokesSchurComplementPreconditioner *const	navier_stokes_schur_complement_preconditioner_pt ()
	Access function to the Navier Stokes Schur complement preconditioner. More...

void	enable_navier_stokes_schur_complement_preconditioner ()

void	disable_navier_stokes_schur_complement_preconditioner ()

Public Member Functions inherited from oomph::BlockPreconditioner< CRDoubleMatrix >
	BlockPreconditioner ()
	Constructor. More...

	BlockPreconditioner (const BlockPreconditioner &)=delete
	Broken copy constructor. More...

virtual	~BlockPreconditioner ()
	Destructor. More...

void	operator= (const BlockPreconditioner &)=delete
	Broken assignment operator. More...

CRDoubleMatrix *	matrix_pt () const

void	turn_on_recursive_debug_flag ()

void	turn_off_recursive_debug_flag ()

void	turn_on_debug_flag ()
	Toggles on the debug flag. More...

void	turn_off_debug_flag ()
	Toggles off the debug flag. More...

void	turn_into_subsidiary_block_preconditioner (BlockPreconditioner< CRDoubleMatrix > *master_block_prec_pt, const Vector< unsigned > &doftype_in_master_preconditioner_coarse)

void	turn_into_subsidiary_block_preconditioner (BlockPreconditioner< CRDoubleMatrix > *master_block_prec_pt, const Vector< unsigned > &doftype_in_master_preconditioner_coarse, const Vector< Vector< unsigned >> &doftype_coarsen_map_coarse)

virtual void	block_setup ()

void	block_setup (const Vector< unsigned > &dof_to_block_map)

void	get_block (const unsigned &i, const unsigned &j, CRDoubleMatrix &output_matrix, const bool &ignore_replacement_block=false) const

CRDoubleMatrix	get_block (const unsigned &i, const unsigned &j, const bool &ignore_replacement_block=false) const

void	set_master_matrix_pt (CRDoubleMatrix *in_matrix_pt)
	Set the matrix_pt in the upper-most master preconditioner. More...

void	get_block_other_matrix (const unsigned &i, const unsigned &j, CRDoubleMatrix *in_matrix_pt, CRDoubleMatrix &output_matrix)

void	get_blocks (DenseMatrix< bool > &required_blocks, DenseMatrix< CRDoubleMatrix * > &block_matrix_pt) const

void	get_dof_level_block (const unsigned &i, const unsigned &j, CRDoubleMatrix &output_block, const bool &ignore_replacement_block=false) const

void	get_dof_level_block (const unsigned &block_i, const unsigned &block_j, CRDoubleMatrix &output_block, const bool &ignore_replacement_block) const

CRDoubleMatrix	get_concatenated_block (const VectorMatrix< BlockSelector > &selected_block)

void	get_concatenated_block_vector (const Vector< unsigned > &block_vec_number, const DoubleVector &v, DoubleVector &b)

void	return_concatenated_block_vector (const Vector< unsigned > &block_vec_number, const DoubleVector &b, DoubleVector &v) const
	Takes concatenated block ordered vector, b, and copies its. More...

void	get_block_vectors (const Vector< unsigned > &block_vec_number, const DoubleVector &v, Vector< DoubleVector > &s) const

void	get_block_vectors (const DoubleVector &v, Vector< DoubleVector > &s) const

void	return_block_vectors (const Vector< unsigned > &block_vec_number, const Vector< DoubleVector > &s, DoubleVector &v) const

void	return_block_vectors (const Vector< DoubleVector > &s, DoubleVector &v) const

void	get_block_vector (const unsigned &n, const DoubleVector &v, DoubleVector &b) const

void	return_block_vector (const unsigned &n, const DoubleVector &b, DoubleVector &v) const

void	get_block_ordered_preconditioner_vector (const DoubleVector &v, DoubleVector &w)

void	return_block_ordered_preconditioner_vector (const DoubleVector &w, DoubleVector &v) const

unsigned	nblock_types () const
	Return the number of block types. More...

unsigned	ndof_types () const
	Return the total number of DOF types. More...

const Mesh *	mesh_pt (const unsigned &i) const

unsigned	nmesh () const

int	block_number (const unsigned &i_dof) const
	Return the block number corresponding to a global index i_dof. More...

int	index_in_block (const unsigned &i_dof) const

const LinearAlgebraDistribution *	block_distribution_pt (const unsigned &b) const
	Access function to the block distributions (const version). More...

LinearAlgebraDistribution *	block_distribution_pt (const unsigned b)
	Access function to the block distributions (non-const version). More...

LinearAlgebraDistribution *	dof_block_distribution_pt (const unsigned &b)
	Access function to the dof-level block distributions. More...

const LinearAlgebraDistribution *	master_distribution_pt () const

unsigned	ndof_types_in_mesh (const unsigned &i) const

bool	is_subsidiary_block_preconditioner () const

bool	is_master_block_preconditioner () const

void	set_block_output_to_files (const std::string &basefilename)

void	disable_block_output_to_files ()
	Turn off output of blocks (by clearing the basefilename string). More...

bool	block_output_on () const
	Test if output of blocks is on or not. More...

void	output_blocks_to_files (const std::string &basefilename, const unsigned &precision=8) const

void	post_block_matrix_assembly_partial_clear ()

BlockPreconditioner< CRDoubleMatrix > *	master_block_preconditioner_pt () const
	Access function to the master block preconditioner pt. More...

void	clear_block_preconditioner_base ()

void	document ()

Vector< Vector< unsigned > >	doftype_coarsen_map_fine () const

Vector< unsigned >	get_fine_grain_dof_types_in (const unsigned &i) const

unsigned	nfine_grain_dof_types_in (const unsigned &i) const

MapMatrix< unsigned, CRDoubleMatrix * >	replacement_dof_block_pt () const
	Access function to the replaced dof-level blocks. More...

void	setup_matrix_vector_product (MatrixVectorProduct matvec_prod_pt, CRDoubleMatrix block_pt, const Vector< unsigned > &block_col_indices)

void	setup_matrix_vector_product (MatrixVectorProduct matvec_prod_pt, CRDoubleMatrix block_pt, const unsigned &block_col_index)

void	internal_get_block_ordered_preconditioner_vector (const DoubleVector &v, DoubleVector &w) const

void	internal_return_block_ordered_preconditioner_vector (const DoubleVector &w, DoubleVector &v) const

unsigned	internal_nblock_types () const

unsigned	internal_ndof_types () const

void	internal_return_block_vector (const unsigned &n, const DoubleVector &b, DoubleVector &v) const

void	internal_get_block_vector (const unsigned &n, const DoubleVector &v, DoubleVector &b) const

void	internal_get_block_vectors (const Vector< unsigned > &block_vec_number, const DoubleVector &v, Vector< DoubleVector > &s) const

void	internal_get_block_vectors (const DoubleVector &v, Vector< DoubleVector > &s) const

void	internal_return_block_vectors (const Vector< unsigned > &block_vec_number, const Vector< DoubleVector > &s, DoubleVector &v) const

void	internal_return_block_vectors (const Vector< DoubleVector > &s, DoubleVector &v) const

void	internal_get_block (const unsigned &i, const unsigned &j, CRDoubleMatrix &output_block) const

void	internal_get_block (const unsigned &block_i, const unsigned &block_j, CRDoubleMatrix &output_block) const

int	internal_block_number (const unsigned &i_dof) const

int	internal_index_in_block (const unsigned &i_dof) const

const LinearAlgebraDistribution *	internal_block_distribution_pt (const unsigned &b) const
	Access function to the internal block distributions. More...

void	insert_auxiliary_block_distribution (const Vector< unsigned > &block_vec_number, LinearAlgebraDistribution *dist_pt)

void	block_matrix_test (const unsigned &i, const unsigned &j, const CRDoubleMatrix *block_matrix_pt) const

int	get_index_of_value (const Vector< myType > &vec, const myType val, const bool sorted=false) const

Public Member Functions inherited from oomph::Preconditioner
	Preconditioner ()
	Constructor. More...

	Preconditioner (const Preconditioner &)=delete
	Broken copy constructor. More...

void	operator= (const Preconditioner &)=delete
	Broken assignment operator. More...

virtual	~Preconditioner ()
	Destructor (empty) More...

virtual void	preconditioner_solve_transpose (const DoubleVector &r, DoubleVector &z)

void	setup (DoubleMatrixBase *matrix_pt)

void	setup (const Problem problem_pt, DoubleMatrixBase matrix_pt)

void	enable_silent_preconditioner_setup ()
	Set up the block preconditioner quietly! More...

void	disable_silent_preconditioner_setup ()
	Be verbose in the block preconditioner setup. More...

virtual void	set_matrix_pt (DoubleMatrixBase *matrix_pt)
	Set the matrix pointer. More...

virtual const OomphCommunicator *	comm_pt () const
	Get function for comm pointer. More...

virtual void	set_comm_pt (const OomphCommunicator *const comm_pt)
	Set the communicator pointer. More...

double	setup_time () const
	Returns the time to setup the preconditioner. More...

virtual void	turn_into_subsidiary_block_preconditioner (BlockPreconditioner< CRDoubleMatrix > *master_block_prec_pt, const Vector< unsigned > &doftype_in_master_preconditioner_coarse)

virtual void	turn_into_subsidiary_block_preconditioner (BlockPreconditioner< CRDoubleMatrix > *master_block_prec_pt, const Vector< unsigned > &doftype_in_master_preconditioner_coarse, const Vector< Vector< unsigned >> &doftype_coarsen_map_coarse)

Public Member Functions inherited from oomph::DistributableLinearAlgebraObject
	DistributableLinearAlgebraObject ()
	Default constructor - create a distribution. More...

	DistributableLinearAlgebraObject (const DistributableLinearAlgebraObject &matrix)=delete
	Broken copy constructor. More...

void	operator= (const DistributableLinearAlgebraObject &)=delete
	Broken assignment operator. More...

virtual	~DistributableLinearAlgebraObject ()
	Destructor. More...

LinearAlgebraDistribution *	distribution_pt () const
	access to the LinearAlgebraDistribution More...

unsigned	nrow () const
	access function to the number of global rows. More...

unsigned	nrow_local () const
	access function for the num of local rows on this processor. More...

unsigned	nrow_local (const unsigned &p) const
	access function for the num of local rows on this processor. More...

unsigned	first_row () const
	access function for the first row on this processor More...

unsigned	first_row (const unsigned &p) const
	access function for the first row on this processor More...

bool	distributed () const
	distribution is serial or distributed More...

bool	distribution_built () const

void	build_distribution (const LinearAlgebraDistribution *const dist_pt)

void	build_distribution (const LinearAlgebraDistribution &dist)

Private Attributes
PseudoElasticPreconditioner *	Pseudo_elastic_preconditioner_pt
	pointer to the pseudo solid preconditioner More...

Preconditioner *	Navier_stokes_preconditioner_pt
	pointer to the navier stokes precondtioner More...

NavierStokesSchurComplementPreconditioner *	Navier_stokes_schur_complement_preconditioner_pt
	Navier Stokes Schur complement preconditioner. More...

Preconditioner *	Solid_preconditioner_pt
	pointer to the solid preconditioner More...

bool	Using_default_solid_preconditioner

bool	Solid_preconditioner_is_block_preconditioner

MatrixVectorProduct *	Fluid_pseudo_elastic_matvec_pt
	fluid onto pseudosolid matrix vector operator More...

MatrixVectorProduct *	Solid_fluid_matvec_pt
	solid onto fluid matrix vector operatio More...

MatrixVectorProduct *	Solid_pseudo_elastic_matvec_pt
	solid onto pseudo solid matrix vector operatio More...

MatrixVectorProduct *	Lagrange_solid_matvec_pt

Mesh *	Fluid_and_pseudo_elastic_mesh_pt
	Mesh containing the combined fluid and pseudo solid element. More...

Mesh *	Solid_mesh_pt
	Mesh containing the solid elements. More...

Mesh *	Lagrange_multiplier_mesh_pt
	Mesh containing the lagrange multiplier elements. More...

unsigned	Dim
	the dimension of the fluid More...

bool	Use_navier_stokes_schur_complement_preconditioner

Additional Inherited Members
Protected Member Functions inherited from oomph::BlockPreconditioner< CRDoubleMatrix >
void	set_nmesh (const unsigned &n)

void	set_mesh (const unsigned &i, const Mesh *const mesh_pt, const bool &allow_multiple_element_type_in_mesh=false)

void	set_replacement_dof_block (const unsigned &block_i, const unsigned &block_j, CRDoubleMatrix *replacement_dof_block_pt)

bool	any_mesh_distributed () const

int	internal_dof_number (const unsigned &i_dof) const

unsigned	internal_index_in_dof (const unsigned &i_dof) const

unsigned	internal_block_dimension (const unsigned &b) const

unsigned	internal_dof_block_dimension (const unsigned &i) const

unsigned	master_nrow () const

unsigned	internal_master_dof_number (const unsigned &b) const

const LinearAlgebraDistribution *	internal_preconditioner_matrix_distribution_pt () const

const LinearAlgebraDistribution *	preconditioner_matrix_distribution_pt () const

Protected Member Functions inherited from oomph::DistributableLinearAlgebraObject
void	clear_distribution ()

Protected Attributes inherited from oomph::BlockPreconditioner< CRDoubleMatrix >
MapMatrix< unsigned, CRDoubleMatrix * >	Replacement_dof_block_pt
	The replacement dof-level blocks. More...

Vector< LinearAlgebraDistribution * >	Block_distribution_pt
	The distribution for the blocks. More...

Vector< Vector< unsigned > >	Block_to_dof_map_coarse

Vector< Vector< unsigned > >	Block_to_dof_map_fine
	Mapping for the block types to the most fine grain dof types. More...

Vector< Vector< unsigned > >	Doftype_coarsen_map_coarse

Vector< Vector< unsigned > >	Doftype_coarsen_map_fine

Vector< LinearAlgebraDistribution * >	Internal_block_distribution_pt
	Storage for the default distribution for each internal block. More...

Vector< LinearAlgebraDistribution * >	Dof_block_distribution_pt

Vector< unsigned >	Allow_multiple_element_type_in_mesh

Vector< const Mesh * >	Mesh_pt

Vector< unsigned >	Ndof_types_in_mesh

unsigned	Internal_nblock_types

unsigned	Internal_ndof_types

Protected Attributes inherited from oomph::Preconditioner
bool	Silent_preconditioner_setup
	Boolean to indicate whether or not the build should be done silently. More...

std::ostream *	Stream_pt
	Pointer to the output stream – defaults to std::cout. More...

Detailed Description

Preconditioner for FSI problems with pseudo-elastic fluid node updates. Note: NavierStokesSchurComplementPreconditioner is applied to the Navier Stokes subsidiary system. Default solid preconditioner is SuperLUPreconditioner. Enumeration of Elastic DOF types in the Pseudo-Elastic Elements The method get_dof_types_for_unknowns() must be implemented such that DOFs subject be Lagrange multiplier and DOFs NOT subject to Lagrange multiplier have different labels. For example in a 3D problem there are 6 DOF types and the following labelling must be implemented: 0 - x displacement (without lagr mult traction) 1 - y displacement (without lagr mult traction) 2 - z displacement (without lagr mult traction) 3 - x displacement (with lagr mult traction) 4 - y displacement (with lagr mult traction) 5 - z displacement (with lagr mult traction)

Constructor & Destructor Documentation

◆ PseudoElasticFSIPreconditioner() [1/2]

oomph::PseudoElasticFSIPreconditioner::PseudoElasticFSIPreconditioner	(	const unsigned &	dim,
		Problem *	problem_pt
	)

inline

constructor - just set defaults. Specify the spatial dimension of the fluid and a (non-const) problem pointer needed for the underlying NavierStokesSchurComplementPreconditioner.

       : Dim(dim)
     {
       Use_navier_stokes_schur_complement_preconditioner = true;
  
       // set the number of meshes
       this->set_nmesh(3);
  
       // null pointers
       Fluid_and_pseudo_elastic_mesh_pt = 0;
       Solid_mesh_pt = 0;
       Lagrange_multiplier_mesh_pt = 0;
  
       // create the pseudo solid preconditioner
       Pseudo_elastic_preconditioner_pt = new PseudoElasticPreconditioner();
  
       // using Schur complement preconditioner for NS
       Navier_stokes_preconditioner_pt = new SuperLUPreconditioner;
       Navier_stokes_schur_complement_preconditioner_pt =
         new NavierStokesSchurComplementPreconditioner(problem_pt);
  
       // set defaults
       Using_default_solid_preconditioner = true;
  
       // default super lu
       Solid_preconditioner_pt = new SuperLUPreconditioner;
  
       // create the matrix vector product operatrs
       Solid_fluid_matvec_pt = new MatrixVectorProduct;
       Solid_pseudo_elastic_matvec_pt = new MatrixVectorProduct;
       Fluid_pseudo_elastic_matvec_pt = new MatrixVectorProduct;
       Lagrange_solid_matvec_pt = new MatrixVectorProduct;
     }

References Fluid_and_pseudo_elastic_mesh_pt, Fluid_pseudo_elastic_matvec_pt, Lagrange_multiplier_mesh_pt, Lagrange_solid_matvec_pt, Navier_stokes_preconditioner_pt, Navier_stokes_schur_complement_preconditioner_pt, Pseudo_elastic_preconditioner_pt, oomph::BlockPreconditioner< CRDoubleMatrix >::set_nmesh(), Solid_fluid_matvec_pt, Solid_mesh_pt, Solid_preconditioner_pt, Solid_pseudo_elastic_matvec_pt, Use_navier_stokes_schur_complement_preconditioner, and Using_default_solid_preconditioner.

◆ ~PseudoElasticFSIPreconditioner()

virtual oomph::PseudoElasticFSIPreconditioner::~PseudoElasticFSIPreconditioner ( )

inlinevirtual

     {
       // clean the memory
       this->clean_up_memory();
  
       // delete the pseudo solid preconditioner
       delete Pseudo_elastic_preconditioner_pt;
  
       // delete the navier stokes preconditioner
       delete Navier_stokes_preconditioner_pt;
       delete Navier_stokes_schur_complement_preconditioner_pt;
  
       // delete the solid preconditioner if default
       if (Using_default_solid_preconditioner)
       {
         delete Solid_preconditioner_pt;
       }
  
       // delete the matrix vector product operators
       delete Fluid_pseudo_elastic_matvec_pt;
       delete Solid_fluid_matvec_pt;
       delete Solid_pseudo_elastic_matvec_pt;
       delete Lagrange_solid_matvec_pt;
     }

References clean_up_memory(), Fluid_pseudo_elastic_matvec_pt, Lagrange_solid_matvec_pt, Navier_stokes_preconditioner_pt, Navier_stokes_schur_complement_preconditioner_pt, Pseudo_elastic_preconditioner_pt, Solid_fluid_matvec_pt, Solid_preconditioner_pt, Solid_pseudo_elastic_matvec_pt, and Using_default_solid_preconditioner.

◆ PseudoElasticFSIPreconditioner() [2/2]

oomph::PseudoElasticFSIPreconditioner::PseudoElasticFSIPreconditioner ( const PseudoElasticFSIPreconditioner & )

delete

Broken copy constructor.

Member Function Documentation

◆ clean_up_memory()

void oomph::PseudoElasticFSIPreconditioner::clean_up_memory ( )

virtual

Broken assignment operator.

clean up memory method

Reimplemented from oomph::Preconditioner.

   {
     // wipe the preconditioner
     Pseudo_elastic_preconditioner_pt->clean_up_memory();
     Navier_stokes_preconditioner_pt->clean_up_memory();
     Navier_stokes_schur_complement_preconditioner_pt->clean_up_memory();
     Solid_preconditioner_pt->clean_up_memory();
  
     // clean the subsidiary matvec operators
     Fluid_pseudo_elastic_matvec_pt->clean_up_memory();
     Solid_fluid_matvec_pt->clean_up_memory();
     Solid_pseudo_elastic_matvec_pt->clean_up_memory();
     Lagrange_solid_matvec_pt->clean_up_memory();
   }

References oomph::MatrixVectorProduct::clean_up_memory(), oomph::Preconditioner::clean_up_memory(), oomph::PseudoElasticPreconditioner::clean_up_memory(), oomph::NavierStokesSchurComplementPreconditioner::clean_up_memory(), Fluid_pseudo_elastic_matvec_pt, Lagrange_solid_matvec_pt, Navier_stokes_preconditioner_pt, Navier_stokes_schur_complement_preconditioner_pt, Pseudo_elastic_preconditioner_pt, Solid_fluid_matvec_pt, Solid_preconditioner_pt, and Solid_pseudo_elastic_matvec_pt.

Referenced by setup(), and ~PseudoElasticFSIPreconditioner().

◆ disable_navier_stokes_schur_complement_preconditioner()

void oomph::PseudoElasticFSIPreconditioner::disable_navier_stokes_schur_complement_preconditioner ( )

inline

Call to use the SuperLUPreconditioner is used for the Navier Stokes subsidiary system.

     {
       Use_navier_stokes_schur_complement_preconditioner = false;
     }

References Use_navier_stokes_schur_complement_preconditioner.

Referenced by FSIChannelWithLeafletProblem< ELEMENT >::set_iterative_solver().

◆ enable_navier_stokes_schur_complement_preconditioner()

void oomph::PseudoElasticFSIPreconditioner::enable_navier_stokes_schur_complement_preconditioner ( )

inline

Call to use the Navier Stokes Schur complement preconditioner.

     {
       Use_navier_stokes_schur_complement_preconditioner = true;
     }

References Use_navier_stokes_schur_complement_preconditioner.

Referenced by FSIChannelWithLeafletProblem< ELEMENT >::set_iterative_solver(), PseudoElasticCollapsibleChannelProblem< FLUID_ELEMENT, SOLID_ELEMENT >::set_pseudo_elastic_fsi_solver(), and set_pseudo_elastic_fsi_solver().

◆ navier_stokes_schur_complement_preconditioner_pt()

NavierStokesSchurComplementPreconditioner* const oomph::PseudoElasticFSIPreconditioner::navier_stokes_schur_complement_preconditioner_pt ( )

inline

Access function to the Navier Stokes Schur complement preconditioner.

     {
       return Navier_stokes_schur_complement_preconditioner_pt;
     }

References Navier_stokes_schur_complement_preconditioner_pt.

Referenced by FSIChannelWithLeafletProblem< ELEMENT >::set_iterative_solver(), PseudoElasticCollapsibleChannelProblem< FLUID_ELEMENT, SOLID_ELEMENT >::set_pseudo_elastic_fsi_solver(), and set_pseudo_elastic_fsi_solver().

◆ preconditioner_solve()

void oomph::PseudoElasticFSIPreconditioner::preconditioner_solve	(	const DoubleVector &	r,
		DoubleVector &	z
	)

virtual

Apply the preconditioner.

Implements oomph::Preconditioner.

   {
     // apply the "pseudo solid" component of the pseudo solid preconditioner
     Pseudo_elastic_preconditioner_pt->elastic_preconditioner_solve(r, z);
  
     // apply the fluid on pseudo solid matrix vector product operator
     DoubleVector x;
     this->get_block_vector(2, z, x);
     DoubleVector y;
     Fluid_pseudo_elastic_matvec_pt->multiply(x, y);
     DoubleVector w;
     Solid_pseudo_elastic_matvec_pt->multiply(x, w);
     x.clear();
     this->get_block_vector(0, r, x);
     x -= y;
     y.clear();
  
     // storage for a copy of z
     DoubleVector z_copy;
  
     // apply the ns preconditioner
     if (Use_navier_stokes_schur_complement_preconditioner)
     {
       z_copy.build(z);
       this->return_block_vector(0, x, z_copy);
       x.clear();
       Navier_stokes_schur_complement_preconditioner_pt->preconditioner_solve(
         z_copy, z);
       z_copy.clear();
     }
     else
     {
       Navier_stokes_preconditioner_pt->preconditioner_solve(x, y);
       x.clear();
       this->return_block_vector(0, y, z);
       y.clear();
     }
  
     // apply the solid onto fluid matrix vector product operator
     this->get_block_vector(0, z, x);
     Solid_fluid_matvec_pt->multiply(x, y);
     x.clear();
     this->get_block_vector(1, r, x);
     x -= y;
     y.clear();
  
     // apply the result of the solid onto pseudo solid matrix vector product
     x -= w;
     w.clear();
  
     // apply the solid preconditioner
     if (Solid_preconditioner_is_block_preconditioner)
     {
       DoubleVector z_copy(z);
       this->return_block_vector(1, x, z_copy);
       x.clear();
       (dynamic_cast<GeneralPurposeBlockPreconditioner<CRDoubleMatrix>*>(
          Solid_preconditioner_pt))
         ->preconditioner_solve(z_copy, z);
       this->get_block_vector(1, z, y);
     }
     else
     {
       Solid_preconditioner_pt->preconditioner_solve(x, y);
       x.clear();
       this->return_block_vector(1, y, z);
     }
  
     // apply the lagrange multiplier solid matrix vector product operator
     Lagrange_solid_matvec_pt->multiply(y, x);
     y.clear();
     this->get_block_vector(3, r, y);
     y -= x;
     x.clear();
     z_copy.build(z);
     this->return_block_vector(3, y, z_copy);
  
     // apply the lagrange multiplier compenent of the pseudo solid
     // preconditioner
     Pseudo_elastic_preconditioner_pt->lagrange_multiplier_preconditioner_solve(
       z_copy, z);
     z_copy.clear();
   }

References oomph::DoubleVector::build(), oomph::DoubleVector::clear(), oomph::PseudoElasticPreconditioner::elastic_preconditioner_solve(), Fluid_pseudo_elastic_matvec_pt, oomph::BlockPreconditioner< CRDoubleMatrix >::get_block_vector(), oomph::PseudoElasticPreconditioner::lagrange_multiplier_preconditioner_solve(), Lagrange_solid_matvec_pt, oomph::MatrixVectorProduct::multiply(), Navier_stokes_preconditioner_pt, Navier_stokes_schur_complement_preconditioner_pt, oomph::NavierStokesSchurComplementPreconditioner::preconditioner_solve(), oomph::Preconditioner::preconditioner_solve(), Pseudo_elastic_preconditioner_pt, UniformPSDSelfTest::r, oomph::BlockPreconditioner< CRDoubleMatrix >::return_block_vector(), Solid_fluid_matvec_pt, Solid_preconditioner_is_block_preconditioner, Solid_preconditioner_pt, Solid_pseudo_elastic_matvec_pt, Use_navier_stokes_schur_complement_preconditioner, w, plotDoE::x, and y.

◆ pseudo_elastic_preconditioner_pt()

PseudoElasticPreconditioner* const oomph::PseudoElasticFSIPreconditioner::pseudo_elastic_preconditioner_pt ( )

inline

Access function to the pseudo elastic subsidiary preconditioner.

     {
       return Pseudo_elastic_preconditioner_pt;
     }

References Pseudo_elastic_preconditioner_pt.

Referenced by FSIChannelWithLeafletProblem< ELEMENT >::set_iterative_solver(), PseudoElasticCollapsibleChannelProblem< FLUID_ELEMENT, SOLID_ELEMENT >::set_pseudo_elastic_fsi_solver(), and set_pseudo_elastic_fsi_solver().

◆ set_fluid_and_pseudo_elastic_mesh_pt()

void oomph::PseudoElasticFSIPreconditioner::set_fluid_and_pseudo_elastic_mesh_pt ( Mesh * mesh_pt )

inline

specify the mesh containing the combined fluid/pseudo solid elements

     {
       Fluid_and_pseudo_elastic_mesh_pt = mesh_pt;
     }

References Fluid_and_pseudo_elastic_mesh_pt, and oomph::BlockPreconditioner< CRDoubleMatrix >::mesh_pt().

Referenced by FSIChannelWithLeafletProblem< ELEMENT >::set_iterative_solver(), PseudoElasticCollapsibleChannelProblem< FLUID_ELEMENT, SOLID_ELEMENT >::set_pseudo_elastic_fsi_solver(), and set_pseudo_elastic_fsi_solver().

◆ set_lagrange_multiplier_mesh_pt()

void oomph::PseudoElasticFSIPreconditioner::set_lagrange_multiplier_mesh_pt ( Mesh * mesh_pt )

inline

specify the mesh containing the lagrange multiplier elements

     {
       Lagrange_multiplier_mesh_pt = mesh_pt;
     }

References Lagrange_multiplier_mesh_pt, and oomph::BlockPreconditioner< CRDoubleMatrix >::mesh_pt().

Referenced by FSIChannelWithLeafletProblem< ELEMENT >::set_iterative_solver(), PseudoElasticCollapsibleChannelProblem< FLUID_ELEMENT, SOLID_ELEMENT >::set_pseudo_elastic_fsi_solver(), and set_pseudo_elastic_fsi_solver().

◆ set_solid_mesh_pt()

void oomph::PseudoElasticFSIPreconditioner::set_solid_mesh_pt ( Mesh * mesh_pt )

inline

specify the mesh containing the solid elements

     {
       Solid_mesh_pt = mesh_pt;
     }

References oomph::BlockPreconditioner< CRDoubleMatrix >::mesh_pt(), and Solid_mesh_pt.

Referenced by FSIChannelWithLeafletProblem< ELEMENT >::set_iterative_solver(), PseudoElasticCollapsibleChannelProblem< FLUID_ELEMENT, SOLID_ELEMENT >::set_pseudo_elastic_fsi_solver(), and set_pseudo_elastic_fsi_solver().

◆ set_solid_preconditioner()

void oomph::PseudoElasticFSIPreconditioner::set_solid_preconditioner ( Preconditioner * prec_pt )

inline

speicify a non default solid preconditioner. This preconditioner will not delete it

     {
       if (Using_default_solid_preconditioner)
       {
         delete Solid_preconditioner_pt;
       }
       Solid_preconditioner_pt = prec_pt;
       Using_default_solid_preconditioner = false;
     }

References Solid_preconditioner_pt, and Using_default_solid_preconditioner.

Referenced by PseudoElasticCollapsibleChannelProblem< FLUID_ELEMENT, SOLID_ELEMENT >::set_pseudo_elastic_fsi_solver(), and set_pseudo_elastic_fsi_solver().

◆ setup()

void oomph::PseudoElasticFSIPreconditioner::setup ( )

virtual

Setup the precoonditioner.

Implements oomph::Preconditioner.

   {
     // clean the memory
     this->clean_up_memory();
  
 #ifdef PARANOID
     // paranoid check that the meshes have been set
     if (Fluid_and_pseudo_elastic_mesh_pt == 0)
     {
       std::ostringstream error_message;
       error_message << "The fluid and pseudo elastic mesh must be set.";
       throw OomphLibError(
         error_message.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
     }
     if (Solid_mesh_pt == 0)
     {
       std::ostringstream error_message;
       error_message << "The solid mesh must be set.";
       throw OomphLibError(
         error_message.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
     }
     if (Lagrange_multiplier_mesh_pt == 0)
     {
       std::ostringstream error_message;
       error_message << "The Lagrange multiplier mesh must be set.";
       throw OomphLibError(
         error_message.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
     }
 #endif
  
     // add the meshes
     this->set_mesh(0, Fluid_and_pseudo_elastic_mesh_pt);
     this->set_mesh(1, Solid_mesh_pt);
     this->set_mesh(2, Lagrange_multiplier_mesh_pt);
  
     // determine the number of fluid dofs
     unsigned nfluid_dof = Dim + 1;
  
     // determine the number of pseudo solid dofs
     unsigned npseudo_elastic_dof = this->ndof_types_in_mesh(0) - nfluid_dof;
  
     // determine the number of solid dofs
     unsigned nsolid_dof = this->ndof_types_in_mesh(1);
  
     // determine the number of lagrange multiplier dofs
     unsigned nlagr_mult_dof = this->ndof_types_in_mesh(2);
  
     // total number of dof types
     unsigned ntotal_dof =
       nfluid_dof + npseudo_elastic_dof + nsolid_dof + nlagr_mult_dof;
  
     // setup the block lookup scheme
     // block 0 - fluid
     //       1 - solid
     //       2 - pseudo solid inc. lagrange mult
     Vector<unsigned> dof_to_block_map(ntotal_dof, 0);
     int c = nfluid_dof;
     for (unsigned i = 0; i < npseudo_elastic_dof; i++)
     {
       dof_to_block_map[c] = 2;
       c++;
     }
     for (unsigned i = 0; i < nsolid_dof; i++)
     {
       dof_to_block_map[c] = 1;
       c++;
     }
     for (unsigned i = 0; i < nlagr_mult_dof; i++)
     {
       dof_to_block_map[c] = 3;
       c++;
     }
  
     // Recast Jacobian matrix to CRDoubleMatrix
     CRDoubleMatrix* cr_matrix_pt = dynamic_cast<CRDoubleMatrix*>(matrix_pt());
 #ifdef PARANOID
     if (cr_matrix_pt == 0)
     {
       std::ostringstream error_message;
       error_message << "FSIPreconditioner only works with"
                     << " CRDoubleMatrix matrices" << std::endl;
       throw OomphLibError(
         error_message.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
     }
 #endif
  
     // Call block setup for this preconditioner
     this->block_setup(dof_to_block_map);
  
     // SETUP THE PRECONDITIONERS
     // =========================
  
     // setup the navier stokes preconditioner
     if (Use_navier_stokes_schur_complement_preconditioner)
     {
       Navier_stokes_schur_complement_preconditioner_pt->set_navier_stokes_mesh(
         Fluid_and_pseudo_elastic_mesh_pt);
  
       Vector<unsigned> ns_dof_list(nfluid_dof, 0);
       for (unsigned i = 0; i < nfluid_dof; i++)
       {
         ns_dof_list[i] = i;
       }
  
       Navier_stokes_schur_complement_preconditioner_pt
         ->turn_into_subsidiary_block_preconditioner(this, ns_dof_list);
  
       Navier_stokes_schur_complement_preconditioner_pt->Preconditioner::setup(
         matrix_pt());
     }
     else
     {
       CRDoubleMatrix* ns_matrix_pt = new CRDoubleMatrix;
       this->get_block(0, 0, *ns_matrix_pt);
  
       Navier_stokes_preconditioner_pt->setup(ns_matrix_pt);
       delete ns_matrix_pt;
       ns_matrix_pt = 0;
     }
  
     // next the solid preconditioner
     if (dynamic_cast<BlockPreconditioner<CRDoubleMatrix>*>(
           Solid_preconditioner_pt) != 0)
     {
       Solid_preconditioner_is_block_preconditioner = true;
       GeneralPurposeBlockPreconditioner<CRDoubleMatrix>*
         solid_block_preconditioner_pt =
           dynamic_cast<GeneralPurposeBlockPreconditioner<CRDoubleMatrix>*>(
             Solid_preconditioner_pt);
  
       if (solid_block_preconditioner_pt != 0)
       {
         unsigned offset = nfluid_dof + npseudo_elastic_dof;
         Vector<unsigned> solid_prec_dof_list(nsolid_dof);
         for (unsigned i = 0; i < nsolid_dof; i++)
         {
           solid_prec_dof_list[i] = offset + i;
         }
         solid_block_preconditioner_pt
           ->turn_into_subsidiary_block_preconditioner(this,
                                                       solid_prec_dof_list);
         solid_block_preconditioner_pt->setup(cr_matrix_pt);
       }
       else
       {
         std::ostringstream error_message;
         error_message << "If the (real) solid preconditioner is a "
                       << "BlockPreconditioner then is must be a "
                       << "GeneralPurposeBlockPreconditioner";
         throw OomphLibError(error_message.str(),
                             OOMPH_CURRENT_FUNCTION,
                             OOMPH_EXCEPTION_LOCATION);
       }
     }
     // otherwise it is a solid preconditioner
     else
     {
       Solid_preconditioner_is_block_preconditioner = false;
       CRDoubleMatrix* s_matrix_pt = new CRDoubleMatrix;
       this->get_block(1, 1, *s_matrix_pt);
       Solid_preconditioner_pt->setup(s_matrix_pt);
       delete s_matrix_pt;
       s_matrix_pt = 0;
     }
  
     // next the pseudo solid preconditioner
     unsigned ndof_for_pseudo_elastic_prec = Dim * 3;
     Vector<unsigned> pseudo_elastic_prec_dof_list(ndof_for_pseudo_elastic_prec,
                                                   0);
     for (unsigned i = 0; i < Dim * 2; i++)
     {
       pseudo_elastic_prec_dof_list[i] = nfluid_dof + i;
     }
     for (unsigned i = 0; i < Dim; i++)
     {
       pseudo_elastic_prec_dof_list[i + Dim * 2] =
         nfluid_dof + npseudo_elastic_dof + nsolid_dof + i;
     }
     Pseudo_elastic_preconditioner_pt->turn_into_subsidiary_block_preconditioner(
       this, pseudo_elastic_prec_dof_list);
     Pseudo_elastic_preconditioner_pt->set_elastic_mesh(
       this->Fluid_and_pseudo_elastic_mesh_pt);
     Pseudo_elastic_preconditioner_pt->set_lagrange_multiplier_mesh(
       this->Lagrange_multiplier_mesh_pt);
     Pseudo_elastic_preconditioner_pt->Preconditioner::setup(matrix_pt());
  
     // SETUP THE MATRIX VECTOR PRODUCT OPERATORS
     // =========================================
  
     // setup the fluid pseudo-solid matvec operator
     CRDoubleMatrix* fp_matrix_pt = new CRDoubleMatrix;
     get_block(0, 2, *fp_matrix_pt);
     //  Fluid_pseudo_elastic_matvec_pt->setup(fp_matrix_pt);
     this->setup_matrix_vector_product(
       Fluid_pseudo_elastic_matvec_pt, fp_matrix_pt, 2);
     delete fp_matrix_pt;
     fp_matrix_pt = 0;
  
     // setup the solid fluid matvec operator
     CRDoubleMatrix* sf_matrix_pt = new CRDoubleMatrix;
     get_block(1, 0, *sf_matrix_pt);
     //  Solid_fluid_matvec_pt->setup(sf_matrix_pt);
     this->setup_matrix_vector_product(Solid_fluid_matvec_pt, sf_matrix_pt, 0);
     delete sf_matrix_pt;
     sf_matrix_pt = 0;
  
     // setup the solid pseudo-solid matvec operator
     CRDoubleMatrix* sp_matrix_pt = new CRDoubleMatrix;
     get_block(1, 2, *sp_matrix_pt);
     //  Solid_pseudo_elastic_matvec_pt->setup(sp_matrix_pt);
     this->setup_matrix_vector_product(
       Solid_pseudo_elastic_matvec_pt, sp_matrix_pt, 2);
     delete sp_matrix_pt;
     sp_matrix_pt = 0;
  
     // build the lagrange solid matvec operator
     CRDoubleMatrix* ls_matrix_pt = new CRDoubleMatrix;
     get_block(3, 1, *ls_matrix_pt);
     //  Lagrange_solid_matvec_pt->setup(ls_matrix_pt);
     this->setup_matrix_vector_product(
       Lagrange_solid_matvec_pt, ls_matrix_pt, 1);
     delete ls_matrix_pt;
     ls_matrix_pt = 0;
   }

References oomph::BlockPreconditioner< CRDoubleMatrix >::block_setup(), calibrate::c, clean_up_memory(), Dim, Fluid_and_pseudo_elastic_mesh_pt, Fluid_pseudo_elastic_matvec_pt, oomph::BlockPreconditioner< CRDoubleMatrix >::get_block(), i, Lagrange_multiplier_mesh_pt, Lagrange_solid_matvec_pt, oomph::BlockPreconditioner< CRDoubleMatrix >::matrix_pt(), Navier_stokes_preconditioner_pt, Navier_stokes_schur_complement_preconditioner_pt, oomph::BlockPreconditioner< CRDoubleMatrix >::ndof_types_in_mesh(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, Pseudo_elastic_preconditioner_pt, oomph::PseudoElasticPreconditioner::set_elastic_mesh(), oomph::PseudoElasticPreconditioner::set_lagrange_multiplier_mesh(), oomph::BlockPreconditioner< CRDoubleMatrix >::set_mesh(), oomph::NavierStokesSchurComplementPreconditioner::set_navier_stokes_mesh(), oomph::Preconditioner::setup(), oomph::BlockPreconditioner< CRDoubleMatrix >::setup_matrix_vector_product(), Solid_fluid_matvec_pt, Solid_mesh_pt, Solid_preconditioner_is_block_preconditioner, Solid_preconditioner_pt, Solid_pseudo_elastic_matvec_pt, oomph::BlockPreconditioner< MATRIX >::turn_into_subsidiary_block_preconditioner(), and Use_navier_stokes_schur_complement_preconditioner.

Member Data Documentation

◆ Dim

unsigned oomph::PseudoElasticFSIPreconditioner::Dim

private

the dimension of the fluid

Referenced by setup().

◆ Fluid_and_pseudo_elastic_mesh_pt

Mesh* oomph::PseudoElasticFSIPreconditioner::Fluid_and_pseudo_elastic_mesh_pt

private

Mesh containing the combined fluid and pseudo solid element.

Referenced by PseudoElasticFSIPreconditioner(), set_fluid_and_pseudo_elastic_mesh_pt(), and setup().

◆ Fluid_pseudo_elastic_matvec_pt

MatrixVectorProduct* oomph::PseudoElasticFSIPreconditioner::Fluid_pseudo_elastic_matvec_pt

private

fluid onto pseudosolid matrix vector operator

Referenced by clean_up_memory(), preconditioner_solve(), PseudoElasticFSIPreconditioner(), setup(), and ~PseudoElasticFSIPreconditioner().

◆ Lagrange_multiplier_mesh_pt

Mesh* oomph::PseudoElasticFSIPreconditioner::Lagrange_multiplier_mesh_pt

private

Mesh containing the lagrange multiplier elements.

Referenced by PseudoElasticFSIPreconditioner(), set_lagrange_multiplier_mesh_pt(), and setup().

◆ Lagrange_solid_matvec_pt

MatrixVectorProduct* oomph::PseudoElasticFSIPreconditioner::Lagrange_solid_matvec_pt

private

Referenced by clean_up_memory(), preconditioner_solve(), PseudoElasticFSIPreconditioner(), setup(), and ~PseudoElasticFSIPreconditioner().

◆ Navier_stokes_preconditioner_pt

Preconditioner* oomph::PseudoElasticFSIPreconditioner::Navier_stokes_preconditioner_pt

private

pointer to the navier stokes precondtioner

Referenced by clean_up_memory(), preconditioner_solve(), PseudoElasticFSIPreconditioner(), setup(), and ~PseudoElasticFSIPreconditioner().

◆ Navier_stokes_schur_complement_preconditioner_pt

NavierStokesSchurComplementPreconditioner* oomph::PseudoElasticFSIPreconditioner::Navier_stokes_schur_complement_preconditioner_pt

private

Navier Stokes Schur complement preconditioner.

Referenced by clean_up_memory(), navier_stokes_schur_complement_preconditioner_pt(), preconditioner_solve(), PseudoElasticFSIPreconditioner(), setup(), and ~PseudoElasticFSIPreconditioner().

◆ Pseudo_elastic_preconditioner_pt

PseudoElasticPreconditioner* oomph::PseudoElasticFSIPreconditioner::Pseudo_elastic_preconditioner_pt

private

pointer to the pseudo solid preconditioner

Referenced by clean_up_memory(), preconditioner_solve(), pseudo_elastic_preconditioner_pt(), PseudoElasticFSIPreconditioner(), setup(), and ~PseudoElasticFSIPreconditioner().

◆ Solid_fluid_matvec_pt

MatrixVectorProduct* oomph::PseudoElasticFSIPreconditioner::Solid_fluid_matvec_pt

private

solid onto fluid matrix vector operatio

Referenced by clean_up_memory(), preconditioner_solve(), PseudoElasticFSIPreconditioner(), setup(), and ~PseudoElasticFSIPreconditioner().

◆ Solid_mesh_pt

Mesh* oomph::PseudoElasticFSIPreconditioner::Solid_mesh_pt

private

Mesh containing the solid elements.

Referenced by PseudoElasticFSIPreconditioner(), set_solid_mesh_pt(), and setup().

◆ Solid_preconditioner_is_block_preconditioner

bool oomph::PseudoElasticFSIPreconditioner::Solid_preconditioner_is_block_preconditioner

private

boolean flag to indicate whether the Solid preconditioner is a block preconditioner

Referenced by preconditioner_solve(), and setup().

◆ Solid_preconditioner_pt

Preconditioner* oomph::PseudoElasticFSIPreconditioner::Solid_preconditioner_pt

private

pointer to the solid preconditioner

Referenced by clean_up_memory(), preconditioner_solve(), PseudoElasticFSIPreconditioner(), set_solid_preconditioner(), setup(), and ~PseudoElasticFSIPreconditioner().

◆ Solid_pseudo_elastic_matvec_pt

MatrixVectorProduct* oomph::PseudoElasticFSIPreconditioner::Solid_pseudo_elastic_matvec_pt

private

solid onto pseudo solid matrix vector operatio

Referenced by clean_up_memory(), preconditioner_solve(), PseudoElasticFSIPreconditioner(), setup(), and ~PseudoElasticFSIPreconditioner().

◆ Use_navier_stokes_schur_complement_preconditioner

bool oomph::PseudoElasticFSIPreconditioner::Use_navier_stokes_schur_complement_preconditioner

private

If true the Navier Stokes Schur complement preconditioner is used. Otherwise SuperLUPreconditioner is used for the Navier Stokes subsidiary system.

Referenced by disable_navier_stokes_schur_complement_preconditioner(), enable_navier_stokes_schur_complement_preconditioner(), preconditioner_solve(), PseudoElasticFSIPreconditioner(), and setup().

◆ Using_default_solid_preconditioner

bool oomph::PseudoElasticFSIPreconditioner::Using_default_solid_preconditioner

private

boolean flag to indicate whether default Solid preconditioner is used

Referenced by PseudoElasticFSIPreconditioner(), set_solid_preconditioner(), and ~PseudoElasticFSIPreconditioner().

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

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ PseudoElasticFSIPreconditioner() [1/2]

◆ ~PseudoElasticFSIPreconditioner()

◆ PseudoElasticFSIPreconditioner() [2/2]

Member Function Documentation

◆ clean_up_memory()

◆ disable_navier_stokes_schur_complement_preconditioner()

◆ enable_navier_stokes_schur_complement_preconditioner()

◆ navier_stokes_schur_complement_preconditioner_pt()

◆ preconditioner_solve()

◆ pseudo_elastic_preconditioner_pt()

◆ set_fluid_and_pseudo_elastic_mesh_pt()

◆ set_lagrange_multiplier_mesh_pt()

◆ set_solid_mesh_pt()

◆ set_solid_preconditioner()

◆ setup()

Member Data Documentation

◆ Dim

◆ Fluid_and_pseudo_elastic_mesh_pt

◆ Fluid_pseudo_elastic_matvec_pt

◆ Lagrange_multiplier_mesh_pt

◆ Lagrange_solid_matvec_pt

◆ Navier_stokes_preconditioner_pt

◆ Navier_stokes_schur_complement_preconditioner_pt

◆ Pseudo_elastic_preconditioner_pt

◆ Solid_fluid_matvec_pt

◆ Solid_mesh_pt

◆ Solid_preconditioner_is_block_preconditioner

◆ Solid_preconditioner_pt

◆ Solid_pseudo_elastic_matvec_pt

◆ Use_navier_stokes_schur_complement_preconditioner

◆ Using_default_solid_preconditioner