Public API Reference

csutil/csendian.h

Go to the documentation of this file.

/*
    Copyright (C) 1998 by Jorrit Tyberghein

    This library is free software; you can redistribute it and/or
    modify it under the terms of the GNU Library General Public
    License as published by the Free Software Foundation; either
    version 2 of the License, or (at your option) any later version.

    This library is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
    Library General Public License for more details.

    You should have received a copy of the GNU Library General Public
    License along with this library; if not, write to the Free
    Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
*/

#ifndef __CS_CSENDIAN_H__
#define __CS_CSENDIAN_H__

#include <math.h>
#include "cstypes.h"
#include "csgeom/math.h"
#include "csutil/bitops.h"
#if defined(CS_HAVE_BYTESWAP_H)
#include <byteswap.h>
#endif

#define csQroundSure(x) (int ((x) + ((x < 0) ? -0.5 : +0.5)))

struct csSwapBytes
{
public:

  static CS_FORCEINLINE uint16 Swap (uint16 s) 
  { 
  #if defined(CS_COMPILER_MSVC) && (_MSC_VER >= 1300)
    return _byteswap_ushort (s);
  #elif defined(CS_HAVE_BYTESWAP_H)
    return bswap_16 (s);
  #else
    return (s >> 8) | (s << 8); 
  #endif
  }
  static CS_FORCEINLINE int16  Swap (int16 s)
  { return (int16)Swap ((uint16)s); }
  static CS_FORCEINLINE uint32 Swap (uint32 l)
  { 
  #if defined(CS_COMPILER_MSVC) && (_MSC_VER >= 1300)
    return _byteswap_ulong (l);
  #elif defined(CS_HAVE_BYTESWAP_H)
    return bswap_32 (l);
  #else
    return (l >> 24) | ((l >> 8) & 0xff00) | ((l << 8) & 0xff0000) | (l << 24); 
  #endif
  }
  static CS_FORCEINLINE int32  Swap (int32 l)
  { return (int32)Swap ((uint32)l); }
  static CS_FORCEINLINE uint64 Swap (uint64 l)
  {
  #if defined(CS_COMPILER_MSVC) && (_MSC_VER >= 1300)
    return _byteswap_uint64 (l);
  #elif defined(CS_HAVE_BYTESWAP_H) && !defined(__STRICT_ANSI__)
    return bswap_64 (l);
  #else
    union
    {
      uint64 ui64;
      uint32 ui32[2];
    } u1, u2;
    u1.ui64 = l;
    u2.ui32[0] = Swap (u1.ui32[1]);
    u2.ui32[1] = Swap (u1.ui32[0]);
    return u2.ui64;
  #endif
  }
  static CS_FORCEINLINE int64  Swap (int64 l)
  { return (int64)Swap ((uint64)l); }
  
  static CS_FORCEINLINE uint16 UInt16 (uint16 x) { return Swap (x); }
  static CS_FORCEINLINE int16  Int16  (int16 x)  { return Swap (x); }
  static CS_FORCEINLINE uint32 UInt32 (uint32 x) { return Swap (x); }
  static CS_FORCEINLINE int32  Int32  (int32 x)  { return Swap (x); }
  static CS_FORCEINLINE uint64 UInt64 (uint64 x) { return Swap (x); }
  static CS_FORCEINLINE int64  Int64  (int64 x)  { return Swap (x); }
};

#ifdef CS_BIG_ENDIAN
struct csBigEndian
#else
struct csLittleEndian
#endif
{

  static CS_FORCEINLINE uint16 Convert (uint16 x) { return x; }
  static CS_FORCEINLINE int16  Convert (int16 x)  { return x; }
  static CS_FORCEINLINE uint32 Convert (uint32 x) { return x; }
  static CS_FORCEINLINE int32  Convert (int32 x)  { return x; }
  static CS_FORCEINLINE uint64 Convert (uint64 x) { return x; }
  static CS_FORCEINLINE int64  Convert (int64 x)  { return x; }
  
  static CS_FORCEINLINE uint16 UInt16 (uint16 x) { return Convert (x); }
  static CS_FORCEINLINE int16  Int16  (int16 x)  { return Convert (x); }
  static CS_FORCEINLINE uint32 UInt32 (uint32 x) { return Convert (x); }
  static CS_FORCEINLINE int32  Int32  (int32 x)  { return Convert (x); }
  static CS_FORCEINLINE uint64 UInt64 (uint64 x) { return Convert (x); }
  static CS_FORCEINLINE int64  Int64  (int64 x)  { return Convert (x); }
};

