| 263 |
|
} |
| 264 |
|
} |
| 265 |
|
|
| 266 |
< |
void GB::calcForce(InteractionData idat) { |
| 266 |
> |
void GB::calcForce(InteractionData &idat) { |
| 267 |
|
|
| 268 |
|
if (!initialized_) initialize(); |
| 269 |
|
|
| 357 |
|
Vector3d rxu2 = cross(idat.d, ul2); |
| 358 |
|
Vector3d uxu = cross(ul1, ul2); |
| 359 |
|
|
| 360 |
< |
idat.pot += U*idat.sw; |
| 360 |
> |
idat.pot[0] += U*idat.sw; |
| 361 |
|
idat.f1 += dUdr * rhat + dUda * ul1 + dUdb * ul2; |
| 362 |
|
idat.t1 += dUda * rxu1 - dUdg * uxu; |
| 363 |
|
idat.t2 += dUdb * rxu2 - dUdg * uxu; |
| 364 |
< |
idat.vpair += U*idat.sw; |
| 364 |
> |
idat.vpair[0] += U*idat.sw; |
| 365 |
|
|
| 366 |
|
return; |
| 367 |
|
|
