summaryrefslogtreecommitdiffstats
path: root/src/vZ/BS23.hpp
diff options
context:
space:
mode:
authorTavian Barnes <tavianator@gmail.com>2010-10-06 23:05:48 -0400
committerTavian Barnes <tavianator@gmail.com>2010-10-06 23:05:48 -0400
commite4857d7abd64993c18986b7b5f7e3b0c14eb27e0 (patch)
tree7c772c730fa8e7dd37d97a417f83bfa75fe547cb /src/vZ/BS23.hpp
parent486a188b074c0c4f3ad771f536159e1bd2045fe2 (diff)
downloadvz-e4857d7abd64993c18986b7b5f7e3b0c14eb27e0.tar.xz
Clean up integrator comments.
Diffstat (limited to 'src/vZ/BS23.hpp')
-rw-r--r--src/vZ/BS23.hpp20
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>
{