src/gmxlib/xdrd.c

   1 /*
   2  * This file is part of the GROMACS molecular simulation package.
   3  *
   4  * Copyright (c) 1991-2000, University of Groningen, The Netherlands.
   5  * Copyright (c) 2001-2004, The GROMACS development team,
   6  * check out http://www.gromacs.org for more information.
   7  * Copyright (c) 2012, by the GROMACS development team, led by
   8  * David van der Spoel, Berk Hess, Erik Lindahl, and including many
   9  * others, as listed in the AUTHORS file in the top-level source
  10  * directory and at http://www.gromacs.org.
  11  *
  12  * GROMACS is free software; you can redistribute it and/or
  13  * modify it under the terms of the GNU Lesser General Public License
  14  * as published by the Free Software Foundation; either version 2.1
  15  * of the License, or (at your option) any later version.
  16  *
  17  * GROMACS is distributed in the hope that it will be useful,
  18  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  19  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  20  * Lesser General Public License for more details.
  21  *
  22  * You should have received a copy of the GNU Lesser General Public
  23  * License along with GROMACS; if not, see
  24  * http://www.gnu.org/licenses, or write to the Free Software Foundation,
  25  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA.
  26  *
  27  * If you want to redistribute modifications to GROMACS, please
  28  * consider that scientific software is very special. Version
  29  * control is crucial - bugs must be traceable. We will be happy to
  30  * consider code for inclusion in the official distribution, but
  31  * derived work must not be called official GROMACS. Details are found
  32  * in the README & COPYING files - if they are missing, get the
  33  * official version at http://www.gromacs.org.
  34  *
  35  * To help us fund GROMACS development, we humbly ask that you cite
  36  * the research papers on the package. Check out http://www.gromacs.org.
  37  */
  38 #ifdef HAVE_CONFIG_H
  39 #include <config.h>
  40 #endif
  41
  42 #include "typedefs.h"
  43 #include "xdrf.h"
  44 #include "gmx_fatal.h"
  45 #include "smalloc.h"
  46
  47 int xdr_real(XDR *xdrs,real *r)
  48 {
  49 #ifdef GMX_DOUBLE
  50   float f;
  51   int   ret;
  52
  53   f=*r;
  54   ret=xdr_float(xdrs,&f);
  55   *r=f;
  56
  57   return ret;
  58 #else
  59   return xdr_float(xdrs,(float *)r);
  60 #endif
  61 }
  62
  63 int xdr3drcoord(XDR *xdrs, real *fp, int *size, real *precision)
  64 {
  65 #ifdef GMX_DOUBLE
  66   float *ffp;
  67   float  fprec;
  68   int    i,ret,isize;
  69
  70   isize=*size*DIM;
  71   if (isize <= 0)
  72     gmx_fatal(FARGS,"Don't know what to malloc for ffp, isize = %d",isize);
  73
  74   snew(ffp,isize);
  75
  76   for(i=0; (i<isize); i++)
  77     ffp[i]=fp[i];
  78   fprec=*precision;
  79   ret=xdr3dfcoord(xdrs,ffp,size,&fprec);
  80
  81   *precision=fprec;
  82   for(i=0; (i<isize); i++)
  83     fp[i]=ffp[i];
  84
  85   sfree(ffp);
  86   return ret;
  87 #else
  88   return xdr3dfcoord(xdrs,(float *)fp,size,(float *)precision);
  89 #endif
  90 }
  91
  92 int xdr_gmx_large_int(XDR *xdrs,gmx_large_int_t *i,const char *warn)
  93 {
  94   /* This routine stores values compatible with xdr_int64_t */
  95
  96   int imaj,imin;
  97   int ret;
  98
  99 #if ((defined SIZEOF_GMX_LARGE_INT) && SIZEOF_GMX_LARGE_INT == 8)
 100   static const gmx_large_int_t two_p32_m1 = 0xFFFFFFFF;
 101   gmx_large_int_t imaj64,imin64;
 102
 103   imaj64 = ((*i)>>32) & two_p32_m1;
 104   imin64 = (*i) & two_p32_m1;
 105   imaj = (int)imaj64;
 106   imin = (int)imin64;
 107 #else
 108   /* Our code has 4 bytes, but we should make sure that this value
 109    * will be correctly read by 8 byte code.
 110    */
 111   if (*i >= 0) {
 112     imaj = 0;
 113   } else {
 114     imaj = -1;
 115   }
 116   imin = *i;
 117 #endif
 118   ret = xdr_int(xdrs,&imaj);
 119   ret = xdr_int(xdrs,&imin);
 120
 121 #if ((defined SIZEOF_GMX_LARGE_INT) && SIZEOF_GMX_LARGE_INT == 8)
 122   *i = (((gmx_large_int_t)imaj << 32) | ((gmx_large_int_t)imin & two_p32_m1));
 123 #else
 124   *i = imin;
 125
 126   if (warn != NULL && (imaj < -1 || imaj > 0)) {
 127     fprintf(stderr,"\nWARNING during %s:\n",warn);
 128     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);
 129   }
 130 #endif
 131
 132   return ret;
 133 }