[gmx-users] Re: RE: problem with g_bar

Berk Hess gmx3 at hotmail.com
Tue Mar 20 10:03:33 CET 2012


Hi,

Apparently I didn't try to compile it myself, there is a pointer dereference missing.
Pleae try the fix below.

Cheers,

Berk

--- a/src/tools/gmx_bar.c
+++ b/src/tools/gmx_bar.c
@@ -2049,7 +2049,7 @@ static void read_bar_xvg(char *fn, real *temp, lambda_t *lambda_head)
         gmx_fatal(FARGS,"File '%s' contains fewer than two columns", fn);
     }
 
-    if ( ( *temp != barsim->temp) && (*temp > 0) )
+    if ( !gmx_within_tol(*temp,barsim->temp,GMX_FLOAT_EPS) && (*temp > 0) )
     {
         gmx_fatal(FARGS,"Temperature in file %s different from earlier files or setting\n", fn);
     }


> Date: Tue, 20 Mar 2012 09:23:38 +0100
> From: tomkirch at mis.mpg.de
> To: gmx-users at gromacs.org
> Subject: [gmx-users] Re: RE: problem with g_bar
> 
> Hi Berk,
> 
> Thanks a lot for your answer.
> 
> I changed the line and recompiled Gromacs, sadly I got this error:
> 
> gmx_bar.c: In function "read_bar_xvg":
> gmx_bar.c:2052:5: error: incompatible type for argument 1 of 
> "gmx_within_tol"
> ../../include/maths.h:125:1: note: expected "double" but argument is of 
> type "real *"
> make[3]: *** [gmx_bar.lo] Error 1
> 
> Best,
> Tom
> 
> On 03/20/2012 03:15 AM, gmx-users-request at gromacs.org wrote:
> > ------------------------------ Message: 3 Date: Mon, 19 Mar 2012 
> > 18:34:31 +0100 From: Berk Hess <gmx3 at hotmail.com> Subject: RE: 
> > [gmx-users] problem with g_bar To: Discussion list for GROMACS users 
> > <gmx-users at gromacs.org> Message-ID: 
> > <COL113-W32C019BA58913FEC173E168E420 at phx.gbl> Content-Type: 
> > text/plain; charset="iso-8859-1" Hi, Yes, there is a problem with 
> > different temperature variables being single and double precision. 
> > Does the one line change below fix the problem? Cheers, Berk - if ( ( 
> > *temp != barsim->temp) && (*temp > 0) ) + if ( 
> > !gmx_within_tol(temp,barsim->temp,GMX_FLOAT_EPS) && (*temp > 0) )
> >> >  Date: Mon, 19 Mar 2012 13:57:51 +0100
> >> >  From:tomkirch at mis.mpg.de
> >> >  To:gmx-users at gromacs.org
> >> >  Subject: [gmx-users] problem with g_bar
> >> >  
> >> >  Dear all,
> >> >  
> >> >  Recently I posted my problem with g_bar on this mailing list but got no
> >> >  solution for it. After additional testing I decided to send an enhanced
> >> >  description of my problem.
> >> >  
> >> >  I am doing Free Energy calculation. The dhdl values are put in separate
> >> >  .xvg files. The routine g_bar is then started using:
> >> >  
> >> >  g_bar -b 100 -f md_*.xvg -o -oi -oh
> >> >  
> >> >  With the result:
> >> >  
> >> >  Program g_bar, VERSION 4.5.5
> >> >  Source code file: .../src/gromacs-4.5.5/src/tools/gmx_bar.c, line: 2054
> >> >  
> >> >  Fatal error:
> >> >  Temperature in file md_0.05.xvg different from earlier files or setting
> >> >  
> >> >  
> >> >  Checking the header of all .xvg files showed the same line for the
> >> >  temperature for all files:
> >> >  
> >> >  @ subtitle "T = 298.15 (K), \xl\f{} = 0.025"
> >> >  
> >> >  Changing the temperature manually to any integer (but still with .0) let
> >> >  g_bar run fine.
> >> >  
> >> >  sed -i 's/T\ \=\ 298.15/T\ \=\ 298.0/g' md*
> >> >  
> >> >  Do you have any suggestions what's wrong here?
> >> >  
> >> >  I'm not that familiar with C, but my suggestion would be an error with
> >> >  floating point arithmetic in line 2052 of gmx_bar.c
> >> >  
> >> >        if ( ( *temp != barsim->temp)&&  (*temp>  0) )
> >> >  
> >> >  where *temp and barsim->temp are compared. Maybe it has something to do
> >> >  with the input as long double (%lf) and statement "real *temp" in the
> >> >  function header?
> >> >  
> >> >  I appreciate any hints.
> >> >  
> >> >  Best regards
> >> >  Tom
> 
> -- 
> gmx-users mailing list    gmx-users at 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-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20120320/2ab8bd63/attachment.html>


More information about the gromacs.org_gmx-users mailing list