| 2 | 
  | 
 | 
| 3 | 
  | 
#include "Atom.hpp" | 
| 4 | 
  | 
 | 
| 5 | 
– | 
double* Atom::pos; // the position array | 
| 6 | 
– | 
double* Atom::vel; // the velocity array | 
| 7 | 
– | 
double* Atom::frc; // the forc array | 
| 8 | 
– | 
double* Atom::trq; // the torque vector  ( space fixed ) | 
| 9 | 
– | 
double* Atom::Amat; // the rotation matrix | 
| 10 | 
– | 
double* Atom::mu;   // the array of dipole moments | 
| 11 | 
– | 
double* Atom::ul;   // the lab frame unit directional vector | 
| 5 | 
  | 
 | 
| 6 | 
+ | 
 | 
| 7 | 
  | 
void DirectionalAtom::setA( double the_A[3][3] ){ | 
| 8 | 
  | 
   | 
| 9 | 
  | 
  Amat[Axx] = the_A[0][0]; Amat[Axy] = the_A[0][1]; Amat[Axz] = the_A[0][2]; | 
| 60 | 
  | 
  the_A[2][2] = Amat[Azz]; | 
| 61 | 
  | 
} | 
| 62 | 
  | 
 | 
| 63 | 
+ | 
void DirectionalAtom::printAmatIndex( void ){ | 
| 64 | 
  | 
 | 
| 65 | 
+ | 
  std::cerr << "Atom[" << index << "] index =>\n"  | 
| 66 | 
+ | 
            << "[ " << Axx << ", " << Axy << ", " << Axz << " ]\n" | 
| 67 | 
+ | 
            << "[ " << Ayx << ", " << Ayy << ", " << Ayz << " ]\n" | 
| 68 | 
+ | 
            << "[ " << Azx << ", " << Azy << ", " << Azz << " ]\n"; | 
| 69 | 
+ | 
} | 
| 70 | 
+ | 
 | 
| 71 | 
+ | 
 | 
| 72 | 
  | 
void DirectionalAtom::getU( double the_u[3] ){ | 
| 73 | 
  | 
   | 
| 74 | 
  | 
  the_u[0] = sux; |