rt_nonfinite.c

使用c语言及c++来实现matlab编程

/*
 * rt_nonfinite.c
 *
 * Embedded MATLAB Coder code generation for M-function 'PointingAngleCalc2'
 *
 * C source code generated on: Thu May 22 09:59:50 2008
 *
 */

/*
 * Abstract:
 *      Real-Time Workshop function to intialize non-finites,
 *      (Inf, NaN and -Inf).
 */
#include "rt_nonfinite.h"
#define NumBitsPerChar	8

real_T rtInf;
real_T rtMinusInf;
real_T rtNaN;
real32_T rtInfF;
real32_T rtMinusInfF;
real32_T rtNaNF;

/* Function: rt_InitInfAndNaN ==================================================
 * Abstract:
 *	Initialize the rtInf, rtMinusInf, and rtNaN needed by the
 *	generated code. NaN is initialized as non-signaling. Assumes IEEE.
 */
void rt_InitInfAndNaN(size_t realSize)
{
  size_t bitsPerReal = realSize * (NumBitsPerChar);
  uint16_T one = 1;
  enum {
    LittleEndian,
    BigEndian
  } machByteOrder = (*((uint8_T *) &one) == 1) ? LittleEndian : BigEndian;
  switch (machByteOrder) {
   case LittleEndian:
    {
      typedef struct {
        uint32_T fraction : 23;
        uint32_T exponent : 8;
        uint32_T sign : 1;
      } LittleEndianIEEESingle;

      union {
        LittleEndianIEEESingle bitVal;
        real32_T fltVal;
      } tmpValF;

      typedef struct {
        struct {
          uint32_T fraction2;
        } wordH;

        struct {
          uint32_T fraction1 : 20;
          uint32_T exponent : 11;
          uint32_T sign : 1;
        } wordL;
      } LittleEndianIEEEDouble;

      union {
        LittleEndianIEEEDouble bitVal;
        real_T fltVal;
      } tmpVal;

      tmpValF.bitVal.sign = 0;
      tmpValF.bitVal.exponent = 0xFF;
      tmpValF.bitVal.fraction = 0;
      rtInfF = tmpValF.fltVal;
      tmpValF.bitVal.sign = 1;
      rtMinusInfF = tmpValF.fltVal;
      tmpValF.bitVal.fraction = 0x400000;
      rtNaNF = tmpValF.fltVal;
      if (bitsPerReal == 32) {
        rtInf = rtInfF;
        rtMinusInf = rtMinusInfF;
        rtNaN = rtNaNF;
      } else {
        tmpVal.bitVal.wordL.sign = 0;
        tmpVal.bitVal.wordL.exponent = 0x7FF;
        tmpVal.bitVal.wordL.fraction1 = 0;
        tmpVal.bitVal.wordH.fraction2 = 0;
        rtInf = tmpVal.fltVal;
        tmpVal.bitVal.wordL.sign = 1;
        rtMinusInf = tmpVal.fltVal;
        tmpVal.bitVal.wordL.fraction1 = 0x80000;
        rtNaN = tmpVal.fltVal;
      }
    }
    break;

   case BigEndian:
    {
      typedef struct {
        uint32_T sign : 1;
        uint32_T exponent : 8;
        uint32_T fraction : 23;
      } BigEndianIEEESingle;

      union {
        BigEndianIEEESingle bitVal;
        real32_T fltVal;
      } tmpValF;

      typedef struct {
        struct {
          uint32_T sign : 1;
          uint32_T exponent : 11;
          uint32_T fraction1 : 20;
        } wordL;

        struct {
          uint32_T fraction2;
        } wordH;
      } BigEndianIEEEDouble;

      union {
        BigEndianIEEEDouble bitVal;
        real_T fltVal;
      } tmpVal;

      tmpValF.bitVal.sign = 0;
      tmpValF.bitVal.exponent = 0xFF;
      tmpValF.bitVal.fraction = 0;
      rtInfF = tmpValF.fltVal;
      tmpValF.bitVal.sign = 1;
      rtMinusInfF = tmpValF.fltVal;
      tmpValF.bitVal.fraction = 0x400000;
      rtNaNF = tmpValF.fltVal;
      tmpValF.bitVal.sign = 0;
      tmpValF.bitVal.fraction = 0x7fffff;
      rtNaNF = tmpValF.fltVal;
      if (bitsPerReal == 32) {
        rtInf = rtInfF;
        rtMinusInf = rtMinusInfF;
        rtNaN = rtNaNF;
      } else {
        tmpVal.bitVal.wordL.sign = 0;
        tmpVal.bitVal.wordL.exponent = 0x7FF;
        tmpVal.bitVal.wordL.fraction1 = 0;
        tmpVal.bitVal.wordH.fraction2 = 0;
        rtInf = tmpVal.fltVal;
        tmpVal.bitVal.wordL.sign = 1;
        rtMinusInf = tmpVal.fltVal;
        tmpVal.bitVal.wordL.sign = 0;
        tmpVal.bitVal.wordL.fraction1 = 0xFFFFF;
        tmpVal.bitVal.wordH.fraction2 = 0xFFFFFFFF;
        rtNaN = tmpVal.fltVal;
      }
    }
    break;
  }
}

/* Function: rtIsInf ==================================================
 * Abstract:
 *	Test if value is infinite
 */
boolean_T rtIsInf(real_T value)
{
  return ((value==rtInf || value==rtMinusInf) ? 1U : 0U);
}

/* Function: rtIsInfF =================================================
 * Abstract:
 *	Test if single-precision value is infinite
 */
boolean_T rtIsInfF(real32_T value)
{
  return(((value)==rtInfF || (value)==rtMinusInfF) ? 1U : 0U);
}

/* Function: rtIsNaN ==================================================
 * Abstract:
 *	Test if value is not a number
 */
boolean_T rtIsNaN(real_T value)
{
#ifdef _MSC_VER
  return _isnan(value)? TRUE:FALSE;
#else
  return (value!=value)? 1U:0U;
#endif
}

/* Function: rtIsNaNF =================================================
 * Abstract:
 *	Test if single-precision value is not a number
 */
boolean_T rtIsNaNF(real32_T value)
{
#ifdef _MSC_VER
  return _isnan((real_T)value)? true:false;
#else
  return (value!=value)? 1U:0U;
#endif
}

#undef NumBitsPerChar

/* End of Embedded MATLAB Coder code generation (rt_nonfinite.c) */