我的班级NRRanNormal
代表一个正态分布的随机变量。默认情况下,实例通常以均值0和stdev 1(即标准正态随机变量)分布。
有时当我复制NRRanNormal
个对象时,复制到(或通过复制构造函数构造)的对象的均值和stdev是乱码和无意义的。我很难找到这种杂乱的原因。
出于测试目的,以下函数显示给定NRRanNormal
对象的mean和stdev:
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 */
}
};
答案 0 :(得分:3)
这是你的问题
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);
}
您的复制构造函数无法复制mean,stddev和seed。您的赋值运算符应该具有相同的问题
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;
}
我猜你们非常专注于课堂上那些你忘了基本内容的棘手指针。
BTW代码
MyClass A;
MyClass B = A;
实际上是由编译器转换为
MyClass A;
MyClass B(A);
换句话说MyClass B = A;
调用复制构造函数(假设A和B的类型相同)。