<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
</style>
</head>
<body class='hmmessage'>
Hi,<br><br>The error should not be in the tables, if you are using Gromacs 4.<br>We made sure that the accuracy of the tabulation and interpolation<br>is nearly the full precision (both in single and double precision).<br><br>For these kind of testing purposes you should be using double precision<br>anyhow, which gives you far more accuracy than 10^-4.<br><br>Berk<br><br>&gt; Date: Tue, 6 Jan 2009 14:24:17 +0100<br>&gt; From: dommert@fias.uni-frankfurt.de<br>&gt; To: gmx-users@gromacs.org<br>&gt; Subject: Re: [gmx-users] force evaluation<br>&gt; <br>&gt; * Mark Abraham &lt;Mark.Abraham@anu.edu.au&gt; [2009-01-07 00:13:00 +1100]:<br>&gt; <br>&gt; &gt; Suman Chakrabarty wrote:<br>&gt; &gt;&gt; On Tue, Jan 6, 2009 at 2:29 PM, Florian Dommert<br>&gt; &gt;&gt; &lt;dommert@fias.uni-frankfurt.de&gt; wrote:<br>&gt; &gt;&gt;&gt;  Finally it would be very interested in a possiblity to switch off the<br>&gt; &gt;&gt;&gt;  vdw interactions easily. So far I used an appropiate topology files with <br>&gt; &gt;&gt;&gt; C6<br>&gt; &gt;&gt;&gt; and<br>&gt; &gt;&gt;&gt;  C12 set to zero. Is there an easier solution as using a table filled<br>&gt; &gt;&gt;&gt;  with zeros ?<br>&gt; &gt;<br>&gt; &gt; Zeroed tables will be considerably slower than zeroed parameters.<br>&gt; &gt;<br>&gt; &gt; Depending on the way the force field is organized, you might create a <br>&gt; &gt; version of the force field files that zeroes out these parameters. Another <br>&gt; &gt; possibility might be to using the preprocessor mechanism to define the VDW <br>&gt; &gt; parameters only under some circumstances.<br>&gt; <br>&gt; Thank you very much for the help. It seems I went the a good way with<br>&gt; using a preprocessing command to define the zero VDW IA.<br>&gt; <br>&gt; Now I have to find where the difference in the calculation of the<br>&gt; electrostatic forces. Perhaps somebody can tell me, where the forces get<br>&gt; tabulated, that I can take a look at the code there ?<br>&gt; <br>&gt; Cheers,<br>&gt; <br>&gt; Flo<br>&gt; <br>&gt; <br>&gt; &gt;<br>&gt; &gt; Mark<br>&gt; &gt; _______________________________________________<br>&gt; &gt; gmx-users mailing list    gmx-users@gromacs.org<br>&gt; &gt; http://www.gromacs.org/mailman/listinfo/gmx-users<br>&gt; &gt; Please search the archive at http://www.gromacs.org/search before posting!<br>&gt; &gt; Please don't post (un)subscribe requests to the list. Use the www interface <br>&gt; &gt; or send it to gmx-users-request@gromacs.org.<br>&gt; &gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php<br>&gt; <br>&gt; -- <br>&gt; Florian Dommert<br>&gt; Dipl.-Phys.<br>&gt; <br>&gt; Computational and Theoretical Softmatter &amp; Biophysics group<br>&gt; <br>&gt; Frankfurt Institute for Advanced Studies<br>&gt; Johann-Wolfgang-Goethe University<br>&gt; <br>&gt; Ruth-Moufang-Str. 1<br>&gt; 60438 Frankfurt am Main<br>&gt; <br>&gt; Phone: +49(0)69 / 798 - 47529<br>&gt; Fax:   +49(0)69 / 798 - 47611<br>&gt; <br>&gt; EMail: dommert@fias.uni-frankfurt.de<br>&gt; Home: http://fias.uni-frankfurt.de/~simbio/Florian_Dommert<br><br /><hr />Express yourself instantly with MSN Messenger! <a href='http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/' target='_new'>MSN Messenger</a></body>
</html>