c2_fibonacci_determinant_ml71.c

m files for SIAM book, Numerical Computing with MATLAB

/* Include files */
#include "Fibonacci_determinant_ML71_sfun.h"
#include "c2_Fibonacci_determinant_ML71.h"
#define CHARTINSTANCE_CHARTNUMBER       (chartInstance.chartNumber)
#define CHARTINSTANCE_INSTANCENUMBER    (chartInstance.instanceNumber)
#include "Fibonacci_determinant_ML71_sfun_debug_macros.h"

/* Type Definitions */

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

/* Variable Declarations */

/* Variable Definitions */
static SFc2_Fibonacci_determinant_ML71InstanceStruct chartInstance;

/* Function Declarations */
static void initialize_c2_Fibonacci_determinant_ML71(void);
static void initialize_params_c2_Fibonacci_determinant_ML71(void);
static void enable_c2_Fibonacci_determinant_ML71(void);
static void disable_c2_Fibonacci_determinant_ML71(void);
static void finalize_c2_Fibonacci_determinant_ML71(void);
static void sf_c2_Fibonacci_determinant_ML71(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 void init_dsm_address_info(void);

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

static void initialize_params_c2_Fibonacci_determinant_ML71(void)
{
}

static void enable_c2_Fibonacci_determinant_ML71(void)
{
}

static void disable_c2_Fibonacci_determinant_ML71(void)
{
}

static void finalize_c2_Fibonacci_determinant_ML71(void)
{
}

static void sf_c2_Fibonacci_determinant_ML71(void)
{
  int32_T c2_i0;
  int32_T c2_i1;
  uint8_T c2_previousEvent;
  int32_T c2_i2;
  int32_T c2_i3;
  real_T c2_u[4];
  real_T c2_d;
  real_T *c2_b_d;
  real_T (*c2_b_u)[4];
  c2_b_u = (real_T (*)[4])ssGetInputPortSignal(chartInstance.S, 0);
  c2_b_d = (real_T *)ssGetOutputPortSignal(chartInstance.S, 1);
  _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_u)[c2_i1 + 2 * c2_i0], 0U);
    }
  }
  _SFD_DATA_RANGE_CHECK(*c2_b_d, 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_u[c2_i3 + 2 * c2_i2] = (*c2_b_u)[c2_i3 + 2 * c2_i2];
    }
  }
  sf_debug_symbol_scope_push(2U, 0U);
  sf_debug_symbol_scope_add("d", &c2_d, c2_b_sf_marshall);
  sf_debug_symbol_scope_add("u", &c2_u, c2_sf_marshall);
  CV_EML_FCN(0, 0);
  /*  This block supports an embeddable subset of the MATLAB language. */
  /*  See the help menu for details.  */
  _SFD_EML_CALL(STATE_DURING_DURING_ACTION_TAG,0,5);
  c2_d = c2_u[0] * c2_u[3] - c2_u[2] * c2_u[1];
  _SFD_EML_CALL(STATE_DURING_DURING_ACTION_TAG,0,-5);
  sf_debug_symbol_scope_pop();
  *c2_b_d = c2_d;
  _SFD_CC_CALL(EXIT_OUT_OF_FUNCTION_TAG,0);
  _sfEvent_ = c2_previousEvent;
  sf_debug_check_for_state_inconsistency(_Fibonacci_determinant_ML71MachineNumber_,
   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_i4;
  int32_T c2_i5;
  real_T c2_b_u[4];
  const mxArray *c2_b_y = NULL;
  c2_y = NULL;
  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_b_u[c2_i5 + 2 * c2_i4] = (*((real_T (*)[4])c2_u))[c2_i5 + 2 * c2_i4];
    }
  }
  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;
}

