Skip to main content
Computational Science

Return to Question

edited body
Source Link
GertVdE
GertVdE
  • 6.2k
  • 1
  • 21
  • 36

I want soto solve a relatively small system of stiff ODEs (< 10 first-order equations) using high precision floating point arithmetic (using MPFR or alike). What would be the easiest algorithm to implement/port? I'm not only interested in the "stepper" algorithm, but in the combination of stepper/error estimation/step size control.

I want so solve a relatively small system of stiff ODEs (< 10 first-order equations) using high precision floating point arithmetic (using MPFR or alike). What would be the easiest algorithm to implement/port? I'm not only interested in the "stepper" algorithm, but in the combination of stepper/error estimation/step size control.

I want to solve a relatively small system of stiff ODEs (< 10 first-order equations) using high precision floating point arithmetic (using MPFR or alike). What would be the easiest algorithm to implement/port? I'm not only interested in the "stepper" algorithm, but in the combination of stepper/error estimation/step size control.

Tweeted twitter.com/#!/StackSciComp/status/214531846485192705
Source Link
GertVdE
GertVdE
  • 6.2k
  • 1
  • 21
  • 36

What algorithm for solving a set of stiff ODEs would be easiest to port to high precision floating point arithmetic?

I want so solve a relatively small system of stiff ODEs (< 10 first-order equations) using high precision floating point arithmetic (using MPFR or alike). What would be the easiest algorithm to implement/port? I'm not only interested in the "stepper" algorithm, but in the combination of stepper/error estimation/step size control.