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