77 |
|
for(i=0; i < 3; i++) |
78 |
|
for (j=0; j < 3; j++) Hmat[i][j] = theBox[i][j]; |
79 |
|
|
80 |
< |
cerr |
81 |
< |
<< "setting Hmat ->\n" |
82 |
< |
<< "[ " << Hmat[0][0] << ", " << Hmat[0][1] << ", " << Hmat[0][2] << " ]\n" |
83 |
< |
<< "[ " << Hmat[1][0] << ", " << Hmat[1][1] << ", " << Hmat[1][2] << " ]\n" |
84 |
< |
<< "[ " << Hmat[2][0] << ", " << Hmat[2][1] << ", " << Hmat[2][2] << " ]\n"; |
80 |
> |
// cerr |
81 |
> |
// << "setting Hmat ->\n" |
82 |
> |
// << "[ " << Hmat[0][0] << ", " << Hmat[0][1] << ", " << Hmat[0][2] << " ]\n" |
83 |
> |
// << "[ " << Hmat[1][0] << ", " << Hmat[1][1] << ", " << Hmat[1][2] << " ]\n" |
84 |
> |
// << "[ " << Hmat[2][0] << ", " << Hmat[2][1] << ", " << Hmat[2][2] << " ]\n"; |
85 |
|
|
86 |
|
calcBoxL(); |
87 |
|
calcHmatInv(); |
153 |
|
double theBox[3][3]; |
154 |
|
int i, j; |
155 |
|
|
156 |
< |
cerr << "Scaling box by " << scale << "\n"; |
156 |
> |
// cerr << "Scaling box by " << scale << "\n"; |
157 |
|
|
158 |
|
for(i=0; i<3; i++) |
159 |
|
for (j=0; j<3; j++) theBox[i][j] = Hmat[i][j]*scale; |