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

/* Include files */
#include "rt_nonfinite.h"
#include "Select_Ground_Motions.h"
#include "std.h"
#include "sqrt.h"
#include "Select_Ground_Motions_emxutil.h"
#include "isequal.h"
#include "Select_Ground_Motions_mexutil.h"
#include "Select_Ground_Motions_data.h"

/* Variable Definitions */
static emlrtRSInfo sm_emlrtRSI = { 12, "std",
  "/usr/local/MATLAB/R2013a/toolbox/eml/lib/matlab/datafun/std.m" };

static emlrtRSInfo tm_emlrtRSI = { 17, "var",
  "/usr/local/MATLAB/R2013a/toolbox/eml/lib/matlab/datafun/var.m" };

static emlrtRSInfo um_emlrtRSI = { 20, "var",
  "/usr/local/MATLAB/R2013a/toolbox/eml/lib/matlab/datafun/var.m" };

static emlrtRSInfo vm_emlrtRSI = { 97, "var",
  "/usr/local/MATLAB/R2013a/toolbox/eml/lib/matlab/datafun/var.m" };

static emlrtMCInfo mc_emlrtMCI = { 18, 9, "var",
  "/usr/local/MATLAB/R2013a/toolbox/eml/lib/matlab/datafun/var.m" };

static emlrtMCInfo nc_emlrtMCI = { 17, 19, "var",
  "/usr/local/MATLAB/R2013a/toolbox/eml/lib/matlab/datafun/var.m" };

static emlrtMCInfo oc_emlrtMCI = { 23, 9, "var",
  "/usr/local/MATLAB/R2013a/toolbox/eml/lib/matlab/datafun/var.m" };

static emlrtMCInfo pc_emlrtMCI = { 20, 19, "var",
  "/usr/local/MATLAB/R2013a/toolbox/eml/lib/matlab/datafun/var.m" };

static emlrtRTEInfo yc_emlrtRTEI = { 90, 1, "var",
  "/usr/local/MATLAB/R2013a/toolbox/eml/lib/matlab/datafun/var.m" };

/* Function Definitions */
void b_std(const emxArray_real_T *varargin_1, emxArray_real_T *y)
{
  boolean_T overflow;
  const mxArray *b_y;
  static const int32_T iv63[2] = { 1, 30 };

  const mxArray *m22;
  char_T cv85[30];
  int32_T i;
  static const char_T cv86[30] = { 'C', 'o', 'd', 'e', 'r', ':', 't', 'o', 'o',
    'l', 'b', 'o', 'x', ':', 'v', 'a', 'r', '_', 's', 'p', 'e', 'c', 'i', 'a',
    'l', 'E', 'm', 'p', 't', 'y' };

  const mxArray *c_y;
  static const int32_T iv64[2] = { 1, 36 };

  char_T cv87[36];
  static const char_T cv88[36] = { 'C', 'o', 'd', 'e', 'r', ':', 't', 'o', 'o',
    'l', 'b', 'o', 'x', ':', 'a', 'u', 't', 'o', 'D', 'i', 'm', 'I', 'n', 'c',
    'o', 'm', 'p', 'a', 't', 'i', 'b', 'i', 'l', 'i', 't', 'y' };

  int32_T d;
  uint32_T sz[2];
  int32_T ix;
  int32_T iy;
  real_T d_y;
  int32_T b_ix;
  real_T xbar;
  int32_T k;
  real_T r;
  emlrtPushRtStackR2012b(&sm_emlrtRSI, emlrtRootTLSGlobal);
  overflow = !b_isequal(varargin_1);
  if (overflow) {
  } else {
    emlrtPushRtStackR2012b(&tm_emlrtRSI, emlrtRootTLSGlobal);
    b_y = NULL;
    m22 = mxCreateCharArray(2, iv63);
    for (i = 0; i < 30; i++) {
      cv85[i] = cv86[i];
    }

    emlrtInitCharArrayR2013a(emlrtRootTLSGlobal, 30, m22, cv85);
    emlrtAssign(&b_y, m22);
    error(message(b_y, &mc_emlrtMCI), &nc_emlrtMCI);
    emlrtPopRtStackR2012b(&tm_emlrtRSI, emlrtRootTLSGlobal);
  }

  if (((varargin_1->size[0] == 1) && (varargin_1->size[1] == 1)) ||
      (varargin_1->size[0] != 1)) {
    overflow = TRUE;
  } else {
    overflow = FALSE;
  }

  if (overflow) {
  } else {
    emlrtPushRtStackR2012b(&um_emlrtRSI, emlrtRootTLSGlobal);
    c_y = NULL;
    m22 = mxCreateCharArray(2, iv64);
    for (i = 0; i < 36; i++) {
      cv87[i] = cv88[i];
    }

    emlrtInitCharArrayR2013a(emlrtRootTLSGlobal, 36, m22, cv87);
    emlrtAssign(&c_y, m22);
    error(message(c_y, &oc_emlrtMCI), &pc_emlrtMCI);
    emlrtPopRtStackR2012b(&um_emlrtRSI, emlrtRootTLSGlobal);
  }

  if (varargin_1->size[0] > 1) {
    d = varargin_1->size[0] - 1;
  } else {
    d = varargin_1->size[0];
  }

  for (i = 0; i < 2; i++) {
    sz[i] = (uint32_T)varargin_1->size[i];
  }

  i = y->size[0] * y->size[1];
  y->size[0] = 1;
  y->size[1] = (int32_T)sz[1];
  emxEnsureCapacity((emxArray__common *)y, i, (int32_T)sizeof(real_T),
                    &yc_emlrtRTEI);
  ix = 0;
  iy = -1;
  emlrtPushRtStackR2012b(&vm_emlrtRSI, emlrtRootTLSGlobal);
  if (1 > varargin_1->size[1]) {
    overflow = FALSE;
  } else {
    overflow = (varargin_1->size[1] > 2147483646);
  }

  if (overflow) {
    emlrtPushRtStackR2012b(&lc_emlrtRSI, emlrtRootTLSGlobal);
    check_forloop_overflow_error();
    emlrtPopRtStackR2012b(&lc_emlrtRSI, emlrtRootTLSGlobal);
  }

  emlrtPopRtStackR2012b(&vm_emlrtRSI, emlrtRootTLSGlobal);
  for (i = 1; i <= varargin_1->size[1]; i++) {
    iy++;
    if ((varargin_1->size[0] == 0) || (varargin_1->size[1] == 0)) {
      d_y = rtNaN;
    } else {
      b_ix = ix;
      xbar = varargin_1->data[ix];
      for (k = 0; k <= varargin_1->size[0] - 2; k++) {
        b_ix++;
        xbar += varargin_1->data[b_ix];
      }

      xbar /= (real_T)varargin_1->size[0];
      b_ix = ix;
      r = varargin_1->data[ix] - xbar;
      d_y = r * r;
      for (k = 0; k <= varargin_1->size[0] - 2; k++) {
        b_ix++;
        r = varargin_1->data[b_ix] - xbar;
        d_y += r * r;
      }

      d_y /= (real_T)d;
    }

    y->data[iy] = d_y;
    ix += varargin_1->size[0];
  }

  c_sqrt(y);
  emlrtPopRtStackR2012b(&sm_emlrtRSI, emlrtRootTLSGlobal);
}

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