A class that is used to define the functions used when assembling the derivatives of the residuals with respect to a parameter. The idea is to replace get_residuals with get_dresiduals_dparameter with a particular parameter and assembly handler that are passed on assembly. More...

#include <assembly_handler.h>

Inheritance diagram for oomph::ParameterDerivativeHandler:

Public Member Functions
	ParameterDerivativeHandler (AssemblyHandler const &assembly_handler_pt, double const &parameter_pt)
	Store the original assembly handler and parameter. More...

unsigned	ndof (GeneralisedElement *const &elem_pt)
	Return the number of degrees of freedom in the element elem_pt Pass through to the original assembly handler. More...

unsigned long	eqn_number (GeneralisedElement *const &elem_pt, const unsigned &ieqn_local)
	Return the global equation number of the local unknown ieqn_local in elem_pt.Pass through to the original assembly handler. More...

void	get_residuals (GeneralisedElement *const &elem_pt, Vector< double > &residuals)
	Return the contribution to the residuals of the element elem_pt by using the derivatives. More...

void	get_jacobian (GeneralisedElement *const &elem_pt, Vector< double > &residuals, DenseMatrix< double > &jacobian)
	Calculate the elemental Jacobian matrix "d equation / d variable" for elem_pt. Overloaded to return the derivatives wrt the parameter. More...

Public Member Functions inherited from oomph::AssemblyHandler
	AssemblyHandler ()
	Empty constructor. More...

virtual void	dof_vector (GeneralisedElement *const &elem_pt, const unsigned &t, Vector< double > &dof)
	Return vector of dofs at time level t in the element elem_pt. More...

virtual void	dof_pt_vector (GeneralisedElement const &elem_pt, Vector< double > &dof_pt)
	Return vector of pointers to dofs in the element elem_pt. More...

virtual double &	local_problem_dof (Problem *const &problem_pt, const unsigned &t, const unsigned &i)
	Return the t-th level of storage associated with the i-th (local) dof stored in the problem. More...

virtual void	get_all_vectors_and_matrices (GeneralisedElement *const &elem_pt, Vector< Vector< double >> &vec, Vector< DenseMatrix< double >> &matrix)
	Calculate all desired vectors and matrices provided by the element elem_pt. More...

virtual void	get_dresiduals_dparameter (GeneralisedElement const &elem_pt, double const &parameter_pt, Vector< double > &dres_dparam)
	Calculate the derivative of the residuals with respect to a parameter. More...

virtual void	get_djacobian_dparameter (GeneralisedElement const &elem_pt, double const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam)
	Calculate the derivative of the residuals and jacobian with respect to a parameter. More...

virtual void	get_hessian_vector_products (GeneralisedElement *const &elem_pt, Vector< double > const &Y, DenseMatrix< double > const &C, DenseMatrix< double > &product)
	Calculate the product of the Hessian (derivative of Jacobian with respect to all variables) an eigenvector, Y, and other specified vectors, C (d(J_{ij})/d u_{k}) Y_{j} C_{k}. More...

virtual int	bifurcation_type () const
	Return an unsigned integer to indicate whether the handler is a bifurcation tracking handler. The default is zero (not) More...

virtual double *	bifurcation_parameter_pt () const
	Return a pointer to the bifurcation parameter in bifurcation tracking problems. More...

virtual void	get_eigenfunction (Vector< DoubleVector > &eigenfunction)
	Return the eigenfunction(s) associated with the bifurcation that has been detected in bifurcation tracking problems. More...

virtual void	get_inner_products (GeneralisedElement *const &elem_pt, Vector< std::pair< unsigned, unsigned >> const &history_index, Vector< double > &inner_product)
	Compute the inner products of the given vector of pairs of history values over the element. More...

virtual void	get_inner_product_vectors (GeneralisedElement *const &elem_pt, Vector< unsigned > const &history_index, Vector< Vector< double >> &inner_product_vector)
	Compute the vectors that when taken as a dot product with other history values give the inner product over the element. More...

virtual void	synchronise ()
	Function that is used to perform any synchronisation required during the solution. More...

virtual	~AssemblyHandler ()
	Empty virtual destructor. More...

Private Attributes
double *	Parameter_pt
	The value of the parameter. More...

AssemblyHandler *	Assembly_handler_pt
	The original assembly handler. More...

Detailed Description

A class that is used to define the functions used when assembling the derivatives of the residuals with respect to a parameter. The idea is to replace get_residuals with get_dresiduals_dparameter with a particular parameter and assembly handler that are passed on assembly.

Definition at line 338 of file assembly_handler.h.

Constructor & Destructor Documentation

◆ ParameterDerivativeHandler()

oomph::ParameterDerivativeHandler::ParameterDerivativeHandler	(	AssemblyHandler *const &	assembly_handler_pt,
		double *const &	parameter_pt
	)

inline

Store the original assembly handler and parameter.

Definition at line 348 of file assembly_handler.h.

Member Function Documentation

◆ eqn_number()

unsigned long oomph::ParameterDerivativeHandler::eqn_number	(	GeneralisedElement *const &	elem_pt,
		const unsigned &	ieqn_local
	)

inlinevirtual

Return the global equation number of the local unknown ieqn_local in elem_pt.Pass through to the original assembly handler.

Reimplemented from oomph::AssemblyHandler.

Definition at line 363 of file assembly_handler.h.

References Assembly_handler_pt, and oomph::AssemblyHandler::eqn_number().

◆ get_jacobian()

void oomph::ParameterDerivativeHandler::get_jacobian	(	GeneralisedElement *const &	elem_pt,
		Vector< double > &	residuals,
		DenseMatrix< double > &	jacobian
	)