#include <periodic_orbit_handler.h>

Inheritance diagram for oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >:

Public Member Functions
	PeriodicOrbitAssemblyHandler (Problem *const &problem_pt, const unsigned &n_element_in_period, const DenseMatrix< double > &initial_guess, const double &omega)
	Constructor, initialises values and constructs mesh of elements. More...

void	set_previous_dofs_to_current_dofs ()
	Update the previous dofs. More...

unsigned	ndof (GeneralisedElement *const &elem_pt)
	Return the number of degrees of freedom in the element elem_pt. 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. More...

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

void	get_dofs_for_element (GeneralisedElement *const elem_pt, Vector< double > &dofs)

void	get_previous_dofs_for_element (GeneralisedElement *const elem_pt, Vector< double > &dofs)

void	set_dofs_for_element (GeneralisedElement *const elem_pt, Vector< double > const &dofs)

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. More...

void	orbit_output (std::ostream &outfile, const unsigned &n_plot)
	Calculate all desired vectors and matrices provided by the element elem_pt. More...

void	discrete_times (Vector< double > &t)
	Tell me the times at which you want the solution. More...

void	adapt_temporal_mesh ()
	Adapt the time mesh. More...

	~PeriodicOrbitAssemblyHandler ()
	Destructor, destroy the time mesh. More...

Public Member Functions inherited from oomph::PeriodicOrbitAssemblyHandlerBase
	PeriodicOrbitAssemblyHandlerBase ()

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
PeriodicOrbitTimeDiscretisation *	Time_stepper_pt
	Pointer to the timestepper. More...

Problem *	Problem_pt
	Pointer to the problem. More...

PeriodicOrbitTemporalMesh< SpectralPeriodicOrbitElement< NNODE_1D > > *	Time_mesh_pt
	Storage for mesh of temporal elements. More...

Mesh *	Basic_time_mesh_pt
	Storage for the mesh of temporal elements with a simple mesh pointer. More...

Vector< double >	Previous_dofs
	Storage for the previous solution. More...

unsigned	Ndof
	Store number of degrees of freedom in the original problem. More...

unsigned	N_element_in_period
	Storage for number of elements in the period. More...

unsigned	N_tstorage
	Storage for the number of unknown time values. More...

double	Omega
	Storage for the frequency of the orbit (scaled by 2pi) More...

Detailed Description

template<unsigned NNODE_1D>
class oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >

A class that is used to assemble and solve the augmented system of equations associated with calculating periodic orbits directly

Definition at line 567 of file periodic_orbit_handler.h.

Constructor & Destructor Documentation

◆ PeriodicOrbitAssemblyHandler()

template<unsigned NNODE_1D>

oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::PeriodicOrbitAssemblyHandler	(	Problem *const &	problem_pt,
		const unsigned &	n_element_in_period,
		const DenseMatrix< double > &	initial_guess,
		const double &	omega
	)

inline

Constructor, initialises values and constructs mesh of elements.

Definition at line 711 of file periodic_orbit_handler.h.

References e, i, oomph::Data::is_pinned(), oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::ndof(), oomph::Data::ntstorage(), oomph::Data::nvalue(), oomph::PeriodicOrbitEquations::omega_pt(), oomph::oomph_info, oomph::PeriodicOrbitTimeDiscretisation::PeriodicOrbitTimeDiscretisation(), oomph::PeriodicOrbitEquations::set_ntstorage(), oomph::Data::set_time_stepper(), oomph::Data::set_value(), t, oomph::PeriodicOrbitEquations::time_pt(), oomph::GeomObject::time_stepper_pt(), oomph::Data::value(), and oomph::Node::value().

◆ ~PeriodicOrbitAssemblyHandler()

template<unsigned NNODE_1D>

oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::~PeriodicOrbitAssemblyHandler ( )

inline

Destructor, destroy the time mesh.

Definition at line 1430 of file periodic_orbit_handler.h.

Member Function Documentation

◆ adapt_temporal_mesh()

template<unsigned NNODE_1D>

void oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::adapt_temporal_mesh ( )

inline

Adapt the time mesh.

