static char mc_version[] = "MATLAB Compiler 1.2.1 Jan 21 1999 infun"; /* * MATLAB Compiler: 1.2.1 * Date: Jan 21 1999 * Arguments: -Z -i -r -w simulate intact_init_cond hlu_init_cond sc_init_cond lv_init_cond init_cond nac_init_cond third_init_cond third_nac_init_cond apr_init_cond a_init_cond rk4 intact_eval_deriv hlu_eval_deriv sc_eval_deriv lv_eval_deriv eval_deriv third_eval_deriv nac_eval_deriv third_nac_eval_deriv apr_eval_deriv a_eval_deriv var_cap var_vcap var_acap vent_vol resp_act rand_int_breath abreflex dncm_filt ilv_dec cpreflex ans_filter oneoverf_filt bl_filt sinc conserve_vol read_param param_change */ #ifndef ARRAY_ACCESS_INLINING #error You must use the -inline option when compiling MATLAB compiler generated code with MEX or MBUILD #endif #ifndef MATLAB_COMPILER_GENERATED_CODE #define MATLAB_COMPILER_GENERATED_CODE #endif #include #include "mex.h" #include "mcc.h" static void rand_int_breath_(mxArray *); static void resp_act_(mxArray *, mxArray *, mxArray *, double, int, double); static void dncm_filt_(mxArray *, mxArray *, mxArray *); static void ilv_dec_(mxArray *, mxArray *, double, double); static void bl_filt_(mxArray *, double, double); static void oneoverf_filt_(mxArray *, mxArray *, double, int, int); static void intact_init_cond_(mxArray *,mxArray *,mxArray *,mxArray *, mxArray *); static void hlu_init_cond_(mxArray *,mxArray *,mxArray *,mxArray *, mxArray *); static void sc_init_cond_(mxArray *,mxArray *,mxArray *,mxArray *, mxArray *); static void lv_init_cond_(mxArray *,mxArray *,mxArray *,mxArray *, mxArray *); static void init_cond_(mxArray *,mxArray *,mxArray *,mxArray *, mxArray *); static void third_init_cond_(mxArray *,mxArray *,mxArray *,mxArray *, mxArray *); static void nac_init_cond_(mxArray *,mxArray *,mxArray *,mxArray *, mxArray *); static void third_nac_init_cond_(mxArray *,mxArray *,mxArray *,mxArray *, mxArray *); static void apr_init_cond_(mxArray *,mxArray *,mxArray *,mxArray *, mxArray *); static void a_init_cond_(mxArray *,mxArray *,mxArray *,mxArray *, mxArray *); static void rk4_(mxArray *, mxArray *, mxArray *, mxArray *, mxArray *, double, mxArray *, double, double); static void abreflex_(mxArray *, mxArray *, mxArray *, mxArray *, mxArray *, mxArray *); static void abreflex__1(mxArray *, mxArray *, mxArray *, mxArray *, mxArray *, mxArray *); static void cpreflex_(mxArray *, mxArray *, mxArray *, mxArray *, mxArray *); static void vent_vol_(double *, double, double, double); static void var_cap_(double *,double *, double, double, double, mxArray *); static void var_vcap_(double *,double *, double, double, double, mxArray *); static void var_acap_(double *,double *, double, double, double, mxArray *); static void read_param_(mxArray *,mxArray *, int); static void param_change_(mxArray *, mxArray *, mxArray *, mxArray *); static void resp_act__1(mxArray *, mxArray *, mxArray *, double, int, double); static void resp_act__2(mxArray *, mxArray *, mxArray *, double, double, double); static void resp_act__3(mxArray *, mxArray *, mxArray *, double, double, double); static void conserve_vol_(mxArray *, mxArray *, mxArray *, mxArray *, double); /* static array S0_ (1 x 6) int, line 1096 */ static int S0r_[] = { 1, 0, 0, 1, 1, 1, }; static mxArray S0_ = mccCINIT( mccINT, 1, 6, S0r_, 0 ); /* static array S1_ (1 x 6) int, line 1098 */ static int S1r_[] = { 1, 0, 0, 1, 1, 1, }; static mxArray S1_ = mccCINIT( mccINT, 1, 6, S1r_, 0 ); /* static array S2_ (1 x 2) int, line 1099 */ static int S2r_[] = { 1, 1, }; static mxArray S2_ = mccCINIT( mccINT, 1, 2, S2r_, 0 ); /* static array S3_ (1 x 32) text, line 130: 'flag is not of correct di...' */ static unsigned short S3__r_[] = { 102, 108, 97, 103, 32, 105, 115, 32, 110, 111, 116, 32, 111, 102, 32, 99, 111, 114, 114, 101, 99, 116, 32, 100, 105, 109, 101, 110, 115, 105, 111, 110, }; static mxArray S3_ = mccCINIT( mccTEXT, 1, 32, S3__r_, 0); /* static array S4_ (1 x 31) text, line 240: 'breathing is not a valid ...' */ static unsigned short S4__r_[] = { 98, 114, 101, 97, 116, 104, 105, 110, 103, 32, 105, 115, 32, 110, 111, 116, 32, 97, 32, 118, 97, 108, 105, 100, 32, 111, 112, 116, 105, 111, 110, }; static mxArray S4_ = mccCINIT( mccTEXT, 1, 31, S4__r_, 0); /* static array S5_ (1 x 26) text, line 274: 'dncm is not a valid optio...' */ static unsigned short S5__r_[] = { 100, 110, 99, 109, 32, 105, 115, 32, 110, 111, 116, 32, 97, 32, 118, 97, 108, 105, 100, 32, 111, 112, 116, 105, 111, 110, }; static mxArray S5_ = mccCINIT( mccTEXT, 1, 26, S5__r_, 0); /* static array S6_ (1 x 25) text, line 326: 'dra is not a valid option...' */ static unsigned short S6__r_[] = { 100, 114, 97, 32, 105, 115, 32, 110, 111, 116, 32, 97, 32, 118, 97, 108, 105, 100, 32, 111, 112, 116, 105, 111, 110, }; static mxArray S6_ = mccCINIT( mccTEXT, 1, 25, S6__r_, 0); /* static array S7_ (1 x 18) text, line 354: 'dpasp is not valid' */ static unsigned short S7__r_[] = { 100, 112, 97, 115, 112, 32, 105, 115, 32, 110, 111, 116, 32, 118, 97, 108, 105, 100, }; static mxArray S7_ = mccCINIT( mccTEXT, 1, 18, S7__r_, 0); /* static array S8_ (1 x 24) text, line 371: 'df is not a valid option' */ static unsigned short S8__r_[] = { 100, 102, 32, 105, 115, 32, 110, 111, 116, 32, 97, 32, 118, 97, 108, 105, 100, 32, 111, 112, 116, 105, 111, 110, }; static mxArray S8_ = mccCINIT( mccTEXT, 1, 24, S8__r_, 0); /* static array S9_ (1 x 26) text, line 389: 'dqvo is not a valid optio...' */ static unsigned short S9__r_[] = { 100, 113, 118, 111, 32, 105, 115, 32, 110, 111, 116, 32, 97, 32, 118, 97, 108, 105, 100, 32, 111, 112, 116, 105, 111, 110, }; static mxArray S9_ = mccCINIT( mccTEXT, 1, 26, S9__r_, 0); /* static array S10_ (1 x 33) text, line 556: 'preparation is not a vali...' */ static unsigned short S10__r_[] = { 112, 114, 101, 112, 97, 114, 97, 116, 105, 111, 110, 32, 105, 115, 32, 110, 111, 116, 32, 97, 32, 118, 97, 108, 105, 100, 32, 111, 112, 116, 105, 111, 110, }; static mxArray S10_ = mccCINIT( mccTEXT, 1, 33, S10__r_, 0); /* static array S11_ (1 x 26) text, line 623: 'baro is not a valid optio...' */ static unsigned short S11__r_[] = { 98, 97, 114, 111, 32, 105, 115, 32, 110, 111, 116, 32, 97, 32, 118, 97, 108, 105, 100, 32, 111, 112, 116, 105, 111, 110, }; static mxArray S11_ = mccCINIT( mccTEXT, 1, 26, S11__r_, 0); /* static array S12_ (1 x 2) text, line 636: '-1' */ static unsigned short S12__r_[] = { 45, 49, }; static mxArray S12_ = mccCINIT( mccTEXT, 1, 2, S12__r_, 0); /* static array S13_ (1 x 4) text, line 641: '.dat' */ static unsigned short S13__r_[] = { 46, 100, 97, 116, }; static mxArray S13_ = mccCINIT( mccTEXT, 1, 4, S13__r_, 0); /* static array S14_ (1 x 2) text, line 641: 'wa' */ static unsigned short S14__r_[] = { 119, 97, }; static mxArray S14_ = mccCINIT( mccTEXT, 1, 2, S14__r_, 0); /* static array S15_ (1 x 4) text, line 644: '.qrs' */ static unsigned short S15__r_[] = { 46, 113, 114, 115, }; static mxArray S15_ = mccCINIT( mccTEXT, 1, 4, S15__r_, 0); /* static array S16_ (1 x 2) text, line 644: 'wa' */ static unsigned short S16__r_[] = { 119, 97, }; static mxArray S16_ = mccCINIT( mccTEXT, 1, 2, S16__r_, 0); /* static array S17_ (1 x 5) text, line 645: 'bit10' */ static unsigned short S17__r_[] = { 98, 105, 116, 49, 48, }; static mxArray S17_ = mccCINIT( mccTEXT, 1, 5, S17__r_, 0); /* static array S18_ (1 x 5) text, line 646: 'ubit6' */ static unsigned short S18__r_[] = { 117, 98, 105, 116, 54, }; static mxArray S18_ = mccCINIT( mccTEXT, 1, 5, S18__r_, 0); /* static array S19_ (1 x 5) text, line 647: 'bit10' */ static unsigned short S19__r_[] = { 98, 105, 116, 49, 48, }; static mxArray S19_ = mccCINIT( mccTEXT, 1, 5, S19__r_, 0); /* static array S20_ (1 x 5) text, line 648: 'ubit6' */ static unsigned short S20__r_[] = { 117, 98, 105, 116, 54, }; static mxArray S20_ = mccCINIT( mccTEXT, 1, 5, S20__r_, 0); /* static array S21_ (1 x 5) text, line 650: 'bit10' */ static unsigned short S21__r_[] = { 98, 105, 116, 49, 48, }; static mxArray S21_ = mccCINIT( mccTEXT, 1, 5, S21__r_, 0); /* static array S22_ (1 x 5) text, line 651: 'ubit6' */ static unsigned short S22__r_[] = { 117, 98, 105, 116, 54, }; static mxArray S22_ = mccCINIT( mccTEXT, 1, 5, S22__r_, 0); /* static array S23_ (1 x 1) text, line 652: '.' */ static unsigned short S23__r_[] = { 46, }; static mxArray S23_ = mccCINIT( mccTEXT, 1, 1, S23__r_, 0); /* static array S24_ (1 x 6) text, line 652: '*uchar' */ static unsigned short S24__r_[] = { 42, 117, 99, 104, 97, 114, }; static mxArray S24_ = mccCINIT( mccTEXT, 1, 6, S24__r_, 0); /* static array S25_ (1 x 2) text, line 670: './' */ static unsigned short S25__r_[] = { 46, 47, }; static mxArray S25_ = mccCINIT( mccTEXT, 1, 2, S25__r_, 0); /* static array S26_ (1 x 2) text, line 1165: '-1' */ static unsigned short S26__r_[] = { 45, 49, }; static mxArray S26_ = mccCINIT( mccTEXT, 1, 2, S26__r_, 0); /* static array S27_ (1 x 1) text, line 1174: 'r' */ static unsigned short S27__r_[] = { 114, }; static mxArray S27_ = mccCINIT( mccTEXT, 1, 1, S27__r_, 0); /* static array S28_ (1 x 2) text, line 1567: '-1' */ static unsigned short S28__r_[] = { 45, 49, }; static mxArray S28_ = mccCINIT( mccTEXT, 1, 2, S28__r_, 0); /* static array S29_ (1 x 5) text, line 1571: 'short' */ static unsigned short S29__r_[] = { 115, 104, 111, 114, 116, }; static mxArray S29_ = mccCINIT( mccTEXT, 1, 5, S29__r_, 0); /* static array S30_ (1 x 5) text, line 1581: 'bit10' */ static unsigned short S30__r_[] = { 98, 105, 116, 49, 48, }; static mxArray S30_ = mccCINIT( mccTEXT, 1, 5, S30__r_, 0); /* static array S31_ (1 x 5) text, line 1582: 'ubit6' */ static unsigned short S31__r_[] = { 117, 98, 105, 116, 54, }; static mxArray S31_ = mccCINIT( mccTEXT, 1, 5, S31__r_, 0); /* static array S32_ (1 x 5) text, line 1589: 'bit10' */ static unsigned short S32__r_[] = { 98, 105, 116, 49, 48, }; static mxArray S32_ = mccCINIT( mccTEXT, 1, 5, S32__r_, 0); /* static array S33_ (1 x 5) text, line 1590: 'ubit6' */ static unsigned short S33__r_[] = { 117, 98, 105, 116, 54, }; static mxArray S33_ = mccCINIT( mccTEXT, 1, 5, S33__r_, 0); /* static array S34_ (1 x 1) text, line 1591: '.' */ static unsigned short S34__r_[] = { 46, }; static mxArray S34_ = mccCINIT( mccTEXT, 1, 1, S34__r_, 0); /* static array S35_ (1 x 6) text, line 1591: '*uchar' */ static unsigned short S35__r_[] = { 42, 117, 99, 104, 97, 114, }; static mxArray S35_ = mccCINIT( mccTEXT, 1, 6, S35__r_, 0); /* static array S36_ (1 x 2) text, line 1624: '-1' */ static unsigned short S36__r_[] = { 45, 49, }; static mxArray S36_ = mccCINIT( mccTEXT, 1, 2, S36__r_, 0); /* static array S37_ (1 x 5) text, line 1628: 'short' */ static unsigned short S37__r_[] = { 115, 104, 111, 114, 116, }; static mxArray S37_ = mccCINIT( mccTEXT, 1, 5, S37__r_, 0); /* static array S38_ (1 x 5) text, line 1635: 'bit10' */ static unsigned short S38__r_[] = { 98, 105, 116, 49, 48, }; static mxArray S38_ = mccCINIT( mccTEXT, 1, 5, S38__r_, 0); /* static array S39_ (1 x 5) text, line 1636: 'ubit6' */ static unsigned short S39__r_[] = { 117, 98, 105, 116, 54, }; static mxArray S39_ = mccCINIT( mccTEXT, 1, 5, S39__r_, 0); /* static array S40_ (1 x 5) text, line 1644: 'bit10' */ static unsigned short S40__r_[] = { 98, 105, 116, 49, 48, }; static mxArray S40_ = mccCINIT( mccTEXT, 1, 5, S40__r_, 0); /* static array S41_ (1 x 5) text, line 1645: 'ubit6' */ static unsigned short S41__r_[] = { 117, 98, 105, 116, 54, }; static mxArray S41_ = mccCINIT( mccTEXT, 1, 5, S41__r_, 0); /* static array S42_ (1 x 1) text, line 1646: '.' */ static unsigned short S42__r_[] = { 46, }; static mxArray S42_ = mccCINIT( mccTEXT, 1, 1, S42__r_, 0); /* static array S43_ (1 x 6) text, line 1646: '*uchar' */ static unsigned short S43__r_[] = { 42, 117, 99, 104, 97, 114, }; static mxArray S43_ = mccCINIT( mccTEXT, 1, 6, S43__r_, 0); /* static array S44_ (1 x 5) text, line 1651: 'bit10' */ static unsigned short S44__r_[] = { 98, 105, 116, 49, 48, }; static mxArray S44_ = mccCINIT( mccTEXT, 1, 5, S44__r_, 0); /* static array S45_ (1 x 5) text, line 1652: 'ubit6' */ static unsigned short S45__r_[] = { 117, 98, 105, 116, 54, }; static mxArray S45_ = mccCINIT( mccTEXT, 1, 5, S45__r_, 0); /* static array S46_ (1 x 6) int, line 98 */ static int S46r_[] = { 1, 0, 0, 1, 1, 1, }; static mxArray S46_ = mccCINIT( mccINT, 1, 6, S46r_, 0 ); /* static array S47_ (1 x 6) int, line 92 */ static int S47r_[] = { 1, 0, 0, 1, 1, 1, }; static mxArray S47_ = mccCINIT( mccINT, 1, 6, S47r_, 0 ); /* static array S48_ (1 x 2) int, line 44 */ static int S48r_[] = { 0, 0, }; static mxArray S48_ = mccCINIT( mccINT, 1, 2, S48r_, 0 ); /* static array S49_ (1 x 6) int, line 88 */ static int S49r_[] = { 1, 0, 0, 1, 1, 1, }; static mxArray S49_ = mccCINIT( mccINT, 1, 6, S49r_, 0 ); /* static array S50_ (1 x 6) int, line 85 */ static int S50r_[] = { 1, 0, 0, 1, 1, 1, }; static mxArray S50_ = mccCINIT( mccINT, 1, 6, S50r_, 0 ); /* static array S51_ (1 x 6) int, line 91 */ static int S51r_[] = { 1, 0, 0, 1, 1, 1, }; static mxArray S51_ = mccCINIT( mccINT, 1, 6, S51r_, 0 ); /* static array S52_ (1 x 6) int, line 98 */ static int S52r_[] = { 1, 0, 0, 1, 1, 1, }; static mxArray S52_ = mccCINIT( mccINT, 1, 6, S52r_, 0 ); /* static array S53_ (1 x 7) int, line 96 */ static int S53r_[] = { 1, 0, 0, 1, 1, 1, 0, }; static mxArray S53_ = mccCINIT( mccINT, 1, 7, S53r_, 0 ); /* static array S54_ (1 x 7) int, line 97 */ static int S54r_[] = { 1, 0, 0, 1, 1, 1, 0, }; static mxArray S54_ = mccCINIT( mccINT, 1, 7, S54r_, 0 ); /* static array S55_ (1 x 6) int, line 96 */ static int S55r_[] = { 1, 0, 0, 1, 1, 1, }; static mxArray S55_ = mccCINIT( mccINT, 1, 6, S55r_, 0 ); /* static array S56_ (1 x 8) int, line 104 */ static int S56r_[] = { 1, 0, 0, 1, 1, 1, 1, 1, }; static mxArray S56_ = mccCINIT( mccINT, 1, 8, S56r_, 0 ); static void intact_eval_deriv_(mxArray *, mxArray *, mxArray *, mxArray *, mxArray *, mxArray *); static void hlu_eval_deriv_(mxArray *, mxArray *, mxArray *, mxArray *, mxArray *, mxArray *); static void sc_eval_deriv_(mxArray *, mxArray *, mxArray *, mxArray *, mxArray *, mxArray *); static void lv_eval_deriv_(mxArray *, mxArray *, mxArray *, mxArray *, mxArray *, mxArray *); static void eval_deriv_(mxArray *, mxArray *, mxArray *, mxArray *, mxArray *, mxArray *); static void third_eval_deriv_(mxArray *, mxArray *, mxArray *, mxArray *, mxArray *, mxArray *); static void nac_eval_deriv_(mxArray *, mxArray *, mxArray *, mxArray *, mxArray *, mxArray *); static void third_nac_eval_deriv_(mxArray *, mxArray *, mxArray *, mxArray *, mxArray *, mxArray *); static void apr_eval_deriv_(mxArray *, mxArray *, mxArray *, mxArray *, mxArray *, mxArray *); static void a_eval_deriv_(mxArray *, mxArray *, mxArray *, mxArray *, mxArray *, mxArray *); static void vent_vol__1(double *, mxArray *, double, double); /* static array S57_ (1 x 16) text, line 64: 'No root is found' */ static unsigned short S57__r_[] = { 78, 111, 32, 114, 111, 111, 116, 32, 105, 115, 32, 102, 111, 117, 110, 100, }; static mxArray S57_ = mccCINIT( mccTEXT, 1, 16, S57__r_, 0); /* static array S58_ (1 x 16) text, line 64: 'No root is found' */ static unsigned short S58__r_[] = { 78, 111, 32, 114, 111, 111, 116, 32, 105, 115, 32, 102, 111, 117, 110, 100, }; static mxArray S58_ = mccCINIT( mccTEXT, 1, 16, S58__r_, 0); static void sinc_(mxArray *, mxArray *); /* static array S59_ (1 x 4) int, line 24 */ static int S59r_[] = { 2, 3, 5, 6, }; static mxArray S59_ = mccCINIT( mccINT, 1, 4, S59r_, 0 ); /* static array S60_ (1 x 4) int, line 24 */ static int S60r_[] = { 2, 3, 5, 6, }; static mxArray S60_ = mccCINIT( mccINT, 1, 4, S60r_, 0 ); /* static array S61_ (1 x 4) int, line 24 */ static int S61r_[] = { 10, 11, 13, 14, }; static mxArray S61_ = mccCINIT( mccINT, 1, 4, S61r_, 0 ); /* static array S62_ (1 x 4) int, line 24 */ static int S62r_[] = { 3, 4, 7, 8, }; static mxArray S62_ = mccCINIT( mccINT, 1, 4, S62r_, 0 ); /* static array S63_ (4 x 1) int, line 24 */ static int S63r_[] = { 0, 0, 1, 1, }; static mxArray S63_ = mccCINIT( mccINT, 4, 1, S63r_, 0 ); /* static array S64_ (1 x 17) text, line 94: '%%%END OF FILE%%%' */ static unsigned short S64__r_[] = { 37, 37, 37, 69, 78, 68, 32, 79, 70, 32, 70, 73, 76, 69, 37, 37, 37, }; static mxArray S64_ = mccCINIT( mccTEXT, 1, 17, S64__r_, 0); /* static array S65_ (1 x 1) text, line 96: '%' */ static unsigned short S65__r_[] = { 37, }; static mxArray S65_ = mccCINIT( mccTEXT, 1, 1, S65__r_, 0); /* static array S66_ (1 x 1) text, line 100: ':' */ static unsigned short S66__r_[] = { 58, }; static mxArray S66_ = mccCINIT( mccTEXT, 1, 1, S66__r_, 0); /* static array S67_ (1 x 1) text, line 104: ':' */ static unsigned short S67__r_[] = { 58, }; static mxArray S67_ = mccCINIT( mccTEXT, 1, 1, S67__r_, 0); /* static array S68_ (1 x 3) text, line 125: 'Cls' */ static unsigned short S68__r_[] = { 67, 108, 115, }; static mxArray S68_ = mccCINIT( mccTEXT, 1, 3, S68__r_, 0); /* static array S69_ (1 x 5) text, line 126: '%s %f' */ static unsigned short S69__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S69_ = mccCINIT( mccTEXT, 1, 5, S69__r_, 0); /* static array S70_ (1 x 3) text, line 128: 'Cld' */ static unsigned short S70__r_[] = { 67, 108, 100, }; static mxArray S70_ = mccCINIT( mccTEXT, 1, 3, S70__r_, 0); /* static array S71_ (1 x 5) text, line 129: '%s %f' */ static unsigned short S71__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S71_ = mccCINIT( mccTEXT, 1, 5, S71__r_, 0); /* static array S72_ (1 x 2) text, line 131: 'Ca' */ static unsigned short S72__r_[] = { 67, 97, }; static mxArray S72_ = mccCINIT( mccTEXT, 1, 2, S72__r_, 0); /* static array S73_ (1 x 5) text, line 132: '%s %f' */ static unsigned short S73__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S73_ = mccCINIT( mccTEXT, 1, 5, S73__r_, 0); /* static array S74_ (1 x 2) text, line 134: 'Cv' */ static unsigned short S74__r_[] = { 67, 118, }; static mxArray S74_ = mccCINIT( mccTEXT, 1, 2, S74__r_, 0); /* static array S75_ (1 x 5) text, line 135: '%s %f' */ static unsigned short S75__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S75_ = mccCINIT( mccTEXT, 1, 5, S75__r_, 0); /* static array S76_ (1 x 3) text, line 137: 'Crs' */ static unsigned short S76__r_[] = { 67, 114, 115, }; static mxArray S76_ = mccCINIT( mccTEXT, 1, 3, S76__r_, 0); /* static array S77_ (1 x 5) text, line 138: '%s %f' */ static unsigned short S77__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S77_ = mccCINIT( mccTEXT, 1, 5, S77__r_, 0); /* static array S78_ (1 x 3) text, line 140: 'Crd' */ static unsigned short S78__r_[] = { 67, 114, 100, }; static mxArray S78_ = mccCINIT( mccTEXT, 1, 3, S78__r_, 0); /* static array S79_ (1 x 5) text, line 141: '%s %f' */ static unsigned short S79__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S79_ = mccCINIT( mccTEXT, 1, 5, S79__r_, 0); /* static array S80_ (1 x 3) text, line 143: 'Cpa' */ static unsigned short S80__r_[] = { 67, 112, 97, }; static mxArray S80_ = mccCINIT( mccTEXT, 1, 3, S80__r_, 0); /* static array S81_ (1 x 5) text, line 144: '%s %f' */ static unsigned short S81__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S81_ = mccCINIT( mccTEXT, 1, 5, S81__r_, 0); /* static array S82_ (1 x 3) text, line 146: 'Cpv' */ static unsigned short S82__r_[] = { 67, 112, 118, }; static mxArray S82_ = mccCINIT( mccTEXT, 1, 3, S82__r_, 0); /* static array S83_ (1 x 5) text, line 147: '%s %f' */ static unsigned short S83__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S83_ = mccCINIT( mccTEXT, 1, 5, S83__r_, 0); /* static array S84_ (1 x 3) text, line 149: 'Qlo' */ static unsigned short S84__r_[] = { 81, 108, 111, }; static mxArray S84_ = mccCINIT( mccTEXT, 1, 3, S84__r_, 0); /* static array S85_ (1 x 5) text, line 150: '%s %f' */ static unsigned short S85__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S85_ = mccCINIT( mccTEXT, 1, 5, S85__r_, 0); /* static array S86_ (1 x 3) text, line 152: 'Qao' */ static unsigned short S86__r_[] = { 81, 97, 111, }; static mxArray S86_ = mccCINIT( mccTEXT, 1, 3, S86__r_, 0); /* static array S87_ (1 x 5) text, line 153: '%s %f' */ static unsigned short S87__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S87_ = mccCINIT( mccTEXT, 1, 5, S87__r_, 0); /* static array S88_ (1 x 3) text, line 155: 'Qvo' */ static unsigned short S88__r_[] = { 81, 118, 111, }; static mxArray S88_ = mccCINIT( mccTEXT, 1, 3, S88__r_, 0); /* static array S89_ (1 x 5) text, line 156: '%s %f' */ static unsigned short S89__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S89_ = mccCINIT( mccTEXT, 1, 5, S89__r_, 0); /* static array S90_ (1 x 3) text, line 158: 'Qro' */ static unsigned short S90__r_[] = { 81, 114, 111, }; static mxArray S90_ = mccCINIT( mccTEXT, 1, 3, S90__r_, 0); /* static array S91_ (1 x 5) text, line 159: '%s %f' */ static unsigned short S91__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S91_ = mccCINIT( mccTEXT, 1, 5, S91__r_, 0); /* static array S92_ (1 x 4) text, line 161: 'Qpao' */ static unsigned short S92__r_[] = { 81, 112, 97, 111, }; static mxArray S92_ = mccCINIT( mccTEXT, 1, 4, S92__r_, 0); /* static array S93_ (1 x 5) text, line 162: '%s %f' */ static unsigned short S93__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S93_ = mccCINIT( mccTEXT, 1, 5, S93__r_, 0); /* static array S94_ (1 x 4) text, line 164: 'Qpvo' */ static unsigned short S94__r_[] = { 81, 112, 118, 111, }; static mxArray S94_ = mccCINIT( mccTEXT, 1, 4, S94__r_, 0); /* static array S95_ (1 x 5) text, line 165: '%s %f' */ static unsigned short S95__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S95_ = mccCINIT( mccTEXT, 1, 5, S95__r_, 0); /* static array S96_ (1 x 2) text, line 167: 'Rl' */ static unsigned short S96__r_[] = { 82, 108, }; static mxArray S96_ = mccCINIT( mccTEXT, 1, 2, S96__r_, 0); /* static array S97_ (1 x 5) text, line 168: '%s %f' */ static unsigned short S97__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S97_ = mccCINIT( mccTEXT, 1, 5, S97__r_, 0); /* static array S98_ (1 x 2) text, line 170: 'Ra' */ static unsigned short S98__r_[] = { 82, 97, }; static mxArray S98_ = mccCINIT( mccTEXT, 1, 2, S98__r_, 0); /* static array S99_ (1 x 5) text, line 171: '%s %f' */ static unsigned short S99__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S99_ = mccCINIT( mccTEXT, 1, 5, S99__r_, 0); /* static array S100_ (1 x 2) text, line 173: 'Rv' */ static unsigned short S100__r_[] = { 82, 118, }; static mxArray S100_ = mccCINIT( mccTEXT, 1, 2, S100__r_, 0); /* static array S101_ (1 x 5) text, line 174: '%s %f' */ static unsigned short S101__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S101_ = mccCINIT( mccTEXT, 1, 5, S101__r_, 0); /* static array S102_ (1 x 2) text, line 176: 'Rr' */ static unsigned short S102__r_[] = { 82, 114, }; static mxArray S102_ = mccCINIT( mccTEXT, 1, 2, S102__r_, 0); /* static array S103_ (1 x 5) text, line 177: '%s %f' */ static unsigned short S103__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S103_ = mccCINIT( mccTEXT, 1, 5, S103__r_, 0); /* static array S104_ (1 x 3) text, line 179: 'Rpa' */ static unsigned short S104__r_[] = { 82, 112, 97, }; static mxArray S104_ = mccCINIT( mccTEXT, 1, 3, S104__r_, 0); /* static array S105_ (1 x 5) text, line 180: '%s %f' */ static unsigned short S105__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S105_ = mccCINIT( mccTEXT, 1, 5, S105__r_, 0); /* static array S106_ (1 x 3) text, line 182: 'Rpv' */ static unsigned short S106__r_[] = { 82, 112, 118, }; static mxArray S106_ = mccCINIT( mccTEXT, 1, 3, S106__r_, 0); /* static array S107_ (1 x 5) text, line 183: '%s %f' */ static unsigned short S107__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S107_ = mccCINIT( mccTEXT, 1, 5, S107__r_, 0); /* static array S108_ (1 x 4) text, line 185: 'Qtot' */ static unsigned short S108__r_[] = { 81, 116, 111, 116, }; static mxArray S108_ = mccCINIT( mccTEXT, 1, 4, S108__r_, 0); /* static array S109_ (1 x 5) text, line 186: '%s %f' */ static unsigned short S109__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S109_ = mccCINIT( mccTEXT, 1, 5, S109__r_, 0); /* static array S110_ (1 x 1) text, line 188: 'F' */ static unsigned short S110__r_[] = { 70, }; static mxArray S110_ = mccCINIT( mccTEXT, 1, 1, S110__r_, 0); /* static array S111_ (1 x 5) text, line 189: '%s %f' */ static unsigned short S111__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S111_ = mccCINIT( mccTEXT, 1, 5, S111__r_, 0); /* static array S112_ (1 x 5) text, line 197: 'Qlmax' */ static unsigned short S112__r_[] = { 81, 108, 109, 97, 120, }; static mxArray S112_ = mccCINIT( mccTEXT, 1, 5, S112__r_, 0); /* static array S113_ (1 x 5) text, line 198: '%s %f' */ static unsigned short S113__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S113_ = mccCINIT( mccTEXT, 1, 5, S113__r_, 0); /* static array S114_ (1 x 5) text, line 200: 'Qrmax' */ static unsigned short S114__r_[] = { 81, 114, 109, 97, 120, }; static mxArray S114_ = mccCINIT( mccTEXT, 1, 5, S114__r_, 0); /* static array S115_ (1 x 5) text, line 201: '%s %f' */ static unsigned short S115__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S115_ = mccCINIT( mccTEXT, 1, 5, S115__r_, 0); /* static array S116_ (1 x 3) text, line 203: 'Ppa' */ static unsigned short S116__r_[] = { 80, 112, 97, }; static mxArray S116_ = mccCINIT( mccTEXT, 1, 3, S116__r_, 0); /* static array S117_ (1 x 5) text, line 204: '%s %f' */ static unsigned short S117__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S117_ = mccCINIT( mccTEXT, 1, 5, S117__r_, 0); /* static array S118_ (1 x 2) text, line 206: 'Pa' */ static unsigned short S118__r_[] = { 80, 97, }; static mxArray S118_ = mccCINIT( mccTEXT, 1, 2, S118__r_, 0); /* static array S119_ (1 x 5) text, line 207: '%s %f' */ static unsigned short S119__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S119_ = mccCINIT( mccTEXT, 1, 5, S119__r_, 0); /* static array S120_ (1 x 2) text, line 209: 'Pv' */ static unsigned short S120__r_[] = { 80, 118, }; static mxArray S120_ = mccCINIT( mccTEXT, 1, 2, S120__r_, 0); /* static array S121_ (1 x 5) text, line 210: '%s %f' */ static unsigned short S121__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S121_ = mccCINIT( mccTEXT, 1, 5, S121__r_, 0); /* static array S122_ (1 x 5) text, line 212: 'Plmax' */ static unsigned short S122__r_[] = { 80, 108, 109, 97, 120, }; static mxArray S122_ = mccCINIT( mccTEXT, 1, 5, S122__r_, 0); /* static array S123_ (1 x 5) text, line 213: '%s %f' */ static unsigned short S123__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S123_ = mccCINIT( mccTEXT, 1, 5, S123__r_, 0); /* static array S124_ (1 x 5) text, line 215: 'Prmax' */ static unsigned short S124__r_[] = { 80, 114, 109, 97, 120, }; static mxArray S124_ = mccCINIT( mccTEXT, 1, 5, S124__r_, 0); /* static array S125_ (1 x 5) text, line 216: '%s %f' */ static unsigned short S125__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S125_ = mccCINIT( mccTEXT, 1, 5, S125__r_, 0); /* static array S126_ (1 x 3) text, line 218: 'Pms' */ static unsigned short S126__r_[] = { 80, 109, 115, }; static mxArray S126_ = mccCINIT( mccTEXT, 1, 3, S126__r_, 0); /* static array S127_ (1 x 5) text, line 219: '%s %f' */ static unsigned short S127__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S127_ = mccCINIT( mccTEXT, 1, 5, S127__r_, 0); /* static array S128_ (1 x 4) text, line 221: 'Qfrs' */ static unsigned short S128__r_[] = { 81, 102, 114, 115, }; static mxArray S128_ = mccCINIT( mccTEXT, 1, 4, S128__r_, 0); /* static array S129_ (1 x 5) text, line 222: '%s %f' */ static unsigned short S129__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S129_ = mccCINIT( mccTEXT, 1, 5, S129__r_, 0); /* static array S130_ (1 x 6) text, line 224: 'againr' */ static unsigned short S130__r_[] = { 97, 103, 97, 105, 110, 114, }; static mxArray S130_ = mccCINIT( mccTEXT, 1, 6, S130__r_, 0); /* static array S131_ (1 x 5) text, line 225: '%s %f' */ static unsigned short S131__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S131_ = mccCINIT( mccTEXT, 1, 5, S131__r_, 0); /* static array S132_ (1 x 5) text, line 227: 'bgain' */ static unsigned short S132__r_[] = { 98, 103, 97, 105, 110, }; static mxArray S132_ = mccCINIT( mccTEXT, 1, 5, S132__r_, 0); /* static array S133_ (1 x 5) text, line 228: '%s %f' */ static unsigned short S133__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S133_ = mccCINIT( mccTEXT, 1, 5, S133__r_, 0); /* static array S134_ (1 x 5) text, line 230: 'pgain' */ static unsigned short S134__r_[] = { 112, 103, 97, 105, 110, }; static mxArray S134_ = mccCINIT( mccTEXT, 1, 5, S134__r_, 0); /* static array S135_ (1 x 5) text, line 231: '%s %f' */ static unsigned short S135__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S135_ = mccCINIT( mccTEXT, 1, 5, S135__r_, 0); /* static array S136_ (1 x 4) text, line 233: 'Pasp' */ static unsigned short S136__r_[] = { 80, 97, 115, 112, }; static mxArray S136_ = mccCINIT( mccTEXT, 1, 4, S136__r_, 0); /* static array S137_ (1 x 5) text, line 234: '%s %f' */ static unsigned short S137__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S137_ = mccCINIT( mccTEXT, 1, 5, S137__r_, 0); /* static array S138_ (1 x 7) text, line 236: 'Pratrsp' */ static unsigned short S138__r_[] = { 80, 114, 97, 116, 114, 115, 112, }; static mxArray S138_ = mccCINIT( mccTEXT, 1, 7, S138__r_, 0); /* static array S139_ (1 x 5) text, line 237: '%s %f' */ static unsigned short S139__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S139_ = mccCINIT( mccTEXT, 1, 5, S139__r_, 0); /* static array S140_ (1 x 2) text, line 239: 'Tr' */ static unsigned short S140__r_[] = { 84, 114, }; static mxArray S140_ = mccCINIT( mccTEXT, 1, 2, S140__r_, 0); /* static array S141_ (1 x 5) text, line 240: '%s %f' */ static unsigned short S141__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S141_ = mccCINIT( mccTEXT, 1, 5, S141__r_, 0); /* static array S142_ (1 x 3) text, line 242: 'fco' */ static unsigned short S142__r_[] = { 102, 99, 111, }; static mxArray S142_ = mccCINIT( mccTEXT, 1, 3, S142__r_, 0); /* static array S143_ (1 x 5) text, line 243: '%s %f' */ static unsigned short S143__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S143_ = mccCINIT( mccTEXT, 1, 5, S143__r_, 0); /* static array S144_ (1 x 2) text, line 245: 'Sa' */ static unsigned short S144__r_[] = { 83, 97, }; static mxArray S144_ = mccCINIT( mccTEXT, 1, 2, S144__r_, 0); /* static array S145_ (1 x 5) text, line 246: '%s %f' */ static unsigned short S145__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S145_ = mccCINIT( mccTEXT, 1, 5, S145__r_, 0); /* static array S146_ (1 x 5) text, line 248: 'stdwr' */ static unsigned short S146__r_[] = { 115, 116, 100, 119, 114, }; static mxArray S146_ = mccCINIT( mccTEXT, 1, 5, S146__r_, 0); /* static array S147_ (1 x 5) text, line 249: '%s %f' */ static unsigned short S147__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S147_ = mccCINIT( mccTEXT, 1, 5, S147__r_, 0); /* static array S148_ (1 x 5) text, line 251: 'stdwf' */ static unsigned short S148__r_[] = { 115, 116, 100, 119, 102, }; static mxArray S148_ = mccCINIT( mccTEXT, 1, 5, S148__r_, 0); /* static array S149_ (1 x 5) text, line 252: '%s %f' */ static unsigned short S149__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S149_ = mccCINIT( mccTEXT, 1, 5, S149__r_, 0); /* static array S150_ (1 x 6) text, line 254: 'dgainp' */ static unsigned short S150__r_[] = { 100, 103, 97, 105, 110, 112, }; static mxArray S150_ = mccCINIT( mccTEXT, 1, 6, S150__r_, 0); /* static array S151_ (1 x 5) text, line 255: '%s %f' */ static unsigned short S151__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S151_ = mccCINIT( mccTEXT, 1, 5, S151__r_, 0); /* static array S152_ (1 x 6) text, line 257: 'dgains' */ static unsigned short S152__r_[] = { 100, 103, 97, 105, 110, 115, }; static mxArray S152_ = mccCINIT( mccTEXT, 1, 6, S152__r_, 0); /* static array S153_ (1 x 5) text, line 258: '%s %f' */ static unsigned short S153__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S153_ = mccCINIT( mccTEXT, 1, 5, S153__r_, 0); /* static array S154_ (1 x 6) text, line 260: 'againq' */ static unsigned short S154__r_[] = { 97, 103, 97, 105, 110, 113, }; static mxArray S154_ = mccCINIT( mccTEXT, 1, 6, S154__r_, 0); /* static array S155_ (1 x 5) text, line 261: '%s %f' */ static unsigned short S155__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S155_ = mccCINIT( mccTEXT, 1, 5, S155__r_, 0); /* static array S156_ (1 x 6) text, line 263: 'cgainr' */ static unsigned short S156__r_[] = { 99, 103, 97, 105, 110, 114, }; static mxArray S156_ = mccCINIT( mccTEXT, 1, 6, S156__r_, 0); /* static array S157_ (1 x 5) text, line 264: '%s %f' */ static unsigned short S157__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S157_ = mccCINIT( mccTEXT, 1, 5, S157__r_, 0); /* static array S158_ (1 x 6) text, line 266: 'cgainq' */ static unsigned short S158__r_[] = { 99, 103, 97, 105, 110, 113, }; static mxArray S158_ = mccCINIT( mccTEXT, 1, 6, S158__r_, 0); /* static array S159_ (1 x 5) text, line 267: '%s %f' */ static unsigned short S159__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S159_ = mccCINIT( mccTEXT, 1, 5, S159__r_, 0); /* static array S160_ (1 x 7) text, line 269: 'cgainfp' */ static unsigned short S160__r_[] = { 99, 103, 97, 105, 110, 102, 112, }; static mxArray S160_ = mccCINIT( mccTEXT, 1, 7, S160__r_, 0); /* static array S161_ (1 x 5) text, line 270: '%s %f' */ static unsigned short S161__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S161_ = mccCINIT( mccTEXT, 1, 5, S161__r_, 0); /* static array S162_ (1 x 2) text, line 272: 'Sc' */ static unsigned short S162__r_[] = { 83, 99, }; static mxArray S162_ = mccCINIT( mccTEXT, 1, 2, S162__r_, 0); /* static array S163_ (1 x 5) text, line 273: '%s %f' */ static unsigned short S163__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S163_ = mccCINIT( mccTEXT, 1, 5, S163__r_, 0); /* static array S164_ (1 x 2) text, line 275: 'Qt' */ static unsigned short S164__r_[] = { 81, 116, }; static mxArray S164_ = mccCINIT( mccTEXT, 1, 2, S164__r_, 0); /* static array S165_ (1 x 5) text, line 276: '%s %f' */ static unsigned short S165__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S165_ = mccCINIT( mccTEXT, 1, 5, S165__r_, 0); /* static array S166_ (1 x 3) text, line 278: 'Pvs' */ static unsigned short S166__r_[] = { 80, 118, 115, }; static mxArray S166_ = mccCINIT( mccTEXT, 1, 3, S166__r_, 0); /* static array S167_ (1 x 5) text, line 279: '%s %f' */ static unsigned short S167__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S167_ = mccCINIT( mccTEXT, 1, 5, S167__r_, 0); /* static array S168_ (1 x 4) text, line 281: 'Crds' */ static unsigned short S168__r_[] = { 67, 114, 100, 115, }; static mxArray S168_ = mccCINIT( mccTEXT, 1, 4, S168__r_, 0); /* static array S169_ (1 x 5) text, line 282: '%s %f' */ static unsigned short S169__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S169_ = mccCINIT( mccTEXT, 1, 5, S169__r_, 0); /* static array S170_ (1 x 3) text, line 284: 'Pas' */ static unsigned short S170__r_[] = { 80, 97, 115, }; static mxArray S170_ = mccCINIT( mccTEXT, 1, 3, S170__r_, 0); /* static array S171_ (1 x 5) text, line 285: '%s %f' */ static unsigned short S171__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S171_ = mccCINIT( mccTEXT, 1, 5, S171__r_, 0); /* static array S172_ (1 x 3) text, line 287: 'Pvc' */ static unsigned short S172__r_[] = { 80, 118, 99, }; static mxArray S172_ = mccCINIT( mccTEXT, 1, 3, S172__r_, 0); /* static array S173_ (1 x 5) text, line 288: '%s %f' */ static unsigned short S173__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S173_ = mccCINIT( mccTEXT, 1, 5, S173__r_, 0); /* static array S174_ (1 x 5) text, line 290: 'again' */ static unsigned short S174__r_[] = { 97, 103, 97, 105, 110, }; static mxArray S174_ = mccCINIT( mccTEXT, 1, 5, S174__r_, 0); /* static array S175_ (1 x 5) text, line 291: '%s %f' */ static unsigned short S175__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S175_ = mccCINIT( mccTEXT, 1, 5, S175__r_, 0); /* static array S176_ (1 x 6) text, line 293: 'againc' */ static unsigned short S176__r_[] = { 97, 103, 97, 105, 110, 99, }; static mxArray S176_ = mccCINIT( mccTEXT, 1, 6, S176__r_, 0); /* static array S177_ (1 x 5) text, line 294: '%s %f' */ static unsigned short S177__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S177_ = mccCINIT( mccTEXT, 1, 5, S177__r_, 0); /* static array S178_ (1 x 6) text, line 296: 'cgainc' */ static unsigned short S178__r_[] = { 99, 103, 97, 105, 110, 99, }; static mxArray S178_ = mccCINIT( mccTEXT, 1, 6, S178__r_, 0); /* static array S179_ (1 x 5) text, line 297: '%s %f' */ static unsigned short S179__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S179_ = mccCINIT( mccTEXT, 1, 5, S179__r_, 0); /* static array S180_ (1 x 7) text, line 299: 'againfs' */ static unsigned short S180__r_[] = { 97, 103, 97, 105, 110, 102, 115, }; static mxArray S180_ = mccCINIT( mccTEXT, 1, 7, S180__r_, 0); /* static array S181_ (1 x 5) text, line 300: '%s %f' */ static unsigned short S181__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S181_ = mccCINIT( mccTEXT, 1, 5, S181__r_, 0); /* static array S182_ (1 x 7) text, line 302: 'cgainfs' */ static unsigned short S182__r_[] = { 99, 103, 97, 105, 110, 102, 115, }; static mxArray S182_ = mccCINIT( mccTEXT, 1, 7, S182__r_, 0); /* static array S183_ (1 x 5) text, line 303: '%s %f' */ static unsigned short S183__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S183_ = mccCINIT( mccTEXT, 1, 5, S183__r_, 0); /* static array S184_ (1 x 7) text, line 305: 'againfp' */ static unsigned short S184__r_[] = { 97, 103, 97, 105, 110, 102, 112, }; static mxArray S184_ = mccCINIT( mccTEXT, 1, 7, S184__r_, 0); /* static array S185_ (1 x 5) text, line 306: '%s %f' */ static unsigned short S185__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S185_ = mccCINIT( mccTEXT, 1, 5, S185__r_, 0); /* static array S186_ (1 x 4) text, line 308: 'time' */ static unsigned short S186__r_[] = { 116, 105, 109, 101, }; static mxArray S186_ = mccCINIT( mccTEXT, 1, 4, S186__r_, 0); /* static array S187_ (1 x 5) text, line 309: '%s %f' */ static unsigned short S187__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S187_ = mccCINIT( mccTEXT, 1, 5, S187__r_, 0); /* static array S188_ (1 x 11) text, line 311: 'preparation' */ static unsigned short S188__r_[] = { 112, 114, 101, 112, 97, 114, 97, 116, 105, 111, 110, }; static mxArray S188_ = mccCINIT( mccTEXT, 1, 11, S188__r_, 0); /* static array S189_ (1 x 5) text, line 312: '%s %f' */ static unsigned short S189__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S189_ = mccCINIT( mccTEXT, 1, 5, S189__r_, 0); /* static array S190_ (1 x 9) text, line 314: 'breathing' */ static unsigned short S190__r_[] = { 98, 114, 101, 97, 116, 104, 105, 110, 103, }; static mxArray S190_ = mccCINIT( mccTEXT, 1, 9, S190__r_, 0); /* static array S191_ (1 x 5) text, line 315: '%s %f' */ static unsigned short S191__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S191_ = mccCINIT( mccTEXT, 1, 5, S191__r_, 0); /* static array S192_ (1 x 4) text, line 317: 'dncm' */ static unsigned short S192__r_[] = { 100, 110, 99, 109, }; static mxArray S192_ = mccCINIT( mccTEXT, 1, 4, S192__r_, 0); /* static array S193_ (1 x 5) text, line 318: '%s %f' */ static unsigned short S193__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S193_ = mccCINIT( mccTEXT, 1, 5, S193__r_, 0); /* static array S194_ (1 x 4) text, line 320: 'baro' */ static unsigned short S194__r_[] = { 98, 97, 114, 111, }; static mxArray S194_ = mccCINIT( mccTEXT, 1, 4, S194__r_, 0); /* static array S195_ (1 x 5) text, line 321: '%s %f' */ static unsigned short S195__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S195_ = mccCINIT( mccTEXT, 1, 5, S195__r_, 0); /* static array S196_ (1 x 3) text, line 323: 'dra' */ static unsigned short S196__r_[] = { 100, 114, 97, }; static mxArray S196_ = mccCINIT( mccTEXT, 1, 3, S196__r_, 0); /* static array S197_ (1 x 5) text, line 324: '%s %f' */ static unsigned short S197__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S197_ = mccCINIT( mccTEXT, 1, 5, S197__r_, 0); /* static array S198_ (1 x 2) text, line 326: 'df' */ static unsigned short S198__r_[] = { 100, 102, }; static mxArray S198_ = mccCINIT( mccTEXT, 1, 2, S198__r_, 0); /* static array S199_ (1 x 5) text, line 327: '%s %f' */ static unsigned short S199__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S199_ = mccCINIT( mccTEXT, 1, 5, S199__r_, 0); /* static array S200_ (1 x 2) text, line 329: 'Fs' */ static unsigned short S200__r_[] = { 70, 115, }; static mxArray S200_ = mccCINIT( mccTEXT, 1, 2, S200__r_, 0); /* static array S201_ (1 x 5) text, line 330: '%s %f' */ static unsigned short S201__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S201_ = mccCINIT( mccTEXT, 1, 5, S201__r_, 0); /* static array S202_ (1 x 3) text, line 332: 'Qfr' */ static unsigned short S202__r_[] = { 81, 102, 114, }; static mxArray S202_ = mccCINIT( mccTEXT, 1, 3, S202__r_, 0); /* static array S203_ (1 x 5) text, line 333: '%s %f' */ static unsigned short S203__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S203_ = mccCINIT( mccTEXT, 1, 5, S203__r_, 0); /* static array S204_ (1 x 3) text, line 335: 'Qds' */ static unsigned short S204__r_[] = { 81, 100, 115, }; static mxArray S204_ = mccCINIT( mccTEXT, 1, 3, S204__r_, 0); /* static array S205_ (1 x 5) text, line 336: '%s %f' */ static unsigned short S205__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S205_ = mccCINIT( mccTEXT, 1, 5, S205__r_, 0); /* static array S206_ (1 x 4) text, line 338: 'Rair' */ static unsigned short S206__r_[] = { 82, 97, 105, 114, }; static mxArray S206_ = mccCINIT( mccTEXT, 1, 4, S206__r_, 0); /* static array S207_ (1 x 5) text, line 339: '%s %f' */ static unsigned short S207__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S207_ = mccCINIT( mccTEXT, 1, 5, S207__r_, 0); /* static array S208_ (1 x 3) text, line 341: 'Clu' */ static unsigned short S208__r_[] = { 67, 108, 117, }; static mxArray S208_ = mccCINIT( mccTEXT, 1, 3, S208__r_, 0); /* static array S209_ (1 x 5) text, line 342: '%s %f' */ static unsigned short S209__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S209_ = mccCINIT( mccTEXT, 1, 5, S209__r_, 0); /* static array S210_ (1 x 4) text, line 344: 'Qluo' */ static unsigned short S210__r_[] = { 81, 108, 117, 111, }; static mxArray S210_ = mccCINIT( mccTEXT, 1, 4, S210__r_, 0); /* static array S211_ (1 x 5) text, line 345: '%s %f' */ static unsigned short S211__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S211_ = mccCINIT( mccTEXT, 1, 5, S211__r_, 0); /* static array S212_ (1 x 6) text, line 350: 'window' */ static unsigned short S212__r_[] = { 119, 105, 110, 100, 111, 119, }; static mxArray S212_ = mccCINIT( mccTEXT, 1, 6, S212__r_, 0); /* static array S213_ (1 x 5) text, line 351: '%s %f' */ static unsigned short S213__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S213_ = mccCINIT( mccTEXT, 1, 5, S213__r_, 0); /* static array S214_ (1 x 11) text, line 353: 'annotations' */ static unsigned short S214__r_[] = { 97, 110, 110, 111, 116, 97, 116, 105, 111, 110, 115, }; static mxArray S214_ = mccCINIT( mccTEXT, 1, 11, S214__r_, 0); /* static array S215_ (1 x 5) text, line 354: '%s %f' */ static unsigned short S215__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S215_ = mccCINIT( mccTEXT, 1, 5, S215__r_, 0); /* static array S216_ (1 x 8) text, line 356: 'numerics' */ static unsigned short S216__r_[] = { 110, 117, 109, 101, 114, 105, 99, 115, }; static mxArray S216_ = mccCINIT( mccTEXT, 1, 8, S216__r_, 0); /* static array S217_ (1 x 5) text, line 357: '%s %f' */ static unsigned short S217__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S217_ = mccCINIT( mccTEXT, 1, 5, S217__r_, 0); /* static array S218_ (1 x 8) text, line 359: 'waveform' */ static unsigned short S218__r_[] = { 119, 97, 118, 101, 102, 111, 114, 109, }; static mxArray S218_ = mccCINIT( mccTEXT, 1, 8, S218__r_, 0); /* static array S219_ (1 x 1) text, line 360: ':' */ static unsigned short S219__r_[] = { 58, }; static mxArray S219_ = mccCINIT( mccTEXT, 1, 1, S219__r_, 0); /* static array S220_ (1 x 5) text, line 362: 'alpha' */ static unsigned short S220__r_[] = { 97, 108, 112, 104, 97, }; static mxArray S220_ = mccCINIT( mccTEXT, 1, 5, S220__r_, 0); /* static array S221_ (1 x 5) text, line 363: '%s %f' */ static unsigned short S221__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S221_ = mccCINIT( mccTEXT, 1, 5, S221__r_, 0); /* static array S222_ (1 x 4) text, line 365: 'step' */ static unsigned short S222__r_[] = { 115, 116, 101, 112, }; static mxArray S222_ = mccCINIT( mccTEXT, 1, 4, S222__r_, 0); /* static array S223_ (1 x 5) text, line 366: '%s %f' */ static unsigned short S223__r_[] = { 37, 115, 32, 37, 102, }; static mxArray S223_ = mccCINIT( mccTEXT, 1, 5, S223__r_, 0); /* static array S224_ (1 x 3) int, line 54 */ static int S224r_[] = { 36, 42, 77, }; static mxArray S224_ = mccCINIT( mccINT, 1, 3, S224r_, 0 ); /* static array S225_ (1 x 3) int, line 54 */ static int S225r_[] = { 36, 42, 77, }; static mxArray S225_ = mccCINIT( mccINT, 1, 3, S225r_, 0 ); /* static array S226_ (1 x 3) int, line 58 */ static int S226r_[] = { 36, 42, 77, }; static mxArray S226_ = mccCINIT( mccINT, 1, 3, S226r_, 0 ); /* static array S227_ (1 x 3) int, line 58 */ static int S227r_[] = { 36, 42, 77, }; static mxArray S227_ = mccCINIT( mccINT, 1, 3, S227r_, 0 ); /* static array S228_ (1 x 2) int, line 99 */ static int S228r_[] = { 44, 47, }; static mxArray S228_ = mccCINIT( mccINT, 1, 2, S228r_, 0 ); /* static array S229_ (1 x 2) int, line 99 */ static int S229r_[] = { 44, 47, }; static mxArray S229_ = mccCINIT( mccINT, 1, 2, S229r_, 0 ); /***************** Compiler Assumptions **************** * M-File: /data/rmukkama/cvsim/rcvsim/src/simulate.m * * B0_ boolean scalar temporary * B1_ boolean scalar temporary * B2_ boolean scalar temporary * B3_ boolean scalar temporary * B4_ boolean scalar temporary * B5_ boolean scalar temporary * B6_ boolean scalar temporary * B7_ boolean scalar temporary * Clssp real scalar * Crd real scalar * Crdo real scalar * Crs real scalar * Crssp real scalar * El real scalar * Ela real scalar * Er real scalar * Era real scalar * Fsp real scalar * I0_ integer scalar temporary * I1_ integer scalar temporary * I2_ integer scalar temporary * I3_ integer scalar temporary * I4_ integer scalar temporary * IM0_ integer vector/matrix temporary * IM10_ integer vector/matrix temporary * IM11_ integer vector/matrix temporary * IM12_ integer vector/matrix temporary * IM13_ integer vector/matrix temporary * IM14_ integer vector/matrix temporary * IM15_ integer vector/matrix temporary * IM16_ integer vector/matrix temporary * IM1_ integer vector/matrix temporary * IM2_ integer vector/matrix temporary * IM3_ integer vector/matrix temporary * IM4_ integer vector/matrix temporary * IM5_ integer vector/matrix temporary * IM6_ integer vector/matrix temporary * IM7_ integer vector/matrix temporary * IM8_ integer vector/matrix temporary * IM9_ integer vector/matrix temporary * Nbytes integer scalar * P real vector/matrix * Pasp real scalar * Pratrsp real scalar * Ptemp real vector/matrix * Ptempn real vector/matrix * Q real vector/matrix * Qamax real scalar * Qao real scalar * Qtot real vector/matrix * Qvosp real scalar * R0_ real scalar temporary * R1_ real scalar temporary * R2_ real scalar temporary * R3_ real scalar temporary * R4_ real scalar temporary * R5_ real scalar temporary * RM0_ real vector/matrix temporary * RM10_ real vector/matrix temporary * RM11_ real vector/matrix temporary * RM12_ real vector/matrix temporary * RM13_ real vector/matrix temporary * RM14_ real vector/matrix temporary * RM15_ real vector/matrix temporary * RM16_ real vector/matrix temporary * RM1_ real vector/matrix temporary * RM2_ real vector/matrix temporary * RM3_ real vector/matrix temporary * RM4_ real vector/matrix temporary * RM5_ real vector/matrix temporary * RM6_ real vector/matrix temporary * RM7_ real vector/matrix temporary * RM8_ real vector/matrix temporary * RM9_ real vector/matrix temporary * Ralvr real vector/matrix * Ralvrp real vector/matrix * Ralvrw real vector/matrix * Rasp real scalar * S0_ * S1_ * S2_ * Sgran real scalar * Slength real scalar * Sratio integer scalar * TM0_ string vector/matrix temporary * TM1_ string vector/matrix temporary * TM2_ string vector/matrix temporary * TM3_ string vector/matrix temporary * TM4_ string vector/matrix temporary * Vvo real vector/matrix * Vvow real vector/matrix * X real vector/matrix * a_init_cond * abreflex * alpha real scalar * alphatau real scalar * amp real scalar * annotator string vector/matrix * ans_filt * ap real vector/matrix * apr_init_cond * at real vector/matrix * atan * atcount integer scalar * avebP real vector/matrix * avebPratr real vector/matrix * bP real vector/matrix * bPratr real vector/matrix * bPratrtemp real vector/matrix * bPtemp real vector/matrix * baro real scalar * beta real scalar * bl_filt * breathing real scalar * c real scalar * c1 real scalar * ceil * ceostep integer scalar * conserve_vol * conv * cos * count_pc integer scalar * countj integer scalar * countn integer scalar * cpreflex * d real scalar * d1 real scalar * dEl real scalar * dEla real scalar * dEr real scalar * dEra real scalar * deltaPsp real vector/matrix * deltaT real scalar * denb real scalar * denb1 real scalar * df real scalar * diff * dmax integer scalar * dmin integer scalar * dncm real scalar * dncm_filt * dpasp real scalar * dqvo real scalar * dra real scalar * dthbreathe real vector/matrix * end * end_time real scalar * error * exp * fb_int integer scalar * fclose * fid integer scalar * fid2 integer scalar * fidp integer scalar * final real scalar * flag real vector/matrix * flagnew real vector/matrix * floor * fopen * frac real scalar * freq real scalar * fwrite * h1 real scalar * h2 real scalar * hans real vector/matrix * hf real vector/matrix * hilvhr real vector/matrix * hlu_init_cond * honeoverf real vector/matrix * honeoverfra real vector/matrix * hra real vector/matrix * hvvo real vector/matrix * i integer scalar * i real vector/matrix => i_1 * i integer vector/matrix => i_2 * ilv real vector/matrix * ilv_dec * init_cond * intact_init_cond * ipfm real scalar * ipfmcond boolean scalar * k integer scalar * kscale integer scalar * kstart integer scalar * kt integer vector/matrix * kt integer scalar => kt_1 * length * lengthra integer scalar * lengthsa integer scalar * lengtht integer scalar * lv_init_cond * lvrstep integer scalar * mbintscalar * mco real vector/matrix * mrap real vector/matrix * n real scalar * n integer scalar => n_1 * nac_init_cond * nargin * num real vector/matrix * num2str * nvvostep integer scalar * oneoverf_filt * ones * outputfile real vector/matrix * outputfile2 string vector/matrix * param_change * parameter_change integer scalar * parameterfile real vector/matrix * pi * preparation real scalar * q real vector/matrix * qpv real scalar * qrs real vector/matrix * qrs_index integer scalar * qrs_index_start integer scalar * qrso real vector/matrix * qv real scalar * rand_int_breath * randn * rarange real vector/matrix * rarangei real vector/matrix * read_key * read_param * reallog * realsqrt * resp_act * rk4 * round * rsastep integer scalar * sc_init_cond * signals real vector/matrix * simulate * start_index real scalar * start_time real scalar * strcmp * suffix string vector/matrix * sum * sumP real vector/matrix * sumPratr real vector/matrix * sumdelta integer scalar * sumdeltaT real vector/matrix * summco real vector/matrix * summrap real vector/matrix * sumt real scalar * t real vector/matrix * tempPsp real vector/matrix * tempilv real vector/matrix * tempp real vector/matrix * tempv real vector/matrix * th real vector/matrix * tha real vector/matrix * thbreathe real vector/matrix * thbreathes real vector/matrix * thc real vector/matrix * third_init_cond * third_nac_init_cond * thn real vector/matrix * thold real vector/matrix * thp real vector/matrix * thsp real vector/matrix * thspold real vector/matrix * tics real scalar * tmp integer scalar * tn real scalar * tp real scalar * tr real vector/matrix * var_acap * var_cap * var_vcap * ve real vector/matrix * vent_vol * wave_remote * wf real vector/matrix * whitenoise real vector/matrix * wra real vector/matrix * wrflag integer scalar * write_param * x boolean scalar * x0 real scalar * xl real scalar * xr real scalar * xvol real scalar * y boolean scalar * yl real scalar * ypres real scalar * yr real scalar * z boolean scalar * zeros *******************************************************/ void mexFunction( int nlhs_, mxArray *plhs_[], int nrhs_, const mxArray *prhs_[] ) { mxArray *Mplhs_[1]; mxArray *Mprhs_[4]; if (nrhs_ > 5 ) { mexErrMsgTxt( "Too many input arguments." ); } if (nlhs_ > 8 ) { mexErrMsgTxt( "Too many output arguments." ); } mcmSetLineNumber(0); { mxArray P; mxArray Q; mxArray q; mxArray ap; mxArray ve; mxArray qrs; mxArray t; mxArray num; mxArray flag; mxArray th; mxArray outputfile; mxArray parameterfile; mxArray signals; double preparation = 0.0; double breathing = 0.0; double dncm = 0.0; double baro = 0.0; double dra = 0.0; double dpasp = 0.0; double df = 0.0; double dqvo = 0.0; double tics = 0.0; double deltaT = 0.0; double ipfm = 0.0; int qrs_index = 0; mxArray sumdeltaT; double Sgran = 0.0; int Sratio = 0; double Slength = 0.0; mxArray at; int i = 0; unsigned short x = 0; mxArray thbreathe; mxArray thbreathes; mxArray dthbreathe; mxArray hilvhr; mxArray ilv; int rsastep = 0; double amp = 0.0; double freq = 0.0; double n = 0.0; mxArray i_1; mxArray tr; mxArray Ralvrp; mxArray Ralvr; mxArray kt; int lvrstep = 0; mxArray hra; int n_1 = 0; mxArray whitenoise; int dmin = 0; int dmax = 0; mxArray honeoverfra; mxArray Ralvrw; mxArray wra; mxArray tempPsp; mxArray deltaPsp; int kt_1 = 0; int ceostep = 0; mxArray hans; mxArray honeoverf; mxArray hf; mxArray wf; mxArray hvvo; mxArray Vvow; mxArray Vvo; int nvvostep = 0; double Crd = 0.0; double final = 0.0; mxArray rarange; int lengthra = 0; int lengthsa = 0; int lengtht = 0; int countn = 0; int countj = 0; int tmp = 0; mxArray summrap; mxArray summco; int fb_int = 0; double Crs = 0.0; mxArray Ptemp; mxArray tempp; mxArray tempv; double Pasp = 0.0; double Pratrsp = 0.0; double Fsp = 0.0; double Rasp = 0.0; double Clssp = 0.0; double Crssp = 0.0; double Qvosp = 0.0; mxArray thsp; mxArray bP; mxArray bPratr; mxArray i_2; mxArray bPtemp; mxArray bPratrtemp; mxArray sumP; mxArray sumPratr; double sumt = 0.0; mxArray thn; mxArray thp; double tn = 0.0; double tp = 0.0; int kstart = 0; int qrs_index_start = 0; int fid = 0; int fid2 = 0; int Nbytes = 0; mxArray annotator; int count_pc = 0; int wrflag = 0; mxArray outputfile2; int k = 0; mxArray Ptempn; double frac = 0.0; mxArray thc; mxArray avebP; mxArray avebPratr; mxArray tempilv; mxArray tha; double qpv = 0.0; unsigned short y = 0; double qv = 0.0; int sumdelta = 0; unsigned short ipfmcond = 0; mxArray mrap; mxArray mco; mxArray rarangei; double Crdo = 0.0; double ypres = 0.0; double xvol = 0.0; double alpha = 0.0; double x0 = 0.0; int kscale = 0; double c1 = 0.0; double d1 = 0.0; double denb1 = 0.0; double beta = 0.0; double c = 0.0; double d = 0.0; double denb = 0.0; double El = 0.0; double dEl = 0.0; double Er = 0.0; double dEr = 0.0; double Ela = 0.0; double dEla = 0.0; double Era = 0.0; double dEra = 0.0; double alphatau = 0.0; double Qamax = 0.0; double Qao = 0.0; double yl = 0.0; double xl = 0.0; double yr = 0.0; double xr = 0.0; mxArray Qtot; mxArray flagnew; mxArray thspold; mxArray thold; int fidp = 0; int parameter_change = 0; int atcount = 0; double start_index = 0.0; double start_time = 0.0; double end_time = 0.0; double h1 = 0.0; double h2 = 0.0; unsigned short z = 0; mxArray X; mxArray qrso; mxArray suffix; int I0_ = 0; unsigned short B0_ = 0; double R0_ = 0.0; mxArray RM0_; unsigned short B1_ = 0; unsigned short B2_ = 0; mxArray IM0_; double R1_ = 0.0; int I1_ = 0; mxArray IM1_; mxArray RM1_; mxArray RM2_; int I2_ = 0; mxArray IM2_; mxArray RM3_; double R2_ = 0.0; double R3_ = 0.0; double R4_ = 0.0; double R5_ = 0.0; mxArray RM4_; mxArray RM5_; mxArray TM0_; mxArray TM1_; mxArray TM2_; mxArray TM3_; mxArray TM4_; unsigned short B3_ = 0; unsigned short B4_ = 0; unsigned short B5_ = 0; unsigned short B6_ = 0; mxArray IM3_; int I3_ = 0; int I4_ = 0; mxArray IM4_; mxArray IM5_; mxArray IM6_; mxArray IM7_; mxArray IM8_; mxArray IM9_; mxArray IM10_; mxArray IM11_; mxArray IM12_; mxArray IM13_; mxArray IM14_; mxArray IM15_; mxArray IM16_; unsigned short B7_ = 0; mxArray RM6_; mxArray RM7_; mxArray RM8_; mxArray RM9_; mxArray RM10_; mxArray RM11_; mxArray RM12_; mxArray RM13_; mxArray RM14_; mxArray RM15_; mxArray RM16_; mccRealInit(flag); mccImport(&flag, ((nrhs_>0) ? prhs_[0] : 0), 0, 0); mccRealInit(th); mccImportCopy(&th, ((nrhs_>1) ? prhs_[1] : 0), 0, 0); mccRealInit(outputfile); mccImport(&outputfile, ((nrhs_>2) ? prhs_[2] : 0), 0, 0); mccRealInit(parameterfile); mccImport(¶meterfile, ((nrhs_>3) ? prhs_[3] : 0), 0, 0); mccRealInit(signals); mccImportCopy(&signals, ((nrhs_>4) ? prhs_[4] : 0), 0, 0); mccRealInit(P); mccRealInit(Q); mccRealInit(q); mccRealInit(ap); mccRealInit(ve); mccRealInit(qrs); mccRealInit(t); mccRealInit(num); mccRealInit(sumdeltaT); mccRealInit(at); mccRealInit(thbreathe); mccRealInit(thbreathes); mccRealInit(dthbreathe); mccRealInit(hilvhr); mccRealInit(ilv); mccRealInit(i_1); mccRealInit(tr); mccRealInit(Ralvrp); mccRealInit(Ralvr); mccIntInit(kt); mccRealInit(hra); mccRealInit(whitenoise); mccRealInit(honeoverfra); mccRealInit(Ralvrw); mccRealInit(wra); mccRealInit(tempPsp); mccRealInit(deltaPsp); mccRealInit(hans); mccRealInit(honeoverf); mccRealInit(hf); mccRealInit(wf); mccRealInit(hvvo); mccRealInit(Vvow); mccRealInit(Vvo); mccRealInit(rarange); mccRealInit(summrap); mccRealInit(summco); mccRealInit(Ptemp); mccRealInit(tempp); mccRealInit(tempv); mccRealInit(thsp); mccRealInit(bP); mccRealInit(bPratr); mccIntInit(i_2); mccRealInit(bPtemp); mccRealInit(bPratrtemp); mccRealInit(sumP); mccRealInit(sumPratr); mccRealInit(thn); mccRealInit(thp); mccTextInit(annotator); mccTextInit(outputfile2); mccRealInit(Ptempn); mccRealInit(thc); mccRealInit(avebP); mccRealInit(avebPratr); mccRealInit(tempilv); mccRealInit(tha); mccRealInit(mrap); mccRealInit(mco); mccRealInit(rarangei); mccRealInit(Qtot); mccRealInit(flagnew); mccRealInit(thspold); mccRealInit(thold); mccRealInit(X); mccRealInit(qrso); mccTextInit(suffix); mccRealInit(RM0_); mccIntInit(IM0_); mccIntInit(IM1_); mccRealInit(RM1_); mccRealInit(RM2_); mccIntInit(IM2_); mccRealInit(RM3_); mccRealInit(RM4_); mccRealInit(RM5_); mccTextInit(TM0_); mccTextInit(TM1_); mccTextInit(TM2_); mccTextInit(TM3_); mccTextInit(TM4_); mccIntInit(IM3_); mccIntInit(IM4_); mccIntInit(IM5_); mccIntInit(IM6_); mccIntInit(IM7_); mccIntInit(IM8_); mccIntInit(IM9_); mccIntInit(IM10_); mccIntInit(IM11_); mccIntInit(IM12_); mccIntInit(IM13_); mccIntInit(IM14_); mccIntInit(IM15_); mccIntInit(IM16_); mccRealInit(RM6_); mccRealInit(RM7_); mccRealInit(RM8_); mccRealInit(RM9_); mccRealInit(RM10_); mccRealInit(RM11_); mccRealInit(RM12_); mccRealInit(RM13_); mccRealInit(RM14_); mccRealInit(RM15_); mccRealInit(RM16_); /* % Declaring C functions. */ /* % #function wave_remote read_key write_param */ /* % Checking number of arguments. If there are only two arguments, */ /* % this indicates MATLAB execution and thus on-line viewing is */ /* % turned off. (If it were left on, MATLAB would produce an error.) */ /* if (nargin == 2) */ I0_ = mccNargin(); B0_ = (I0_ == 2); if ((double)B0_) { /* signals = '-1'; */ mccCreateString(&signals, "-1"); /* end */ } /* % Checking dimensionality of the argument flag -- status parameters. */ /* % If flag is of incorrect dimension, the program will exit in error. */ /* if (length(flag) ~= 8) */ if(mccNOTSET(&flag)) { mexErrMsgTxt( "variable flag undefined, line 129" ); } I0_ = mccGetLength(&flag); B0_ = (I0_ != 8); if ((double)B0_) { /* error('flag is not of correct dimension'); */ mccError(&S3_); /* end */ } /* % Assigning status parameters to new variable names. */ /* preparation = flag(1); */ preparation = (mccPR(&flag)[(1-1)]); /* breathing = flag(2); */ breathing = (mccPR(&flag)[(2-1)]); /* dncm = flag(3); */ dncm = (mccPR(&flag)[(3-1)]); /* baro = flag(4); */ baro = (mccPR(&flag)[(4-1)]); /* dra = flag(5); */ dra = (mccPR(&flag)[(5-1)]); /* dpasp = flag(6); */ dpasp = (mccPR(&flag)[(6-1)]); /* df = flag(7); */ df = (mccPR(&flag)[(7-1)]); /* dqvo = flag(8); */ dqvo = (mccPR(&flag)[(8-1)]); /* % Assigning numerical integration variables. */ /* % Number of integration steps. */ /* tics = th(106); */ if(mccNOTSET(&th)) { mexErrMsgTxt( "variable th undefined, line 147" ); } tics = (mccPR(&th)[(106-1)]); /* % Integration step size (s). */ /* deltaT = 1/th(113); */ deltaT = (1 / (double) (mccPR(&th)[(113-1)])); /* % Assigning integrate and fire variables. */ /* % Fraction of current cardiac cycle that integration step is at. */ /* ipfm = 0; */ ipfm = 0; /* % Current number of ventricular contractions. */ /* qrs_index = 1; */ qrs_index = 1; /* % Time in current cardiac cycle that integration is at. */ /* sumdeltaT = 0; */ { double tr_ = 0; mccAllocateMatrix(&sumdeltaT, 1, 1); *mccPR(&sumdeltaT) = tr_; } /* % Short-term regulatory systems and resting physiologic */ /* % perturbation variables. */ /* % Note there are restrictions on altering Sgran and Sratio. */ /* % 1) Sgran*Sratio <= 1.0. */ /* % 2) (1.5/Sgran) must be an integer. */ /* % 3) ((2/Sgran)-((Sratio*Sgran/2)/Sgran)) must be an integer. */ /* % 4) (0.5-(Sratio*Sgran/2))/Sgran must be an integer. */ /* % Granularity (s). */ /* Sgran = 0.0625; */ Sgran = 0.0625; /* % Ratio of length of averaging window to Sgran */ /* % (must be an integer). */ /* Sratio = 4; */ Sratio = 4; /* % Duration (samples). */ /* Slength = 40/Sgran; */ Slength = (40 / (double) Sgran); /* % Generation of a vector of arrival times for each respiratory cycle. */ /* % at(1) = 0 & length(at) ~= 1: Qlu periodic */ /* if (breathing == 1) */ if (( (breathing == 1) && !mccREL_NAN(breathing) )) { /* at = zeros(ceil((deltaT*tics)+1)+1,1); */ R0_ = ceil(((deltaT * (double) tics) + 1)); mccZerosMN(&at, (((int)R0_)+1), 1); /* i = 2; */ i = 2; /* x = (at(i) <= deltaT*tics+1); */ x = ( (((int)mccPR(&at)[(i-1)]) <= ((deltaT * (double) tics) + 1)) && !mccREL_NAN(((deltaT * (double) tics) + 1)) ); /* mbintscalar(x); */ /* while (x) */ while (1) { if (!((double)x)) break; /* i = i+1; */ i = (i + 1); /* at(i) = at(i-1)+th(42); */ mccPR(&at)[(i-1)] = ((mccPR(&at)[((i-1)-1)]) + (mccPR(&th)[(42-1)])); /* x = (at(i) <= deltaT*tics+1); */ x = ( ((mccPR(&at)[(i-1)]) <= ((deltaT * (double) tics) + 1)) && !mccREL_NAN((mccPR(&at)[(i-1)])) && !mccREL_NAN(((deltaT * (double) tics) + 1)) ); /* mbintscalar(x); */ /* end */ } /* at = at(1:i); */ { int i_, j_; int m_=1, n_=1, cx_ = 0; double *p_RM0_; int I_RM0_=1; double *p_at; int I_at=1; m_ = mccCalcSubscriptDimensions(m_, &n_, 1, ((int)(i - 1) + 1), &at); mccAllocateMatrix(&RM0_, m_, n_); I_RM0_ = (mccM(&RM0_) != 1 || mccN(&RM0_) != 1); p_RM0_ = mccPR(&RM0_); I_at = (mccM(&at) != 1 || mccN(&at) != 1); p_at = mccPR(&at) + ((int)(1 - .5)); if (m_ != 0) { for (j_=0; j_ 1 & lengthsa > 1) */ B2_ = (lengthra > 1); if ((double)B2_) { B2_ = (lengthsa > 1); } if ((double)B2_) { /* countn = 0; */ countn = 0; /* countj = 0; */ countj = 0; /* th(30) = final; */ mccPR(&th)[(30-1)] = final; /* th(29) = 30; */ mccPR(&th)[(29-1)] = 30; /* tmp = 0; */ tmp = 0; /* % Initializing variables for a single cardiac function curve. */ } else { /* elseif (lengthra > 1 & lengthsa == 1) */ B2_ = (lengthra > 1); if ((double)B2_) { B2_ = (lengthsa == 1); } if ((double)B2_) { /* countn = 0; */ countn = 0; /* countj = 0; */ countj = 0; /* th(30) = final; */ mccPR(&th)[(30-1)] = final; /* tmp = 0; */ tmp = 0; /* % Initializing variables for curve of average ql versus Pa. */ } else { /* elseif (lengthra == 1 & lengthsa > 1) */ B2_ = (lengthra == 1); if ((double)B2_) { B2_ = (lengthsa > 1); } if ((double)B2_) { /* final = th(30); */ final = (mccPR(&th)[(30-1)]); /* rarange = final; */ { double tr_ = final; mccAllocateMatrix(&rarange, 1, 1); *mccPR(&rarange) = tr_; } /* th(29) = 30; */ mccPR(&th)[(29-1)] = 30; /* tmp = 0; */ tmp = 0; /* end */ } } } /* % Initializing variables for averaging Pra and ql over each cardiac cycle. */ /* summrap = 0; */ { double tr_ = 0; mccAllocateMatrix(&summrap, 1, 1); *mccPR(&summrap) = tr_; } /* summco = 0; */ { double tr_ = 0; mccAllocateMatrix(&summco, 1, 1); *mccPR(&summco) = tr_; } /* % Initializing variable representing number of beats until voltage */ /* % source(s) is varied. */ /* fb_int = 10; */ fb_int = 10; } else { /* elseif (preparation == 2) */ if (( (preparation == 2) && !mccREL_NAN(preparation) )) { /* % Pre-allocating memory. */ /* Crs = ((th(27)-th(12))/th(32))*th(5); */ Crs = ((((mccPR(&th)[(27-1)]) - (mccPR(&th)[(12-1)])) / (double) (mccPR(&th)[(32-1)])) * (double) (mccPR(&th)[(5-1)])); /* lengthra = length([Crs:th(89):60]); */ mccColon(&RM3_, Crs, (mccPR(&th)[(89-1)]), (double)60); lengthra = mccGetLength(&RM3_); /* num = zeros(2,lengthra); */ mccZerosMN(&num, 2, lengthra); mccLOG(&num) = 0; mccSTRING(&num) = 0; /* if (length(num(1,:)) > 1) */ { int i_, j_; int m_=1, n_=1, cx_ = 0; int *p_IM1_; int I_IM1_=1; double *p_num; int I_num=1, J_num; m_ = mcmCalcResultSize(m_, &n_, 1, mccN(&num)); mccAllocateMatrix(&IM1_, m_, n_); I_IM1_ = (mccM(&IM1_) != 1 || mccN(&IM1_) != 1); p_IM1_ = mccIPR(&IM1_); if (mccN(&num) == 1) { I_num = J_num = 0;} else { I_num = 1; J_num=mccM(&num)-m_; } p_num = mccPR(&num) + (1-1) + mccM(&num) * 0; if (m_ != 0) { for (j_=0; j_ 1); if ((double)B2_) { /* countn = 0; */ countn = 0; /* Crd = Crs; */ Crd = Crs; /* th(6) = (th(32)/(th(27)-th(12)))*Crd; */ mccPR(&th)[(6-1)] = (((mccPR(&th)[(32-1)]) / (double) ((mccPR(&th)[(27-1)]) - (mccPR(&th)[(12-1)]))) * (double) Crd); /* end */ } /* % Initializing variables for integrating Pra and qv over each cardiac cycle. */ /* summrap = 0; */ { double tr_ = 0; mccAllocateMatrix(&summrap, 1, 1); *mccPR(&summrap) = tr_; } /* summco = 0; */ { double tr_ = 0; mccAllocateMatrix(&summco, 1, 1); *mccPR(&summco) = tr_; } /* % Initializing variable representing number of beats until Crd */ /* % is varied. */ /* fb_int = 10; */ fb_int = 10; } else { /* else */ /* num = zeros(2,1); */ mccZerosMN(&num, 2, 1); mccLOG(&num) = 0; mccSTRING(&num) = 0; /* end */ } } /* % Resetting ventricular compliance values for the intact circulatory */ /* % preparation with only linear elements. */ /* if (preparation == 5) */ if (( (preparation == 5) && !mccREL_NAN(preparation) )) { /* th(1) = th(1)*((th(26)-th(9))/th(31)); */ mccPR(&th)[(1-1)] = ((mccPR(&th)[(1-1)]) * (double) (((mccPR(&th)[(26-1)]) - (mccPR(&th)[(9-1)])) / (double) (mccPR(&th)[(31-1)]))); /* th(2) = th(2)*((th(26)-th(9))/th(31)); */ mccPR(&th)[(2-1)] = ((mccPR(&th)[(2-1)]) * (double) (((mccPR(&th)[(26-1)]) - (mccPR(&th)[(9-1)])) / (double) (mccPR(&th)[(31-1)]))); /* th(5) = th(5)*((th(27)-th(12))/th(32)); */ mccPR(&th)[(5-1)] = ((mccPR(&th)[(5-1)]) * (double) (((mccPR(&th)[(27-1)]) - (mccPR(&th)[(12-1)])) / (double) (mccPR(&th)[(32-1)]))); /* th(6) = th(6)*((th(27)-th(12))/th(32)); */ mccPR(&th)[(6-1)] = ((mccPR(&th)[(6-1)]) * (double) (((mccPR(&th)[(27-1)]) - (mccPR(&th)[(12-1)])) / (double) (mccPR(&th)[(32-1)]))); /* th(31) = th(26)-th(9); */ mccPR(&th)[(31-1)] = ((mccPR(&th)[(26-1)]) - (mccPR(&th)[(9-1)])); /* th(32) = th(27)-th(12); */ mccPR(&th)[(32-1)] = ((mccPR(&th)[(27-1)]) - (mccPR(&th)[(12-1)])); /* end */ } /* % Initializing the function output waveforms. */ /* % Adjustable parameters. */ /* ap(1:6,1) = [th(1)*((th(26)-th(9))/th(31)) th(5)*((th(27)-th(12))/th(32)) th(11) th(16) th(22) .3*sqrt(th(24))]'; */ mccCatenateColumns(&RM3_, mccTempMatrix(((mccPR(&th)[(1-1)]) * (double) (((mccPR(&th)[(26-1)]) - (mccPR(&th)[(9-1)])) / (double) (mccPR(&th)[(31-1)]))), 0., mccREAL, 0 ), mccTempMatrix(((mccPR(&th)[(5-1)]) * (double) (((mccPR(&th)[(27-1)]) - (mccPR(&th)[(12-1)])) / (double) (mccPR(&th)[(32-1)]))), 0., mccREAL, 0 )); mccCatenateColumns(&RM0_, &RM3_, mccTempVectorElement(&th, 11)); mccCatenateColumns(&RM2_, &RM0_, mccTempVectorElement(&th, 16)); mccCatenateColumns(&RM1_, &RM2_, mccTempVectorElement(&th, 22)); R1_ = sqrt((mccPR(&th)[(24-1)])); mccCatenateColumns(&RM4_, &RM1_, mccTempMatrix((.3 * (double) R1_), 0., mccREAL, 0 )); mccConjTrans(&RM5_, &RM4_); { int i_, j_; int m_=1, n_=1, cx_ = 0; double *p_ap; int I_ap=1, J_ap; double *p_RM5_; int I_RM5_=1; m_ = mcmCalcResultSize(m_, &n_, ((int)(6 - 1) + 1), 1); m_ = mcmCalcResultSize(m_, &n_, mccM(&RM5_), mccN(&RM5_)); if (mccM(&ap) == 1) { I_ap = J_ap = 0;} else { I_ap = 1; J_ap=mccM(&ap)-m_; } p_ap = mccPR(&ap) + ((int)(1 - .5)) + mccM(&ap) * (1-1); I_RM5_ = (mccM(&RM5_) != 1 || mccN(&RM5_) != 1); p_RM5_ = mccPR(&RM5_); if (m_ != 0) { for (j_=0; j_