| 374 |
|
Vector3d uxu = cross(ul1, ul2); |
| 375 |
|
|
| 376 |
|
(*(idat.pot))[VANDERWAALS_FAMILY] += U * *(idat.sw); |
| 377 |
< |
*(idat.f1) += dUdr * rhat + dUda * ul1 + dUdb * ul2; |
| 378 |
< |
*(idat.t1) += dUda * rxu1 - dUdg * uxu; |
| 379 |
< |
*(idat.t2) += dUdb * rxu2 + dUdg * uxu; |
| 380 |
< |
*(idat.vpair) += U * *(idat.sw); |
| 377 |
> |
*(idat.f1) += (dUdr * rhat + dUda * ul1 + dUdb * ul2) * *(idat.sw); |
| 378 |
> |
*(idat.t1) += (dUda * rxu1 - dUdg * uxu) * *(idat.sw); |
| 379 |
> |
*(idat.t2) += (dUdb * rxu2 + dUdg * uxu) * *(idat.sw); |
| 380 |
> |
*(idat.vpair) += U; |
| 381 |
|
|
| 382 |
|
return; |
| 383 |
|
|
