From 2e25da27f14566000fb34d3859bfb470bf5fd1da Mon Sep 17 00:00:00 2001
From: Tavian Barnes <tavianator@gmail.com>
Date: Wed, 6 Oct 2010 14:38:22 -0400
Subject: Add a new RKIntegrator base class.

---
 src/vZ/RK.hpp | 48 ++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 48 insertions(+)
 create mode 100644 src/vZ/RK.hpp

(limited to 'src/vZ/RK.hpp')

diff --git a/src/vZ/RK.hpp b/src/vZ/RK.hpp
new file mode 100644
index 0000000..3ba4b8a
--- /dev/null
+++ b/src/vZ/RK.hpp
@@ -0,0 +1,48 @@
+/*************************************************************************
+ * Copyright (C) 2009-2010 Tavian Barnes <tavianator@gmail.com>          *
+ *                                                                       *
+ * This file is part of The vZ Library.                                  *
+ *                                                                       *
+ * The vZ Library is free software; you can redistribute it and/or       *
+ * modify it under the terms of the GNU Lesser General Public License as *
+ * published by the Free Software Foundation; either version 3 of the    *
+ * License, or (at your option) any later version.                       *
+ *                                                                       *
+ * The vZ Library is distributed in the hope that it will be useful, but *
+ * WITHOUT ANY WARRANTY; without even the implied warranty of            *
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU     *
+ * Lesser General Public License for more details.                       *
+ *                                                                       *
+ * You should have received a copy of the GNU Lesser General Public      *
+ * License along with this program.  If not, see                         *
+ * <http://www.gnu.org/licenses/>.                                       *
+ *************************************************************************/
+
+#ifndef VZ_RK_HPP
+#define VZ_RK_HPP
+
+#include <vector>
+
+namespace vZ
+{
+  // Base class for Runge-Kutta type algorithms
+  template <typename T>
+  class GenericRKIntegrator : public GenericIntegrator<T>
+  {
+  public:
+    typedef typename GenericIntegrator<T>::Function Function;
+
+  protected:
+    // Coefficients in the tableau representation of the RK algorithm
+    typedef std::vector<std::vector<T> > ACoefficients;
+    typedef std::vector<T>               BCoefficients;
+
+    GenericRKIntegrator(Function f, T dt) : GenericIntegrator<T>(f, dt) { }
+    virtual ~GenericRKIntegrator() { }
+  };
+
+  // Type alias
+  typedef GenericRKIntegrator<double> RKIntegrator;
+}
+
+#endif // VZ_RK_HPP
-- 
cgit v1.2.3