Virtual base class for ODE/DAE right-hand-side function evaluators. More...

#include <FuncEval.h>

Inheritance diagram for FuncEval:

Detailed Description

Virtual base class for ODE/DAE right-hand-side function evaluators.

Classes derived from FuncEval evaluate the right-hand-side function \( \vec{F}(t,\vec{y}) \) in

\[ \dot{\vec{y}} = \vec{F}(t,\vec{y}). \]

Definition at line 31 of file FuncEval.h.

Public Member Functions
virtual void	eval (double t, double y, double ydot, double *p)
	Evaluate the right-hand-side ODE function. More...

virtual void	evalDae (double t, double y, double ydot, double p, double residual)
	Evaluate the right-hand-side DAE function. More...

virtual void	getConstraints (double *constraints)
	Given a vector of length neq(), mark which variables should be considered algebraic constraints. More...

int	evalNoThrow (double t, double y, double ydot)
	Evaluate the right-hand side using return code to indicate status. More...

int	evalDaeNoThrow (double t, double y, double ydot, double *residual)
	Evaluate the right-hand side using return code to indicate status. More...

virtual void	preconditionerSetup (double t, double *y, double gamma)
	Evaluate the setup processes for the Jacobian preconditioner. More...

virtual void	preconditionerSolve (double rhs, double output)
	Evaluate the linear system Ax=b where A is the preconditioner. More...

virtual void	updatePreconditioner (double gamma)
	Update the preconditioner based on already computed jacobian values. More...

int	preconditioner_setup_nothrow (double t, double *y, double gamma)
	Preconditioner setup that doesn't throw an error but returns a CVODES flag. More...

int	preconditioner_solve_nothrow (double rhs, double output)
	Preconditioner solve that doesn't throw an error but returns a CVODES flag. More...

virtual void	getState (double *y)
	Fill in the vector y with the current state of the system. More...

virtual void	getStateDae (double y, double ydot)
	Fill in the vectors y and ydot with the current state of the system. More...

virtual size_t	neq () const =0
	Number of equations. More...

virtual size_t	nparams () const
	Number of sensitivity parameters. More...

void	suppressErrors (bool suppress)
	Enable or disable suppression of errors when calling eval() More...

bool	suppressErrors () const
	Get current state of error suppression. More...

string	getErrors () const
	Return a string containing the text of any suppressed errors. More...

void	clearErrors ()
	Clear any previously-stored suppressed errors. More...

Public Attributes
vector< double >	m_sens_params
	Values for the problem parameters for which sensitivities are computed This is the array which is perturbed and passed back as the fourth argument to eval(). More...

vector< double >	m_paramScales
	Scaling factors for each sensitivity parameter. More...

Protected Attributes
bool	m_suppress_errors = false

vector< string >	m_errors
	Errors occurring during function evaluations. More...

Member Function Documentation

◆ eval()

virtual void eval	(	double	t,
		double *	y,
		double *	ydot,
		double *	p
	)

inlinevirtual

Evaluate the right-hand-side ODE function.

Called by the integrator.

Parameters

[in]	t	time.
[in]	y	solution vector, length neq()
[out]	ydot	rate of change of solution vector, length neq()
[in]	p	sensitivity parameter vector, length nparams()

Reimplemented in ReactorNet, and ImplicitSurfChem.

Definition at line 44 of file FuncEval.h.

◆ evalDae()

virtual void evalDae	(	double	t,
		double *	y,
		double *	ydot,
		double *	p,
		double *	residual
	)

inlinevirtual

Evaluate the right-hand-side DAE function.

Called by the integrator.

Parameters

[in]	t	time.
[in]	y	solution vector, length neq()
[out]	ydot	rate of change of solution vector, length neq()
[in]	p	sensitivity parameter vector, length nparams()
[out]	residual	the DAE residuals, length nparams()

Reimplemented in ReactorNet.

Definition at line 56 of file FuncEval.h.

◆ getConstraints()

virtual void getConstraints ( double * constraints )

inlinevirtual

Given a vector of length neq(), mark which variables should be considered algebraic constraints.

Reimplemented in ReactorNet.

Definition at line 63 of file FuncEval.h.

◆ evalNoThrow()

int evalNoThrow	(	double	t,
		double *	y,
		double *	ydot
	)

Evaluate the right-hand side using return code to indicate status.

Errors are indicated using the return value, rather than by throwing exceptions. This method is used when calling from a C-based integrator such as CVODES. Exceptions may either be stored or printed, based on the setting of suppressErrors().

Returns: 0 for a successful evaluation; 1 after a potentially- recoverable error; -1 after an unrecoverable error.

Definition at line 7 of file FuncEval.cpp.

◆ evalDaeNoThrow()

int evalDaeNoThrow	(	double	t,
		double *	y,
		double *	ydot,
		double *	residual
	)

Evaluate the right-hand side using return code to indicate status.

Errors are indicated using the return value, rather than by throwing exceptions. This method is used when calling from a C-based integrator such as CVODES. Exceptions may either be stored or printed, based on the setting of suppressErrors().

Returns: 0 for a successful evaluation; 1 after a potentially- recoverable error; -1 after an unrecoverable error.

