c6_lmdap_library_ncs.c

m files for SIAM book, Numerical Computing with MATLAB

  _SFD_EML_CALL(STATE_DURING_DURING_ACTION_TAG,0,39);
  c6_k_x = c6_u * (real_T)(c6_u >= 0.0);
  c6_tcalc = c6_k_x;
  if(c6_tcalc < 0.0) {
    for(c6_i2 = 0; c6_i2 < 31; c6_i2 = c6_i2 + 1) {
      c6_b_x1[c6_i2] = c6_cv0[c6_i2];
    }
    for(c6_i3 = 0; c6_i3 < 77; c6_i3 = c6_i3 + 1) {
      c6_b_x2[c6_i3] = c6_cv1[c6_i3];
    }
    for(c6_i4 = 0; c6_i4 < 31; c6_i4 = c6_i4 + 1) {
      c6_b_u[c6_i4] = c6_b_x1[c6_i4];
    }
    c6_p_y = NULL;
    sf_mex_assign(&c6_p_y, sf_mex_create(&c6_b_u, "y", 8, 0U, 1U, 2, 1, 31));
    for(c6_i5 = 0; c6_i5 < 77; c6_i5 = c6_i5 + 1) {
      c6_c_u[c6_i5] = c6_b_x2[c6_i5];
    }
    c6_q_y = NULL;
    sf_mex_assign(&c6_q_y, sf_mex_create(&c6_c_u, "y", 8, 0U, 1U, 2, 1, 77));
    sf_mex_call("error", 0U, 2U, 14, c6_p_y, 14, c6_q_y);
  }
  c6_tcalc = sqrt(c6_tcalc);
  _SFD_EML_CALL(STATE_DURING_DURING_ACTION_TAG,0,-39);
  sf_debug_symbol_scope_pop();
  *c6_b_Firefct1 = c6_Firefct1;
  *c6_b_Coastfct1 = c6_Coastfct1;
  *c6_b_Firefct2 = c6_Firefct2;
  *c6_b_Coastfct2 = c6_Coastfct2;
  *c6_b_tcalc1 = c6_tcalc1;
  *c6_b_tcalc = c6_tcalc;
  _SFD_CC_CALL(EXIT_OUT_OF_FUNCTION_TAG,1);
  _sfEvent_ = c6_previousEvent;
  sf_debug_check_for_state_inconsistency(_LMdap_Library_NCSMachineNumber_,
   chartInstance->chartNumber, chartInstance->instanceNumber);
}

static void init_script_number_translation(uint32_T c6_machineNumber, uint32_T
 c6_chartNumber)
{
}

