| 38 | 
  | 
 | 
| 39 | 
  | 
void DirectionalAtom::setPrevA(const RotMat3x3d& a) { | 
| 40 | 
  | 
    ((snapshotMan_->getPrevSnapshot())->*storage_).aMat[localIndex_] = a; | 
| 41 | 
< | 
    ((snapshotMan_->getPrevSnapshot())->*storage_).unitVector[localIndex_] = a.inverse() * sU_.getColum(2); | 
| 41 | 
> | 
    ((snapshotMan_->getPrevSnapshot())->*storage_).unitVector[localIndex_] = a.inverse() * sU_.getColumn(2); | 
| 42 | 
  | 
} | 
| 43 | 
  | 
 | 
| 44 | 
  | 
       | 
| 45 | 
  | 
void DirectionalAtom::setA(const RotMat3x3d& a) { | 
| 46 | 
  | 
    ((snapshotMan_->getCurrentSnapshot())->*storage_).aMat[localIndex_] = a; | 
| 47 | 
< | 
    ((snapshotMan_->getCurrentSnapshot())->*storage_).unitVector[localIndex_] = a.inverse() * sU_.getColum(2); | 
| 47 | 
> | 
    ((snapshotMan_->getCurrentSnapshot())->*storage_).unitVector[localIndex_] = a.inverse() * sU_.getColumn(2); | 
| 48 | 
  | 
}     | 
| 49 | 
  | 
     | 
| 50 | 
  | 
void DirectionalAtom::setA(const RotMat3x3d& a, int snapshotNo) { | 
| 51 | 
  | 
    ((snapshotMan_->getSnapshot(snapshotNo))->*storage_).aMat[localIndex_] = a; | 
| 52 | 
< | 
    ((snapshotMan_->getSnapshot(snapshotNo))->*storage_).unitVector[localIndex_] = a.inverse() * sU_.getColum(2);     | 
| 52 | 
> | 
    ((snapshotMan_->getSnapshot(snapshotNo))->*storage_).unitVector[localIndex_] = a.inverse() * sU_.getColumn(2);     | 
| 53 | 
  | 
}     | 
| 54 | 
  | 
 | 
| 55 | 
  | 
void DirectionalAtom::rotateBy(const RotMat3x3d& m) { |