c2_lyap_uisng_simulink.c

m files for SIAM book, Numerical Computing with MATLAB

/* Include files */
#include "Lyap_uisng_Simulink_sfun.h"
#include "c2_Lyap_uisng_Simulink.h"
#define CHARTINSTANCE_CHARTNUMBER       (chartInstance.chartNumber)
#define CHARTINSTANCE_INSTANCENUMBER    (chartInstance.instanceNumber)
#include "Lyap_uisng_Simulink_sfun_debug_macros.h"

/* Type Definitions */

/* Named Constants */
#define c2_IN_NO_ACTIVE_CHILD           (0)

/* Variable Declarations */

/* Variable Definitions */
static SFc2_Lyap_uisng_SimulinkInstanceStruct chartInstance;

/* Function Declarations */
static void initialize_c2_Lyap_uisng_Simulink(void);
static void initialize_params_c2_Lyap_uisng_Simulink(void);
static void enable_c2_Lyap_uisng_Simulink(void);
static void disable_c2_Lyap_uisng_Simulink(void);
static void finalize_c2_Lyap_uisng_Simulink(void);
static void sf_c2_Lyap_uisng_Simulink(void);
static void init_script_number_translation(uint32_T c2_machineNumber, uint32_T
 c2_chartNumber);
static const mxArray *c2_sf_marshall(void *c2_chartInstance, void *c2_u);
static const mxArray *c2_b_sf_marshall(void *c2_chartInstance, void *c2_u);
static const mxArray *c2_c_sf_marshall(void *c2_chartInstance, void *c2_u);
static void c2_info_helper(c2_ResolvedFunctionInfo c2_info[53]);
static void init_dsm_address_info(void);

/* Function Definitions */
static void initialize_c2_Lyap_uisng_Simulink(void)
{
  _sfTime_ = (real_T)ssGetT(chartInstance.S);
  chartInstance.c2_is_active_c2_Lyap_uisng_Simulink = 0U;
}

static void initialize_params_c2_Lyap_uisng_Simulink(void)
{
}

static void enable_c2_Lyap_uisng_Simulink(void)
{
}

static void disable_c2_Lyap_uisng_Simulink(void)
{
}

static void finalize_c2_Lyap_uisng_Simulink(void)
{
}

