一個函數中的局部變量變化導致一個不相關的函數，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中聲明的，如果這樣使用函數rand（）函數randNorm（const double mean，const double stddev）（來自這裏http://www-personal.umich.edu/~wagnerr/MersenneTwister.h）。 爲了在每次使用rand和srand時都有不同的值。看這裏：http://www.cplusplus.com/reference/clibrary/cstdlib/rand/。