c2_fib_determinant.c

m files for SIAM book, Numerical Computing with MATLAB

  c2_info[52].fileLength = 0U;
  c2_info[52].fileTime1 = 0U;
  c2_info[52].fileTime2 = 0U;
  c2_info[53].context =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/elmat/isfinite.m";
  c2_info[53].name = "isinf";
  c2_info[53].dominantType = "double";
  c2_info[53].resolved = "[I]$matlabroot$/toolbox/eml/lib/matlab/elmat/isinf.m";
  c2_info[53].fileLength = 514U;
  c2_info[53].fileTime1 = 1160446702U;
  c2_info[53].fileTime2 = 0U;
  c2_info[54].context =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/ops/colon.m/unrounded_npoints";
  c2_info[54].name = "mrdivide";
  c2_info[54].dominantType = "double";
  c2_info[54].resolved = "[I]$matlabroot$/toolbox/eml/lib/matlab/ops/mrdivide.m";
  c2_info[54].fileLength = 653U;
  c2_info[54].fileTime1 = 1160446802U;
  c2_info[54].fileTime2 = 0U;
  c2_info[55].context = "[I]$matlabroot$/toolbox/eml/lib/matlab/datafun/max.m";
  c2_info[55].name = "eml_min_or_max";
  c2_info[55].dominantType = "char";
  c2_info[55].resolved =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/eml/eml_min_or_max.m";
  c2_info[55].fileLength = 7178U;
  c2_info[55].fileTime1 = 1163970800U;
  c2_info[55].fileTime2 = 0U;
  c2_info[56].context = "[I]$matlabroot$/toolbox/eml/lib/matlab/elmat/realmin.m";
  c2_info[56].name = "uminus";
  c2_info[56].dominantType = "double";
  c2_info[56].resolved = "[B]uminus";
  c2_info[56].fileLength = 0U;
  c2_info[56].fileTime1 = 0U;
  c2_info[56].fileTime2 = 0U;
  c2_info[57].context = "[I]$matlabroot$/toolbox/eml/lib/matlab/ops/colon.m";
  c2_info[57].name = "isnumeric";
  c2_info[57].dominantType = "double";
  c2_info[57].resolved =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/datatypes/isnumeric.m";
  c2_info[57].fileLength = 250U;
  c2_info[57].fileTime1 = 1160446622U;
  c2_info[57].fileTime2 = 0U;
  c2_info[58].context =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/elmat/isfinite.m";
  c2_info[58].name = "and";
  c2_info[58].dominantType = "logical";
  c2_info[58].resolved = "[B]and";
  c2_info[58].fileLength = 0U;
  c2_info[58].fileTime1 = 0U;
  c2_info[58].fileTime2 = 0U;
  c2_info[59].context = "[I]$matlabroot$/toolbox/eml/lib/matlab/ops/mrdivide.m";
  c2_info[59].name = "rdivide";
  c2_info[59].dominantType = "double";
  c2_info[59].resolved = "[I]$matlabroot$/toolbox/eml/lib/matlab/ops/rdivide.m";
  c2_info[59].fileLength = 4533U;
  c2_info[59].fileTime1 = 1177099748U;
  c2_info[59].fileTime2 = 0U;
  c2_info[60].context =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/ops/colon.m/checkrange";
  c2_info[60].name = "realmax";
  c2_info[60].dominantType = "char";
  c2_info[60].resolved =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/elmat/realmax.m";
  c2_info[60].fileLength = 672U;
  c2_info[60].fileTime1 = 1160446716U;
  c2_info[60].fileTime2 = 0U;
  c2_info[61].context = "[I]$matlabroot$/toolbox/eml/lib/matlab/eml/eml_alloc.m";
  c2_info[61].name = "plus";
  c2_info[61].dominantType = "double";
  c2_info[61].resolved = "[B]plus";
  c2_info[61].fileLength = 0U;
  c2_info[61].fileTime1 = 0U;
  c2_info[61].fileTime2 = 0U;
  c2_info[62].context =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/eml/eml_is_float_class.m";
  c2_info[62].name = "ischar";
  c2_info[62].dominantType = "char";
  c2_info[62].resolved =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/strfun/ischar.m";
  c2_info[62].fileLength = 244U;
  c2_info[62].fileTime1 = 1160446834U;
  c2_info[62].fileTime2 = 0U;
  c2_info[63].context =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/eml/eml_min_or_max.m/eml_bin_extremum";
  c2_info[63].name = "eml_alloc";
  c2_info[63].dominantType = "double";
  c2_info[63].resolved =
    "[I]$matlabroot$/toolbox/eml/lib/matlab/eml/eml_alloc.m";
  c2_info[63].fileLength = 2151U;
  c2_info[63].fileTime1 = 1179464470U;
  c2_info[63].fileTime2 = 0U;
}