Definition at line 1022 of file periodic_orbit_handler.h.

References oomph::Problem::assign_eqn_numbers(), oomph::LinearAlgebraDistribution::build(), oomph::Problem::communicator_pt(), oomph::Mesh::Default_TimeStepper, oomph::Problem::Dof_distribution_pt, oomph::Problem::Dof_pt, e, oomph::Mesh::element_pt(), oomph::Problem::global_data_pt(), i, oomph::GeneralisedElement::internal_data_pt(), oomph::Problem::mesh_pt(), oomph::Mesh::nelement(), oomph::Problem::nglobal_data(), oomph::GeneralisedElement::ninternal_data(), oomph::Mesh::nnode(), oomph::Mesh::node_pt(), oomph::Problem::nsub_mesh(), oomph::Data::nvalue(), oomph::PeriodicOrbitEquations::omega_pt(), oomph::oomph_info, oomph::PeriodicOrbitTimeDiscretisation::PeriodicOrbitTimeDiscretisation(), oomph::PeriodicOrbitEquations::set_ntstorage(), oomph::Data::set_time_stepper(), oomph::Data::set_value(), t, oomph::PeriodicOrbitEquations::time_pt(), oomph::Problem::time_pt(), oomph::GeomObject::time_stepper_pt(), oomph::Data::value(), oomph::Node::value(), and oomph::SolidNode::variable_position_pt().

◆ discrete_times()

template<unsigned NNODE_1D>

void oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::discrete_times ( Vector< double > & t )

inline

Tell me the times at which you want the solution.

Definition at line 1014 of file periodic_orbit_handler.h.

◆ eqn_number()

template<unsigned NNODE_1D>

unsigned long oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::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.

Reimplemented from oomph::AssemblyHandler.

Definition at line 903 of file periodic_orbit_handler.h.

References oomph::GeneralisedElement::eqn_number(), oomph::GeneralisedElement::ndof(), and t.

◆ get_dofs_for_element()

template<unsigned NNODE_1D>

void oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::get_dofs_for_element	(	GeneralisedElement *const	elem_pt,
		Vector< double > &	dofs
	)

inlinevirtual

Implements oomph::PeriodicOrbitAssemblyHandlerBase.

Definition at line 930 of file periodic_orbit_handler.h.

References oomph::Problem::dof(), and i.

◆ get_jacobian()

template<unsigned NNODE_1D>

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

inlinevirtual

Calculate the elemental Jacobian matrix "d equation / d variable" for elem_pt.

Reimplemented from oomph::AssemblyHandler.

Definition at line 975 of file periodic_orbit_handler.h.

◆ get_previous_dofs_for_element()

template<unsigned NNODE_1D>

void oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::get_previous_dofs_for_element	(	GeneralisedElement *const	elem_pt,
		Vector< double > &	dofs
	)

inlinevirtual

Implements oomph::PeriodicOrbitAssemblyHandlerBase.

Definition at line 944 of file periodic_orbit_handler.h.

References i.

◆ get_residuals()

template<unsigned NNODE_1D>

void oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::get_residuals	(	GeneralisedElement *const &	elem_pt,
		Vector< double > &	residuals
	)

inlinevirtual

Return the contribution to the residuals of the element elem_pt.

Reimplemented from oomph::AssemblyHandler.

Definition at line 923 of file periodic_orbit_handler.h.

◆ ndof()

template<unsigned NNODE_1D>

unsigned oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::ndof ( GeneralisedElement *const & elem_pt )

inlinevirtual

Return the number of degrees of freedom in the element elem_pt.

Reimplemented from oomph::AssemblyHandler.

Definition at line 898 of file periodic_orbit_handler.h.

References oomph::GeneralisedElement::ndof().

Referenced by oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::PeriodicOrbitAssemblyHandler().

◆ orbit_output()

template<unsigned NNODE_1D>

void oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::orbit_output	(	std::ostream &	outfile,
		const unsigned &	n_plot
	)

inline

Calculate all desired vectors and matrices provided by the element elem_pt.

