[gmx-developers] pme in release-4-6 with domain decomposition

[Dommert Florian]

Hi,

today I tested the new release-4-6 branch, but inconsistent
electrostatic forces arise.

As test system I used the spc216 water box and calculated the forces
with ewald and pme. If only a single processor is used the forces are
consistent. Also for multiple processes and particle decomposition the
forces are correct. However if domain decomposition used, the forces are
sometimes in the order of 10^102 (correct average maximum force 2402
kJ/mol/nm). 

If the test systems consists of 500 positive and negative atomic ions or
molecular ions, everything is fine, too, indepenent of the decomposition
scheme. 

Only the model systems with water make trouble (I tried SPC/E and
TIP-4P). 
Does anyone has an idea, what the reason for this problem might be ?

Cheers,
Flo



-- 
Florian Dommert
Dipl. - Phys.

Institute for Computational Physics
University Stuttgart

Pfaffenwaldring 27
70569 Stuttgart

EMail: dommert at icp.uni-stuttgart.de
Homepage: http://www.icp.uni-stuttgart.de/~icp/Florian_Dommert

Tel.: +49 - (0)711 - 68563613
Fax.: +49 - (0)711 - 68563658
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 198 bytes
Desc: This is a digitally signed message part
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-developers/attachments/20120409/b1f689bd/attachment.sig>