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

/* Include files */
#include "rt_nonfinite.h"
#include "Select_Ground_Motions.h"
#include "power.h"
#include "Select_Ground_Motions_emxutil.h"
#include "Select_Ground_Motions_data.h"

/* Function Definitions */
void b_power(const emxArray_real_T *a, emxArray_real_T *y)
{
  uint32_T uv4[2];
  int32_T k;
  for (k = 0; k < 2; k++) {
    uv4[k] = (uint32_T)a->size[k];
  }

  k = y->size[0] * y->size[1];
  y->size[0] = 1;
  y->size[1] = (int32_T)uv4[1];
  emxEnsureCapacity((emxArray__common *)y, k, (int32_T)sizeof(real_T),
                    &qc_emlrtRTEI);
  for (k = 0; k < (int32_T)uv4[1]; k++) {
    y->data[k] = a->data[k] * a->data[k];
  }
}

void power(const emxArray_real_T *a, emxArray_real_T *y)
{
  uint32_T uv3[2];
  int32_T i14;
  int32_T k;
  for (i14 = 0; i14 < 2; i14++) {
    uv3[i14] = (uint32_T)a->size[i14];
  }

  i14 = y->size[0] * y->size[1];
  y->size[0] = (int32_T)uv3[0];
  y->size[1] = (int32_T)uv3[1];
  emxEnsureCapacity((emxArray__common *)y, i14, (int32_T)sizeof(real_T),
                    &qc_emlrtRTEI);
  i14 = (int32_T)uv3[0] * (int32_T)uv3[1];
  for (k = 0; k < i14; k++) {
    y->data[(int32_T)(1.0 + (real_T)k) - 1] = a->data[(int32_T)(1.0 + (real_T)k)
      - 1] * a->data[(int32_T)(1.0 + (real_T)k) - 1];
  }
}

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