CVODES

CVODES is a solver for stiff and nonstiff ODE systems (initial value problem) given in explicit form y’ = f(t,y,p) with sensitivity analysis capabilities (both forward and adjoint modes). CVODES is a superset of CVODE and hence all options available to CVODE (with the exception of the FCVODE interface module) are also available for CVODES. Both integration methods (Adams-Moulton and BDF) and the corresponding nonlinear iteration methods, as well as all linear solver and preconditioner modules are available for the integration of the original ODEs, the sensitivity systems, or the adjoint system. Depending on the number of model parameters and the number of functional outputs, one of two sensitivity methods is more appropriate. The forward sensitivity analysis (FSA) method is mostly suitable when the gradients of many outputs (for example the entire solution vector) with respect to relatively few parameters are needed. In this approach, the model is differentiated with respect to each parameter in turn to yield an additional system of the same size as the original one, the result of which is the solution sensitivity. The gradient of any output function depending on the solution can then be directly obtained from these sensitivities by applying the chain rule of differentiation. The adjoint sensitivity analysis (ASA) method is more practical than the forward approach when the number of parameters is large and the gradients of only few output functionals are needed. In this approach, the solution sensitivities need not be computed explicitly. Instead, for each output functional of interest, an additional system, adjoint to the original one, is formed and solved. The solution of the adjoint system can then be used to evaluate the gradient of the output functional with respect to any set of model parameters. The FSA module in CVODES implements a simultaneous corrector method as well as two flavors of staggered corrector methods -- for the case when sensitivity right hand sides are generated all at once or separated for each model parameter. The ASA module provides the infrastructure required for the backward integration in time of systems of differential equations dependent on the solution of the original ODEs. It employs a checkpointing scheme for efficient interpolation of forward solutions during the backward integration.


References in zbMATH (referenced in 28 articles )

Showing results 1 to 20 of 28.
Sorted by year (citations)

1 2 next

  1. Gul, R.; Bernhard, S.: Parametric uncertainty and global sensitivity analysis in a model of the carotid bifurcation: identification and ranking of most sensitive model parameters (2015)
  2. Janka, Dennis: Sequential quadratic programming with indefinite Hessian approximations for nonlinear optimum experimental design for parameter estimation in differential-algebraic equations (2015)
  3. Quirynen, Rien; Vukov, Milan; Diehl, Moritz: Multiple shooting in a microsecond (2015)
  4. Schmitt, Bernhard A.: Peer methods with improved embedded sensitivities for parameter-dependent ODEs (2014)
  5. Schmitt, Bernhard A.: Reprint of: Peer methods with improved embedded sensitivities for parameter-dependent ODEs (2014)
  6. Zeng, X.; Anitescu, M.: Sequential Monte Carlo sampling in hidden Markov models of nonlinear dynamical systems (2014)
  7. Hug, S.; Raue, A.; Hasenauer, J.; Bachmann, J.; Klingmüller, U.; Timmer, J.; Theis, F.J.: High-dimensional Bayesian parameter estimation: case study for a model of JAK2/STAT5 signaling (2013)
  8. Zivari-Piran, Hossein; Enright, Wayne H.: Accurate first-order sensitivity analysis for delay differential equations (2012)
  9. Lunacek, Monte; Nag, Ambarish; Alber, David M.; Gruchalla, Kenny; Chang, Christopher H.; Graf, Peter A.: Simulation, characterization, and optimization of metabolic models with the high performance systems biology toolkit (2011)
  10. Koutsawa, Yao; Belouettar, Salim; Makradi, Ahmed; Nasser, Houssein: Sensitivities of effective properties computed using micromechanics differential schemes and high-order Taylor series: application to piezo-polymer composites (2010)
  11. Stumm, Philipp; Walther, Andrea: New algorithms for optimal online checkpointing (2010)
  12. Alexe, Mihai; Sandu, Adrian: On the discrete adjoints of adaptive time stepping algorithms (2009)
  13. Alexe, Mihai; Sandu, Adrian: Forward and adjoint sensitivity analysis with continuous explicit Runge-Kutta schemes (2009)
  14. Donzé, Alexandre; Krogh, Bruce; Rajhans, Akshay: Parameter synthesis for hybrid systems with an application to Simulink models (2009)
  15. Duran, Ahmet: Sensitivity analysis of asset flow differential equations and volatility comparison of two related variables (2009)
  16. Serban, Radu: A parallel computational model for sensitivity analysis in optimization for robustness (2009)
  17. Wilkins, A.Katharina; Tidor, Bruce; White, Jacob; Barton, Paul I.: Sensitivity analysis for oscillating dynamical systems (2009)
  18. Duran, Ahmet; Caginalp, Gunduz: Parameter optimization for differential equations in asset price forecasting (2008)
  19. Griewank, Andreas; Walther, Andrea: Evaluating derivatives. Principles and techniques of algorithmic differentiation (2008)
  20. Aittokallio, T.; Gyllenberg, M.; Polo, O.; Virkki, A.: Parameter estimation of a respiratory control model from noninvasive carbon dioxide measurements during sleep (2007)

1 2 next