Wrapper class for 'cvodes' integrator from LLNL. More...

#include <CVodesIntegrator.h>

Inheritance diagram for CVodesIntegrator:

Collaboration diagram for CVodesIntegrator:

Public Member Functions
	CVodesIntegrator ()
	Constructor. More...

virtual void	setTolerances (double reltol, size_t n, double *abstol)
	Set error tolerances. More...

virtual void	setTolerances (double reltol, double abstol)
	Set error tolerances. More...

virtual void	setSensitivityTolerances (double reltol, double abstol)
	Set the sensitivity error tolerances. More...

virtual void	setProblemType (int probtype)
	Set the problem type. More...

virtual void	initialize (double t0, FuncEval &func)
	Initialize the integrator for a new problem. More...

virtual void	reinitialize (double t0, FuncEval &func)

virtual void	integrate (double tout)
	Integrate the system of equations. More...

virtual doublereal	step (double tout)
	Integrate the system of equations. More...

virtual double &	solution (size_t k)
	The current value of the solution of equation k. More...

virtual double *	solution ()
	The current value of the solution of the system of equations. More...

virtual int	nEquations () const
	The number of equations. More...

virtual int	nEvals () const
	The number of function evaluations. More...

virtual void	setMaxOrder (int n)
	Set the maximum integration order that will be used. More...

virtual void	setMethod (MethodType t)
	Set the solution method. More...

virtual void	setIterator (IterType t)
	Set the linear iterator. More...

virtual void	setMaxStepSize (double hmax)
	Set the maximum step size. More...

virtual void	setMinStepSize (double hmin)
	Set the minimum step size. More...

virtual void	setMaxSteps (int nmax)

virtual void	setMaxErrTestFails (int n)
	Set the maximum permissible number of error test failures. More...

virtual void	setBandwidth (int N_Upper, int N_Lower)

virtual int	nSensParams ()

virtual double	sensitivity (size_t k, size_t p)

virtual std::string	getErrorInfo (int N)
	Returns a string listing the weighted error estimates associated with each solution component. More...

Public Member Functions inherited from Integrator
	Integrator ()
	Default Constructor. More...

virtual	~Integrator ()
	Destructor. More...

Public Attributes
std::string	m_error_message
	Error message information provide by CVodes. More...

Protected Member Functions
void	applyOptions ()
	Applies user-specified options to the underlying CVODES solver. More...

Private Member Functions
void	sensInit (double t0, FuncEval &func)

Private Attributes
size_t	m_neq

void *	m_cvode_mem

void *	m_linsol
	Sundials linear solver object. More...

void *	m_linsol_matrix
	matrix used by Sundials More...

FuncEval *	m_func

double	m_t0

double	m_time
	The current integrator time. More...

N_Vector	m_y

N_Vector	m_abstol

int	m_type

int	m_itol

int	m_method

int	m_iter

int	m_maxord

double	m_reltol

double	m_abstols

double	m_reltolsens

double	m_abstolsens

size_t	m_nabs

double	m_hmax

double	m_hmin

int	m_maxsteps

int	m_maxErrTestFails

N_Vector *	m_yS

size_t	m_np

int	m_mupper

int	m_mlower

bool	m_sens_ok
	Indicates whether the sensitivities stored in m_yS have been updated for at the current integrator time. More...

Detailed Description

Wrapper class for 'cvodes' integrator from LLNL.

See also: FuncEval.h. Classes that use CVodesIntegrator: ImplicitChem, ImplicitSurfChem, Reactor

Definition at line 25 of file CVodesIntegrator.h.

Constructor & Destructor Documentation

◆ CVodesIntegrator()

CVodesIntegrator ( )

Constructor.

Default settings: dense Jacobian, no user-supplied Jacobian function, Newton iteration.

Definition at line 79 of file CVodesIntegrator.cpp.

Member Function Documentation

◆ setTolerances() [1/2]

void setTolerances	(	double	reltol,
		size_t	n,
		double *	abstol
	)

virtual

Set error tolerances.

Parameters

reltol	scalar relative tolerance
n	Number of equations
abstol	array of N absolute tolerance values

Reimplemented from Integrator.

Definition at line 144 of file CVodesIntegrator.cpp.

◆ setTolerances() [2/2]

void setTolerances	(	double	reltol,
		double	abstol
	)

virtual

Set error tolerances.

Parameters

reltol	scalar relative tolerance
abstol	scalar absolute tolerance

Reimplemented from Integrator.

Definition at line 160 of file CVodesIntegrator.cpp.

◆ setSensitivityTolerances()

void setSensitivityTolerances	(	double	reltol,
		double	abstol
	)

virtual

Set the sensitivity error tolerances.

Parameters

reltol	scalar relative tolerance
abstol	scalar absolute tolerance

Reimplemented from Integrator.

Definition at line 167 of file CVodesIntegrator.cpp.

◆ setProblemType()

void setProblemType ( int probtype )

virtual

Set the problem type.

Parameters

probtype Type of the problem

Reimplemented from Integrator.

Definition at line 173 of file CVodesIntegrator.cpp.

◆ initialize()

void initialize	(	double	t0,
		FuncEval &	func
	)

virtual

Initialize the integrator for a new problem.

Call after all options have been set.

Parameters

t0	initial time
func	RHS evaluator object for system of equations.

Reimplemented from Integrator.

Definition at line 259 of file CVodesIntegrator.cpp.

References FuncEval::clearErrors(), CVodesIntegrator::m_time, and FuncEval::neq().

◆ integrate()

void integrate ( double tout )

virtual

Integrate the system of equations.

Parameters