static void sf_c2_Lyap_uisng_Simulink(void)
{
  int32_T c2_i0;
  int32_T c2_i1;
  uint8_T c2_previousEvent;
  int32_T c2_i2;
  int32_T c2_i3;
  real_T c2_Pdot[4];
  real_T c2_trP;
  real_T c2_stop;
  int32_T c2_i4;
  int32_T c2_i5;
  real_T c2_a[4];
  real_T c2_x;
  real_T c2_b_x;
  real_T c2_r;
  real_T *c2_b_stop;
  real_T (*c2_b_Pdot)[4];
  c2_b_stop = (real_T *)ssGetOutputPortSignal(chartInstance.S, 1);
  c2_b_Pdot = (real_T (*)[4])ssGetInputPortSignal(chartInstance.S, 0);
  _sfTime_ = (real_T)ssGetT(chartInstance.S);
  for(c2_i0 = 0; c2_i0 < 2; c2_i0 = c2_i0 + 1) {
    for(c2_i1 = 0; c2_i1 < 2; c2_i1 = c2_i1 + 1) {
      _SFD_DATA_RANGE_CHECK((*c2_b_Pdot)[c2_i1 + 2 * c2_i0], 0U);
    }
  }
  _SFD_DATA_RANGE_CHECK(*c2_b_stop, 1U);
  c2_previousEvent = _sfEvent_;
  _sfEvent_ = CALL_EVENT;
  _SFD_CC_CALL(CHART_ENTER_DURING_FUNCTION_TAG,0);
  for(c2_i2 = 0; c2_i2 < 2; c2_i2 = c2_i2 + 1) {
    for(c2_i3 = 0; c2_i3 < 2; c2_i3 = c2_i3 + 1) {
      c2_Pdot[c2_i3 + 2 * c2_i2] = (*c2_b_Pdot)[c2_i3 + 2 * c2_i2];
    }
  }
  sf_debug_symbol_scope_push(3U, 0U);
  sf_debug_symbol_scope_add("trP", &c2_trP, c2_c_sf_marshall);
  sf_debug_symbol_scope_add("stop", &c2_stop, c2_b_sf_marshall);
  sf_debug_symbol_scope_add("Pdot", &c2_Pdot, c2_sf_marshall);
  CV_EML_FCN(0, 0);
  /*  This block supports the Embedded MATLAB subset. */
  /*  See the help menu for details.  */
  _SFD_EML_CALL(STATE_DURING_DURING_ACTION_TAG,0,5);
  c2_stop = 0.0;
  _SFD_EML_CALL(STATE_DURING_DURING_ACTION_TAG,0,6);
  for(c2_i4 = 0; c2_i4 < 2; c2_i4 = c2_i4 + 1) {
    for(c2_i5 = 0; c2_i5 < 2; c2_i5 = c2_i5 + 1) {
      c2_a[c2_i5 + 2 * c2_i4] = c2_Pdot[c2_i5 + 2 * c2_i4];
    }
  }
  c2_x = c2_a[0];
  c2_trP = c2_x;
  c2_b_x = c2_a[3];
  c2_r = c2_b_x;
  c2_trP = c2_trP + c2_r;
  _SFD_EML_CALL(STATE_DURING_DURING_ACTION_TAG,0,7);
  if(CV_EML_IF(0, 0, c2_trP < 1.0E-006)) {
    _SFD_EML_CALL(STATE_DURING_DURING_ACTION_TAG,0,8);
    c2_stop = 1.0;
  }
  _SFD_EML_CALL(STATE_DURING_DURING_ACTION_TAG,0,-8);
  sf_debug_symbol_scope_pop();
  *c2_b_stop = c2_stop;
  _SFD_CC_CALL(EXIT_OUT_OF_FUNCTION_TAG,0);
  _sfEvent_ = c2_previousEvent;
  sf_debug_check_for_state_inconsistency(_Lyap_uisng_SimulinkMachineNumber_,
   chartInstance.chartNumber, chartInstance.instanceNumber);
}

static void init_script_number_translation(uint32_T c2_machineNumber, uint32_T
 c2_chartNumber)
{
}

static const mxArray *c2_sf_marshall(void *c2_chartInstance, void *c2_u)
{
  const mxArray *c2_y = NULL;
  int32_T c2_i6;
  int32_T c2_i7;
  real_T c2_b_u[4];
  const mxArray *c2_b_y = NULL;
  c2_y = NULL;
  for(c2_i6 = 0; c2_i6 < 2; c2_i6 = c2_i6 + 1) {
    for(c2_i7 = 0; c2_i7 < 2; c2_i7 = c2_i7 + 1) {
      c2_b_u[c2_i7 + 2 * c2_i6] = (*((real_T (*)[4])c2_u))[c2_i7 + 2 * c2_i6];
    }
  }
  c2_b_y = NULL;
  sf_mex_assign(&c2_b_y, sf_mex_create(&c2_b_u, "y", 0, 0U, 1U, 2, 2, 2));
  sf_mex_assign(&c2_y, c2_b_y);
  return c2_y;
}

static const mxArray *c2_b_sf_marshall(void *c2_chartInstance, void *c2_u)
{
  const mxArray *c2_y = NULL;
  real_T c2_b_u;
  const mxArray *c2_b_y = NULL;
  c2_y = NULL;
  c2_b_u = *((real_T *)c2_u);
  c2_b_y = NULL;
  sf_mex_assign(&c2_b_y, sf_mex_create(&c2_b_u, "y", 0, 0U, 0U, 0));
  sf_mex_assign(&c2_y, c2_b_y);
  return c2_y;
}

static const mxArray *c2_c_sf_marshall(void *c2_chartInstance, void *c2_u)
{
  const mxArray *c2_y = NULL;
  real_T c2_b_u;
  const mxArray *c2_b_y = NULL;
  c2_y = NULL;
  c2_b_u = *((real_T *)c2_u);
  c2_b_y = NULL;
  sf_mex_assign(&c2_b_y, sf_mex_create(&c2_b_u, "y", 0, 0U, 0U, 0));
  sf_mex_assign(&c2_y, c2_b_y);
  return c2_y;
}

