<br /><br /><span>On 01/26/11, <b class="name">TJ Mustard </b> &lt;mustardt@onid.orst.edu&gt; wrote:</span><blockquote cite="mid: &lt;546997728.25225.1295998899004.JavaMail.open-xchange@oxusgw01.lxa.perfora.net" class="iwcQuote" style="border-left: 1px solid rgb(0, 0, 255); padding-left: 13px; margin-left: 0pt;" type="cite"><div class="mimepart text html">

<span xmlns="http://www.w3.org/1999/xhtml"><p>
  
    <meta content="text/html; charset=UTF-8" http-equiv="Content-Type" />
    <title></title>
  <table><tbody><tr><td><p>
    </p><p style="margin: 0px;"></p>

    <div style="margin: 5px 0px;">
      On January 25, 2011 at 3:24 PM &quot;Justin A. Lemkul&quot; &lt;jalemkul@vt.edu&gt; wrote:<br />
      <br />
      &gt;<br />
      &gt;<br />
      &gt; TJ Mustard wrote:<br />
      &gt; &gt;<br />
      &gt; &gt;<br />
      &gt; &gt; <br />
      &gt; &gt;<br />
      &gt; &gt;<br />
      &gt; &gt; On January 25, 2011 at 2:08 PM Mark Abraham &lt;Mark.Abraham@anu.edu.au&gt; wrote:<br />
      &gt; &gt;<br />
      &gt; &gt;&gt; On 26/01/2011 5:50 AM, TJ Mustard wrote:<br />
      &gt; &gt;&gt;&gt;<br />
      &gt; &gt;&gt;&gt; Hi all,<br />
      &gt; &gt;&gt;&gt;<br />
      &gt; &gt;&gt;&gt; <br />
      &gt; &gt;&gt;&gt;<br />
      &gt; &gt;&gt;&gt; I am running MD/FEP on a protein-ligand system with gromacs 4.5.3 and<br />
      &gt; &gt;&gt;&gt; FFTW 3.2.2.<br />
      &gt; &gt;&gt;&gt;<br />
      &gt; &gt;&gt;&gt; <br />
      &gt; &gt;&gt;&gt;<br />
      &gt; &gt;&gt;&gt; My iMac will run the job (over 4000 steps, till I killed it) at 4fs<br />
      &gt; &gt;&gt;&gt; steps. (I am using heavy H)<br />
      &gt; &gt;&gt;&gt;<br />
      &gt; &gt;&gt;&gt; <br />
      &gt; &gt;&gt;&gt;<br />
      &gt; &gt;&gt;&gt; Once I put this on our groups AMD Cluster the jobs fail even with 2fs<br />
      &gt; &gt;&gt;&gt; steps. (with thousands of lincs errors)<br />
      &gt; &gt;&gt;&gt;<br />
      &gt; &gt;&gt;&gt; <br />
      &gt; &gt;&gt;&gt;<br />
      &gt; &gt;&gt;&gt; We have recompiled the clusters gromacs 4.5.3 build, with no change.<br />
      &gt; &gt;&gt;&gt; I know the system is the same since I copied the job from the server<br />
      &gt; &gt;&gt;&gt; to my machine, to rerun it.<br />
      &gt; &gt;&gt;&gt;<br />
      &gt; &gt;&gt;&gt; <br />
      &gt; &gt;&gt;&gt;<br />
      &gt; &gt;&gt;&gt; What is going on? Why can one machine run a job perfectly and the<br />
      &gt; &gt;&gt;&gt; other cannot? I also know there is adequate memory on both machines.<br />
      &gt; &gt;&gt;&gt;<br />
      &gt; &gt;&gt;<br />
      &gt; &gt;&gt; You've posted this before, and I made a number of diagnostic<br />
      &gt; &gt;&gt; suggestions. What did you learn?<br />
      &gt; &gt;&gt;<br />
      &gt; &gt;&gt; Mark<br />
      &gt; &gt;<br />
      &gt; &gt; Mark and all,<br />
      &gt; &gt;<br />
      &gt; &gt; <br />
      &gt; &gt;<br />
      &gt; &gt; First thank you for all our help. What you suggested last time helped<br />
      &gt; &gt; considerably with our jobs/calculations. I have learned that using the<br />
      &gt; &gt; standard mdp settings allow my heavyh 4fs jobs to run on my iMac (intel)<br />
      &gt; &gt; and have made these my new standard for future jobs. We chose to use the<br />
      &gt; &gt; smaller 0.8nm PME/Cutoff due to others papers/tutorials, but now we<br />
      &gt; &gt; understand why we need these standard settings. Now what I see to be our<br />
      &gt; &gt; problem is that our machines have some sort of variable we cannot<br />
      &gt; &gt; account for. If I am blind to my error, please show me. I just don't<br />
      &gt; &gt; understand why one computer works while the other does not. We have<br />
      &gt; &gt; recompiled gromacs 4.5.3 single precission on our cluster, and still<br />
      &gt; &gt; have this problem.<br />
      &gt; &gt;<br />
      &gt;<br />
      &gt; I know the feeling all too well.  PowerPC jobs crash instantly, on our cluster,<br />
      &gt; despite working beautifully on our lab machines.  There's a bug report about<br />
      &gt; that one, but I haven't heard anything about AMD failures.  It remains a<br />
      &gt; possibility that something beyond your control is going on.  To explore a bit<br />
      &gt; further:<br />
      &gt;<br />
      &gt; 1. Do the systems in question crash immediately (i.e., step zero) or do they run<br />
      &gt; for some time?<br />
      &gt;
    </div>

    <p style="margin: 0px;">Step 0, every time.</p>

    <p style="margin: 0px;"> </p>

    <div style="margin: 5px 0px;">
      &gt; 2. If they give you even a little bit of output, you can analyze which energy<br />
      &gt; terms, etc go haywire with the tips listed here:<br />
      &gt;
    </div>

    <p style="margin: 0px;">All I have seen on these is LINCS Errors and Water molecules unable to be settled.</p>

    <p style="margin: 0px;"> </p>

    <p style="margin: 0px;">But I will check this out right now, and email if I smell trouble.</p>

    <p style="margin: 0px;"> </p>

    <div style="margin: 5px 0px;">
      &gt; http://www.gromacs.org/Documentation/Terminology/Blowing_Up#Diagnosing_an_Unstable_System<br />
      &gt;<br />
      &gt; That would help in tracking down any potential bug or error.<br />
      &gt;<br />
      &gt; 3. Is it just the production runs that are crashing, or everything?  If EM isn't<br />
      &gt; even working, that smells even buggier.
    </div>

    <p style="margin: 0px;">Awesome question here, we have seen some weird stuff. Sometimes the cluster will give us segmentation faults, then it will fail on our machines or sometimes not on our iMacs. I know weird! If EM starts on the cluster it will finish. Where we have issues is in positional restraint (PR) and MD and MD/FEP. It doesn't matter if FEP is on or off in a MD (although we are using SD for these MD/FEP runs).</p>

    </td></tr></tbody></table></p></span></div></blockquote>Good. That rules out FEP as the source of the problem, like I asked in your previous thread.<br /> <blockquote cite="mid: &lt;546997728.25225.1295998899004.JavaMail.open-xchange@oxusgw01.lxa.perfora.net" class="iwcQuote" style="border-left: 1px solid rgb(0, 0, 255); padding-left: 13px; margin-left: 0pt;" type="cite"><div class="mimepart text html"><span><p><table><tbody><tr><td>

    <div style="margin: 5px 0px;">
      &gt;<br />
      &gt; 4. Are the compilers the same on the iMac vs. AMD cluster?
    </div>

    <p style="margin: 0px;">No I am using x86_64-apple-darwin10 GCC 4.4.4 and the cluster is using x86_64-redhat-linux 4.1.2 GCC.</p>

    <p style="margin: 0px;">I just did a quick yum search and there 
doesn't seem to be a newer GCC. We know you are going to cmake but we 
have yet to get it implemented on our cluster successfully.</p></td></tr></tbody></table></p></span></div></blockquote><br />There
 have been doubts about the 4.1.x series of GCC compilers for GROMACS - 
and IIRC 4.1.2 in particular (do search the archives yourself). Some 
time back, Berk solicited actual accounts of problems and nobody 
presented one. So we no longer have an official warning against using 
it. However I'd say this is a candidate for the source of your problems.
 I would ask your cluster admins to get and compile a source-code version of GCC for you 
to try.<br /><br />Mark<br /><blockquote cite="mid: &lt;546997728.25225.1295998899004.JavaMail.open-xchange@oxusgw01.lxa.perfora.net" class="iwcQuote" style="border-left: 1px solid rgb(0, 0, 255); padding-left: 13px; margin-left: 0pt;" type="cite"><div class="mimepart text html"><span xmlns="http://www.w3.org/1999/xhtml"><p><table ><tbody><tr><td><p style="margin: 0px;"></p></td></tr></tbody></table>
</p></span>
</div></blockquote>