| 1292 |
|
eCos[l][i].z() = eCos[l-1][i].z()*eCos[2][i].z() - eSin[l-1][i].z()*eSin[2][i].z(); |
| 1293 |
|
|
| 1294 |
|
eSin[l][i].x() = eSin[l-1][i].x()*eCos[2][i].x() + eCos[l-1][i].x()*eSin[2][i].x(); |
| 1295 |
< |
eSin[l][i].y() = eSin[l-1][i].y()*eCos[2][i].y() + eCos[l-1][i].y()*eSin[2][i].z(); |
| 1296 |
< |
eSin[l][i].z() = eSin[l-1][i].z()*eCos[2][i].z() + eCos[l-1][i].z()*eSin[2][i].y(); |
| 1295 |
> |
eSin[l][i].y() = eSin[l-1][i].y()*eCos[2][i].y() + eCos[l-1][i].y()*eSin[2][i].y(); |
| 1296 |
> |
eSin[l][i].z() = eSin[l-1][i].z()*eCos[2][i].z() + eCos[l-1][i].z()*eSin[2][i].z(); |
| 1297 |
|
|
| 1298 |
|
|
| 1299 |
|
// a.Vmul(eCos[l-1][i], u); |
| 1462 |
|
MPI::SUM); |
| 1463 |
|
#endif |
| 1464 |
|
|
| 1465 |
< |
// Accumulate potential energy and virial contribution: |
| 1465 |
> |
// Accumulate potential energy and virial contribution: |
| 1466 |
|
|
| 1467 |
|
kPot += 2.0 * rvol * AK[kk]*((ckss+dkcs-qkss)*(ckss+dkcs-qkss) |
| 1468 |
|
+ (ckcs-dkss-qkcs)*(ckcs-dkss-qkss)); |
| 1491 |
|
-ckr[i]*(ckss+dkcs-qkss)); |
| 1492 |
|
RealType qtrq2 = 2.0*AK[kk]*(ckr[i]*(ckcs-dkss-qkcs)+ |
| 1493 |
|
skr[i]*(ckss+dkcs-qkss)); |
| 1494 |
< |
|
| 1494 |
> |
|
| 1495 |
|
atom->addFrc( 4.0 * rvol * qfrc * kVec ); |
| 1496 |
|
|
| 1497 |
|
if (data.is_Dipole) { |
