[gmx-developers] npme broken in latest cvs?
hessb at mpip-mainz.mpg.de
hessb at mpip-mainz.mpg.de
Wed Jun 18 21:02:46 CEST 2008
Hi,
I fixed this bug in CVS.
This was introduced a week ago on June 10.
Berk.
> Hi,
>
> It is a bug I put in recently.
>
> A fast fix is:
> diff -r1.119 pme.c
> 608c608,609
> < rvec_inc(f[i],atc->f[local_pos]);
> ---
> > //rvec_inc(f[i],atc->f[local_pos]);
> > copy_rvec(atc->f[local_pos],f[i]);
> 612c613,614
> < rvec_inc(f[i],pme->bufv[buf_index[node]]);
> ---
> > //rvec_inc(f[i],pme->bufv[buf_index[node]]);
> > copy_rvec(pme->bufv[buf_index[node]],f[i]);
>
> I have to think a bit to understand if this is a proper solution (also
> working
> single processor) or not.
>
> Berk.
>
>
> Atte Sillanpää wrote:
>> Hi,
>>
>> I've tried compiling the latest cvs-versions (mdrun version 55) on two
>> separate Opteron clusters using gnu and pgi-compilers but I get
>> crashing code. The cluster that I refer to below has two dual core
>> cpu:s on each node. I've managed to reproduce this with any input
>> files I've tried. See below, for a gmxcheck between a two step dd and
>> dd:npme run.
>>
>> 1) If I request only domain decomposition, the job runs ok.
>> 2) If I specifically ask for npme with the mdrun flag the jobs either
>> stall or crashe (see below).
>> 3) If I let mdrun to decide division to dd/npme (no dd/npme flags),
>> the code does not stall, but crashes within some tens of steps as
>> waters can't be settled etc.
>>
>> With the same configure options (nothing fancy, fftw 3.1.2) I get
>> working code (the source was checked out on 11th April). Any ideas
>> what could be the problem?
>>
>> Cheers,
>> Atte
>>
>> gmxcheck output between two runs with "-dd 2 2 1" and "-dd 2 1 1 -npme
>> 2":
>> ...
>> comparing energy file e_pme.edr and e_dd.edr
>>
>> Opened e_pme.edr as single precision energy file
>> Opened e_dd.edr as single precision energy file
>> There are 162 terms in the energy files
>>
>> There are 162 terms to compare in the energy files
>>
>> Reading energy frame 1 time 0.000 Vir-ZY
>> step 0: 1058.05, step 0: 1056.45
>> Pres-ZY (bar) step 0: -95.3608, step 0: -95.2099
>> Reading energy frame 3 time 0.004 Kinetic En.
>> step 2: 58262.7, step 2: 57779.1
>> Total Energy step 2: -362651, step 2: -363088
>> Temperature step 2: 322.793, step 2: 320.114
>> Pressure (bar) step 2: 489.052, step 2: 148.929
>> pV step 2: 10406.1, step 2: 3168.91
>> Vir-XX step 2: 14573.3, step 2: 17932.7
>> Vir-XY step 2: -2524.57, step 2: -2370.18
>> Vir-XZ step 2: 2051.88, step 2: 1733.56
>> Vir-YX step 2: -2632.26, step 2: -2475.74
>> Vir-YY step 2: 16702.3, step 2: 19691.3
>> Vir-YZ step 2: 1213.29, step 2: 675.571
>> Vir-ZX step 2: 1978.46, step 2: 1733.89
>> Vir-ZY step 2: 1226.37, step 2: 612.046
>> Vir-ZZ step 2: 11378, step 2: 15401.7
>> Pres-XX (bar) step 2: 497.122, step 2: 167.411
>> Pres-XY (bar) step 2: 246.06, step 2: 232.524
>> Pres-XZ (bar) step 2: -179.978, step 2: -150.756
>> Pres-YX (bar) step 2: 256.182, step 2: 242.446
>> Pres-YY (bar) step 2: 263.995, step 2: -31.6074
>> Pres-YZ (bar) step 2: -112.03, step 2: -61.5908
>> Pres-ZX (bar) step 2: -173.078, step 2: -150.787
>> Pres-ZY (bar) step 2: -113.259, step 2: -55.6198
>> Pres-ZZ (bar) step 2: 706.039, step 2: 310.982
>> #Surf*SurfTen step 2: 2056.54, step 2: 1535.89
>> Coul-SR:CHOL-CHOL step 2: -0.00815666, step 2: -0.00819171
>> Coul-SR:CHOL-SOL step 2: -361.949, step 2: -362.405
>> Coul-SR:CHOL-Na step 2: 3.20663, step 2: 3.20127
>> Coul-SR:CHOL-Cl step 2: -0.0473165, step 2: -0.0468598
>> T-DPPC_POPC_POPG_CHOL step 2: 314.493, step 2: 313.836
>> T-SOL_Na_Cl step 2: 329.66, step 2: 325.308
>> Last energy frame read 3 time 0.004 Files read succesfully
>>
>> _______________________________________________
>> gmx-developers mailing list
>> gmx-developers at gromacs.org
>> http://www.gromacs.org/mailman/listinfo/gmx-developers
>> Please don't post (un)subscribe requests to the list. Use the www
>> interface or send it to gmx-developers-request at gromacs.org.
>>
>> This email was Anti Virus checked by Astaro Security Gateway.
>> http://www.astaro.com
>>
>
> _______________________________________________
> gmx-developers mailing list
> gmx-developers at gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-developers
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-developers-request at gromacs.org.
>
> This email was Anti Virus checked by Astaro Security Gateway.
> http://www.astaro.com
>
>
More information about the gromacs.org_gmx-developers
mailing list