static void init_dsm_address_info(void)
{
}

/* SFunction Glue Code */
void sf_c2_Fib_determinant_get_check_sum(mxArray *plhs[])
{
  ((real_T *)mxGetPr((plhs[0])))[0] = (real_T)(858259405U);
  ((real_T *)mxGetPr((plhs[0])))[1] = (real_T)(2501850434U);
  ((real_T *)mxGetPr((plhs[0])))[2] = (real_T)(2529172946U);
  ((real_T *)mxGetPr((plhs[0])))[3] = (real_T)(4114925903U);
}

mxArray *sf_c2_Fib_determinant_get_autoinheritance_info(void)
{
  const char *autoinheritanceFields[] =
  {"checksum","inputs","parameters","outputs"};
  mxArray *mxAutoinheritanceInfo =
  mxCreateStructMatrix(1,1,4,autoinheritanceFields);
  {
    mxArray *mxChecksum = mxCreateDoubleMatrix(4,1,mxREAL);
    double *pr = mxGetPr(mxChecksum);
    pr[0] = (double)(589133400U);
    pr[1] = (double)(1083622324U);
    pr[2] = (double)(1033889713U);
    pr[3] = (double)(530004206U);
    mxSetField(mxAutoinheritanceInfo,0,"checksum",mxChecksum);
  }
  {
    const char *dataFields[] = {"size","type","complexity"};
    mxArray *mxData = mxCreateStructMatrix(1,1,3,dataFields);
    {
      mxArray *mxSize = mxCreateDoubleMatrix(1,2,mxREAL);
      double *pr = mxGetPr(mxSize);
      pr[0] = (double)(2);
      pr[1] = (double)(2);
      mxSetField(mxData,0,"size",mxSize);
    }
    {
      const char *typeFields[] = {"base","aliasId","fixpt"};
      mxArray *mxType = mxCreateStructMatrix(1,1,3,typeFields);
      mxSetField(mxType,0,"base",mxCreateDoubleScalar(10));
      mxSetField(mxType,0,"aliasId",mxCreateDoubleScalar(0));
      mxSetField(mxType,0,"fixpt",mxCreateDoubleMatrix(0,0,mxREAL));
      mxSetField(mxData,0,"type",mxType);
    }
    mxSetField(mxData,0,"complexity",mxCreateDoubleScalar(0));
    mxSetField(mxAutoinheritanceInfo,0,"inputs",mxData);
  }
  {
    mxSetField(mxAutoinheritanceInfo,0,"parameters",mxCreateDoubleMatrix(0,0,mxREAL));
  }
  {
    const char *dataFields[] = {"size","type","complexity"};
    mxArray *mxData = mxCreateStructMatrix(1,1,3,dataFields);
    {
      mxArray *mxSize = mxCreateDoubleMatrix(1,2,mxREAL);
      double *pr = mxGetPr(mxSize);
      pr[0] = (double)(1);
      pr[1] = (double)(1);
      mxSetField(mxData,0,"size",mxSize);
    }
    {
      const char *typeFields[] = {"base","aliasId","fixpt"};
      mxArray *mxType = mxCreateStructMatrix(1,1,3,typeFields);
      mxSetField(mxType,0,"base",mxCreateDoubleScalar(10));
      mxSetField(mxType,0,"aliasId",mxCreateDoubleScalar(0));
      mxSetField(mxType,0,"fixpt",mxCreateDoubleMatrix(0,0,mxREAL));
      mxSetField(mxData,0,"type",mxType);
    }
    mxSetField(mxData,0,"complexity",mxCreateDoubleScalar(0));
    mxSetField(mxAutoinheritanceInfo,0,"outputs",mxData);
  }
  return(mxAutoinheritanceInfo);
}