#ifdef CS_LITTLE_ENDIAN

struct csBigEndian
#else
struct csLittleEndian
#endif
{
public:

  static CS_FORCEINLINE uint16 Convert (uint16 s) 
  { return csSwapBytes::Swap (s); }
  static CS_FORCEINLINE int16  Convert (int16 s)
  { return csSwapBytes::Swap (s); }
  static CS_FORCEINLINE uint32 Convert (uint32 l)
  { return csSwapBytes::Swap (l); }
  static CS_FORCEINLINE int32  Convert (int32 l)
  { return csSwapBytes::Swap (l); }
  static CS_FORCEINLINE uint64 Convert (uint64 l)
  { return csSwapBytes::Swap (l); }
  static CS_FORCEINLINE int64  Convert (int64 l)
  { return csSwapBytes::Swap (l); }
  
  static CS_FORCEINLINE uint16 UInt16 (uint16 x) { return Convert (x); }
  static CS_FORCEINLINE int16  Int16  (int16 x)  { return Convert (x); }
  static CS_FORCEINLINE uint32 UInt32 (uint32 x) { return Convert (x); }
  static CS_FORCEINLINE int32  Int32  (int32 x)  { return Convert (x); }
  static CS_FORCEINLINE uint64 UInt64 (uint64 x) { return Convert (x); }
  static CS_FORCEINLINE int64  Int64  (int64 x)  { return Convert (x); }
};

struct csIEEEfloat
{
  /* \todo It would be even better if we also check for sizeof (float)
   * in configure or so. */
#ifdef CS_IEEE_DOUBLE_FORMAT


  static CS_FORCEINLINE uint32 FromNative (float f)
  { 
    union
    {
      float f;
      uint32 ui32;
    } u;
    u.f = f;
    return u.ui32; 
  }
  static CS_FORCEINLINE uint64 FromNative (double f)
  { 
    union
    {
      double f;
      uint64 ui64;
    } u;
    u.f = f;
    return u.ui64; 
  }


  static CS_FORCEINLINE float ToNative (uint32 f)
  { 
    union
    {
      float f;
      uint32 ui32;
    } u;
    u.ui32 = f;
    return u.f; 
  }
  static CS_FORCEINLINE double ToNative (uint64 f)
  { 
    union
    {
      double f;
      uint64 ui64;
    } u;
    u.ui64 = f;
    return u.f; 
  }
#else
  #error Do not know how to convert to IEEE floats
#endif

  static CS_FORCEINLINE float ToNative (uint16 half)
  {
    union
    {
      uint32 u;
      float f;
    } u2f;

    uint32 sign = (half & 0x8000) << 16;
    int32 exponent = (half & 0x7C00) >> 10;
    uint32 mantissa = (half & 0x03ff) << 13;
   
    // Check for INF or NaN.
    if (exponent == 0x1F)
    {
      u2f.u = sign | mantissa;

      if (mantissa != 0)
      {
        // NaN
        u2f.u |= 0x7FC00000;
      }
      else
      {
        // INF
        u2f.u |= 0x7f800000;
      }

      return u2f.f;
    }

    // Check for a denorm.
    if(exponent == 0)
    {
      unsigned long index;
      if (CS::Utility::BitOps::ScanBitReverse (mantissa, index))
      {
        exponent -= (index - 9);
        mantissa <<= (index - 8);
        mantissa &= 0x007FFFFF;
      }
      else
      {
        // 0 mantissa -> force float to be 0
        return (sign != 0) ? -0.0f : 0.0f;
      }
    }

    // Convert the exponent...
    exponent += 112;
    exponent <<= 23;

    // And create the float...
    u2f.u = sign | exponent | mantissa;

    return u2f.f;
  }