tout	Integrate to this time. Note that this is the absolute time value, not a time interval.

Reimplemented from Integrator.

Definition at line 434 of file CVodesIntegrator.cpp.

References CVodesIntegrator::getErrorInfo(), FuncEval::getErrors(), CVodesIntegrator::m_error_message, CVodesIntegrator::m_sens_ok, and CVodesIntegrator::m_time.

◆ step()

double step ( double tout )

virtual

Integrate the system of equations.

Parameters

tout	integrate to this time. Note that this is the absolute time value, not a time interval.

Reimplemented from Integrator.

Definition at line 454 of file CVodesIntegrator.cpp.

References CVodesIntegrator::getErrorInfo(), FuncEval::getErrors(), CVodesIntegrator::m_error_message, CVodesIntegrator::m_sens_ok, and CVodesIntegrator::m_time.

◆ solution() [1/2]

double & solution ( size_t k )

virtual

The current value of the solution of equation k.

Reimplemented from Integrator.

Definition at line 134 of file CVodesIntegrator.cpp.

◆ solution() [2/2]

double * solution ( )

virtual

The current value of the solution of the system of equations.

Reimplemented from Integrator.

Definition at line 139 of file CVodesIntegrator.cpp.

◆ nEquations()

virtual int nEquations ( ) const

inlinevirtual

The number of equations.

Reimplemented from Integrator.

Definition at line 44 of file CVodesIntegrator.h.

◆ nEvals()

int nEvals ( ) const

virtual

The number of function evaluations.

Reimplemented from Integrator.

Definition at line 473 of file CVodesIntegrator.cpp.

◆ setMaxOrder()

virtual void setMaxOrder ( int n )

inlinevirtual

Set the maximum integration order that will be used.

Reimplemented from Integrator.

Definition at line 48 of file CVodesIntegrator.h.

◆ setMethod()

void setMethod ( MethodType t )

virtual

Set the solution method.

Reimplemented from Integrator.

Definition at line 178 of file CVodesIntegrator.cpp.

References Cantera::Adams_Method, and Cantera::BDF_Method.

◆ setIterator()

void setIterator ( IterType t )

virtual

Set the linear iterator.

Reimplemented from Integrator.

Definition at line 221 of file CVodesIntegrator.cpp.

References Cantera::Functional_Iter, and Cantera::Newton_Iter.

◆ setMaxStepSize()

void setMaxStepSize ( double hmax )

virtual

Set the maximum step size.

Reimplemented from Integrator.

Definition at line 189 of file CVodesIntegrator.cpp.

◆ setMinStepSize()

void setMinStepSize ( double hmin )

virtual

Set the minimum step size.

Reimplemented from Integrator.

Definition at line 197 of file CVodesIntegrator.cpp.

◆ setMaxErrTestFails()

void setMaxErrTestFails ( int n )

virtual

Set the maximum permissible number of error test failures.

Reimplemented from Integrator.

Definition at line 213 of file CVodesIntegrator.cpp.

◆ getErrorInfo()

string getErrorInfo ( int N )

virtual

Returns a string listing the weighted error estimates associated with each solution component.

This information can be used to identify which variables are responsible for integrator failures or unexpected small timesteps.

Definition at line 506 of file CVodesIntegrator.cpp.

Referenced by CVodesIntegrator::integrate(), and CVodesIntegrator::step().

◆ applyOptions()

void applyOptions ( )

protected

Applies user-specified options to the underlying CVODES solver.

Called during integrator initialization or reinitialization.

Definition at line 359 of file CVodesIntegrator.cpp.

References CVodesIntegrator::m_linsol, and CVodesIntegrator::m_linsol_matrix.

Member Data Documentation

◆ m_error_message

std::string m_error_message

Error message information provide by CVodes.

Definition at line 73 of file CVodesIntegrator.h.

Referenced by Cantera::cvodes_err(), CVodesIntegrator::integrate(), and CVodesIntegrator::step().

◆ m_linsol

void* m_linsol

private

Sundials linear solver object.

Definition at line 85 of file CVodesIntegrator.h.

Referenced by CVodesIntegrator::applyOptions().

◆ m_linsol_matrix

void* m_linsol_matrix

private

matrix used by Sundials

Definition at line 86 of file CVodesIntegrator.h.

Referenced by CVodesIntegrator::applyOptions().

◆ m_time

double m_time

private

The current integrator time.

Definition at line 89 of file CVodesIntegrator.h.

Referenced by CVodesIntegrator::initialize(), CVodesIntegrator::integrate(), and CVodesIntegrator::step().

◆ m_sens_ok

bool m_sens_ok

private

Indicates whether the sensitivities stored in m_yS have been updated for at the current integrator time.

Definition at line 109 of file CVodesIntegrator.h.

Referenced by CVodesIntegrator::integrate(), and CVodesIntegrator::step().

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

Public Member Functions

Public Attributes

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ CVodesIntegrator()

Member Function Documentation

◆ setTolerances() [1/2]

◆ setTolerances() [2/2]

◆ setSensitivityTolerances()

◆ setProblemType()

◆ initialize()

◆ integrate()

◆ step()

◆ solution() [1/2]

◆ solution() [2/2]

◆ nEquations()

◆ nEvals()

◆ setMaxOrder()

◆ setMethod()

◆ setIterator()

◆ setMaxStepSize()

◆ setMinStepSize()

◆ setMaxErrTestFails()

◆ getErrorInfo()

◆ applyOptions()

Member Data Documentation

◆ m_error_message

◆ m_linsol

◆ m_linsol_matrix

◆ m_time

◆ m_sens_ok