static void chart_debug_initialization(SimStruct *S, unsigned int
 fullDebuggerInitialization)
{
  if(ssIsFirstInitCond(S) && fullDebuggerInitialization==1) {
    /* do this only if simulation is starting */
    if(!sim_mode_is_rtw_gen(S)) {
      {
        unsigned int chartAlreadyPresent;
        chartAlreadyPresent =
          sf_debug_initialize_chart(_Fib_determinantMachineNumber_,
         2,
         1,
         1,
         2,
         0,
         0,
         0,
         0,
         0,
         &(chartInstance.chartNumber),
         &(chartInstance.instanceNumber),
         ssGetPath(S),
         (void *)S);
        if(chartAlreadyPresent==0) {
          /* this is the first instance */
          init_script_number_translation(_Fib_determinantMachineNumber_,chartInstance.chartNumber);
          sf_debug_set_chart_disable_implicit_casting(_Fib_determinantMachineNumber_,chartInstance.chartNumber,1);
          sf_debug_set_chart_event_thresholds(_Fib_determinantMachineNumber_,
           chartInstance.chartNumber,
           0,
           0,
           0);

          {
            unsigned int dimVector[2];
            dimVector[0]= 2;
            dimVector[1]= 2;
            _SFD_SET_DATA_PROPS(0,1,1,0,SF_DOUBLE,2,&(dimVector[0]),0,0,0,0.0,1.0,0,"u",0,c2_sf_marshall);
          }
          _SFD_SET_DATA_PROPS(1,2,0,1,SF_DOUBLE,0,NULL,0,0,0,0.0,1.0,0,"d",0,c2_b_sf_marshall);
          _SFD_STATE_INFO(0,0,2);
          _SFD_CH_SUBSTATE_COUNT(0);
          _SFD_CH_SUBSTATE_DECOMP(0);
        }
        _SFD_CV_INIT_CHART(0,0,0,0);
        {
          _SFD_CV_INIT_STATE(0,0,0,0,0,0,NULL,NULL);
        }

        _SFD_CV_INIT_TRANS(0,0,NULL,NULL,0,NULL);

        /* Initialization of EML Model Coverage */
        _SFD_CV_INIT_EML(0,1,0,0,0,0,0,0);
        _SFD_CV_INIT_EML_FCN(0,0,"eML_blk_kernel",0,-1,286);
        _SFD_TRANS_COV_WTS(0,0,0,1,0);
        if(chartAlreadyPresent==0)
        {
          _SFD_TRANS_COV_MAPS(0,
           0,NULL,NULL,
           0,NULL,NULL,
           1,NULL,NULL,
           0,NULL,NULL);
        }

        {
          real_T (*c2_u)[4];
          real_T *c2_d;
          c2_u = (real_T (*)[4])ssGetInputPortSignal(chartInstance.S, 0);
          c2_d = (real_T *)ssGetOutputPortSignal(chartInstance.S, 1);
          _SFD_SET_DATA_VALUE_PTR(0U, c2_u);
          _SFD_SET_DATA_VALUE_PTR(1U, c2_d);
        }
      }
    }
  } else {
    sf_debug_reset_current_state_configuration(_Fib_determinantMachineNumber_,chartInstance.chartNumber,chartInstance.instanceNumber);
  }
}

static void sf_opaque_initialize_c2_Fib_determinant(void *chartInstanceVar)
{
  chart_debug_initialization(chartInstance.S,0);
  initialize_params_c2_Fib_determinant();
  initialize_c2_Fib_determinant();
}

static void sf_opaque_enable_c2_Fib_determinant(void *chartInstanceVar)
{
  enable_c2_Fib_determinant();
}

static void sf_opaque_disable_c2_Fib_determinant(void *chartInstanceVar)
{
  disable_c2_Fib_determinant();
}

static void sf_opaque_gateway_c2_Fib_determinant(void *chartInstanceVar)
{
  sf_c2_Fib_determinant();
}

static void sf_opaque_terminate_c2_Fib_determinant(void *chartInstanceVar)
{
  if (sim_mode_is_rtw_gen(chartInstance.S) ||
   sim_mode_is_external(chartInstance.S)) {
    sf_clear_rtw_identifier(chartInstance.S);
  }
  finalize_c2_Fib_determinant();
}

extern unsigned int sf_machine_global_initializer_called(void);
static void mdlProcessParameters_c2_Fib_determinant(SimStruct *S)
{
  int i;
  for(i=0;i<ssGetNumRunTimeParams(S);i++) {
    if(ssGetSFcnParamTunable(S,i)) {
      ssUpdateDlgParamAsRunTimeParam(S,i);
    }
  }
  if(sf_machine_global_initializer_called()) {
    initialize_params_c2_Fib_determinant();
  }
}

