<br><br><div class="gmail_quote">On Mon, Sep 27, 2010 at 9:58 PM, Mark Abraham <span dir="ltr">&lt;<a href="mailto:mark.abraham@anu.edu.au">mark.abraham@anu.edu.au</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">

<div class="im"><br><br>----- Original Message -----<br>From: &quot;Justin A. Lemkul&quot; &lt;<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>&gt;<br></div><div class="im">Date: Tuesday, September 28, 2010 11:39<br>

Subject: Re: [gmx-users] Hardware-specific crash with 4.5.1<br></div><div><div></div><div class="h5">To: Discussion list for GROMACS users &lt;<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>&gt;<br>

<br>&gt; <br>&gt; <br>&gt; Mark Abraham wrote:<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt;----- Original Message -----<br>&gt; &gt;From: &quot;Justin A. Lemkul&quot; &lt;<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>&gt;<br>

&gt; &gt;Date: Tuesday, September 28, 2010 11:11<br>&gt; &gt;Subject: Re: [gmx-users] Hardware-specific crash with 4.5.1<br>&gt; &gt;To: Gromacs Users&#39; List &lt;<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>&gt;<br>

&gt; &gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Roland Schulz wrote:<br>&gt; &gt; &gt; &gt;Justin,<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt;I think the interaction kernel is not OK on your PowerPC<br>
&gt; &gt; &gt; machine. I assume that from: 1) The force seems to be zero<br>
&gt; &gt; &gt; (minimization output). 2) When you use the all-to-all kernel<br>&gt; &gt; &gt; which is not available for the powerpc kernel, it automatically<br>&gt; &gt; &gt; falls back to the C kernel and then it works.<br>

&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Sounds about right.<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; &gt;What is the kernel you are using? It should say in the log<br>&gt; &gt; &gt; file. Look for: &quot;Configuring single precision IBM Power6-<br>

&gt; &gt; &gt; specific Fortran kernels&quot; or &quot;Testing Altivec/VMX support&quot;<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; I&#39;m not finding either in the config.log - weird?<br>&gt; &gt;<br>&gt; &gt;You were meant to look in the mdrun.log for runtime <br>

&gt; confirmation of what kernels GROMACS has decided to use.<br>&gt; &gt; <br>&gt; <br>&gt; That seems entirely obvious, now that you mention it :)  <br>&gt; Conveniently, I find the following in the md.log file from the <br>

&gt; (failing) autoconf-assembled mdrun:<br>&gt; <br>&gt; Configuring nonbonded kernels...<br>&gt; Configuring standard C nonbonded kernels...<br>&gt; Testing Altivec/VMX support... present.<br>&gt; Configuring PPC/Altivec nonbonded kernels...<br>

&gt; <br>&gt; The (non)MPI CMake build shows the following:<br>&gt; <br>&gt; Configuring nonbonded kernels...<br>&gt; Configuring standard C nonbonded kernels...<br>&gt; <br>&gt; So it seems clear to me that autoconf built faulty nonbonded <br>

&gt; kernels, and CMake didn&#39;t.<br><br></div></div>OK, so assuming that PPC/Altivec kernels are supposed to be good for Mac (as they were in 4.0.x, I believe):<br><br>1) CMake doesn&#39;t detect that it should be using those kernels, and so appears to work, but does an inefficient run. autoconf detects that it should use those kernels, but the mdrun fails for reasons that are not yet clear.<div>

<div></div><div class="h5"><br> <br>&gt; &gt; &gt; &gt;You can also look in the config.h whether  GMX_POWER6<br>&gt; &gt; &gt; and/or GMX_PPC_ALTIVEC is set. I suggest you try to compile with<br>&gt; &gt; &gt; one/both of them deactivated and see whether that solves it.<br>

&gt; &gt; &gt; This will make it slower too. Thus if this is indeed the<br>&gt; &gt; &gt; problem, you will probably want to figure out why the fastest<br>&gt; &gt; &gt; kernel doesn&#39;t work correctly to get good performance.<br>

&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; It looks like GMX_PPC_ALTIVEC is set.  I suppose I <br>&gt; could re-<br>&gt; &gt; &gt; compile with this turned off.<br>&gt; &gt;<br>&gt; &gt;This is supposed to be fine for Mac, as I understand.<br>

&gt; &gt;<br>&gt; &gt; &gt; Here&#39;s what&#39;s even weirder.  The problematic version was<br>&gt; &gt; &gt; compiled using the standard autoconf procedure.  If I <br>&gt; use a<br>&gt; &gt; &gt; CMake-compiled version, the energy minimization runs fine,<br>

&gt; &gt; &gt; giving the same results (energy and force) as the two <br>&gt; systems I<br>&gt; &gt; &gt; know are good.  So I guess there&#39;s something wrong with the<br>&gt; &gt; &gt; way autoconf installed Gromacs.  Perhaps this isn&#39;t of<br>

&gt; &gt; &gt; concern since Gromacs will require CMake in subsequent releases,<br>&gt; &gt; &gt; but I figure I should at least report it in case it affects<br>&gt; &gt; &gt; anyone else.<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; If I may tack one more question on here, I&#39;m wondering why my<br>

