<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">

<html xmlns="http://www.w3.org/1999/xhtml">
  <head>
    <meta content="text/html; charset=UTF-8" http-equiv="Content-Type" />
    <title></title>
  </head>

  <body>
    <p style="margin: 0px;">Found that if I set the setting of &quot;-dlb no&quot; when running the mdrun it would not fail. How could dynamic load balancing do this?</p>

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

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

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

    <div style="margin: 5px 0px 5px 0px;">
      <br />
       On September 28, 2010 at 3:21 AM TJ Mustard &lt;mustardt@onid.orst.edu&gt; wrote:<br />
      <br />

      <blockquote type="cite" style="margin-left: 0px; padding-left: 10px; border-left: solid 1px blue;">
        <p style="margin: 0px;">Hey all,</p>

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

        <p style="margin: 0px;">Here is another error that I keep getting. I am trying to &quot;speed up&quot; my md runs with -heavyh and longer time steps. I don&#39;t get LINCS errors but I do get this...</p>

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

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

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

        <p style="margin: 0px;">Back Off! I just backed up prlog.log to ./#prlog.log.1#<br />
         Getting Loaded...<br />
         Reading file monomer_pr.tpr, VERSION 4.5.1 (single precision)<br />
         Starting 2 threads<br />
         Loaded with Money<br />
        <br />
         Making 1D domain decomposition 2 x 1 x 1<br />
        <br />
         Back Off! I just backed up pr.edr to ./#pr.edr.1#<br />
         starting mdrun &#39;Protein in water&#39;<br />
         25000 steps,&#160;&#160;&#160; 100.0 ps.<br />
         step 900, will finish Mon Sep 27 18:19:07 2010imb F 18%<br />
         NOTE: Turning on dynamic load balancing<br />
        <br />
         step 1600, will finish Mon Sep 27 18:18:42 2010vol 0.92&#160; imb F&#160; 1%<br />
         A list of missing interactions:<br />
         &#160;&#160;&#160;&#160;&#160;&#160;&#160; LJC Pairs NB of&#160;&#160;&#160; 278 missing&#160;&#160;&#160;&#160;&#160; 1<br />
         &#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; exclusions of&#160;&#160; 6966 missing&#160;&#160;&#160;&#160;&#160; 1<br />
        <br />
         -------------------------------------------------------<br />
         Program g4.5.1-mdrun, VERSION 4.5.1<br />
         Source code file: domdec_top.c, line: 173<br />
        <br />
         Software inconsistency error:<br />
         Some interactions seem to be assigned multiple times<br />
         For more information and tips for troubleshooting, please check the GROMACS<br />
         website at http://www.gromacs.org/Documentation/Errors<br />
         -------------------------------------------------------<span></span></p>

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

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

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

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

        <p style="margin: 0px;">define&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = -DPOSRES</p>

        <p style="margin: 0px;">dt&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 0.004<br />
         nsteps&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 25000</p>

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

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

        <p style="margin: 0px;">; OPTIONS FOR ELECTROSTATICS AND VDW<br />
         ; Method for doing electrostatics<br />
         coulombtype&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = PME&#160; ; = Cutoff<br />
         rcoulomb-switch&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 0&#160;&#160;&#160; ; = 0<br />
         rcoulomb&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 0.9&#160; ; = 1<br />
         ; Relative dielectric constant for the medium and the reaction field<br />
         epsilon_r&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 1&#160;&#160;&#160; ; = 1<br />
         epsilon_rf&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 1&#160;&#160;&#160; ; = 1<br />
         ; Method for doing Van der Waals<br />
         vdw-type&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = Cut-off&#160;&#160;&#160;&#160;&#160; ; = Cut-off<br />
         ; cut-off lengths&#160;&#160;&#160;&#160;&#160;&#160;<br />
         rvdw-switch&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 0&#160;&#160;&#160; ; = 0<br />
         rvdw&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 1&#160;&#160;&#160; ; = 1<br />
         ; Apply long range dispersion corrections for Energy and Pressure<br />
         DispCorr&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = No&#160;&#160; ; = No<br />
         ; Extension of the potential lookup tables beyond the cut-off<br />
         table-extension&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 1&#160;&#160;&#160; ; = 1<br />
         ; Seperate tables between energy group pairs<br />
         energygrp_table&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; =&#160;&#160;&#160;&#160;&#160; ; =<br />
         ; Spacing for the PME/PPPM FFT grid<br />
         fourierspacing&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 0.12 ; = 0.12<br />
         ; FFT grid size, when a value is 0 fourierspacing will be used<br />
         fourier_nx&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 0&#160;&#160;&#160; ; = 0<br />
         fourier_ny&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 0&#160;&#160;&#160; ; = 0<br />
         fourier_nz&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 0&#160;&#160;&#160; ; = 0<br />
         ; EWALD/PME/PPPM parameters<br />
         pme_order&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 6&#160;&#160;&#160; ; = 4<br />
         ewald_rtol&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 1e-05&#160;&#160;&#160;&#160;&#160;&#160;&#160; ; = 1e-05<br />
         ewald_geometry&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 3d&#160;&#160; ; = 3d<br />
         epsilon_surface&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 0&#160;&#160;&#160; ; = 0<br />
         optimize_fft&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = yes&#160; ; = no</p>

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

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

        <p style="margin: 0px;">; OPTIONS FOR WEAK COUPLING ALGORITHMS<br />
         ; Temperature coupling&#160;<br />
         tcoupl&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = v-rescale&#160;&#160;&#160; ; = No<br />
         nsttcouple&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = -1&#160;&#160; ; = -1<br />
         nh-chain-length&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 10&#160;&#160; ; = 10<br />
         ; Groups to couple separately<br />
         tc-grps&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = RNA SOL&#160;&#160;&#160;&#160;&#160; ; =<br />
         ; Time constant (ps) and reference temperature (K)<br />
         tau-t&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 0.1 0.1&#160;&#160;&#160;&#160;&#160; ; =<br />
         ref-t&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 300 300&#160;&#160;&#160;&#160;&#160; ; =<br />
         ; Pressure coupling&#160;&#160;&#160;&#160;<br />
         Pcoupl&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = Parrinello-Rahman&#160;&#160;&#160; ; = No<br />
         Pcoupltype&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = Isotropic<br />
         nstpcouple&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = -1&#160;&#160; ; = -1<br />
         ; Time constant (ps), compressibility (1/bar) and reference P (bar)<br />
         tau-p&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 1&#160;&#160;&#160; ; = 1<br />
         compressibility&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 4.5e-5&#160;&#160;&#160;&#160;&#160;&#160; ; =<br />
         ref-p&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 1.0&#160; ; =<br />
         ; Scaling of reference coordinates, No, All or COM<br />
         refcoord_scaling&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = No&#160;&#160; ; = No<br />
         ; Random seed for Andersen thermostat<br />
         andersen_seed&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 815131&#160;&#160;&#160;&#160;&#160;&#160; ; = 815131</p>

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

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

        <p style="margin: 0px;">gen-vel&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = yes&#160; ; = no<br />
         gen-temp&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 300&#160; ; = 300<br />
         gen-seed&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = 173529&#160;&#160;&#160;&#160;&#160;&#160; ; = 173529</p>

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

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

        <p style="margin: 0px;">constraints&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; = all-bonds</p>

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

        <p style="margin: 0px;">Any help would be appreciated. It also seems to be intermittent as I have 21 identical runs (with different lambda values) and some work and some don&#39;t. It also changes every time I run them.</p>
        <br />
         Thank you,<br />

        <p style="font-family: monospace; white-space: nowrap; margin: 5px 0px 5px 0px;">TJ Mustard<br />
         Email: mustardt@onid.orst.edu</p>
      </blockquote>
    </div>

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

    <p style="font-family: monospace; white-space: nowrap; margin: 5px 0px 5px 0px;">TJ Mustard<br />
    Email: mustardt@onid.orst.edu</p>
  </body>
</html>