/* file generated by oo2c -- do not edit */
#include "__oo2c.h"
#include "__libc.h"
#include "LowReal.d"

static _ModId _mid;

void LowReal__DefaultHandler(short int errno) {
  LowReal__err = errno;
}

void LowReal__ClearError(void) {
  LowReal__err = 0;
}

short int LowReal__exponent(float x) {
  register int i0;
  i0 = x == 0.00000000F;
  if (i0) goto l0;
  {
    float _src;
    _src = (float)x;
    i0 = *(int*)&_src;
  }
  i0 = (int)_lshr(i0, 23, (unsigned int));
  i0 &= 0xFFU;
  i0 -= 127;
  goto l1;
l0:
  i0 = 0;
l1:
  return (short int)i0;
}

short int LowReal__exponent10(float x) {
  register int i0, i1, i2, i3;
  register float f0, f1;
  f1 = _abs(x);
  i0 = f1 == 0.00000000F;
  if (i0) goto l5;
  i2 = f1 >= 1.00000000E+1F;
  if (i2) goto l0;
  f0 = f1;
  i0 = 0;
  goto l2;
l0:
  i0 = 0;
  f0 = f1;
l1:
  f0 /= 1.00000000E+1F;
  i0++;
  i2 = f0 >= 1.00000000E+1F;
  if (i2) goto l1;
l2:
  i2 = f0 > 0.00000000F;
  if (!(i2)) goto l4;
  i2 = f0 >= 1.00000000F;
  if (i2) goto l4;
  i3 = i0;
l3:
  f0 *= 1.00000000E+1F;
  i3--;
  i1 = f0 <= 0.00000000F;
  if (i1) goto l6;
  i1 = f0 < 1.00000000F;
  if (i1) goto l3;
  goto l6;
l4:
  i3 = i0;
  goto l6;
l5:
  i3 = 0;
l6:
  return (short int)i3;
}

float LowReal__fraction(float x) {
  register int i0;
  register float f0;
  i0 = x == 0.00000000F;
  if (i0) goto l0;
  {
    float _src;
    _src = (float)x;
    i0 = *(unsigned int*)&_src;
  }
  i0 &= 0x807FFFFFU;
  i0 |= 0x3F000000U;
  {
    unsigned int _src;
    _src = (unsigned int)i0;
    f0 = *(float*)&_src;
  }
  f0 *= 2.00000000F;
  goto l1;
l0:
  f0 = 0.00000000F;
l1:
  return (float)f0;
}

unsigned char LowReal__IsInfinity(float real) {
  register int i0;
  {
    float _src;
    _src = (float)real;
    i0 = *(unsigned int*)&_src;
  }
  i0 &= 0x7FFFFFFFU;
  i0 = i0 == 0x7F800000U;
  return (unsigned char)i0;
}

unsigned char LowReal__IsNaN(float real) {
  register int i0, i1;
  {
    float _src;
    _src = (float)real;
    i0 = *(unsigned int*)&_src;
  }
  i1 = i0 & 0x7F800000U;
  i1 = i1 != 0x7F800000U;
  if (i1) goto l0;
  i0 &= 0x7FFFFFU;
  i0 = i0 != 0x0U;
  if (i0) goto l1;
l0:
  i0 = 0;
  goto l2;
l1:
  i0 = 1;
l2:
  return (unsigned char)i0;
}

float LowReal__sign(float x) {
  register int i0;
  register float f0;
  i0 = x < 0.00000000F;
  if (i0) goto l0;
  f0 = 1.00000000F;
  goto l1;
l0:
  f0 = -1.00000000F;
l1:
  return (float)f0;
}

float LowReal__scale(float x, short int n) {
  register int i0, i1;
  register float f0;
  i0 = x == 0.00000000F;
  if (i0) goto l3;
  i0 = LowReal__exponent((float)x);
  i0 = n + i0;
  i1 = i0 > 127;
  if (i1) goto l1;
  i1 = i0 < -126;
  if (i1) goto l0;
  goto l2;
l0:
  f0 = LowReal__sign((float)x);
  f0 *= 1.17549435E-38F;
  goto l2;
l1:
  f0 = LowReal__sign((float)x);
  f0 *= 3.40282347E+38F;
  i1 = 1;
l2:
  if (i1) goto l4;
  i0 += 127;
  i0 = (int)_lshl(i0, 23, (unsigned int));
  {
    float _src;
    _src = (float)x;
    i1 = *(unsigned int*)&_src;
  }
  i1 &= 0x807FFFFFU;
  i0 |= i1;
  {
    unsigned int _src;
    _src = (unsigned int)i0;
    f0 = *(float*)&_src;
  }
  goto l4;
l3:
  f0 = 0.00000000F;
l4:
  return (float)f0;
}

