Define a residual function using the boundary conditions for the Blasius profile. More...

Inheritance diagram for CppNoddy::Example::FS_residual:

Public Member Functions
	FS_residual ()

	~FS_residual ()

void	residual_fn (const DenseVector< double > &z, DenseVector< double > &f) 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

FS_eqn *	eqn

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 Blasius profile.

Definition at line 46 of file ArcShootFalknerSkan.cpp.

Constructor & Destructor Documentation

◆ FS_residual()

CppNoddy::Example::FS_residual::FS_residual ( )

inline

Definition at line 52 of file ArcShootFalknerSkan.cpp.

                    : Residual<double> ( 1 )
      {
        eqn = new FS_eqn;
        eqn -> beta = 0.0; // init the Hartree parameter
        ode = new ODE_IVP<double> ( eqn, 0.0, 20.0, 2000 );
      }

References eqn, and ode.

◆ ~FS_residual()

CppNoddy::Example::FS_residual::~FS_residual ( )

inline

Definition at line 59 of file ArcShootFalknerSkan.cpp.

      {
        delete ode;
        delete eqn;
      }

References eqn, and ode.

Member Function Documentation

◆ residual_fn()

void CppNoddy::Example::FS_residual::residual_fn	(	const DenseVector< double > &	z,
		DenseVector< double > &	f
	)		const

inlinevirtual

implement a residual function.

Reimplemented from CppNoddy::Residual< double >.

Definition at line 66 of file ArcShootFalknerSkan.cpp.

      {
        DenseVector<double> u( 3, 0.0 ); // 3 variables of FS problem
        u[ 0 ] = 0.0;          // f = 0
        u[ 1 ] = 0.0;          // f' = 0
        u[ 2 ] = z[ 0 ];       // f'' = unknown
        // shoot the ODE
        u = ode -> shoot4( u );
        // return the residual of f'(infty) = 1 condition
        f[ 0 ] = u[ 1 ] - 1.;
      }