The problem I'm having with my code is located at
Code:
void fGR(double x[NSTEPS], double y[NSTEPS], double sig[NSTEPS],
double **GRalpha, double **GRbeta, int ma, double chi,
int nu, double Q, double GRa[GRma])
* * *
fout1 << "GRma:" << ma << " GR[]:" << GRa[0] << " " << GRa[1] << endl;
fout1 << "FF L50: GRalpha[][0]\tGRalpha[][1]\tGRbeta[][0]" << endl;
for(i=1; i<=ma; i++)
{
fout1 << endl << "54GR:" << GRa[0] << " " << GRa[1] << endl;
fout1 << "[" << i << "][] \t ";
for(j=1; j<=ma; j++)
{
GRalpha[i][j] = 0.0;
fout1 << GRalpha[i][j] << "\t\t";
}
fout1 << endl << "61GR:" << GRa[0] << " " << GRa[1] << endl;
GRbeta[i][0] = 0.0;
fout1 << GRbeta[i][0] << endl;
}
As you can see, there doesn't seem to be any trouble; I'm mostly just outputting some data to a file. The problem I'm having is (output; not code):
Code:
GRma:2 GR[]:430.819 6.0904
FF L50: GRalpha[][0] GRalpha[][1] GRbeta[][0]
54GR:430.819 6.0904
[1][] 0 0
61GR:430.819 6.0904
0
54GR:430.819 6.0904
[2][] 0 0
61GR:0 6.0904
0
I input into the function GRa[0] = 430.819, GRa[1] = 6.0904. After the second iteration of the outer for loop, it suddenly changes GRa[0] to 0 for no apparent reason. Why would a program suddenly do this?