const mxArray
*sf_c2_Fibonacci_determinant_ML71_get_eml_resolved_functions_info(void)
{
  const mxArray *c2_nameCaptureInfo = NULL;
  c2_ResolvedFunctionInfo c2_info[2];
  c2_ResolvedFunctionInfo (*c2_b_info)[2];
  const mxArray *c2_m0 = NULL;
  int32_T c2_i6;
  c2_ResolvedFunctionInfo *c2_r0;
  c2_nameCaptureInfo = NULL;
  c2_b_info = &c2_info;
  (*c2_b_info)[0].context = "";
  (*c2_b_info)[0].name = "mtimes";
  (*c2_b_info)[0].dominantType = "double";
  (*c2_b_info)[0].resolved = "[B]mtimes";
  (*c2_b_info)[0].fileLength = 0U;
  (*c2_b_info)[0].fileTime1 = 0U;
  (*c2_b_info)[0].fileTime2 = 0U;
  (*c2_b_info)[1].context = "";
  (*c2_b_info)[1].name = "minus";
  (*c2_b_info)[1].dominantType = "double";
  (*c2_b_info)[1].resolved = "[B]minus";
  (*c2_b_info)[1].fileLength = 0U;
  (*c2_b_info)[1].fileTime1 = 0U;
  (*c2_b_info)[1].fileTime2 = 0U;
  sf_mex_assign(&c2_m0, sf_mex_createstruct("nameCaptureInfo", 1, 2));
  for(c2_i6 = 0; c2_i6 < 2; c2_i6 = c2_i6 + 1) {
    c2_r0 = &c2_info[c2_i6];
    sf_mex_addfield(c2_m0, sf_mex_create(c2_r0->context, "nameCaptureInfo", 13,
      0U, 0U, 2, 1, strlen(c2_r0->context)), "context",
     "nameCaptureInfo", c2_i6);
    sf_mex_addfield(c2_m0, sf_mex_create(c2_r0->name, "nameCaptureInfo", 13, 0U,
      0U, 2, 1, strlen(c2_r0->name)), "name",
     "nameCaptureInfo", c2_i6);
    sf_mex_addfield(c2_m0, sf_mex_create(c2_r0->dominantType, "nameCaptureInfo",
      13, 0U, 0U, 2, 1, strlen(c2_r0->dominantType)),
     "dominantType", "nameCaptureInfo", c2_i6);
    sf_mex_addfield(c2_m0, sf_mex_create(c2_r0->resolved, "nameCaptureInfo", 13,
      0U, 0U, 2, 1, strlen(c2_r0->resolved)), "resolved",
     "nameCaptureInfo", c2_i6);
    sf_mex_addfield(c2_m0, sf_mex_create(&c2_r0->fileLength, "nameCaptureInfo",
      7, 0U, 0U, 0), "fileLength", "nameCaptureInfo", c2_i6);
    sf_mex_addfield(c2_m0, sf_mex_create(&c2_r0->fileTime1, "nameCaptureInfo",
      7, 0U, 0U, 0), "fileTime1", "nameCaptureInfo", c2_i6);
    sf_mex_addfield(c2_m0, sf_mex_create(&c2_r0->fileTime2, "nameCaptureInfo",
      7, 0U, 0U, 0), "fileTime2", "nameCaptureInfo", c2_i6);
  }
  sf_mex_assign(&c2_nameCaptureInfo, c2_m0);
  return c2_nameCaptureInfo;
}

static void init_dsm_address_info(void)
{
}

/* SFunction Glue Code */
void sf_c2_Fibonacci_determinant_ML71_get_check_sum(mxArray *plhs[])
{
  ((real_T *)mxGetPr((plhs[0])))[0] = (real_T)(3321098180U);
  ((real_T *)mxGetPr((plhs[0])))[1] = (real_T)(400844051U);
  ((real_T *)mxGetPr((plhs[0])))[2] = (real_T)(1265011157U);
  ((real_T *)mxGetPr((plhs[0])))[3] = (real_T)(2893462783U);
}

mxArray *sf_c2_Fibonacci_determinant_ML71_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)(4127722219U);
    pr[1] = (double)(1367233816U);
    pr[2] = (double)(3468663682U);
    pr[3] = (double)(2893138099U);
    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);