我从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版本。更新到最新的稳定版本后,所描述的行为消失了。

10-04 23:20