Flipped force sign in polar_real, made sure that multipole_real is true for precompute() to be invoked, ubdirect2b() is segfault and needs work
This commit is contained in:
@ -1609,12 +1609,12 @@ __kernel void k_amoeba_polar(const __global numtyp4 *restrict x_,
|
||||
numtyp vyz = (numtyp)0.5 * (zr*frcy+yr*frcz);
|
||||
numtyp vzz = zr * frcz;
|
||||
|
||||
virial[0] += vxx;
|
||||
virial[1] += vyy;
|
||||
virial[2] += vzz;
|
||||
virial[3] += vxy;
|
||||
virial[4] += vxz;
|
||||
virial[5] += vyz;
|
||||
virial[0] -= vxx;
|
||||
virial[1] -= vyy;
|
||||
virial[2] -= vzz;
|
||||
virial[3] -= vxy;
|
||||
virial[4] -= vxz;
|
||||
virial[5] -= vyz;
|
||||
}
|
||||
} // nbor
|
||||
|
||||
|
||||
Reference in New Issue
Block a user