/*
 * eps.c
 *
 * Code generation for function 'eps'
 *
 * C source code generated on: Wed Aug 26 14:59:33 2015
 *
 */

/* Include files */
#include "rt_nonfinite.h"
#include "Select_Ground_Motions.h"
#include "eps.h"

/* Function Definitions */
void eps(const real_T x_data[1], const int32_T x_size[2], real_T r_data[1],
         int32_T r_size[2])
{
  int32_T k;
  real_T absxk;
  int32_T exponent;
  for (k = 0; k < 2; k++) {
    r_size[k] = x_size[k];
  }

  k = 0;
  while (k <= x_size[1] - 1) {
    absxk = muDoubleScalarAbs(x_data[0]);
    if ((!muDoubleScalarIsInf(absxk)) && (!muDoubleScalarIsNaN(absxk))) {
      if (absxk <= 2.2250738585072014E-308) {
        r_data[0] = 4.94065645841247E-324;
      } else {
        frexp(absxk, &exponent);
        r_data[0] = ldexp(1.0, exponent - 53);
      }
    } else {
      r_data[0] = rtNaN;
    }

    k = 1;
  }
}

/* End of code generation (eps.c) */