[gmx-users] Reading double-precision trr file using xdrfile

[James]

Dear Gromacs developers and users,

I would like to read a .trr trajectory file in double-precision after
running a simulation in double-precision, however reading the atomic
positions etc with the xdrfile appears to yield only single precision.

Are the trajectory position, velocity and forces not stored in
double-precision in the trr file? Certainly the file size is much larger,
so I imagine it's in double precision. xdrfile however is expecting only
floats:

  extern int read_trr(XDRFILE *xd,int natoms,int *step,float *t,float
*lambda,
              matrix box,rvec *x,rvec *v,rvec *f);

(rvec is defined as  typedef float rvec[DIM];)

Is there any alternative to get them in double precision?

Thank you for advice.

Best regards,

James