0

我的班级NRRanNormal代表一个正态分布的随机变量。默认情况下，实例正态分布，均值为 0，标准偏差为 1（即标准正态随机变量）。

有时当我复制NRRanNormal对象时，复制到（或通过复制构造函数构造）的对象的均值和标准差是乱码和胡说八道。我很难找到这种乱码的原因。

NRRanNormal出于测试目的，以下函数显示给定对象的均值和标准差：

void go(NRRanNormal& rv, const string label) {
    std::cout << label     << "\n"
              << "Mean:  " << rv.getMean()  << "\n"
              << "Stdev: " << rv.getStdev() << "\n\n";
}

现在，让我们看看以下 4 种情况会发生什么：

NRRanNormal foo;
go(foo, "foo");

NRRanNormal bar1 = foo;
go(bar1, "bar1");

NRRanNormal bar2;
bar2 = foo;
go(bar2, "bar2");

NRRanNormal bar3(foo);
go(bar3, "bar3");

上述语句的输出如下：

foo
Mean:  0
Stdev: 1

bar1
Mean:  5.55633e-317
Stdev: 6.95332e-310

bar2
Mean:  0
Stdev: 1

bar3
Mean:  0
Stdev: 0

如您所见，只需实例化一个对象 ( foo) 即可按预期工作。

现在，当我这样做时NRRanNormal bar1 = foo;，对象bar1是乱码。但是，当我这样做时，NRRanNormal bar2; bar2 = foo;对象不会乱码。这让我很困惑。我认为一个语句块如bar2

MyClass A;
MyClass B = A;

实际上是由编译器转换为语句块

MyClass A;
MyClass B;
B = A;

因此，除非我刚刚在上面写的内容不正确，否则似乎bar1并且bar2应该具有完全相同的成员值。但正如您从上面粘贴的输出中看到的那样，bar1是乱码，而bar2很好。

怎么会这样？

您还会注意到这bar3是乱码。我不确定这是同一个问题，还是另一个问题。

这是接口和实现的简化版本NRRanNormal：

class NRRanNormal {
public:

    NRRanNormal();

    ~NRRanNormal();

    NRRanNormal(const NRRanNormal& nrran);

    NRRanNormal& operator= (const NRRanNormal& nrran);

    double getMean()  const;    
    double getStdev() const;    
    long   getSeed()  const;

private:    
    double m_mean, m_stdev;     
    long m_seed;    
    Normaldev* stream; // underlying C struct RN generator

};

NRRanNormal::NRRanNormal() { // by default, N(0,1)
    m_mean  = 0.0;
    m_stdev = 1.0;
    m_seed  = 12345L;
    stream  = new Normaldev(m_mean, m_stdev, m_seed);
}

NRRanNormal::~NRRanNormal() { delete stream; }

NRRanNormal::NRRanNormal(const NRRanNormal& nrran) {
    stream = new Normaldev(nrran.getMean(),nrran.getStdev(),nrran.getSeed());
    *stream = *(nrran.stream);
}

NRRanNormal& NRRanNormal::operator= (const NRRanNormal& nrran) {            
    if(this == &nrran)
        return *this;

    delete stream;
    stream = new Normaldev(nrran.getMean(),nrran.getStdev(),nrran.getSeed());
    *stream = *(nrran.stream);

    return *this;
}

double NRRanNormal::getMean()  const { return m_mean; }    
double NRRanNormal::getStdev() const { return m_stdev; }    
long   NRRanNormal::getSeed()  const { return m_seed; }

该Normaldev结构来自 Numerical Recipes 3d Edition。

我的复制赋值运算符或复制构造函数有问题吗？

这是Normaldev，剥离了专有计算。

typedef double Doub;        
typedef unsigned long long int Ullong;
typedef unsigned int Uint;

struct Ranq1 {
    Ullong v;
    Ranq1(Ullong j) : v(/* some long number here */) {
        /* proprietary calculations here */
    }
    inline Ullong int64() {
        /* proprietary calculations here */
    }
    inline Doub doub() { /* proprietary calculations here */ }
    inline Uint int32() { return (Uint)int64(); }
};

struct Normaldev : Ranq1 {
    Doub mu,sig;
    Normaldev(Doub mmu, Doub ssig, Ullong i):
    Ranq1(i), mu(mmu), sig(ssig){}
    Doub dev() {
        /* proprietary calculations here */
    }
};
4

1 回答 1

4

这是你的问题

NRRanNormal::NRRanNormal(const NRRanNormal& nrran) {
    stream = new Normaldev(nrran.getMean(),nrran.getStdev(),nrran.getSeed());
    *stream = *(nrran.stream);
}

应该

NRRanNormal::NRRanNormal(const NRRanNormal& nrran) : 
    m_mean(nrran.m_mean), 
    m_stdev(nrran.m_stdev), 
    m_seed(nrran.m_seed)
{
    stream = new Normaldev(nrran.getMean(),nrran.getStdev(),nrran.getSeed());
    *stream = *(nrran.stream);
}

您的复制构造函数无法复制均值、标准差和种子。您的赋值运算符有同样的问题，它应该是

NRRanNormal& NRRanNormal::operator= (const NRRanNormal& nrran) {            
    if(this == &nrran)
        return *this;

    m_mean  = nrran.m_mean;
    m_stdev = nrran.m_stdev;
    m_seed  = nrran.m_seed;
    delete stream;
    stream = new Normaldev(nrran.getMean(),nrran.getStdev(),nrran.getSeed());
    *stream = *(nrran.stream);

    return *this;
}

我猜你在课堂上太专注于棘手的指针，而忘记了基本的东西。

顺便说一句代码

MyClass A;
MyClass B = A;

实际上是由编译器转换为

MyClass A;
MyClass B(A);

换句话说，MyClass B = A;调用复制构造函数（假设 A 和 B 是相同的类型）。

于 2013-09-27T05:46:13.063 回答