Yep. There are lots of known bugs in old compilers. It's not an accident that the GROMACS installation instructions tell you to get up to date compilers - both for correctness and performance of the resulting code. 4.4.7 might be correct (if you can run the 4.6.3 regression tests on it, it might pass), but I would bet 4.7 or 4.8 will be faster.
Mark On Wed, Sep 11, 2013 at 3:28 PM, Xu, Jianqing <x...@medimmune.com> wrote: > Thanks again for the message! > > I am using the Gromacs 4.6.3. > > Hmm, indeed a lot of information from "mdrun -version"..... > > Both servers used gcc, and they are in different versions. > Server 1 (AMD) has gcc 4.1.2 > Server 2 (Intel) has gcc 4.4.7 > > I know the compilers are a bit old, but can this cause big difference in > results? > > Thanks a lot, > > Jianqing > > > ________________________________________ > From: gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] on behalf > of Mark Abraham [mark.j.abra...@gmail.com] > Sent: 11 September 2013 14:13 > To: Discussion list for GROMACS users > Subject: Re: [gmx-users] Gromacs on Intel Xeon and AMD Opteron > > Hmm, that looks like a problem. Which GROMACS version, and which > compiler version? (See mdrun -version) > > On Wed, Sep 11, 2013 at 3:06 PM, Xu, Jianqing <x...@medimmune.com> wrote: >> >> Thanks for the comments, Mark. >> >> Sorry that I did not explain my question clear enough. >> >> I am actually talking about post-equilibration time average. I understand >> that some instantaneous deviations are expected, but the simulations should >> be correct as long as the ensemble averages are the same (or close enough), >> but the differences that I observed seemed beyond that. >> >> I firstly realized the total energy reported from two servers are obviously >> different (again, two simulations started from the same *.tpr file, and the >> executable on each server were compiled independently using the same >> arguments). I then checked the potential energy and kinetic energy as shown >> the below link: >> >> https://docs.google.com/file/d/0B7MGQ6N_2Ln7emRySENsZU40d2s/edit?usp=sharing >> >> >> Looks like the differences are from the potential energy. I then took a >> quick look at the LJ (SR) and Coulomb (SR) components, and it looks like the >> Coulomb (SR) component has big deviation which beyond the numerical errors. >> There are other differences, but looks like the Coulomb (SR) is the most >> obvious one. I was not able to understand how this can happen. >> >> https://docs.google.com/file/d/0B7MGQ6N_2Ln7bzB4YVhvMjdsNkE/edit?usp=sharing >> >> I hope the figures can clarify my situation better, and I will appreciate if >> you or anyone else on the list could provide more suggestions. >> >> Thanks a lot! >> >> Regards, >> >> Jianqing >> >> >> >> ________________________________________ >> From: gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] on >> behalf of Mark Abraham [mark.j.abra...@gmail.com] >> Sent: 11 September 2013 10:56 >> To: Discussion list for GROMACS users >> Subject: Re: [gmx-users] Gromacs on Intel Xeon and AMD Opteron >> >> It's not clear whether you are reporting single points or >> post-equilibration time averages, but either way you must expect >> differences. MD simulations are chaotic. However, the long-time >> ensemble averages should agree well - that's the point of the >> simulation. >> >> Mark >> >> On Wed, Sep 11, 2013 at 10:35 AM, Xu, Jianqing <x...@medimmune.com> wrote: >>> Dear all, >>> >>> This is my first time here. I apologize if I am not aware any rules for >>> posting a new message. >>> >>> I was testing Gromacs on two servers, one server has AMD Opteron processor >>> (Server 1), and the other has Intel Xeon processor (Server 2). The >>> simulation system that I tried is the lysozyme system from the online >>> tutorial >>> (http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/lysozyme/). >>> Briefly, it has a lysozyme in 10824 water, and 8 Cl-. Box size is around >>> 6.99 nm. >>> >>> I finished all the preparation steps advised by the tutorials, and tested >>> the final NPT production runs on two servers, using the same final input >>> file "md_0_1.tpr". It is strange that the results from Server 1 (AMD) has >>> higher potential energy than the results from Server2 (Intel), >>> specifically, results from Server1 (AMD) has higher "Coulomb-(SR)" than >>> Server2 (Intel) ( ~ -5.5*10^5 .vs. ~ -5.7*10^5 kj/mol ). >>> >>> I was wondering if anyone could give me any suggestions about what could be >>> the reasons causing this discrepancy. >>> >>> Some additional information: >>> 1). Jobs run on my local computer, which also has Intel CPUs (Intel(R) >>> Xeon(R) CPU E5620 @ 2.40GHz), reported the same results as the Server2 >>> (Intel(R) Xeon(R) CPU X7550 @ 2.00GHz). >>> 2). All installations (Gromacs and OpenMPI) on Server 1 and 2 were compiled >>> independently, but with the same arguments. >>> >>> >>> I understand that the information I provided may not be sufficient. If I >>> should provide more specific information, please let me know. >>> >>> Any help will be highly appreciated! >>> >>> Thanks, >>> >>> Jianqing >>> To the extent this electronic communication or any of its attachments >>> contain information that is not in the public domain, such information is >>> considered by MedImmune to be confidential and proprietary. This >>> communication is expected to be read and/or used only by the individual(s) >>> for whom it is intended. If you have received this electronic communication >>> in error, please reply to the sender advising of the error in transmission >>> and delete the original message and any accompanying documents from your >>> system immediately, without copying, reviewing or otherwise using them for >>> any purpose. Thank you for your cooperation. >>> -- >>> gmx-users mailing list gmx-users@gromacs.org >>> http://lists.gromacs.org/mailman/listinfo/gmx-users >>> * Please search the archive at >>> http://www.gromacs.org/Support/Mailing_Lists/Search before posting! >>> * Please don't post (un)subscribe requests to the list. Use the >>> www interface or send it to gmx-users-requ...@gromacs.org. >>> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists >> -- >> gmx-users mailing list gmx-users@gromacs.org >> http://lists.gromacs.org/mailman/listinfo/gmx-users >> * Please search the archive at >> http://www.gromacs.org/Support/Mailing_Lists/Search before posting! >> * Please don't post (un)subscribe requests to the list. Use the >> www interface or send it to gmx-users-requ...@gromacs.org. >> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists >> To the extent this electronic communication or any of its attachments >> contain information that is not in the public domain, such information is >> considered by MedImmune to be confidential and proprietary. This >> communication is expected to be read and/or used only by the individual(s) >> for whom it is intended. If you have received this electronic communication >> in error, please reply to the sender advising of the error in transmission >> and delete the original message and any accompanying documents from your >> system immediately, without copying, reviewing or otherwise using them for >> any purpose. Thank you for your cooperation. >> -- >> gmx-users mailing list gmx-users@gromacs.org >> http://lists.gromacs.org/mailman/listinfo/gmx-users >> * Please search the archive at >> http://www.gromacs.org/Support/Mailing_Lists/Search before posting! >> * Please don't post (un)subscribe requests to the list. Use the >> www interface or send it to gmx-users-requ...@gromacs.org. >> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists > -- > gmx-users mailing list gmx-users@gromacs.org > http://lists.gromacs.org/mailman/listinfo/gmx-users > * Please search the archive at > http://www.gromacs.org/Support/Mailing_Lists/Search before posting! > * Please don't post (un)subscribe requests to the list. Use the > www interface or send it to gmx-users-requ...@gromacs.org. > * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists > To the extent this electronic communication or any of its attachments contain > information that is not in the public domain, such information is considered > by MedImmune to be confidential and proprietary. This communication is > expected to be read and/or used only by the individual(s) for whom it is > intended. If you have received this electronic communication in error, please > reply to the sender advising of the error in transmission and delete the > original message and any accompanying documents from your system immediately, > without copying, reviewing or otherwise using them for any purpose. Thank you > for your cooperation. > -- > gmx-users mailing list gmx-users@gromacs.org > http://lists.gromacs.org/mailman/listinfo/gmx-users > * Please search the archive at > http://www.gromacs.org/Support/Mailing_Lists/Search before posting! > * Please don't post (un)subscribe requests to the list. Use the > www interface or send it to gmx-users-requ...@gromacs.org. > * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists -- gmx-users mailing list gmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! * Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists