framefast/fftype.hh

/*----------------------------------------------------------------------*/
/*                                                                      */
/* Module Name: fftype                                                  */
/*                                                                      */
/* Module Description: Basic data types for framefast and fantom        */
/*                                                                      */
/* Revision History:                                                    */
/* Rel   Date     Programmer    Comments                                */
/* 0.1   26Sep00  D. Sigg       First release                           */
/*                                                                      */
/* Documentation References:                                            */
/*      Man Pages: fftype.html                                          */
/*      References: none                                                */
/*                                                                      */
/* Author Information:                                                  */
/* Name          Telephone       Fax             e-mail                 */
/* Daniel Sigg   (509) 372-8132  (509) 372-8137  sigg_d@ligo.mit.edu    */
/*                                                                      */
/*                                                                      */
/*                      -------------------                             */
/*                                                                      */
/*                             LIGO                                     */
/*                                                                      */
/*        THE LASER INTERFEROMETER GRAVITATIONAL WAVE OBSERVATORY.      */
/*                                                                      */
/*                     (C) The LIGO Project, 1999.                      */
/*                                                                      */
/*                                                                      */
/* Caltech                              MIT                             */
/* LIGO Project MS 51-33                LIGO Project NW-17 161          */
/* Pasadena CA 91125                    Cambridge MA 01239              */
/*                                                                      */
/* LIGO Hanford Observatory             LIGO Livingston Observatory     */
/* P.O. Box 1970 S9-02                  19100 LIGO Lane Rd.             */
/* Richland WA 99352                    Livingston, LA 70754            */
/*                                                                      */
/*----------------------------------------------------------------------*/

#ifndef _LIGO_FFTYPE_H
#define _LIGO_FFTYPE_H

#include "PConfig.h"
#include <complex>
#ifdef HAVE_INT_TYPES
#include <inttypes.h>
#endif


namespace framefast {

   const int kDefaultFrameVersion = 6;

#ifdef HAVE_INT_TYPES
   typedef int8_t int_1s_t;
   typedef uint8_t int_1u_t;
   typedef int16_t int_2s_t;
   typedef uint16_t int_2u_t;
   typedef int32_t int_4s_t;
   typedef uint32_t int_4u_t;
   typedef int64_t int_8s_t;
   typedef uint64_t int_8u_t;
#else // take a guess
   // CHAR
   typedef char int_1s_t;
   // CHAR_U
   typedef unsigned char int_1u_t;
   // INT_2S
   typedef short int_2s_t;
   // INT_2U
   typedef unsigned short int_2u_t;
   // INT_4S
   typedef int int_4s_t;
   // INT_4U
   typedef unsigned int int_4u_t;
   // INT_8S
   typedef long long int_8s_t;
   // INT_8U
   typedef unsigned long long int_8u_t;
#endif

   typedef float real_4_t;
   typedef double real_8_t;
   typedef std::complex<float> complex_8_t;
   typedef std::complex<double> complex_16_t;

   enum frametype {
   NONE = 0,
   FF = 1,
   STF = 2,
   MTF = 3
   };




#ifndef __CINT__
   inline void swap (int_8u_t* ww);
   inline void swap (int_8s_t* ww);
#endif
   inline void swap (int_4u_t* w);
   inline void swap (int_4s_t* ww);
   inline void swap (int_2u_t* w);
   inline void swap (int_2s_t* ww);
   inline void swap (int_1s_t* ww);
   inline void swap (int_1u_t* ww);
   inline void swap (real_8_t* ww);
   inline void swap (real_4_t* w);
   inline void swap (complex_8_t* ww);
   inline void swap (complex_16_t* ww);
   inline bool littleendian ();

}



//                                                                      //
// inlines                                                              //
//                                                                      //
namespace framefast {
#define byteswap2(a,b) ((char *)&(a))[0] = ((char *)&(b))[1]; ((char *)&(a))[1] = ((char *)&(b))[0];
#define byteswap4(a,b) ((char *)&(a))[0] = ((char *)&(b))[3]; ((char *)&(a))[1] = ((char *)&(b))[2];((char *)&(a))[2] = ((char *)&(b))[1];((char *)&(a))[3] = ((char *)&(b))[0];
#define byteswap8(a,b) ((char *)&(a))[0] = ((char *)&(b))[7]; ((char *)&(a))[1] = ((char *)&(b))[6];((char *)&(a))[2] = ((char *)&(b))[5];((char *)&(a))[3] = ((char *)&(b))[4];((char *)&(a))[4] = ((char *)&(b))[3]; ((char *)&(a))[5] = ((char *)&(b))[2]; ((char *)&(a))[6] = ((char *)&(b))[1]; ((char *)&(a))[7] = ((char *)&(b))[0];

#ifndef __CINT__
   inline void swap (int_8u_t* ww)
   {
      int_8u_t tmp;
      byteswap8(tmp, *ww);
      *ww = tmp;
   }

//______________________________________________________________________________
   inline void swap (int_8s_t* ww) 
   {
      swap ((int_8u_t*)ww);
   }
#endif

//______________________________________________________________________________
   inline void swap (int_4u_t* w)
   {
      int_4u_t tmp;
      byteswap4(tmp, *w);
      *w = tmp;
   }

//______________________________________________________________________________
   inline void swap (int_4s_t* ww) 
   {
      swap ((int_4u_t*)ww);
   }

//______________________________________________________________________________
   inline void swap (int_2u_t* w)
   {
      int_2u_t tmp;
      byteswap2(tmp, *w);
      *w = tmp;
   }

//______________________________________________________________________________
   inline void swap (int_2s_t* ww) 
   {
      swap ((int_2u_t*)ww);
   }

//______________________________________________________________________________
   inline void swap (int_1s_t* ww)
   {
   }

//______________________________________________________________________________
   inline void swap (int_1u_t* ww)
   {
   }

//______________________________________________________________________________
   inline void swap (real_8_t* ww)
   {
      real_8_t tmp;
      byteswap8(tmp, *ww);
      *ww = tmp;
   }

//______________________________________________________________________________
   inline void swap (real_4_t* w)
   {
      real_4_t tmp;
      byteswap4(tmp, *w);
      *w = tmp;
   }

//______________________________________________________________________________
   inline void swap (complex_8_t* ww)
   {
      struct complex {
         real_4_t re;
         real_4_t im;
      };
      swap (&((complex*)ww)->re);
      swap (&((complex*)ww)->im);
   }

//______________________________________________________________________________
   inline void swap (complex_16_t* ww)
   { 
      struct complex {
         real_8_t re;
         real_8_t im;
      } ;
      swap (&((complex*)ww)->re);
      swap (&((complex*)ww)->im);
   }

   inline bool littleendian () {
      int_4u_t test = 0; *(int_1u_t*) &test = 1;
      return (test == 1);
   }

#undef byteswap2
#undef byteswap4
#undef byteswap8
}

#endif // _LIGO_FFTYPE_H

---