ODEBND
Real-valued (non-validated) integration of parametric ODEs

Classes

struct  mc::ODESLV_GSL< T >::Options
 Integrator options. More...
class  mc::ODESLV_GSL< T >::Exceptions
 Structure for setting up storing the solver exceptions. More...
struct  mc::ODESLV_GSL< T >::Results
 Store integration bounds at a given time instant. More...

Functions

 mc::ODESLV_GSL< T >::ODESLV_GSL ()
 Default class constructor.
virtual mc::ODESLV_GSL< T >::~ODESLV_GSL ()
 Default destructor.
STATUS mc::ODESLV_GSL< T >::states (const unsigned int ns, const double *tk, const double *p, double **xk, double *f, std::ostream &os=std::cout)
 Integrate trajectory of parametric ODEs.
STATUS mc::ODESLV_GSL< T >::states_ASA (const unsigned int ns, const double *tk, const double *p, double **xk, double *f, double **lk, double *df, std::ostream &os=std::cout)
 Integrate trajectory of parametric ODEs and adjoint ODEs.
STATUS mc::ODESLV_GSL< T >::bounds (const unsigned int ns, const double *tk, const T *Ip, T **Ixk, T *If, const unsigned int nsamp, std::ostream &os=std::cout)
 Compute approximate interval enclosure of reachable set of parametric ODEs using parameter sampling.
STATUS mc::ODESLV_GSL< T >::bounds_ASA (const unsigned int ns, const double *tk, const T *Ip, T **Ixk, T *If, T **Ilk, T *Idf, const unsigned int nsamp, std::ostream &os=std::cout)
 Compute approximate interval enclosure of reachable set of parametric ODEs and adjoint ODEs using parameter sampling.
void mc::ODESLV_GSL< T >::record (std::ofstream &bndrec, const unsigned int iprec=5) const
 Record results in file bndrec, with accuracy of iprec digits.

Variables

Stats mc::ODESLV_GSL< T >::stats_traj
 Statistics for state integration.
Stats mc::ODESLV_GSL< T >::stats_adj
 Statistics for adjoint integration.

Function Documentation

template<typename T >
ODESLV_GSL< T >::STATUS mc::ODESLV_GSL< T >::bounds ( const unsigned int  ns,
const double *  tk,
const T *  Ip,
T **  Ixk,
T *  f,
const unsigned int  nsamp,
std::ostream &  os = std::cout 
) [inline]

This function computes an approximate interval enclosure of the reachable set of the parametric ODEs defined in IVP using equally spaced samples:

  • ns [input] number of time stages
  • tk [input] stage times, including the initial time
  • Ip [input] interval parameter set
  • Ixk [output] approximate interval state enclosures at stage times
  • If [output] approximate function enclosure (only if If != 0)
  • nsamp [input] number of samples for each parameter
  • os [input] output stream

The return value is the status.

template<typename T >
ODESLV_GSL< T >::STATUS mc::ODESLV_GSL< T >::bounds_ASA ( const unsigned int  ns,
const double *  tk,
const T *  Ip,
T **  Ixk,
T *  If,
T **  Ilk,
T *  df,
const unsigned int  nsamp,
std::ostream &  os = std::cout 
) [inline]

This function computes an approximate interval enclosure of the reachable set of the parametric ODEs defined in IVP using equally spaced samples:

  • ns [input] number of time stages
  • tk [input] stage times, including the initial time
  • Ip [input] interval parameter set
  • Ixk [output] approximate interval state enclosures at stage times
  • If [output] approximate function enclosure (only if If != 0)
  • Ilk [output] approximate interval adjoint enclosures at stage times
  • Idf [output] approximate function derivative enclosures (stored row-wise)
  • nsamp [input] number of samples for each parameter
  • os [input] output stream

The return value is the status.

template<typename T >
ODESLV_GSL< T >::STATUS mc::ODESLV_GSL< T >::states ( const unsigned int  ns,
const double *  tk,
const double *  p,
double **  xk,
double *  f,
std::ostream &  os = std::cout 
) [inline]

This function computes the solution of the parametric ODEs defined in IVP:

  • ns [input] number of time stages
  • tk [input] stage times, including the initial time
  • p [input] parameter values
  • xk [output] state values at stage times
  • f [output] function values (only if f != 0)
  • os [input] output stream

The return value is the status.

template<typename T >
ODESLV_GSL< T >::STATUS mc::ODESLV_GSL< T >::states_ASA ( const unsigned int  ns,
const double *  tk,
const double *  p,
double **  xk,
double *  f,
double **  lk,
double *  df,
std::ostream &  os = std::cout 
) [inline]

This function computes the solution of the parametric ODEs defined in IVP:

  • ns [input] number of time stages
  • tk [input] stage times, including the initial time
  • p [input] parameter values
  • xk [output] state values at stage times
  • f [output] function values (only if f != 0)
  • lk [output] adjoint values at stage times
  • f [output] function derivatives (stored row-wise)
  • os [input] output stream

The return value is the status.