Question

我最近开始使用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);
        }

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

Answer 1

最简单的方法是在集成函数内部制作数组的本地副本。

改变你传递的方式＆＃39; n＆＃39;要运行到＆＃39; const int n＆＃39;，这样你就可以在里面创建类似double currentSolution[n];的东西，并将元素从初始数组复制到新数组。这种方法可以保存您的初始数组，除非您意外地＃34;在某处修改它。

为了防止这种事故修改的可能性，我们需要更深入地使用其中一个stl容器。我认为你对std::valarray<T>没问题。

更改传递方式为const std::valarray<double>&并再次制作非常规本地副本。

Answer 2

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

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 ++功能而不是原始指针等更好。

用于解决ode系统的C ++数值积分器

2 个答案: