2

如果这很明显，我很抱歉，但我对来自 Python / MATLAB / Mathematica 背景的 C++ 非常陌生。我使用有限差分空间离散化为经典的一维热方程编写了一个简单的求解器，以便使用Odeint库的功能并将性能与其他库进行比较。代码应该是不言自明的：

#include <iostream>
#include <boost/math/constants/constants.hpp>
#include <boost/array.hpp>
#include <boost/numeric/odeint.hpp>

using namespace std;
using namespace boost::numeric::odeint;

const double a_sq = 1;
const int p = 10;
const double h = 1 / p;

double pi = boost::math::constants::pi<double>();

typedef boost::array<double, p> state_type;


void heat_equation(const state_type &x, state_type &dxdt, double t)
{
    int i;
    for (i=1; i<p; i++)
    {
        dxdt[i] = a_sq * (dxdt[i+1] - 2*dxdt[i] + dxdt[i-1]) / h / h;
    }     
    dxdt[0] = 0;
    dxdt[p] = 0;
}

void initial_conditions(state_type &x)
{
   int i;
   for (i=0; i<=p; i++)
   {
      x[i] = sin(pi*i*h);
   }
}

void write_heat_equation(const state_type &x, const double t)
{
   cout << t << '\t' << x[0] << '\t' << x[p] << '\t' << endl;
}

int main()
{
    state_type x;
    initial_conditions(x);    
    integrate(heat_equation, x, 0.0, 10.0, 0.1, write_heat_equation);
}

这在 Ubuntu 14.04 上使用 g++ 4.8.2 和来自 Ubuntu 存储库的最新 boost 库编译得很好。但是，当我运行生成的可执行文件时，出现以下错误：

***** Internal Program Error - assertion (i < N) failed in T& boost::array<T, N>::operator[](boost::array<T, N>::size_type) [with T = double; long unsigned int N = 10ul; boost::array<T, N>::reference = double&; boost::array<T, N>::size_type = long unsigned int]:
/usr/include/boost/array.hpp(123): out of range
Aborted (core dumped)

不幸的是，这对我的新手大脑并不是特别有帮助，我不知道如何解决这个问题。是什么导致了错误？

4

1 回答 1

0

从零开始计算数组或 N 元素向量的元素。最后一个元素的索引为 N-1。因此，您需要更改 for 循环以从 i 迭代到 p-1 并且您需要修改 dxdt[p] = 0; 行 到 dxdt[p-1] = 0：

for (i=1; i<p-1; i++)
{
    dxdt[i] = a_sq * (dxdt[i+1] - 2*dxdt[i] + dxdt[i-1]) / h / h;
}     
dxdt[0] = 0;
dxdt[p-1] = 0;
于 2014-06-05T18:54:03.567 回答