c++ - C++数值积分器来求解ode的系统

我最近开始使用C ++，我刚刚创建了一个类，该类允许集成用户定义的ode系统。它使用两个不同的积分器来比较其性能。这是代码的一般布局：

    class integrators {
    private:
      double  ti;    // initial time
            double *xi;  // initial solution
            double  tf;    // end time
            double  dt;  // time step
            int      n;  // number of ode's

        public:
            // Function prototypes
            double   f(double, double *, double *);  // function to integrate
            double rk4(int, double, double, double, double *, double *);
            double dp8(int, double, double, double, double *, double *);
        };

        // 4th Order Runge-Kutta function
        double integrators::rk4(int n, double ti, double tf, double dt, double *xi, double *xf) {
            // Function statements
        }

        // 8th Order Dormand-Prince function
        double integrators::dp8(int n, double ti, double tf, double dt, double *xi, double *xf) {
            // Function statements
        }

       // System of first order differential equations
       double integrators::f(double t, double *x, double *dx) {
           // Function statements
       }

        int main() {
            // Initial conditions and time related parameters
            const int n = 4;
            double t0, tmax, dt;
            double x0[n], xf[n];

            x0[0] = 0.0;
            x0[1] = 0.0;
            x0[2] = 1.0;
            x0[3] = 2.0;

            // Calling class integrators
            integrators example01;
            integrators example02;

            // First integrator
            example02.dp8(n, t0, tmax, dt, x0, xf);

            // Second integrator
            example01.rk4(n, t0, tmax, dt, x0, xf);
        }

问题是，包含主条件x0的数组在执行第一个积分器后会发生变化，除非我定义另一个具有相同初始条件（x0_rk4和x0_dp8）的数组，否则我无法对第二个积分器使用相同的初始条件。有没有更专业的方法来保持此数组恒定以便在两个积分器中都使用？

最佳答案

不，不是。但是存在一个更优雅的解决方案：

std::array<double, n> x0_rk4 = { 0.0, 0.0, 1.0, 2.0 };
auto x0_dp8 = x0_rk4; // copy!

您将必须使用x0_rk4.data()访问基础数组。请注意，如果您使用std::array和其他现代C ++功能而不是原始指针等会更好。

关于c++ - C++数值积分器来求解ode的系统，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/42747227/