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

3

您的运营商返回 void

void operator+(double d)                        // These only overload the Vector + int 
{
    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;}
}

所以当你调用AddVec(x, dydx*h)它时AddVec(x, void)

于 2013-09-24T18:03:01.630 回答
1

您重载了与传统运算符不同的运算符，不要那样做！！！ 这样做会给 C++ 一个坏名声。

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

这不会返回任何内容，因此您不能像这样使用它：

 z.fillRow(1,AddVec(x, dydx*h)); 
                       ^^^^^^

operator*应该返回一个新对象，而不是修改其左操作数。如果你想修改左操作数，使用operator*=但你仍然应该返回一些东西。

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

        Vector operator*(double d)
        {
            Vector v(*this); 
            v *= d;
            return v;
        }

作为operator*非成员会更好，并且应该通过引用来获取它的论点：

        Vector operator*(const Vector& v, double d)
        {
            Vector v2(v); 
            v2 *= d;
            return v2;
        }
于 2013-09-24T18:05:18.187 回答