Define a residual function using the boundary conditions for the Berman equation. More...

Inheritance diagram for CppNoddy::Example::Berman_residual:

Public Member Functions
	Berman_residual ()

	~Berman_residual ()

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

Public Member Functions inherited from CppNoddy::Residual< double >
	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< double > &state)
	Update the Residual object for the current set of state variables. More...

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

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

double &	delta ()

const double &	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< double > &state, DenseVector< double > &f) const
	A blank virtual residual function method. More...

Public Attributes
ODE_IVP< double > *	ode

Additional Inherited Members
Protected Member Functions inherited from CppNoddy::Residual< double >
virtual void	jacobian (const DenseVector< double > &state, DenseMatrix< double > &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< double >
DenseMatrix< double >	JAC_AT_LAST_STATE
	Jacobian for the last state vector. More...

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

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

double	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 Berman equation.

Definition at line 52 of file BermanShooting.cpp.

Constructor & Destructor Documentation

◆ Berman_residual()

CppNoddy::Example::Berman_residual::Berman_residual ( )

inline

Definition at line 58 of file BermanShooting.cpp.

                        : Residual<double> ( 2 )
      {
        // instantiate the Berman equation as our problem
        eqn = new Berman_equation;
        // set the Reynolds number
        eqn -> Re = 1.0;
        // instantiate an ODE in [-1.,1] with a maximum of
        // 100 integrator steps, using this problem equation
        ode = new ODE_IVP<double>( eqn, -1.0, 1.0, 500 );
 
      }

References ode, and CppNoddy::Example::Re.

◆ ~Berman_residual()

CppNoddy::Example::Berman_residual::~Berman_residual ( )

inline

Definition at line 70 of file BermanShooting.cpp.

      {
        delete ode;
        delete eqn;
      }

References ode.

Member Function Documentation

◆ residual_fn()

void CppNoddy::Example::Berman_residual::residual_fn	(	const DenseVector< double > &	unknown,
		DenseVector< double > &	BC
	)		const

inlinevirtual

implement a residual function.

Reimplemented from CppNoddy::Residual< double >.

Definition at line 77 of file BermanShooting.cpp.

      {
        DenseVector<double> u( 4, 0.0 ); // 4 variables of Berman problem
 
        u[ 0 ] = -1.0;         // V = 0
        u[ 1 ] = 0.0;          // V' = 0
        u[ 2 ] = unknown[ 0 ]; // V'' = unknown0
        u[ 3 ] = unknown[ 1 ]; // A = unknown1
 
        u = ode -> shoot4( u );
 
        BC[ 0 ] = u[ 0 ] - 1.;
        BC[ 1 ] = u[ 1 ];
      }