Commit dc86a1be authored by Simon Spannagel's avatar Simon Spannagel
Browse files

Add Cash-Karp method to RK tableaus

parent 44299cb6

src/tools/runge_kutta.h

+12 −0
Original line number Diff line number Diff line
@@ -203,6 +203,18 @@ namespace allpix { 
            -8.0/27, 2, -3544.0/2565, 1859.0/4104, -11.0/40, 0,

            16.0/135, 0, 6656.0/12825, 28561.0/56430, -9.0/50, 2.0/55,

            25.0/216, 0, 1408.0/2565, 2197.0/4104, -1.0/5, 0).finished());

        /**

         * @brief Runge-Kutta-Cash-Karp method

         */

        static const auto RKCK((Eigen::Matrix<double, 8, 6>() <<

            0, 0, 0, 0, 0, 0,

            1.0/5, 0, 0, 0, 0, 0,

            3.0/40, 9.0/40, 0, 0, 0, 0,

            3.0/10, -9.0/10, 6.0/5, 0, 0, 0,

            -11.0/54, 5.0/2, -70.0/27, 35.0/27, 0, 0,

            1631.0/55296, 175.0/512, 575.0/13824, 44275.0/110592, 253.0/4096, 0,

            37.0/378, 0, 250.0/621, 125.0/594, 0, 512.0/1771,

            2825.0/27648, 0, 18575.0/48384, 13525.0/55296, 277.0/14336, 1.0/4).finished());

    }

    // clang-format on