Beta release!
Please note that the library has not been "officially" released. While we continue to work on the documentation, these web pages are likely to contain broken links and documents in draft form. Please send an email to

oomph-lib AT maths DOT man DOT ac DOT uk
if you wish to be informed of the library's "official" release.

oomph::AugmentedBlockPitchForkLinearSolver

oomph::AugmentedBlockPitchForkLinearSolver Class Reference

#include <assembly_handler.h>

Inheritance diagram for oomph::AugmentedBlockPitchForkLinearSolver:

List of all members.

Detailed Description

A custom linear solver class that is used to solve a block-factorised version of the PitchFork bifurcation detection problem.

Definition at line 450 of file assembly_handler.h.


Public Member Functions
	AugmentedBlockPitchForkLinearSolver (LinearSolver *const linear_solver_pt)
	Constructor, inherits the original linear solver.
	~AugmentedBlockPitchForkLinearSolver ()
	Clean up the memory that may have been allocated by the solver.
void	solve (Problem *const &problem_pt, DoubleVector &result)
	The solve function uses the block factorisation.
void	solve (DoubleMatrixBase *const &matrix_pt, const DoubleVector &rhs, DoubleVector &result)
	The interface is deliberately broken.
void	solve (DoubleMatrixBase *const &matrix_pt, const Vector< double > &rhs, Vector< double > &result)
	The linear-algebra-type solver does not make sense. The interface is deliberately broken.
void	resolve (const DoubleVector &rhs, DoubleVector &result)
	The resolve function also uses the block factorisation.
LinearSolver *	linear_solver_pt () const
	Access function to the original linear solver.
Private Attributes
LinearSolver *	Linear_solver_pt
	Pointer to the original linear solver.
Problem *	Problem_pt
	Pointer to the problem, used in the resolve.
DoubleVector *	Alpha_pt
	Pointer to the storage for the vector alpha.
DoubleVector *	E_pt
	Pointer to the storage for the vector e.

Constructor & Destructor Documentation

oomph::AugmentedBlockPitchForkLinearSolver::AugmentedBlockPitchForkLinearSolver ( LinearSolver *const linear_solver_pt ) [inline]

Constructor, inherits the original linear solver.

Definition at line 467 of file assembly_handler.h.

oomph::AugmentedBlockPitchForkLinearSolver::~AugmentedBlockPitchForkLinearSolver ( )

Clean up the memory that may have been allocated by the solver.

Definition at line 1796 of file assembly_handler.cc.

References Alpha_pt, and E_pt.

Member Function Documentation

LinearSolver* oomph::AugmentedBlockPitchForkLinearSolver::linear_solver_pt ( ) const [inline]

Access function to the original linear solver.

Definition at line 506 of file assembly_handler.h.

References Linear_solver_pt.

Referenced by oomph::PitchForkHandler::~PitchForkHandler().

void oomph::AugmentedBlockPitchForkLinearSolver::resolve	(	const DoubleVector &	rhs,
		DoubleVector &	result
	)			`[virtual]`

The resolve function also uses the block factorisation.

Reimplemented from oomph::LinearSolver.

Definition at line 2019 of file assembly_handler.cc.

References oomph::Problem::actions_after_change_in_bifurcation_parameter(), Alpha_pt, oomph::Problem::assembly_handler_pt(), oomph::DoubleVector::build(), oomph::Problem::communicator_pt(), oomph::LinearSolver::disable_resolve(), oomph::DistributableLinearAlgebraObject::Distribution_pt, oomph::DistributableLinearAlgebraObject::distribution_setup(), oomph::Problem::dof(), oomph::Mesh::element_pt(), oomph::LinearSolver::enable_resolve(), oomph::PitchForkHandler::eqn_number(), oomph::BlackBoxFDNewtonSolver::FD_step, oomph::PitchForkHandler::get_jacobian(), Linear_solver_pt, oomph::Problem::mesh_pt(), oomph::PitchForkHandler::ndof(), oomph::Problem::ndof(), oomph::Mesh::nelement(), OOMPH_EXCEPTION_LOCATION, Problem_pt, oomph::LinearAlgebraDistribution::rebuild(), oomph::LinearSolver::resolve(), oomph::PitchForkHandler::solve_augmented_block_system(), oomph::PitchForkHandler::solve_full_system(), and oomph::PitchForkHandler::Y.

void oomph::AugmentedBlockPitchForkLinearSolver::solve	(	DoubleMatrixBase *const &	matrix_pt,
		const Vector< double > &	rhs,
		Vector< double > &	result
	)			`[inline, virtual]`

The linear-algebra-type solver does not make sense. The interface is deliberately broken.

Reimplemented from oomph::LinearSolver.

Definition at line 491 of file assembly_handler.h.

References OOMPH_EXCEPTION_LOCATION.

void oomph::AugmentedBlockPitchForkLinearSolver::solve	(	DoubleMatrixBase *const &	matrix_pt,
		const DoubleVector &	rhs,
		DoubleVector &	result
	)			`[inline, virtual]`