diff options
| author | schnetter <schnetter@578cdeb0-5ea1-4b81-8215-5a3b8777ee0b> | 2006-02-19 00:31:54 +0000 |
|---|---|---|
| committer | schnetter <schnetter@578cdeb0-5ea1-4b81-8215-5a3b8777ee0b> | 2006-02-19 00:31:54 +0000 |
| commit | c4896403254a6de126501a42f9cb046538a44727 (patch) | |
| tree | 5044935d014dab7bca78235a52a54f64dad45402 /src/IndexArrays.c | |
| parent | a7134900411d8adb1cb1f4d6962dfcf892983326 (diff) | |
Implement RK45 Cash-Karp integrator in MoL.
This integrator is similar to the existing RK45 integrator. It also
supports adaptive time stepping, but uses slightly different
coefficients. The Numerical Recipes say that it has "slightly better
error properties".
git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@110 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
Diffstat (limited to 'src/IndexArrays.c')
| -rw-r--r-- | src/IndexArrays.c | 6 |
1 files changed, 5 insertions, 1 deletions
diff --git a/src/IndexArrays.c b/src/IndexArrays.c index 505b1ab..5e59734 100644 --- a/src/IndexArrays.c +++ b/src/IndexArrays.c @@ -306,7 +306,11 @@ void MoL_SetupIndexArrays(CCTK_ARGUMENTS) } else if (CCTK_EQUALS(ODE_Method,"RK45")) { - sprintf(infoline, "Runge-Kutta 45"); + sprintf(infoline, "Runge-Kutta 45 (Fehlberg)"); + } + else if (CCTK_EQUALS(ODE_Method,"RK45CK")) + { + sprintf(infoline, "Runge-Kutta 45 (Cash-Karp)"); } else if (CCTK_EQUALS(ODE_Method,"RK65")) { |