src/gmxlib/xdrd.c

   1 /*
   2  *
   3  *                This source code is part of
   4  *
   5  *                 G   R   O   M   A   C   S
   6  *
   7  *          GROningen MAchine for Chemical Simulations
   8  *
   9  *                        VERSION 3.2.0
  10  * Written by David van der Spoel, Erik Lindahl, Berk Hess, and others.
  11  * Copyright (c) 1991-2000, University of Groningen, The Netherlands.
  12  * Copyright (c) 2001-2004, The GROMACS development team,
  13  * check out http://www.gromacs.org for more information.
  14
  15  * This program is free software; you can redistribute it and/or
  16  * modify it under the terms of the GNU General Public License
  17  * as published by the Free Software Foundation; either version 2
  18  * of the License, or (at your option) any later version.
  19  *
  20  * If you want to redistribute modifications, please consider that
  21  * scientific software is very special. Version control is crucial -
  22  * bugs must be traceable. We will be happy to consider code for
  23  * inclusion in the official distribution, but derived work must not
  24  * be called official GROMACS. Details are found in the README & COPYING
  25  * files - if they are missing, get the official version at www.gromacs.org.
  26  *
  27  * To help us fund GROMACS development, we humbly ask that you cite
  28  * the papers on the package - you can find them in the top README file.
  29  *
  30  * For more info, check our website at http://www.gromacs.org
  31  *
  32  * And Hey:
  33  * GROningen Mixture of Alchemy and Childrens' Stories
  34  */
  35 #ifdef HAVE_CONFIG_H
  36 #include <config.h>
  37 #endif
  38
  39 #include "typedefs.h"
  40 #include "xdrf.h"
  41 #include "gmx_fatal.h"
  42 #include "smalloc.h"
  43
  44 int xdr_real(XDR *xdrs,real *r)
  45 {
  46 #ifdef GMX_DOUBLE
  47   float f;
  48   int   ret;
  49
  50   f=*r;
  51   ret=xdr_float(xdrs,&f);
  52   *r=f;
  53
  54   return ret;
  55 #else
  56   return xdr_float(xdrs,(float *)r);
  57 #endif
  58 }
  59
  60 int xdr3drcoord(XDR *xdrs, real *fp, int *size, real *precision)
  61 {
  62 #ifdef GMX_DOUBLE
  63   float *ffp;
  64   float  fprec;
  65   int    i,ret,isize;
  66
  67   isize=*size*DIM;
  68   if (isize <= 0)
  69     gmx_fatal(FARGS,"Don't know what to malloc for ffp, isize = %d",isize);
  70
  71   snew(ffp,isize);
  72
  73   for(i=0; (i<isize); i++)
  74     ffp[i]=fp[i];
  75   fprec=*precision;
  76   ret=xdr3dfcoord(xdrs,ffp,size,&fprec);
  77
  78   *precision=fprec;
  79   for(i=0; (i<isize); i++)
  80     fp[i]=ffp[i];
  81
  82   sfree(ffp);
  83   return ret;
  84 #else
  85   return xdr3dfcoord(xdrs,(float *)fp,size,(float *)precision);
  86 #endif
  87 }
  88
  89 int xdr_gmx_large_int(XDR *xdrs,gmx_large_int_t *i,const char *warn)
  90 {
  91   /* This routine stores values compatible with xdr_int64_t */
  92
  93   int imaj,imin;
  94   int ret;
  95
  96 #if ((defined SIZEOF_GMX_LARGE_INT) && SIZEOF_GMX_LARGE_INT == 8)
  97   static const gmx_large_int_t two_p32_m1 = 0xFFFFFFFF;
  98   gmx_large_int_t imaj64,imin64;
  99
 100   imaj64 = ((*i)>>32) & two_p32_m1;
 101   imin64 = (*i) & two_p32_m1;
 102   imaj = (int)imaj64;
 103   imin = (int)imin64;
 104 #else
 105   /* Our code has 4 bytes, but we should make sure that this value
 106    * will be correctly read by 8 byte code.
 107    */
 108   if (*i >= 0) {
 109     imaj = 0;
 110   } else {
 111     imaj = -1;
 112   }
 113   imin = *i;
 114 #endif
 115   ret = xdr_int(xdrs,&imaj);
 116   ret = xdr_int(xdrs,&imin);
 117
 118 #if ((defined SIZEOF_GMX_LARGE_INT) && SIZEOF_GMX_LARGE_INT == 8)
 119   *i = (((gmx_large_int_t)imaj << 32) | ((gmx_large_int_t)imin & two_p32_m1));
 120 #else
 121   *i = imin;
 122
 123   if (warn != NULL && (imaj < -1 || imaj > 0)) {
 124     fprintf(stderr,"\nWARNING during %s:\n",warn);
 125     fprintf(stderr,"a step value written by code supporting 64bit integers is read by code that only supports 32bit integers, out of range step value has been converted to %d\n\n",*i);
 126   }
 127 #endif
 128
 129   return ret;
 130 }