Definition at line 39 of file FuncEval.cpp.

◆ preconditionerSetup()

virtual void preconditionerSetup	(	double	t,
		double *	y,
		double	gamma
	)

inlinevirtual

Evaluate the setup processes for the Jacobian preconditioner.

Parameters

[in]	t	time.
[in]	y	solution vector, length neq()
	gamma	the gamma in M=I-gamma*J

Warning: This function is an experimental part of the Cantera API and may be changed or removed without notice.

Reimplemented in ReactorNet.

Definition at line 97 of file FuncEval.h.

◆ preconditionerSolve()

virtual void preconditionerSolve	(	double *	rhs,
		double *	output
	)

inlinevirtual

Evaluate the linear system Ax=b where A is the preconditioner.

Parameters

[in]	rhs	right hand side vector used in linear system
[out]	output	output vector for solution

Warning: This function is an experimental part of the Cantera API and may be changed or removed without notice.

Reimplemented in ReactorNet.

Definition at line 108 of file FuncEval.h.

◆ updatePreconditioner()

virtual void updatePreconditioner ( double gamma )

inlinevirtual

Update the preconditioner based on already computed jacobian values.

Reimplemented in ReactorNet.

Definition at line 113 of file FuncEval.h.

◆ preconditioner_setup_nothrow()

int preconditioner_setup_nothrow	(	double	t,
		double *	y,
		double	gamma
	)

Preconditioner setup that doesn't throw an error but returns a CVODES flag.

It also helps as a first level of polymorphism which identifies the specific FuncEval, e.g., ReactorNet.

Parameters

[in]	t	time.
[in]	y	solution vector, length neq()
	gamma	the gamma in M=I-gamma*J

Warning: This function is an experimental part of the Cantera API and may be changed or removed without notice.

Definition at line 80 of file FuncEval.cpp.

◆ preconditioner_solve_nothrow()

int preconditioner_solve_nothrow	(	double *	rhs,
		double *	output
	)

Preconditioner solve that doesn't throw an error but returns a CVODES flag.

It also helps as a first level of polymorphism which identifies the specific FuncEval, e.g., ReactorNet.

Parameters

[in]	rhs	right hand side vector used in linear system
[out]	output	output vector for solution

Warning: This function is an experimental part of the Cantera API and may be changed or removed without notice.

Definition at line 113 of file FuncEval.cpp.

◆ getState()

virtual void getState ( double * y )

inlinevirtual

Fill in the vector y with the current state of the system.

Used for getting the initial state for ODE systems.

Reimplemented in ReactorNet, and ImplicitSurfChem.

Definition at line 142 of file FuncEval.h.

◆ getStateDae()

virtual void getStateDae	(	double *	y,
		double *	ydot
	)

inlinevirtual

Fill in the vectors y and ydot with the current state of the system.

Used for getting the initial state for DAE systems.

Reimplemented in ReactorNet.

Definition at line 148 of file FuncEval.h.

◆ neq()

virtual size_t neq ( ) const

pure virtual

Number of equations.

Implemented in ReactorNet, and ImplicitSurfChem.

◆ nparams()

virtual size_t nparams ( ) const

inlinevirtual

Number of sensitivity parameters.

Reimplemented in ReactorNet.

Definition at line 156 of file FuncEval.h.

◆ suppressErrors() [1/2]

void suppressErrors ( bool suppress )

inline

Enable or disable suppression of errors when calling eval()

Definition at line 161 of file FuncEval.h.

◆ suppressErrors() [2/2]

bool suppressErrors ( ) const

inline

Get current state of error suppression.

Definition at line 166 of file FuncEval.h.

◆ getErrors()

string getErrors ( ) const

Return a string containing the text of any suppressed errors.

Definition at line 71 of file FuncEval.cpp.

◆ clearErrors()

void clearErrors ( )

inline

Clear any previously-stored suppressed errors.

Definition at line 174 of file FuncEval.h.

Member Data Documentation

◆ m_sens_params

vector<double> m_sens_params

Values for the problem parameters for which sensitivities are computed This is the array which is perturbed and passed back as the fourth argument to eval().

Definition at line 181 of file FuncEval.h.

◆ m_paramScales

vector<double> m_paramScales

Scaling factors for each sensitivity parameter.

Definition at line 184 of file FuncEval.h.

◆ m_errors

vector<string> m_errors

protected

Errors occurring during function evaluations.

Definition at line 191 of file FuncEval.h.

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

Detailed Description

Public Member Functions

Public Attributes

Protected Attributes

Member Function Documentation

◆ eval()

◆ evalDae()

◆ getConstraints()

◆ evalNoThrow()

◆ evalDaeNoThrow()

◆ preconditionerSetup()

◆ preconditionerSolve()

◆ updatePreconditioner()

◆ preconditioner_setup_nothrow()

◆ preconditioner_solve_nothrow()

◆ getState()

◆ getStateDae()

◆ neq()

◆ nparams()

◆ suppressErrors() [1/2]

◆ suppressErrors() [2/2]

◆ getErrors()

◆ clearErrors()

Member Data Documentation

◆ m_sens_params

◆ m_paramScales

◆ m_errors