static real_T c6_mpower(SFc6_LMdap_Library_NCSInstanceStruct *chartInstance,
 real_T c6_a)
{
  real_T c6_b_a;
  real_T c6_b;
  real_T c6_ak;
  real_T c6_bk;
  real_T c6_x;
  real_T c6_b_x;
  int32_T c6_i6;
  static char c6_cv2[32] = { 'E', 'm', 'b', 'e', 'd', 'd', 'e', 'd', 'M', 'A',
    'T', 'L', 'A', 'B', ':', 'p', 'o', 'w', 'e', 'r', ':',
    'd', 'o', 'm', 'a', 'i', 'n', 'E', 'r', 'r', 'o', 'r' };
  char c6_x1[32];
  int32_T c6_i7;
  static char c6_cv3[102] = { 'D', 'o', 'm', 'a', 'i', 'n', ' ', 'e', 'r', 'r',
    'o', 'r', '.', ' ', 'T', 'o', ' ', 'c', 'o', 'm', 'p'
    , 'u', 't', 'e', ' ', 'c', 'o', 'm', 'p', 'l', 'e', 'x', ' ', 'r', 'e', 's',
    'u', 'l', 't', 's', ',',
    ' ', 'm', 'a', 'k', 'e', ' ', 'a', 't', ' ', 'l', 'e', 'a', 's', 't', ' ',
    'o', 'n', 'e', ' ', 'i', 'n'
    , 'p', 'u', 't', ' ', 'c', 'o', 'm', 'p', 'l', 'e', 'x', ',', ' ', 'e', '.',
    'g', '.', ' ', '\'', 'p',
    'o', 'w', 'e', 'r', '(', 'c', 'o', 'm', 'p', 'l', 'e', 'x', '(', 'a', ')',
    ',', 'b', ')', '\'', '.' };
  char c6_x2[102];
  int32_T c6_i8;
  char c6_u[32];
  const mxArray *c6_y = NULL;
  int32_T c6_i9;
  char c6_b_u[102];
  const mxArray *c6_b_y = NULL;
  c6_b_a = c6_a;
  c6_b = 2.0;
  c6_ak = c6_b_a;
  c6_bk = c6_b;
  if(c6_ak < 0.0) {
    c6_x = c6_bk;
    c6_b_x = c6_x;
    c6_b_x = floor(c6_b_x);
    if(c6_b_x != c6_bk) {
      for(c6_i6 = 0; c6_i6 < 32; c6_i6 = c6_i6 + 1) {
        c6_x1[c6_i6] = c6_cv2[c6_i6];
      }
      for(c6_i7 = 0; c6_i7 < 102; c6_i7 = c6_i7 + 1) {
        c6_x2[c6_i7] = c6_cv3[c6_i7];
      }
      for(c6_i8 = 0; c6_i8 < 32; c6_i8 = c6_i8 + 1) {
        c6_u[c6_i8] = c6_x1[c6_i8];
      }
      c6_y = NULL;
      sf_mex_assign(&c6_y, sf_mex_create(&c6_u, "y", 8, 0U, 1U, 2, 1, 32));
      for(c6_i9 = 0; c6_i9 < 102; c6_i9 = c6_i9 + 1) {
        c6_b_u[c6_i9] = c6_x2[c6_i9];
      }
      c6_b_y = NULL;
      sf_mex_assign(&c6_b_y, sf_mex_create(&c6_b_u, "y", 8, 0U, 1U, 2, 1, 102));
      sf_mex_call("error", 0U, 2U, 14, c6_y, 14, c6_b_y);
    }
  }
  return pow(c6_ak, c6_bk);
}

static const mxArray *c6_sf_marshall(void *chartInstanceVoid, void *c6_u)
{
  const mxArray *c6_y = NULL;
  int32_T c6_i10;
  real_T c6_b_u[2];
  const mxArray *c6_b_y = NULL;
  SFc6_LMdap_Library_NCSInstanceStruct *chartInstance;
  chartInstance = (SFc6_LMdap_Library_NCSInstanceStruct *)chartInstanceVoid;
  c6_y = NULL;
  for(c6_i10 = 0; c6_i10 < 2; c6_i10 = c6_i10 + 1) {
    c6_b_u[c6_i10] = (*((real_T (*)[2])c6_u))[c6_i10];
  }
  c6_b_y = NULL;
  sf_mex_assign(&c6_b_y, sf_mex_create(&c6_b_u, "y", 0, 0U, 1U, 1, 2));
  sf_mex_assign(&c6_y, c6_b_y);
  return c6_y;
}

static const mxArray *c6_b_sf_marshall(void *chartInstanceVoid, void *c6_u)
{
  const mxArray *c6_y = NULL;
  real_T c6_b_u;
  const mxArray *c6_b_y = NULL;
  SFc6_LMdap_Library_NCSInstanceStruct *chartInstance;
  chartInstance = (SFc6_LMdap_Library_NCSInstanceStruct *)chartInstanceVoid;
  c6_y = NULL;
  c6_b_u = *((real_T *)c6_u);
  c6_b_y = NULL;
  sf_mex_assign(&c6_b_y, sf_mex_create(&c6_b_u, "y", 0, 0U, 0U, 0));
  sf_mex_assign(&c6_y, c6_b_y);
  return c6_y;
}

static const mxArray *c6_c_sf_marshall(void *chartInstanceVoid, void *c6_u)
{
  const mxArray *c6_y = NULL;
  real_T c6_b_u;
  const mxArray *c6_b_y = NULL;
  SFc6_LMdap_Library_NCSInstanceStruct *chartInstance;
  chartInstance = (SFc6_LMdap_Library_NCSInstanceStruct *)chartInstanceVoid;
  c6_y = NULL;
  c6_b_u = *((real_T *)c6_u);
  c6_b_y = NULL;
  sf_mex_assign(&c6_b_y, sf_mex_create(&c6_b_u, "y", 0, 0U, 0U, 0));
  sf_mex_assign(&c6_y, c6_b_y);
  return c6_y;
}