float LowReal__ulp(float x) {
  register int i0;
  register float f0;
  i0 = LowReal__exponent((float)x);
  i0 -= 24;
  i0++;
  f0 = LowReal__scale((float)1.00000000F, (short int)i0);
  return (float)f0;
}

float LowReal__succ(float x) {
  register float f0, f1;
  f0 = LowReal__sign((float)x);
  f1 = LowReal__ulp((float)x);
  f0 = f1 * f0;
  f0 = x + f0;
  return (float)f0;
}

float LowReal__pred(float x) {
  register float f0, f1;
  f0 = LowReal__sign((float)x);
  f1 = LowReal__ulp((float)x);
  f0 = f1 * f0;
  f0 = x - f0;
  return (float)f0;
}

float LowReal__intpart(float x) {
  register int i0, i1;
  register float f0;
  i0 = LowReal__exponent((float)x);
  i0 = 23 - i0;
  i1 = i0 <= 0;
  if (i1) goto l1;
  i1 = i0 <= 23;
  if (i1) goto l0;
  f0 = 0.00000000F;
  goto l2;
l0:
  i1 = _bit_range((unsigned int), i0, 31);
  {
    float _src;
    _src = (float)x;
    i0 = *(unsigned int*)&_src;
  }
  i0 &= i1;
  {
    unsigned int _src;
    _src = (unsigned int)i0;
    f0 = *(float*)&_src;
  }
  goto l2;
l1:
  f0 = x;
l2:
  return (float)f0;
}

float LowReal__fractpart(float x) {
  register float f0;
  f0 = LowReal__intpart((float)x);
  f0 = x - f0;
  return (float)f0;
}

float LowReal__trunc(float x, short int n) {
  register int i0, i1;
  register float f0;
  i1 = n <= 0;
  i0 = 24 - n;
  if (i1) goto l1;
  i1 = i0 <= 0;
  if (i1) goto l0;
  i1 = _bit_range((unsigned int), i0, 31);
  {
    float _src;
    _src = (float)x;
    i0 = *(unsigned int*)&_src;
  }
  i0 &= i1;
  {
    unsigned int _src;
    _src = (unsigned int)i0;
    f0 = *(float*)&_src;
  }
  goto l2;
l0:
  f0 = x;
  goto l2;
l1:
  f0 = 0.00000000F;
l2:
  return (float)f0;
}

float LowReal__round(float x, short int n) {
  register int i0, i1, i2;
  register float f0, f1, f2;
  i1 = 24 - n;
  i0 = n <= 0;
  if (i0) goto l3;
  i0 = i1 <= 0;
  if (i0) goto l2;
  {
    float _src;
    _src = (float)x;
    i0 = *(unsigned int*)&_src;
  }
  i2 = _bit_range((unsigned int), i1, 31);
  i2 = i0 & i2;
  {
    unsigned int _src;
    _src = (unsigned int)i2;
    f1 = *(float*)&_src;
  }
  i1--;
  i1 = (i0 & ((unsigned int)1 << i1)) != 0;
  if (i1) goto l0;
  f2 = f1;
  goto l4;
l0:
  i1 = LowReal__exponent((float)f1);
  i1 -= n;
  i0 = (i0 & ((unsigned int)1 << 31)) != 0;
  i1++;
  f0 = LowReal__scale((float)1.00000000F, (short int)i1);
  if (i0) goto l1;
  f2 = f0 + f1;
  goto l4;
l1:
  f2 = f1 - f0;
  goto l4;
l2:
  f2 = x;
  goto l4;
l3:
  f2 = 0.00000000F;
l4:
  return (float)f2;
}

float LowReal__synthesize(short int expart, float frapart) {
  register float f0;
  f0 = LowReal__scale((float)frapart, (short int)expart);
  return (float)f0;
}

void LowReal__setMode(unsigned int m) {
}

unsigned int LowReal__currentMode(void) {
  return (unsigned int)0x0U;
}

unsigned char LowReal__IsLowException(void) {
  return (unsigned char)0;
}

void LowReal_init(void) {
  _mid = _register_module(&LowReal_md.md, NULL);
  LowReal__ErrorHandler = (void(*)(short int))(int)&LowReal__DefaultHandler;
}