1186 |
|
idat.A1 = &(snap_->atomData.aMat[atom1]); |
1187 |
|
idat.A2 = &(snap_->atomData.aMat[atom2]); |
1188 |
|
} |
1189 |
+ |
|
1190 |
+ |
RealType ct = dot(idat.A1->getColumn(2), idat.A2->getColumn(2)); |
1191 |
|
|
1192 |
|
if (storageLayout_ & DataStorage::dslTorque) { |
1193 |
|
idat.t1 = &(snap_->atomData.torque[atom1]); |
1308 |
|
#endif |
1309 |
|
|
1310 |
|
RealType rList_ = (largestRcut_ + skinThickness_); |
1309 |
– |
RealType rl2 = rList_ * rList_; |
1311 |
|
Snapshot* snap_ = sman_->getCurrentSnapshot(); |
1312 |
|
Mat3x3d Hmat = snap_->getHmat(); |
1313 |
|
Vector3d Hx = Hmat.getColumn(0); |