static const mxArray *c6_d_sf_marshall(void *chartInstanceVoid, void *c6_u)
{
  const mxArray *c6_y = NULL;
  real_T c6_b_u;
  const mxArray *c6_b_y = NULL;
  SFc6_LMdap_Library_NCSInstanceStruct *chartInstance;
  chartInstance = (SFc6_LMdap_Library_NCSInstanceStruct *)chartInstanceVoid;
  c6_y = NULL;
  c6_b_u = *((real_T *)c6_u);
  c6_b_y = NULL;
  sf_mex_assign(&c6_b_y, sf_mex_create(&c6_b_u, "y", 0, 0U, 0U, 0));
  sf_mex_assign(&c6_y, c6_b_y);
  return c6_y;
}

static const mxArray *c6_e_sf_marshall(void *chartInstanceVoid, void *c6_u)
{
  const mxArray *c6_y = NULL;
  real_T c6_b_u;
  const mxArray *c6_b_y = NULL;
  SFc6_LMdap_Library_NCSInstanceStruct *chartInstance;
  chartInstance = (SFc6_LMdap_Library_NCSInstanceStruct *)chartInstanceVoid;
  c6_y = NULL;
  c6_b_u = *((real_T *)c6_u);
  c6_b_y = NULL;
  sf_mex_assign(&c6_b_y, sf_mex_create(&c6_b_u, "y", 0, 0U, 0U, 0));
  sf_mex_assign(&c6_y, c6_b_y);
  return c6_y;
}

static const mxArray *c6_f_sf_marshall(void *chartInstanceVoid, void *c6_u)
{
  const mxArray *c6_y = NULL;
  real_T c6_b_u;
  const mxArray *c6_b_y = NULL;
  SFc6_LMdap_Library_NCSInstanceStruct *chartInstance;
  chartInstance = (SFc6_LMdap_Library_NCSInstanceStruct *)chartInstanceVoid;
  c6_y = NULL;
  c6_b_u = *((real_T *)c6_u);
  c6_b_y = NULL;
  sf_mex_assign(&c6_b_y, sf_mex_create(&c6_b_u, "y", 0, 0U, 0U, 0));
  sf_mex_assign(&c6_y, c6_b_y);
  return c6_y;
}

static const mxArray *c6_g_sf_marshall(void *chartInstanceVoid, void *c6_u)
{
  const mxArray *c6_y = NULL;
  real_T c6_b_u;
  const mxArray *c6_b_y = NULL;
  SFc6_LMdap_Library_NCSInstanceStruct *chartInstance;
  chartInstance = (SFc6_LMdap_Library_NCSInstanceStruct *)chartInstanceVoid;
  c6_y = NULL;
  c6_b_u = *((real_T *)c6_u);
  c6_b_y = NULL;
  sf_mex_assign(&c6_b_y, sf_mex_create(&c6_b_u, "y", 0, 0U, 0U, 0));
  sf_mex_assign(&c6_y, c6_b_y);
  return c6_y;
}

static const mxArray *c6_h_sf_marshall(void *chartInstanceVoid, void *c6_u)
{
  const mxArray *c6_y = NULL;
  real_T c6_b_u;
  const mxArray *c6_b_y = NULL;
  SFc6_LMdap_Library_NCSInstanceStruct *chartInstance;
  chartInstance = (SFc6_LMdap_Library_NCSInstanceStruct *)chartInstanceVoid;
  c6_y = NULL;
  c6_b_u = *((real_T *)c6_u);
  c6_b_y = NULL;
  sf_mex_assign(&c6_b_y, sf_mex_create(&c6_b_u, "y", 0, 0U, 0U, 0));
  sf_mex_assign(&c6_y, c6_b_y);
  return c6_y;
}

static const mxArray *c6_i_sf_marshall(void *chartInstanceVoid, void *c6_u)
{
  const mxArray *c6_y = NULL;
  real_T c6_b_u;
  const mxArray *c6_b_y = NULL;
  SFc6_LMdap_Library_NCSInstanceStruct *chartInstance;
  chartInstance = (SFc6_LMdap_Library_NCSInstanceStruct *)chartInstanceVoid;
  c6_y = NULL;
  c6_b_u = *((real_T *)c6_u);
  c6_b_y = NULL;
  sf_mex_assign(&c6_b_y, sf_mex_create(&c6_b_u, "y", 0, 0U, 0U, 0));
  sf_mex_assign(&c6_y, c6_b_y);
  return c6_y;
}

