Define a residual function using the boundary conditions for the biharmonic_equation. More...

Inheritance diagram for CppNoddy::Example::Biharmonic_residual:

Public Member Functions
	Biharmonic_residual ()

	~Biharmonic_residual ()

void	residual_fn (const DenseVector< D_complex > &unknown, DenseVector< D_complex > &BC) const
	implement a residual function. More...

Public Member Functions inherited from CppNoddy::Residual< D_complex >
	Residual (const unsigned &order)
	Constructor for a 'square' residual object that is, N residuals for N unknowns. More...

	Residual (const unsigned &order, const unsigned &nvars)
	Constructor for a 'non-square' residual object that is, there are less residual constraints than unknowns. More...

virtual	~Residual ()
	An empty destructor, virtual since we have virtual methods. More...

void	update (const DenseVector< D_complex > &state)
	Update the Residual object for the current set of state variables. More...

const DenseVector< D_complex > &	residual () const
	Return a handle to the residuals corresponding to the last update state. More...

const DenseMatrix< D_complex > &	jacobian () const
	Retrun a handle to the Jacobian of the residual corresponding to the last update state. More...

D_complex &	delta ()

const D_complex &	delta () const

unsigned	get_order () const
	Get the order of the residual vector. More...

unsigned	get_number_of_vars () const
	Get the number of variables that this residual condition is defined for. More...

virtual void	residual_fn (const DenseVector< D_complex > &state, DenseVector< D_complex > &f) const
	A blank virtual residual function method. More...

Public Attributes
ODE_IVP< D_complex > *	ode

Additional Inherited Members
Protected Member Functions inherited from CppNoddy::Residual< D_complex >
virtual void	jacobian (const DenseVector< D_complex > &state, DenseMatrix< D_complex > &jac) const
	Because the residual evaluation at the current state is assumed to have already been done by the 'update' method, this routine is protected. More...

Protected Attributes inherited from CppNoddy::Residual< D_complex >
DenseMatrix< D_complex >	JAC_AT_LAST_STATE
	Jacobian for the last state vector. More...

DenseVector< D_complex >	FN_AT_LAST_STATE
	Residual for the last state vector. More...

DenseVector< D_complex >	LAST_STATE
	The last state vector. More...

D_complex	DELTA
	A default step for FD computation of the Jacobian. More...

unsigned	ORDER_OF_SYSTEM
	The order of the system of equations. More...

unsigned	NUMBER_OF_VARS
	The number of elements in the state vector. More...

Detailed Description

Define a residual function using the boundary conditions for the biharmonic_equation.

Definition at line 44 of file EVPShootBiharmonic.cpp.

Constructor & Destructor Documentation

◆ Biharmonic_residual()

CppNoddy::Example::Biharmonic_residual::Biharmonic_residual ( )

inline

Definition at line 49 of file EVPShootBiharmonic.cpp.

                            : Residual<D_complex> ( 2 )
      {
        // instantiate the 1D bi-harmonic eqn
        eqn = new Biharmonic_equation;
        // set up the ODE domain, step etc
        ode = new ODE_IVP<D_complex>( eqn, 0.0, 1.0, 1000 );
      }

References ode.

◆ ~Biharmonic_residual()

CppNoddy::Example::Biharmonic_residual::~Biharmonic_residual ( )

inline

Definition at line 57 of file EVPShootBiharmonic.cpp.

      {
        delete ode;
        delete eqn;
      }

References ode.

Member Function Documentation

◆ residual_fn()

void CppNoddy::Example::Biharmonic_residual::residual_fn	(	const DenseVector< D_complex > &	unknown,
		DenseVector< D_complex > &	BC
	)		const

inlinevirtual

implement a residual function.

Reimplemented from CppNoddy::Residual< D_complex >.

Definition at line 64 of file EVPShootBiharmonic.cpp.

      {
        DenseVector<D_complex> u( 4, 0.0 );        // 4th order problem
        u[ 0 ] = 0.0;                 // f = 0
        u[ 1 ] = 0.0;                 // f' = 0
        u[ 2 ] = 1.0;                 // f'' = 1
        u[ 3 ] = unknown[ 0 ];        // f''' = unknown[0]
        eqn -> lambda = unknown[ 1 ]; // parameter = unknown[1]
        // shoot using R-K fixed step size algorithm
        u = ode -> shoot4( u );
        // return the BC residuals
        BC[ 0 ] = u[ 0 ];
        BC[ 1 ] = u[ 1 ]; // trying to make f = f' = 0 at other boundary
      }