const mxArray *sf_c2_Lyap_uisng_Simulink_get_eml_resolved_functions_info(void)
{
  const mxArray *c2_nameCaptureInfo = NULL;
  c2_ResolvedFunctionInfo c2_info[53];
  const mxArray *c2_m0 = NULL;
  int32_T c2_i8;
  c2_ResolvedFunctionInfo *c2_r0;
  c2_nameCaptureInfo = NULL;
  c2_info_helper(c2_info);
  sf_mex_assign(&c2_m0, sf_mex_createstruct("nameCaptureInfo", 1, 53));
  for(c2_i8 = 0; c2_i8 < 53; c2_i8 = c2_i8 + 1) {
    c2_r0 = &c2_info[c2_i8];
    sf_mex_addfield(c2_m0, sf_mex_create(c2_r0->context, "nameCaptureInfo", 13,
      0U, 0U, 2, 1, strlen(c2_r0->context)), "context",
     "nameCaptureInfo", c2_i8);
    sf_mex_addfield(c2_m0, sf_mex_create(c2_r0->name, "nameCaptureInfo", 13, 0U,
      0U, 2, 1, strlen(c2_r0->name)), "name",
     "nameCaptureInfo", c2_i8);
    sf_mex_addfield(c2_m0, sf_mex_create(c2_r0->dominantType, "nameCaptureInfo",
      13, 0U, 0U, 2, 1, strlen(c2_r0->dominantType)),
     "dominantType", "nameCaptureInfo", c2_i8);
    sf_mex_addfield(c2_m0, sf_mex_create(c2_r0->resolved, "nameCaptureInfo", 13,
      0U, 0U, 2, 1, strlen(c2_r0->resolved)), "resolved",
     "nameCaptureInfo", c2_i8);
    sf_mex_addfield(c2_m0, sf_mex_create(&c2_r0->fileLength, "nameCaptureInfo",
      7, 0U, 0U, 0), "fileLength", "nameCaptureInfo", c2_i8);
    sf_mex_addfield(c2_m0, sf_mex_create(&c2_r0->fileTime1, "nameCaptureInfo",
      7, 0U, 0U, 0), "fileTime1", "nameCaptureInfo", c2_i8);
    sf_mex_addfield(c2_m0, sf_mex_create(&c2_r0->fileTime2, "nameCaptureInfo",
      7, 0U, 0U, 0), "fileTime2", "nameCaptureInfo", c2_i8);
  }
  sf_mex_assign(&c2_nameCaptureInfo, c2_m0);
  return c2_nameCaptureInfo;
}

