| 348 | 
  | 
  // entire system: | 
| 349 | 
  | 
 | 
| 350 | 
  | 
  ndf = ndf - 3 - nZconstraints; | 
| 351 | 
– | 
 | 
| 352 | 
– | 
  std::cerr << "ndf = " << ndf; | 
| 351 | 
  | 
 | 
| 352 | 
  | 
  return ndf; | 
| 353 | 
  | 
} | 
| 510 | 
  | 
               "\tCurrent Value of LJrcut = %G at time %G\n " | 
| 511 | 
  | 
               "\tThis is larger than half of at least one of the\n" | 
| 512 | 
  | 
               "\tperiodic box vectors.  Right now, the Box matrix is:\n" | 
| 513 | 
< | 
               "\n, %G" | 
| 513 | 
> | 
               "\n" | 
| 514 | 
  | 
               "\t[ %G %G %G ]\n" | 
| 515 | 
  | 
               "\t[ %G %G %G ]\n" | 
| 516 | 
  | 
               "\t[ %G %G %G ]\n", | 
| 517 | 
< | 
               rCut, currentTime, maxCutoff, | 
| 517 | 
> | 
               rCut, currentTime, | 
| 518 | 
  | 
               Hmat[0][0], Hmat[0][1], Hmat[0][2], | 
| 519 | 
  | 
               Hmat[1][0], Hmat[1][1], Hmat[1][2], | 
| 520 | 
  | 
               Hmat[2][0], Hmat[2][1], Hmat[2][2]); |