Return an unsigned integer to indicate whether the handler is a bifurcation tracking handler. The default is zero (not) Return a pointer to the bifurcation parameter in bifurcation tracking problems Return the eigenfunction(s) associated with the bifurcation that has been detected in bifurcation tracking problems Return the contribution to the residuals of the element elem_pt

Definition at line 1003 of file periodic_orbit_handler.h.

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

◆ set_dofs_for_element()

template<unsigned NNODE_1D>

void oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::set_dofs_for_element	(	GeneralisedElement *const	elem_pt,
		Vector< double > const &	dofs
	)

inlinevirtual

Implements oomph::PeriodicOrbitAssemblyHandlerBase.

Definition at line 959 of file periodic_orbit_handler.h.

References oomph::Problem::dof(), and i.

◆ set_previous_dofs_to_current_dofs()

template<unsigned NNODE_1D>

void oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::set_previous_dofs_to_current_dofs ( )

inline

Update the previous dofs.

Definition at line 889 of file periodic_orbit_handler.h.

References oomph::Problem::dof().

Member Data Documentation

◆ Basic_time_mesh_pt

template<unsigned NNODE_1D>

Mesh* oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::Basic_time_mesh_pt

private

Storage for the mesh of temporal elements with a simple mesh pointer.

Definition at line 691 of file periodic_orbit_handler.h.

◆ N_element_in_period

template<unsigned NNODE_1D>

unsigned oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::N_element_in_period

private

Storage for number of elements in the period.

Definition at line 700 of file periodic_orbit_handler.h.

◆ N_tstorage

template<unsigned NNODE_1D>

unsigned oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::N_tstorage

private

Storage for the number of unknown time values.

Definition at line 703 of file periodic_orbit_handler.h.

◆ Ndof

template<unsigned NNODE_1D>

unsigned oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::Ndof

private

Store number of degrees of freedom in the original problem.

Definition at line 697 of file periodic_orbit_handler.h.

◆ Omega

template<unsigned NNODE_1D>

double oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::Omega

private

Storage for the frequency of the orbit (scaled by 2pi)

Definition at line 706 of file periodic_orbit_handler.h.

◆ Previous_dofs

template<unsigned NNODE_1D>

Vector<double> oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::Previous_dofs

private

Storage for the previous solution.

Definition at line 694 of file periodic_orbit_handler.h.

◆ Problem_pt

template<unsigned NNODE_1D>

Problem* oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::Problem_pt

private

Pointer to the problem.

Definition at line 684 of file periodic_orbit_handler.h.

◆ Time_mesh_pt

template<unsigned NNODE_1D>

PeriodicOrbitTemporalMesh<SpectralPeriodicOrbitElement<NNODE_1D> >* oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::Time_mesh_pt

private

Storage for mesh of temporal elements.

Definition at line 688 of file periodic_orbit_handler.h.

◆ Time_stepper_pt

template<unsigned NNODE_1D>

PeriodicOrbitTimeDiscretisation* oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >::Time_stepper_pt

private

Pointer to the timestepper.

Definition at line 681 of file periodic_orbit_handler.h.

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

periodic_orbit_handler.h

Public Member Functions

Private Attributes

Detailed Description

template<unsigned NNODE_1D> class oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >

Constructor & Destructor Documentation

◆ PeriodicOrbitAssemblyHandler()

◆ ~PeriodicOrbitAssemblyHandler()

Member Function Documentation

◆ adapt_temporal_mesh()

◆ discrete_times()

◆ eqn_number()

◆ get_dofs_for_element()

◆ get_jacobian()

◆ get_previous_dofs_for_element()

◆ get_residuals()

◆ ndof()

◆ orbit_output()

◆ set_dofs_for_element()

◆ set_previous_dofs_to_current_dofs()

Member Data Documentation

◆ Basic_time_mesh_pt

◆ N_element_in_period

◆ N_tstorage

◆ Ndof

◆ Omega

◆ Previous_dofs

◆ Problem_pt

◆ Time_mesh_pt

◆ Time_stepper_pt

template<unsigned NNODE_1D>
class oomph::PeriodicOrbitAssemblyHandler< NNODE_1D >