static void c2_info_helper(c2_ResolvedFunctionInfo c2_info[53])
{
  c2_info[0].context =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/eml/eml_assert_valid_size_arg.m/check_arg";
  c2_info[0].name = "isfloat";
  c2_info[0].dominantType = "double";
  c2_info[0].resolved =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/datatypes/isfloat.m";
  c2_info[0].fileLength = 182U;
  c2_info[0].fileTime1 = 1160446620U;
  c2_info[0].fileTime2 = 0U;
  c2_info[1].context = "[I]$matlabroot$/toolbox/eml/lib/matlab/eml/eml_alloc.m";
  c2_info[1].name = "eml_assert_valid_size_arg";
  c2_info[1].dominantType = "double";
  c2_info[1].resolved =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/eml/eml_assert_valid_size_arg.m";
  c2_info[1].fileLength = 1097U;
  c2_info[1].fileTime1 = 1163970788U;
  c2_info[1].fileTime2 = 0U;
  c2_info[2].context = "[I]$matlabroot$/toolbox/eml/lib/matlab/datatypes/isa.m";
  c2_info[2].name = "eml_is_integer_class";
  c2_info[2].dominantType = "char";
  c2_info[2].resolved =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/eml/eml_is_integer_class.m";
  c2_info[2].fileLength = 462U;
  c2_info[2].fileTime1 = 1160446744U;
  c2_info[2].fileTime2 = 0U;
  c2_info[3].context = "[I]$matlabroot$/toolbox/eml/lib/matlab/matfun/trace.m";
  c2_info[3].name = "isvector";
  c2_info[3].dominantType = "double";
  c2_info[3].resolved =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/elmat/isvector.m";
  c2_info[3].fileLength = 284U;
  c2_info[3].fileTime1 = 1160446706U;
  c2_info[3].fileTime2 = 0U;
  c2_info[4].context = "[I]$matlabroot$/toolbox/eml/lib/matlab/matfun/trace.m";
  c2_info[4].name = "cast";
  c2_info[4].dominantType = "double";
  c2_info[4].resolved =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/datatypes/cast.m";
  c2_info[4].fileLength = 907U;
  c2_info[4].fileTime1 = 1160446618U;
  c2_info[4].fileTime2 = 0U;
  c2_info[5].context = "[I]$matlabroot$/toolbox/eml/lib/matlab/matfun/trace.m";
  c2_info[5].name = "plus";
  c2_info[5].dominantType = "double";
  c2_info[5].resolved = "[B]plus";
  c2_info[5].fileLength = 0U;
  c2_info[5].fileTime1 = 0U;
  c2_info[5].fileTime2 = 0U;
  c2_info[6].context = "[I]$matlabroot$/toolbox/eml/lib/matlab/matfun/trace.m";
  c2_info[6].name = "isempty";
  c2_info[6].dominantType = "double";
  c2_info[6].resolved = "[B]isempty";
  c2_info[6].fileLength = 0U;
  c2_info[6].fileTime1 = 0U;
  c2_info[6].fileTime2 = 0U;
  c2_info[7].context = "[I]$matlabroot$/toolbox/eml/lib/matlab/datatypes/isa.m";
  c2_info[7].name = "eml_is_float_class";
  c2_info[7].dominantType = "char";
  c2_info[7].resolved =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/eml/eml_is_float_class.m";
  c2_info[7].fileLength = 273U;
  c2_info[7].fileTime1 = 1160446744U;
  c2_info[7].fileTime2 = 0U;
  c2_info[8].context = "[I]$matlabroot$/toolbox/eml/lib/matlab/elmat/intmax.m";
  c2_info[8].name = "int32";
  c2_info[8].dominantType = "double";
  c2_info[8].resolved = "[B]int32";
  c2_info[8].fileLength = 0U;
  c2_info[8].fileTime1 = 0U;
  c2_info[8].fileTime2 = 0U;
  c2_info[9].context =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/eml/eml_check_dim.m";
  c2_info[9].name = "floor";
  c2_info[9].dominantType = "double";
  c2_info[9].resolved = "[I]$matlabroot$/toolbox/eml/lib/matlab/elfun/floor.m";
  c2_info[9].fileLength = 595U;
  c2_info[9].fileTime1 = 1163970736U;
  c2_info[9].fileTime2 = 0U;
  c2_info[10].context = "[I]$matlabroot$/toolbox/eml/lib/matlab/matfun/trace.m";
  c2_info[10].name = "min";
  c2_info[10].dominantType = "double";
  c2_info[10].resolved = "[I]$matlabroot$/toolbox/eml/lib/matlab/datafun/min.m";
  c2_info[10].fileLength = 672U;
  c2_info[10].fileTime1 = 1160446610U;
  c2_info[10].fileTime2 = 0U;
  c2_info[11].context = "[I]$matlabroot$/toolbox/eml/lib/matlab/datatypes/isa.m";
  c2_info[11].name = "eq";
  c2_info[11].dominantType = "double";
  c2_info[11].resolved = "[B]eq";
  c2_info[11].fileLength = 0U;
  c2_info[11].fileTime1 = 0U;
  c2_info[11].fileTime2 = 0U;
  c2_info[12].context =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/eml/eml_min_or_max.m/eml_extremum_sub";
  c2_info[12].name = "eml_index_plus";
  c2_info[12].dominantType = "int32";
  c2_info[12].resolved =