short_float.h

Go to the documentation of this file.

// ----------------------------------------------------------------------------
// openCARP is an open cardiac electrophysiology simulator.
//
// Copyright (C) 2020 openCARP project
//
// This program is licensed under the openCARP Academic Public License (APL)
// v1.0: You can use and redistribute it and/or modify it in non-commercial
// academic environments under the terms of APL as published by the openCARP
// project v1.0, or (at your option) any later version. Commercial use requires
// a commercial license (info@opencarp.org).
//
// This program is distributed without any warranty; see the openCARP APL for
// more details.
//
// You should have received a copy of the openCARP APL along with this program
// and can find it online: http://www.opencarp.org/license
// ----------------------------------------------------------------------------
 
// If automatically generated, keep above
// comment as first line in file.
#ifndef __SHORT_FLOAT_H__
#define __SHORT_FLOAT_H__
 
/*
ftp://ftp.fox-toolkit.org/pub/fasthalffloatconversion.pdf
 
Fast Half Float Conversions
Jeroen van der Zijp
November 2008
(Revised September 2010)
 
Use the dumpTables code to generate the corresponding C file.
*/
 
 
#include <cstdint>
 
namespace opencarp {
 
extern const uint16_t basetable[512];
extern const uint16_t shifttable[512];
extern const uint32_t mantissatable[2048];
extern const uint16_t offsettable[64];
extern const uint32_t exponenttable[64];
 
typedef uint16_t short_float;
 
static inline short_float shortFromFloat(const float external_ff) {
  union {
    float as_float;
    uint32_t as_uint32;
  } c;
  c.as_float = external_ff;
  return basetable[(c.as_uint32>>23)&0x1ff]+((c.as_uint32&0x007fffff)>>shifttable[(c.as_uint32>>23)&0x1ff]);
}
 
static inline float floatFromShort(const short_float h) {
  union {
    float as_float;
    uint32_t as_uint32;
  } c;
  c.as_uint32 = mantissatable[offsettable[h>>10]+(h&0x3ff)]+exponenttable[h>>10];
  return c.as_float;
}
 
#define SHORT_FLOAT_MACHINE_EPS 0.0009765625
 
}  // namespace opencarp
 
#endif