一个函数中的局部变量更改会导致不相关的函数 C++

Question

预先警告大家，我是一名物理学家，正在编写一个模拟，从专业程序员的角度来看，这可能是玩具代码，所以我的编码知识不是很好。不管怎样，解释一下 700 行，我的代码看起来像下面这样：

#define _USE_MATH_DEFINES
#include "MersenneTwister.h"  //Random number generator
#include "Particle.h"         //Class that just keeps track of coordinates
#include <algorithm>
#include <vector>
using namespace std;

const double SomeGlobalParameters;

//Filling an array with gaussian random numbers

void GaussianDisorder(double EnergyArray[Nx][Ny][Nz], double Sigma)
{
    for (int i=0; i<Nx; i++){
        for (int j=0; j<Ny; j++){
            for (int k=0; k<Nz; k++){
                EnergyArray[i][j][k] = rnd.randNorm(0, Sigma);
            }
        }
    }
}

//Using the above array (read only) to do some calculations

void CalcRates(Particle &charge1, const double (&EnergyArray)[Nx][Ny][Nz], double (&RateArray)[12])
{
    int X1 = charge1.xPosition();
    double NearRate, MagSqr, Dist, OldDist, OldEnergy, NewEnergy;    //Declaring stuff for later
    const double Parity=1.0;

    InitializeStuffAbove(SomeGlobalParameters);           // Give stuff values based on parameters

    if (Dist == 0)
        RateArray[0] = NearRate;
    else
        NewEnergy = Parity*EnergyArray2[X1+1][Y1][Z1] - Efield[0] + Coulomb/(Dist);

int main()
{
    double EnergyArray[Nx][Ny][Nz];
    double RateArray[12];

    GaussianDisorder(EnergyArray);
    CalcRates(charge1, EnergyArray, RateArray);

    return 0;
}

我的问题是，当我将

const double Parity=1.0;

，这是一个局部变量，更改为

const double Parity=-1.0;

我的随机数生成器时，

void GaussianDisorder

给出了不同的结果，即使我给了它一个常量种子，并且它给出了如果我多次运行代码而不更改任何内容，则顺序相同。如果我遗漏了一些重要的内容，我提前道歉。

Answer 1

我猜 rnd 是在 MersenneTwister.h 中声明的，如果是的话，函数 randNorm(const double Mean, const double stddev) （从这里 http://www-personal.umich.edu/~wagnerr/MersenneTwister.h）使用函数 rand()。
为了每次都有不同的值，您应该联合使用 rand 和 srand。看这里： http://www.cplusplus.com/reference/clibrary/cstdlib/rand /。