&gt; &gt; &gt; CMake installation  doesn&#39;t actually appear to be using<br>&gt; &gt; &gt; MPI.  I get the right result, but the problem is, I get a<br>&gt; &gt; &gt; .log, .edr, and .trr for every processor that&#39;s being used, as<br>

&gt; &gt; &gt; if each processor is being given its own job and not<br>&gt; &gt; &gt; distributing the work. Here&#39;s how I compiled my MPI mdrun,<br>&gt; &gt; &gt; version 4.5.1:<br>&gt; &gt;<br>&gt; &gt;At the start and end of the .log files you should get <br>

&gt; indicators about how many MPI processes were actually being used.<br>&gt; &gt; <br>&gt; <br>&gt; That explains it (sort of).  It looks like mdrun thinks <br>&gt; it&#39;s only being run over 1 node, just several times over, and a <br>

&gt; bunch of junk that isn&#39;t getting written properly:<br>&gt; <br>&gt; Log file opened on Mon Sep 27 21:36:00 2010<br>&gt; Host: n235  pid: 6857  nodeid: 0  nnodes:  1<br>&gt; The Gromacs distribution was built @TMP_TIME@ by<br>

&gt; <a href="mailto:jalemkul@sysx2.arc-int.vt.edu" target="_blank">jalemkul@sysx2.arc-int.vt.edu</a> [CMAKE] (@TMP_MACHINE@)<br>&gt; <br>&gt; Frustrating.<br><br></div></div>You can set the GMX_NOOPTIMIZEDKERNELS environment variable with your autoconf build to see whether the MPI issue is CMake-dependent. Normally, I&#39;d say your supercomputer MPI environment isn&#39;t being invoked correctly, but presumably you already know how to do that right...<div class="im">

<br><br>&gt; &gt; &gt; cmake ../gromacs-4.5.1 -<br>&gt; DFFTW3F_LIBRARIES=/home/rdiv1001/fftw-<br>&gt; &gt; &gt; 3.0.1-osx/lib/libfftw3f.a -<br>&gt; &gt; &gt; DFFTW3F_INCLUDE_DIR=/home/rdiv1001/fftw-3.0.1-osx/include/ -<br>

&gt; &gt; &gt; DCMAKE_INSTALL_PREFIX=/home/rdiv1001/gromacs-4.5_cmake-osx -<br>&gt; &gt; &gt; DGMX_BINARY_SUFFIX=_4.5_cmake_mpi -DGMX_THREADS=OFF -<br>&gt; &gt; &gt; DBUILD_SHARED_LIBS=OFF -DGMX_X11=OFF -DGMX_MPI=ON -<br>

&gt; &gt; &gt; DMPI_COMPILER=/home/rdiv1001/compilers/openmpi-1.2.3-<br>&gt; &gt; &gt; osx/bin/mpicxx -<br>&gt; &gt; &gt; DMPI_INCLUDE_PATH=/home/rdiv1001/compilers/openmpi-1.2.3-<br>&gt; osx/include&gt; &gt;<br>&gt; &gt; &gt; $ make mdrun<br>

&gt; &gt; &gt;<br>&gt; &gt; &gt; $ make install-mdrun<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Is there anything obviously wrong with those commands?  Is<br>&gt; &gt; &gt; there any way I should know (before actually using mdrun)<br>

&gt; &gt; &gt; whether or not I&#39;ve done things right?<br>&gt; &gt;<br>&gt; &gt;I think there ought to be, but IMO not enough preparation and <br>&gt; testing has gone into the CMake switch for it to be usable.<br>&gt; &gt;<br>

&gt; <br>&gt; I agree.  After hours of hacking CMake to try to make it <br>&gt; work (and thinking I had gotten it squared away), the MPI <br>&gt; doesn&#39;t seem to function.  The &quot;old&quot; way of doing things <br>

&gt; worked flawlessly, except that somewhere between 4.0.7 and <br>&gt; 4.5.1, the nonbonded kernels that used to work on our <br>&gt; architecture somehow got hosed. So now I&#39;m in limbo.<br></div></blockquote><div>
Well that&#39;s why Autoconf is still supported for 4.5 - to allow a smooth transition. It is difficult to get enough feedback for cmake for just the development version. Only know with the cmake support in the released version their is enough feedback to make the uncommon cases work. But you can still use autoconf if cmake doesn&#39;t work. Thus it should work for everyone and we get the required feedback to use cmake for 5.0.</div>

<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Sounds like Bugzilla time.<br></blockquote><div><br></div><div>I agree you should file 3 bugs. </div><div>

1) CMake + MPI -&gt; MPI doesn&#39;t work, </div><div>2) CMake + Altivec -&gt; Altivec isn&#39;t detected</div><div>3) Altivec produces 0 forces </div><div><br></div><div>As far as I know the POWER6 kernel should work too for your hardware and should be a little bit faster than the standard C kernel. You might want to try that until the Altivec kernel is fixed.</div>

<div><br></div><div>Roland</div><div><br></div></div><br clear="all"><br>-- <br>ORNL/UT Center for Molecular Biophysics <a href="http://cmb.ornl.gov">cmb.ornl.gov</a><br>865-241-1537, ORNL PO BOX 2008 MS6309<br>