| 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; |