static const mxArray *c6_j_sf_marshall(void *chartInstanceVoid, void *c6_u)
{
  const mxArray *c6_y = NULL;
  real_T c6_b_u;
  const mxArray *c6_b_y = NULL;
  SFc6_LMdap_Library_NCSInstanceStruct *chartInstance;
  chartInstance = (SFc6_LMdap_Library_NCSInstanceStruct *)chartInstanceVoid;
  c6_y = NULL;
  c6_b_u = *((real_T *)c6_u);
  c6_b_y = NULL;
  sf_mex_assign(&c6_b_y, sf_mex_create(&c6_b_u, "y", 0, 0U, 0U, 0));
  sf_mex_assign(&c6_y, c6_b_y);
  return c6_y;
}

static const mxArray *c6_k_sf_marshall(void *chartInstanceVoid, void *c6_u)
{
  const mxArray *c6_y = NULL;
  real_T c6_b_u;
  const mxArray *c6_b_y = NULL;
  SFc6_LMdap_Library_NCSInstanceStruct *chartInstance;
  chartInstance = (SFc6_LMdap_Library_NCSInstanceStruct *)chartInstanceVoid;
  c6_y = NULL;
  c6_b_u = *((real_T *)c6_u);
  c6_b_y = NULL;
  sf_mex_assign(&c6_b_y, sf_mex_create(&c6_b_u, "y", 0, 0U, 0U, 0));
  sf_mex_assign(&c6_y, c6_b_y);
  return c6_y;
}

static const mxArray *c6_l_sf_marshall(void *chartInstanceVoid, void *c6_u)
{
  const mxArray *c6_y = NULL;
  real_T c6_b_u;
  const mxArray *c6_b_y = NULL;
  SFc6_LMdap_Library_NCSInstanceStruct *chartInstance;
  chartInstance = (SFc6_LMdap_Library_NCSInstanceStruct *)chartInstanceVoid;
  c6_y = NULL;
  c6_b_u = *((real_T *)c6_u);
  c6_b_y = NULL;
  sf_mex_assign(&c6_b_y, sf_mex_create(&c6_b_u, "y", 0, 0U, 0U, 0));
  sf_mex_assign(&c6_y, c6_b_y);
  return c6_y;
}

static const mxArray *c6_m_sf_marshall(void *chartInstanceVoid, void *c6_u)
{
  const mxArray *c6_y = NULL;
  real_T c6_b_u;
  const mxArray *c6_b_y = NULL;
  SFc6_LMdap_Library_NCSInstanceStruct *chartInstance;
  chartInstance = (SFc6_LMdap_Library_NCSInstanceStruct *)chartInstanceVoid;
  c6_y = NULL;
  c6_b_u = *((real_T *)c6_u);
  c6_b_y = NULL;
  sf_mex_assign(&c6_b_y, sf_mex_create(&c6_b_u, "y", 0, 0U, 0U, 0));
  sf_mex_assign(&c6_y, c6_b_y);
  return c6_y;
}

static const mxArray *c6_n_sf_marshall(void *chartInstanceVoid, void *c6_u)
{
  const mxArray *c6_y = NULL;
  real_T c6_b_u;
  const mxArray *c6_b_y = NULL;
  SFc6_LMdap_Library_NCSInstanceStruct *chartInstance;
  chartInstance = (SFc6_LMdap_Library_NCSInstanceStruct *)chartInstanceVoid;
  c6_y = NULL;
  c6_b_u = *((real_T *)c6_u);
  c6_b_y = NULL;
  sf_mex_assign(&c6_b_y, sf_mex_create(&c6_b_u, "y", 0, 0U, 0U, 0));
  sf_mex_assign(&c6_y, c6_b_y);
  return c6_y;
}

static const mxArray *c6_o_sf_marshall(void *chartInstanceVoid, void *c6_u)
{
  const mxArray *c6_y = NULL;
  real_T c6_b_u;
  const mxArray *c6_b_y = NULL;
  SFc6_LMdap_Library_NCSInstanceStruct *chartInstance;
  chartInstance = (SFc6_LMdap_Library_NCSInstanceStruct *)chartInstanceVoid;
  c6_y = NULL;
  c6_b_u = *((real_T *)c6_u);
  c6_b_y = NULL;