我从GNU科学图书馆复数函数中看到一些奇怪的行为。以下代码段演示了该问题:
#include <iostream>
#include <gsl/gsl_complex.h>
#include <gsl/gsl_complex_math.h>
int main() {
double realNumber = 1.0;
gsl_complex complexNumber;
GSL_SET_COMPLEX(&complexNumber, 1.0, 1.0);
printf("%p: ", &realNumber);
printf("%f\n", realNumber);
gsl_complex_add(complexNumber,complexNumber);
printf("%p: ", &realNumber);
printf("%f\n", realNumber);
return 0;
}
输出是
0061FF18: 1.000000
0061FF1C: 0.000000
函数
gsl_complex_add
将指针的值更改为realNumber
,使其不再指向realNumber
,即使realNumber
似乎没有以任何方式链接到该函数。我认为这是由于某种溢出造成的。如果我使用
std::cout
而不是printf
,则不会发生此问题;我想这是由于编译器造成的。如果我将
realNumber
声明为static
,则不会发生此问题。如果我将
realNumber
声明为const
,则printf("%p: ", &realNumber);
返回错误的值,但是printf("%f\n", realNumber);
返回正确的值。这是一个错误还是我在做我不应该做的事情?
我在Windows机器上使用g++编译器(版本6.3.0)。
在此先感谢您的帮助。
最佳答案
经过进一步调查,我意识到我使用的是非常旧的GSL版本。更新到最新的稳定版本后,所描述的行为消失了。