File: power.c

    1   /*
    2    * power.c
    3    *
    4    * Code generation for function 'power'
    5    *
    6    */
    7   
    8   /* Include files */
    9   #include "rt_nonfinite.h"
   10   #include "gjbf.h"
   11   #include "power.h"
   12   #include "error1.h"
   13   #include "gjbf_emxutil.h"
   14   #include "scalexpAlloc.h"
   15   #include "gjbf_data.h"
   16   #include "blas.h"
   17   
   18   /* Variable Definitions */
   19   static emlrtRSInfo yb_emlrtRSI = { 61, "power",
   20     "C:\\Program Files\\MATLAB\\R2016a\\toolbox\\eml\\lib\\matlab\\ops\\power.m" };
   21   
   22   /* Function Declarations */
   23   static boolean_T b_fltpower_domain_error(void);
   24   static boolean_T fltpower_domain_error(void);
   25   
   26   /* Function Definitions */
   27   static boolean_T b_fltpower_domain_error(void)
   28   {
   29     return false;
   30   }
   31   
   32   static boolean_T fltpower_domain_error(void)
   33   {
   34     return false;
   35   }
   36   
   37   void b_power(const emlrtStack *sp, const real_T a_data[], const int32_T a_size[1],
   38                real_T y_data[], int32_T y_size[1])
   39   {
   40     int32_T k;
   41     emlrtStack st;
   42     emlrtStack b_st;
   43     emlrtStack c_st;
   44     st.prev = sp;
   45     st.tls = sp->tls;
   46     st.site = &xb_emlrtRSI;
   47     b_st.prev = &st;
   48     b_st.tls = st.tls;
   49     c_st.prev = &b_st;
   50     c_st.tls = b_st.tls;
   51     b_st.site = &fd_emlrtRSI;
   52     c_st.site = &gd_emlrtRSI;
   53     y_size[0] = (int16_T)a_size[0];
   54     if (dimagree(y_size, a_size)) {
   55     } else {
   56       emlrtErrorWithMessageIdR2012b(&c_st, &u_emlrtRTEI, "MATLAB:dimagree", 0);
   57     }
   58   
   59     c_st.site = &hd_emlrtRSI;
   60     for (k = 0; k + 1 <= a_size[0]; k++) {
   61       y_data[k] = a_data[k] * a_data[k];
   62     }
   63   
   64     if (fltpower_domain_error()) {
   65       b_st.site = &yb_emlrtRSI;
   66       v_error(&b_st);
   67     }
   68   }
   69   
   70   void c_power(const emlrtStack *sp, const emxArray_real_T *a, emxArray_real_T *y)
   71   {
   72     int32_T n;
   73     int32_T k;
   74     emlrtStack st;
   75     emlrtStack b_st;
   76     emlrtStack c_st;
   77     st.prev = sp;
   78     st.tls = sp->tls;
   79     st.site = &xb_emlrtRSI;
   80     b_st.prev = &st;
   81     b_st.tls = st.tls;
   82     c_st.prev = &b_st;
   83     c_st.tls = b_st.tls;
   84     b_st.site = &fd_emlrtRSI;
   85     c_st.site = &gd_emlrtRSI;
   86     n = y->size[0] * y->size[1];
   87     y->size[0] = (int16_T)a->size[0];
   88     y->size[1] = 16;
   89     emxEnsureCapacity(&c_st, (emxArray__common *)y, n, (int32_T)sizeof(real_T),
   90                       &c_emlrtRTEI);
   91     if (b_dimagree(y, a)) {
   92     } else {
   93       emlrtErrorWithMessageIdR2012b(&c_st, &u_emlrtRTEI, "MATLAB:dimagree", 0);
   94     }
   95   
   96     n = a->size[0] << 4;
   97     c_st.site = &hd_emlrtRSI;
   98     for (k = 0; k + 1 <= n; k++) {
   99       y->data[k] = a->data[k] * a->data[k];
  100     }
  101   
  102     if (b_fltpower_domain_error()) {
  103       b_st.site = &yb_emlrtRSI;
  104       v_error(&b_st);
  105     }
  106   }
  107   
  108   void power(const real_T a[48], real_T y[48])
  109   {
  110     int32_T k;
  111     for (k = 0; k < 48; k++) {
  112       y[k] = a[k] * a[k];
  113     }
  114   }
  115   
  116   /* End of code generation (power.c) */
  117