1つの関数内のローカル変数が無関係な関数に変化する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

randは関数rand（）を使って関数randNorm（const double double、const double stddev）（ここではhttp://www-personal.umich.edu/~wagnerr/MersenneTwister.hから）の場合にMersenneTwister.hで宣言されていると思います。 randとsrandを共同で使用するたびに異なる値を使用するためです。ここを見てください：http://www.cplusplus.com/reference/clibrary/cstdlib/rand/。