1

これは、微分方程式系を解く私のコードの一部です。

Vector dydx(Neq);
void DiffEq(Vector x, Vector &dydx)
{
    dydx(0) = x(1);
    dydx(1) = -x(0);
}

double MidPoint(int n, Vector x)
{
    double h=H/n;
    Matrix z(n+1,n+1);
    z.fillRow(0,x);
    DiffEq(x, dydx);
    z.fillRow(1,AddVec(x, dydx*h));      //Error: Invalid use of void expression
    for (int j=1; j<n; j++)
    {
        DiffEq(z.getRow(j), dydx);
        z.fillRow(j+1, AddVec(z.getRow(j-1), dydx*h*2));     //error: void value not ignored as it ought to be
    }   
    DiffEq(z.getRow(n), dydx);
    return 0.5*AddVec(z.getRow(n), z.getRow(n-1), dydx*h);       //Error: Invalid use of void expression
}

クラス Vector と Matrix はカスタムです。これがベクトル

class Vector
{
    public:
    Vector(size_t size): vSize(size), vData(size){}
    int getSize(){return vSize;}
    double& operator()(size_t i){return vData[i];}
    double operator()(size_t i) const {return vData[i];}
    void operator+(double d)                        // These only overload the Vector + int operator
    {
        for (int i=0; i<vSize; i++) {vData[i]=vData[i]+d;}
    }

    void operator*(double d)
    {
        for (int i=0; i<vSize; i++) {vData[i]=vData[i]*d;}
    }


    size_t vSize;
    vector<double> vData;
};

機能付き

Vector AddVec(Vector v1, Vector v2)
{
    Vector totVec(v1.getSize());
    for (int i=0; i<v1.getSize(); i++) {totVec(i) = v1(i) + v2(i);}
    return totVec;
}

そして、3つのベクトルに対して同じ機能。

エラーは、関数に void を渡していることを意味していることに気付きましたが、どこで問題が発生したのかわかりません。テストプログラムを書いてみると、すべてうまくいっているようです。

4

2 に答える 2