static void mdlSetWorkWidths_c2_Fib_determinant(SimStruct *S)
{
  if(sim_mode_is_rtw_gen(S) || sim_mode_is_external(S)) {
    int_T chartIsInlinable =
      (int_T)sf_is_chart_inlinable("Fib_determinant",2);
    ssSetStateflowIsInlinable(S,chartIsInlinable);
    ssSetEnableFcnIsTrivial(S,1);
    ssSetDisableFcnIsTrivial(S,1);
    ssSetNotMultipleInlinable(S,sf_rtw_info_uint_prop("Fib_determinant",2,"gatewayCannotBeInlinedMultipleTimes"));
    if(chartIsInlinable) {
      ssSetInputPortOptimOpts(S, 0, SS_REUSABLE_AND_LOCAL);
      sf_mark_chart_expressionable_inputs(S,"Fib_determinant",2,1);
      sf_mark_chart_reusable_outputs(S,"Fib_determinant",2,1);
    }
    sf_set_rtw_dwork_info(S,"Fib_determinant",2);
    ssSetHasSubFunctions(S,!(chartIsInlinable));
    ssSetOptions(S,ssGetOptions(S)|SS_OPTION_WORKS_WITH_CODE_REUSE);
  }

  ssSetChecksum0(S,(1628595885U));
  ssSetChecksum1(S,(2593892068U));
  ssSetChecksum2(S,(3223411231U));
  ssSetChecksum3(S,(1805122284U));

  ssSetmdlDerivatives(S, NULL);

  ssSetExplicitFCSSCtrl(S,1);
}

static void mdlRTW_c2_Fib_determinant(SimStruct *S)
{
  if(sim_mode_is_rtw_gen(S)) {
    sf_write_symbol_mapping(S, "Fib_determinant", 2);
    ssWriteRTWStrParam(S, "StateflowChartType", "Embedded MATLAB");
  }
}

static void mdlStart_c2_Fib_determinant(SimStruct *S)
{
  chartInstance.chartInfo.chartInstance = NULL;
  chartInstance.chartInfo.isEMLChart = 1;
  chartInstance.chartInfo.chartInitialized = 0;
  chartInstance.chartInfo.sFunctionGateway =
    sf_opaque_gateway_c2_Fib_determinant;
  chartInstance.chartInfo.initializeChart =
    sf_opaque_initialize_c2_Fib_determinant;
  chartInstance.chartInfo.terminateChart =
    sf_opaque_terminate_c2_Fib_determinant;
  chartInstance.chartInfo.enableChart = sf_opaque_enable_c2_Fib_determinant;
  chartInstance.chartInfo.disableChart = sf_opaque_disable_c2_Fib_determinant;
  chartInstance.chartInfo.zeroCrossings = NULL;
  chartInstance.chartInfo.outputs = NULL;
  chartInstance.chartInfo.derivatives = NULL;
  chartInstance.chartInfo.mdlRTW = mdlRTW_c2_Fib_determinant;
  chartInstance.chartInfo.mdlStart = mdlStart_c2_Fib_determinant;
  chartInstance.chartInfo.mdlSetWorkWidths = mdlSetWorkWidths_c2_Fib_determinant;
  chartInstance.chartInfo.extModeExec = NULL;
  chartInstance.chartInfo.restoreLastMajorStepConfiguration = NULL;
  chartInstance.chartInfo.restoreBeforeLastMajorStepConfiguration = NULL;
  chartInstance.chartInfo.storeCurrentConfiguration = NULL;
  chartInstance.S = S;
  ssSetUserData(S,(void *)(&(chartInstance.chartInfo))); /* register the chart instance with simstruct */

  if(!sim_mode_is_rtw_gen(S)) {
    init_dsm_address_info();
  }
  chart_debug_initialization(S,1);
}

void c2_Fib_determinant_method_dispatcher(SimStruct *S, int_T method, void *data)
{
  switch (method) {
   case SS_CALL_MDL_START:
    mdlStart_c2_Fib_determinant(S);
    break;
   case SS_CALL_MDL_SET_WORK_WIDTHS:
    mdlSetWorkWidths_c2_Fib_determinant(S);
    break;
   case SS_CALL_MDL_PROCESS_PARAMETERS:
    mdlProcessParameters_c2_Fib_determinant(S);
    break;
   default:
    /* Unhandled method */
    sf_mex_error_message("Stateflow Internal Error:\n"
     "Error calling c2_Fib_determinant_method_dispatcher.\n"
     "Can't handle method %d.\n", method);
    break;
  }
}