  static CS_FORCEINLINE uint16 FromNativeRTZ (float f)
  {
    union
    {
      float f;
      unsigned int u;
    } f2u;

    f2u.f = f;
    unsigned short sign = 0x8000 & (f2u.u >> 16);

    // Get the absolute value.
    f2u.u &= 0x7FFFFFFF;

    // Check for a NaN
    if(CS::IsNaN (f2u.f))
    {
      // Construct a silent NaN.
      f2u.u >>= 13;
      f2u.u &= 0x7fff;
      f2u.u |= 0x0200;
      return sign | f2u.u;
    }

    // Check for overflow.
    if(f2u.u >= 0x47800000)
    {
      // Check for INF.
      if(f2u.u == 0x7F800000)
        return sign | 0x7C00;

      return sign | 0x7BFF;
    }

    // Check for underflow and denorms (flush to zero).
    if(f2u.u < 0x38800000)
      return sign;

    // Convert the float to a half (rounding to zero).
    f2u.u &= 0xFFFFE000U;
    f2u.u -= 0x38000000U;
    return sign | (f2u.u >> 13);
  }
};

struct csGetFromAddress
{

  static CS_FORCEINLINE uint16 UInt16 (const void *buff)
  {
  #ifdef CS_STRICT_ALIGNMENT
    uint16 s; memcpy (&s, buff, sizeof (s));
    return s;
  #else
    return *(uint16 *)buff;
  #endif
  }
  static CS_FORCEINLINE int16  Int16 (const void *buff)
  { return (int16)UInt16 (buff); }
  static CS_FORCEINLINE uint32 UInt32 (const void *buff)
  {
  #ifdef CS_STRICT_ALIGNMENT
    uint32 s; memcpy (&s, buff, sizeof (s));
    return s;
  #else
    return *(uint32 *)buff;
  #endif
  }
  static CS_FORCEINLINE int32  Int32 (const void *buff)
  { return (int32)UInt32 (buff); }
  static CS_FORCEINLINE uint64 UInt64 (const void *buff)
  {
  #ifdef CS_STRICT_ALIGNMENT
    uint64 s; memcpy (&s, buff, sizeof (s));
    return s;
  #else
    return *(uint64 *)buff;
  #endif
  }
  static CS_FORCEINLINE int64  Int64 (const void *buff)
  { return (int64)UInt64 (buff); }
};

struct csSetToAddress
{

  static CS_FORCEINLINE void UInt16 (void *buff, uint16 s)
  {
  #ifdef CS_STRICT_ALIGNMENT
    memcpy (buff, &s, sizeof (s));
  #else
    *((uint16 *)buff) = s;
  #endif
  }
  static CS_FORCEINLINE void Int16  (void *buff, int16 s)
  { UInt16 (buff, (uint16)s); }
  static CS_FORCEINLINE void UInt32 (void *buff, uint32 s)
  {
  #ifdef CS_STRICT_ALIGNMENT
    memcpy (buff, &s, sizeof (s));
  #else
    *((uint32 *)buff) = s;
  #endif
  }
  static CS_FORCEINLINE void Int32  (void *buff, int32 s)
  { UInt32 (buff, (uint32)s); }
  static CS_FORCEINLINE void UInt64 (void *buff, uint64 s)
  {
  #ifdef CS_STRICT_ALIGNMENT
    memcpy (buff, &s, sizeof (s));
  #else
    *((uint64 *)buff) = s;
  #endif
  }
  static CS_FORCEINLINE void Int64  (void *buff, int64 s)
  { UInt64 (buff, (uint64)s); }
};


/*
    To be able to painlessly transfer files betwen platforms, we should
    avoid using native floating-point format. Here are a couple of routines
    that are guaranteed to work on all platforms.

    The floating point is converted to a fixed 1.7.25 bits format
    (one bit sign, 7 bits exponent, 25 bits mantissa) and back,
    so that we can binary store floating-point number without
    cross-platform problems. If you wonder why 1+7+25 = 33 while we
    only have 32 bits, we'll ommit the most significant bit of mantissa
    since it is always 1 (we use normalized numbers). This increases the
    precision twice.

    For double, we use one bit sign, 15 bits exponent, 49 bits mantissa.
*/

