diff options
| author | Tavian Barnes <tavianator@gmail.com> | 2010-10-06 23:05:48 -0400 |
|---|---|---|
| committer | Tavian Barnes <tavianator@gmail.com> | 2010-10-06 23:05:48 -0400 |
| commit | e4857d7abd64993c18986b7b5f7e3b0c14eb27e0 (patch) | |
| tree | 7c772c730fa8e7dd37d97a417f83bfa75fe547cb /src/vZ/BS23.hpp | |
| parent | 486a188b074c0c4f3ad771f536159e1bd2045fe2 (diff) | |
| download | vz-e4857d7abd64993c18986b7b5f7e3b0c14eb27e0.tar.xz | |
Clean up integrator comments.
Diffstat (limited to 'src/vZ/BS23.hpp')
| -rw-r--r-- | src/vZ/BS23.hpp | 20 |
1 files changed, 9 insertions, 11 deletions
diff --git a/src/vZ/BS23.hpp b/src/vZ/BS23.hpp index 1092d81..f81cad1 100644 --- a/src/vZ/BS23.hpp +++ b/src/vZ/BS23.hpp @@ -23,20 +23,18 @@ namespace vZ { - // Heun-Euler method + // Bogacki-Shampine method // - // Second-order with embedded first-order + // Third-order with embedded second-order // Its tableau is: // - // 0 | - // 1 | 1 - // --+--------- - // | 1/2 1/2 - // | 1 0 - // - // k1 = dt*f(y[n]) - // k2 = dt*f(y[n] + dt*k1) - // y[n + 1] = y[n] + 1/2*(k1 + k2) + // 0 | + // 1/2 | 1/2 + // 3/4 | 0 3/4 + // 1 | 2/9 1/3 4/9 + // ----+----------------- + // b | 2/9 1/3 4/9 0 + // b* | 7/24 1/4 1/3 1/8 template <typename Y> class GenericBS23Integrator : public GenericAdaptiveIntegrator<Y> { |