| 40 |
|
* [5] Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011). |
| 41 |
|
*/ |
| 42 |
|
|
| 43 |
+ |
#ifdef IS_MPI |
| 44 |
+ |
#include <mpi.h> |
| 45 |
+ |
#endif |
| 46 |
+ |
|
| 47 |
|
#include <stdio.h> |
| 48 |
|
#include <string.h> |
| 49 |
|
|
| 61 |
|
#include "math/erfc.hpp" |
| 62 |
|
#include "math/SquareMatrix.hpp" |
| 63 |
|
#include "primitives/Molecule.hpp" |
| 60 |
– |
#ifdef IS_MPI |
| 61 |
– |
#include <mpi.h> |
| 62 |
– |
#endif |
| 64 |
|
|
| 65 |
|
namespace OpenMD { |
| 66 |
|
|
| 1441 |
|
Q = atom->getQuadrupole() * mPoleConverter; |
| 1442 |
|
Qk = Q * kVec; |
| 1443 |
|
qk = dot(kVec, Qk); |
| 1444 |
< |
qxk[i] = cross(kVec, Qk); |
| 1444 |
> |
qxk[i] = -cross(kVec, Qk); |
| 1445 |
|
qkc[i] = qk * ckr[i]; |
| 1446 |
|
qks[i] = qk * skr[i]; |
| 1447 |
|
} |