CS_DEPRECATED_METHOD_MSG("Use csIEEEfloat methods instead")
static inline int32 csFloatToLong (float f)
{
  int exp;
  int32 mant = csQroundSure (frexp (f, &exp) * 0x1000000);
  int32 sign = mant & 0x80000000;
  if (mant < 0) mant = -mant;
  if (exp > 63) exp = 63; else if (exp < -64) exp = -64;
  return sign | ((exp & 0x7f) << 24) | (mant & 0xffffff);
}

CS_DEPRECATED_METHOD_MSG("Use csIEEEfloat methods instead")
static inline float csLongToFloat (int32 l)
{
  int exp = (l >> 24) & 0x7f;
  if (exp & 0x40) exp = exp | ~0x7f;
  float mant = float (l & 0x00ffffff) / 0x1000000;
  if (l & 0x80000000) mant = -mant;
  return (float) ldexp (mant, exp);
}

/* Implementation note: csDoubleToLongLong() and csLongLongToDouble()
 *
 * We avoid use of CONST_INT64() because 64-bit constants are illegal with g++
 * under -ansi -pedantic, and we want this header to be useful to external
 * projects which use -ansi -pedantic.  Instead, we use bit shifts, such as (1
 * << 59), and construct `mask' manually.
 */

CS_DEPRECATED_METHOD_MSG("Use csIEEEfloat methods instead")
static inline int64 csDoubleToLongLong (double d)
{
  int exp;
  int64 mant = (int64) (frexp (d, &exp) * ((int64)1 << 48));
  int64 sign = mant & ((int64)1 << 59);
  if (mant < 0) mant = -mant;
  if (exp > 32767) exp = 32767; else if (exp < -32768) exp = -32768;
  int64 const mask = ((uint64)0xffff << 32) | (uint64)0xffffffff;
  return sign | ((int64 (exp) & 0x7fff) << 48) | (mant & mask);
}

CS_DEPRECATED_METHOD_MSG("Use csIEEEfloat methods instead")
static inline double csLongLongToDouble (int64 i)
{
  int exp = (i >> 48) & 0x7fff;
  if (exp & 0x4000) exp = exp | ~0x7fff;
  int64 const mask = ((uint64)0xffff << 32) | (uint64)0xffffffff;
  double mant = double (i & mask) / ((int64)1 << 48);
  if (i & ((int64)1 << 59)) mant = -mant;
  return ldexp (mant, exp);
}

/* *\name Floating point conversions
 * These routines are used for converting floating-point numbers
 * into 16-bit shorts and back. This is useful for low-precision data.
 * They use the 1.4.12 format. The range of numbers that can be represented
 * in this format is from 2^-8 to 2^7. The precision for numbers near to
 * 2^-8 (0.00390625) is near 0.000001, for numbers near 2^7 (128) is near 0.03.
 * @{ */

CS_DEPRECATED_METHOD_MSG("Use csIEEEfloat methods instead")
static inline short csFloatToShort (float f)
{
  int exp;
  long mant = csQroundSure (frexp (f, &exp) * 0x1000);
  long sign = mant & 0x8000;
  if (mant < 0) mant = -mant;
  if (exp > 7) mant = 0x7ff, exp = 7; else if (exp < -8) mant = 0, exp = -8;
  return short(sign | ((exp & 0xf) << 11) | (mant & 0x7ff));
}

CS_DEPRECATED_METHOD_MSG("Use csIEEEfloat methods instead")
static inline float csShortToFloat (short s)
{
  int exp = (s >> 11) & 0xf;
  if (exp & 0x8) exp = exp | ~0xf;
  float mant = float ((s & 0x07ff) | 0x0800) / 0x1000;
  if (s & 0x8000) mant = -mant;
  return (float) ldexp (mant, exp);
}

#endif // __CS_CSENDIAN_H__

---