Generated Code

The following is c_ida code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
/*
   There are a total of 30 entries in the algebraic variable array.
   There are a total of 8 entries in each of the rate and state variable arrays.
   There are a total of 81 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (minute).
 * STATES[0] is PYR in component PYR (millimolar).
 * ALGEBRAIC[17] is V_GLYC in component V_GLYC (flux).
 * ALGEBRAIC[18] is V_LDH in component V_LDH (flux).
 * ALGEBRAIC[19] is V_PDH in component V_PDH (flux).
 * ALGEBRAIC[24] is V_ALS in component V_ALS (flux).
 * STATES[1] is ACP in component ACP (millimolar).
 * ALGEBRAIC[20] is V_PTA in component V_PTA (flux).
 * ALGEBRAIC[21] is V_ACK in component V_ACK (flux).
 * STATES[2] is ACAL in component ACAL (millimolar).
 * ALGEBRAIC[22] is V_ACALDH in component V_ACALDH (flux).
 * ALGEBRAIC[23] is V_ADH in component V_ADH (flux).
 * STATES[3] is ACLAC in component ACLAC (millimolar).
 * ALGEBRAIC[25] is V_ALDC in component V_ALDC (flux).
 * ALGEBRAIC[16] is V_NEALC in component V_NEALC (flux).
 * STATES[4] is ACET in component ACET (millimolar).
 * ALGEBRAIC[27] is V_ACETDH in component V_ACETDH (flux).
 * ALGEBRAIC[26] is V_ACETEFF in component V_ACETEFF (flux).
 * STATES[5] is ATP in component ATP (millimolar).
 * ALGEBRAIC[28] is V_ATPase in component V_ATPase (flux).
 * ALGEBRAIC[0] is ADP in component ADP (millimolar).
 * CONSTANTS[0] is A_tot in component ADP (millimolar).
 * STATES[6] is NADH in component NADH (millimolar).
 * ALGEBRAIC[29] is V_NOX in component V_NOX (flux).
 * ALGEBRAIC[1] is NAD in component NAD (millimolar).
 * CONSTANTS[1] is NAD_tot in component NAD (millimolar).
 * STATES[7] is ACCOA in component ACCOA (millimolar).
 * ALGEBRAIC[2] is COA in component COA (millimolar).
 * CONSTANTS[2] is C_tot in component COA (millimolar).
 * CONSTANTS[3] is AC in component AC (millimolar).
 * CONSTANTS[4] is BUT in component BUT (millimolar).
 * CONSTANTS[5] is ETOH in component ETOH (millimolar).
 * CONSTANTS[6] is GLC in component GLC (millimolar).
 * CONSTANTS[7] is LAC in component LAC (millimolar).
 * CONSTANTS[8] is O in component O (millimolar).
 * CONSTANTS[9] is P in component P (millimolar).
 * ALGEBRAIC[3] is V_GLYC_temp in component V_GLYC (flux).
 * CONSTANTS[10] is Km_GLC in component V_GLYC (millimolar).
 * CONSTANTS[11] is Km_NAD in component V_GLYC (millimolar).
 * CONSTANTS[12] is Km_ADP in component V_GLYC (millimolar).
 * CONSTANTS[13] is Km_PYR in component V_GLYC (millimolar).
 * CONSTANTS[14] is Km_NADH in component V_GLYC (millimolar).
 * CONSTANTS[15] is Km_ATP in component V_GLYC (millimolar).
 * CONSTANTS[16] is V_GLYC_max in component V_GLYC (flux).
 * ALGEBRAIC[4] is V_LDH_temp in component V_LDH (flux).
 * CONSTANTS[17] is Keq in component V_LDH (dimensionless).
 * CONSTANTS[18] is Km_LAC in component V_LDH (millimolar).
 * CONSTANTS[19] is Km_NAD in component V_LDH (millimolar).
 * CONSTANTS[20] is Km_PYR in component V_LDH (millimolar).
 * CONSTANTS[21] is Km_NADH in component V_LDH (millimolar).
 * CONSTANTS[22] is V_LDH_max in component V_LDH (flux).
 * ALGEBRAIC[5] is V_PDH_temp in component V_PDH (flux).
 * CONSTANTS[23] is Ki in component V_PDH (dimensionless).
 * CONSTANTS[24] is Km_NAD in component V_PDH (millimolar).
 * CONSTANTS[25] is Km_COA in component V_PDH (millimolar).
 * CONSTANTS[26] is Km_PYR in component V_PDH (millimolar).
 * CONSTANTS[27] is Km_NADH in component V_PDH (millimolar).
 * CONSTANTS[28] is Km_ACCOA in component V_PDH (millimolar).
 * CONSTANTS[29] is V_PDH_max in component V_PDH (flux).
 * ALGEBRAIC[6] is V_PTA_temp in component V_PTA (flux).
 * CONSTANTS[30] is Keq in component V_PTA (dimensionless).
 * CONSTANTS[31] is Km_P in component V_PTA (millimolar).
 * CONSTANTS[32] is Ki_P in component V_PTA (millimolar).
 * CONSTANTS[33] is Ki_COA in component V_PTA (millimolar).
 * CONSTANTS[34] is Km_ACP in component V_PTA (millimolar).
 * CONSTANTS[35] is Ki_ACP in component V_PTA (millimolar).
 * CONSTANTS[36] is Ki_ACCOA in component V_PTA (millimolar).
 * CONSTANTS[37] is V_PTA_max in component V_PTA (flux).
 * ALGEBRAIC[7] is V_ACK_temp in component V_ACK (flux).
 * CONSTANTS[38] is Keq in component V_ACK (dimensionless).
 * CONSTANTS[39] is Km_AC in component V_ACK (millimolar).
 * CONSTANTS[40] is Km_ATP in component V_ACK (millimolar).
 * CONSTANTS[41] is Km_ADP in component V_ACK (millimolar).
 * CONSTANTS[42] is Km_ACP in component V_ACK (millimolar).
 * CONSTANTS[43] is V_ACK_max in component V_ACK (flux).
 * ALGEBRAIC[8] is V_ACALDH_temp in component V_ACALDH (flux).
 * CONSTANTS[44] is Keq in component V_ACALDH (millimolar).
 * CONSTANTS[45] is Km_NAD in component V_ACALDH (millimolar).
 * CONSTANTS[46] is Km_NADH in component V_ACALDH (millimolar).
 * CONSTANTS[47] is Km_COA in component V_ACALDH (millimolar).
 * CONSTANTS[48] is Km_ACCOA in component V_ACALDH (millimolar).
 * CONSTANTS[49] is Km_ACAL in component V_ACALDH (millimolar).
 * CONSTANTS[50] is V_ACALDH_max in component V_ACALDH (flux).
 * ALGEBRAIC[9] is V_ADH_temp in component V_ADH (flux).
 * CONSTANTS[51] is Keq in component V_ADH (dimensionless).
 * CONSTANTS[52] is Km_NAD in component V_ADH (millimolar).
 * CONSTANTS[53] is Km_NADH in component V_ADH (millimolar).
 * CONSTANTS[54] is Km_ETOH in component V_ADH (millimolar).
 * CONSTANTS[55] is Km_ACAL in component V_ADH (millimolar).
 * CONSTANTS[56] is V_ADH_max in component V_ADH (flux).
 * ALGEBRAIC[10] is V_ALS_temp in component V_ALS (flux).
 * CONSTANTS[57] is N in component V_ALS (dimensionless).
 * CONSTANTS[58] is Keq in component V_ALS (dimensionless).
 * CONSTANTS[59] is Km_ACLAC in component V_ALS (millimolar).
 * CONSTANTS[60] is Km_PYR in component V_ALS (millimolar).
 * CONSTANTS[61] is V_ALS_max in component V_ALS (flux).
 * ALGEBRAIC[11] is V_ALDC_temp in component V_ALDC (flux).
 * CONSTANTS[62] is Km_ACLAC in component V_ALDC (millimolar).
 * CONSTANTS[63] is Km_ACET in component V_ALDC (millimolar).
 * CONSTANTS[64] is V_ALDC_max in component V_ALDC (flux).
 * ALGEBRAIC[12] is V_ACETEFF_temp in component V_ACETEFF (flux).
 * CONSTANTS[65] is Km_ACET in component V_ACETEFF (millimolar).
 * CONSTANTS[66] is V_ACETEFF_max in component V_ACETEFF (flux).
 * ALGEBRAIC[13] is V_ACETDH_temp in component V_ACETDH (flux).
 * CONSTANTS[67] is Keq in component V_ACETDH (dimensionless).
 * CONSTANTS[68] is Km_NAD in component V_ACETDH (millimolar).
 * CONSTANTS[69] is Km_NADH in component V_ACETDH (millimolar).
 * CONSTANTS[70] is Km_BUT in component V_ACETDH (millimolar).
 * CONSTANTS[71] is Km_ACET in component V_ACETDH (millimolar).
 * CONSTANTS[72] is V_ACETDH_max in component V_ACETDH (flux).
 * ALGEBRAIC[14] is V_ATPase_temp in component V_ATPase (flux).
 * CONSTANTS[73] is N in component V_ATPase (dimensionless).
 * CONSTANTS[74] is Km_ATP in component V_ATPase (dimensionless).
 * CONSTANTS[75] is V_ATPase_max in component V_ATPase (flux).
 * ALGEBRAIC[15] is V_NOX_temp in component V_NOX (flux).
 * CONSTANTS[76] is Km_NAD in component V_NOX (millimolar).
 * CONSTANTS[77] is Km_NADH in component V_NOX (millimolar).
 * CONSTANTS[78] is Km_O in component V_NOX (millimolar).
 * CONSTANTS[79] is V_NOX_max in component V_NOX (flux).
 * CONSTANTS[80] is k in component V_NEALC (first_order_rate_constant).
 * RATES[0] is d/dt PYR in component PYR (millimolar).
 * RATES[1] is d/dt ACP in component ACP (millimolar).
 * RATES[2] is d/dt ACAL in component ACAL (millimolar).
 * RATES[3] is d/dt ACLAC in component ACLAC (millimolar).
 * RATES[4] is d/dt ACET in component ACET (millimolar).
 * RATES[5] is d/dt ATP in component ATP (millimolar).
 * RATES[6] is d/dt NADH in component NADH (millimolar).
 * RATES[7] is d/dt ACCOA in component ACCOA (millimolar).
 * There are a total of 13 condition variables.
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 1;
STATES[1] = 0.03145;
STATES[2] = 0.11;
STATES[3] = 1e-5;
STATES[4] = 1e-5;
STATES[5] = 0.1;
CONSTANTS[0] = 5;
STATES[6] = 3.67;
CONSTANTS[1] = 10;
STATES[7] = 0.11;
CONSTANTS[2] = 1;
CONSTANTS[3] = 0.01;
CONSTANTS[4] = 0.01;
CONSTANTS[5] = 0.1;
CONSTANTS[6] = 15;
CONSTANTS[7] = 0.1;
CONSTANTS[8] = 0.2;
CONSTANTS[9] = 10;
CONSTANTS[10] = 0.1;
CONSTANTS[11] = 0.1412;
CONSTANTS[12] = 0.04699;
CONSTANTS[13] = 2.5;
CONSTANTS[14] = 0.08999;
CONSTANTS[15] = 0.01867;
CONSTANTS[16] = 2397;
CONSTANTS[17] = 21120.69;
CONSTANTS[18] = 100;
CONSTANTS[19] = 2.4;
CONSTANTS[20] = 1.5;
CONSTANTS[21] = 0.08;
CONSTANTS[22] = 5118;
CONSTANTS[23] = 46.4159;
CONSTANTS[24] = 0.4;
CONSTANTS[25] = 0.014;
CONSTANTS[26] = 1;
CONSTANTS[27] = 0.1;
CONSTANTS[28] = 0.008;
CONSTANTS[29] = 259;
CONSTANTS[30] = 0.0065;
CONSTANTS[31] = 2.6;
CONSTANTS[32] = 2.6;
CONSTANTS[33] = 0.029;
CONSTANTS[34] = 0.7;
CONSTANTS[35] = 0.2;
CONSTANTS[36] = 0.2;
CONSTANTS[37] = 42;
CONSTANTS[38] = 174.217;
CONSTANTS[39] = 7;
CONSTANTS[40] = 0.07;
CONSTANTS[41] = 0.5;
CONSTANTS[42] = 0.16;
CONSTANTS[43] = 2700;
CONSTANTS[44] = 1;
CONSTANTS[45] = 0.08;
CONSTANTS[46] = 0.025;
CONSTANTS[47] = 0.008;
CONSTANTS[48] = 0.007;
CONSTANTS[49] = 10;
CONSTANTS[50] = 97;
CONSTANTS[51] = 12354.9;
CONSTANTS[52] = 0.08;
CONSTANTS[53] = 0.05;
CONSTANTS[54] = 1;
CONSTANTS[55] = 0.03;
CONSTANTS[56] = 162;
CONSTANTS[57] = 2.4;
CONSTANTS[58] = 9e12;
CONSTANTS[59] = 100;
CONSTANTS[60] = 50;
CONSTANTS[61] = 600;
CONSTANTS[62] = 10;
CONSTANTS[63] = 100;
CONSTANTS[64] = 106;
CONSTANTS[65] = 5;
CONSTANTS[66] = 200;
CONSTANTS[67] = 1400;
CONSTANTS[68] = 0.16;
CONSTANTS[69] = 0.02;
CONSTANTS[70] = 2.6;
CONSTANTS[71] = 0.06;
CONSTANTS[72] = 105;
CONSTANTS[73] = 2.58;
CONSTANTS[74] = 6.196;
CONSTANTS[75] = 900;
CONSTANTS[76] = 1;
CONSTANTS[77] = 0.041;
CONSTANTS[78] = 0.2;
CONSTANTS[79] = 118;
CONSTANTS[80] = 0.0003;
RATES[0] = 0.1001;
RATES[1] = 0.1001;
RATES[2] = 0.1001;
RATES[3] = 0.1001;
RATES[4] = 0.1001;
RATES[5] = 0.1001;
RATES[6] = 0.1001;
RATES[7] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
                 double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[0] - ALGEBRAIC[17] - (ALGEBRAIC[18]+ALGEBRAIC[19]+ALGEBRAIC[24]);
resid[1] = RATES[1] - ALGEBRAIC[20] - ALGEBRAIC[21];
resid[2] = RATES[2] - ALGEBRAIC[22] - ALGEBRAIC[23];
resid[3] = RATES[3] -  0.500000*ALGEBRAIC[24] - (ALGEBRAIC[25]+ALGEBRAIC[16]);
resid[4] = RATES[4] - (ALGEBRAIC[25]+ALGEBRAIC[16]) - (ALGEBRAIC[27]+ALGEBRAIC[26]);
resid[5] = RATES[5] - (ALGEBRAIC[17]+ALGEBRAIC[21]) - ALGEBRAIC[28];
resid[6] = RATES[6] - (ALGEBRAIC[17]+ALGEBRAIC[19]) - (ALGEBRAIC[18]+ALGEBRAIC[22]+ALGEBRAIC[23]+ALGEBRAIC[27]+ALGEBRAIC[29]);
resid[7] = RATES[7] - ALGEBRAIC[19] - (ALGEBRAIC[22]+ALGEBRAIC[20]);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[16] =  CONSTANTS[80]*STATES[3];
ALGEBRAIC[0] = CONSTANTS[0] - STATES[5];
ALGEBRAIC[1] = CONSTANTS[1] - STATES[6];
ALGEBRAIC[3] = (( (( (( CONSTANTS[16]*CONSTANTS[6])/CONSTANTS[10])*ALGEBRAIC[1])/CONSTANTS[11])*ALGEBRAIC[0])/CONSTANTS[12])/( (1.00000+CONSTANTS[6]/CONSTANTS[10]+STATES[0]/CONSTANTS[13])*(1.00000+ALGEBRAIC[1]/CONSTANTS[11]+STATES[6]/CONSTANTS[14])*(1.00000+ALGEBRAIC[0]/CONSTANTS[12]+STATES[5]/CONSTANTS[15]));
ALGEBRAIC[17] = (CONDVAR[0]>=0.00000 ? ALGEBRAIC[3] : ALGEBRAIC[3]);
ALGEBRAIC[4] = ( (CONSTANTS[22]/( CONSTANTS[20]*CONSTANTS[21]))*( STATES[0]*STATES[6] - ( CONSTANTS[7]*ALGEBRAIC[1])/CONSTANTS[17]))/( (1.00000+STATES[0]/CONSTANTS[20]+CONSTANTS[7]/CONSTANTS[18])*(1.00000+STATES[6]/CONSTANTS[21]+ALGEBRAIC[1]/CONSTANTS[19]));
ALGEBRAIC[18] = (CONDVAR[1]>=0.00000 ? ALGEBRAIC[4] : ALGEBRAIC[4]);
ALGEBRAIC[2] = CONSTANTS[2] - STATES[7];
ALGEBRAIC[5] = (( (( (( (CONSTANTS[29]/(1.00000+( CONSTANTS[23]*STATES[6])/ALGEBRAIC[1]))*STATES[0])/CONSTANTS[26])*ALGEBRAIC[1])/CONSTANTS[24])*ALGEBRAIC[2])/CONSTANTS[25])/( (1.00000+STATES[0]/CONSTANTS[26])*(1.00000+ALGEBRAIC[1]/CONSTANTS[24]+STATES[6]/CONSTANTS[27])*(1.00000+ALGEBRAIC[2]/CONSTANTS[25]+STATES[7]/CONSTANTS[28]));
ALGEBRAIC[19] = (CONDVAR[2]>=0.00000 ? ALGEBRAIC[5] : ALGEBRAIC[5]);
ALGEBRAIC[6] = ( (CONSTANTS[37]/( CONSTANTS[36]*CONSTANTS[31]))*( STATES[7]*CONSTANTS[9] - ( STATES[1]*ALGEBRAIC[2])/CONSTANTS[30]))/(1.00000+STATES[7]/CONSTANTS[36]+CONSTANTS[9]/CONSTANTS[32]+STATES[1]/CONSTANTS[35]+ALGEBRAIC[2]/CONSTANTS[33]+( STATES[7]*CONSTANTS[9])/( CONSTANTS[36]*CONSTANTS[31])+( STATES[1]*ALGEBRAIC[2])/( CONSTANTS[34]*CONSTANTS[33]));
ALGEBRAIC[20] = (CONDVAR[3]>=0.00000 ? ALGEBRAIC[6] : ALGEBRAIC[6]);
ALGEBRAIC[7] = ( (CONSTANTS[43]/( CONSTANTS[41]*CONSTANTS[42]))*( STATES[1]*ALGEBRAIC[0] - ( CONSTANTS[3]*STATES[5])/CONSTANTS[38]))/( (1.00000+STATES[1]/CONSTANTS[42]+CONSTANTS[3]/CONSTANTS[39])*(1.00000+ALGEBRAIC[0]/CONSTANTS[41]+STATES[5]/CONSTANTS[40]));
ALGEBRAIC[21] = (CONDVAR[4]>=0.00000 ? ALGEBRAIC[7] : ALGEBRAIC[7]);
ALGEBRAIC[8] = ( (CONSTANTS[50]/( CONSTANTS[48]*CONSTANTS[46]))*( STATES[7]*STATES[6] - ( ALGEBRAIC[1]*ALGEBRAIC[2]*STATES[2])/CONSTANTS[44]))/( (1.00000+ALGEBRAIC[1]/CONSTANTS[45]+STATES[6]/CONSTANTS[46])*(1.00000+STATES[7]/CONSTANTS[48]+ALGEBRAIC[2]/CONSTANTS[47])*(1.00000+STATES[2]/CONSTANTS[49]));
ALGEBRAIC[22] = (CONDVAR[5]>=0.00000 ? ALGEBRAIC[8] : ALGEBRAIC[8]);
ALGEBRAIC[9] = ( (CONSTANTS[56]/( CONSTANTS[55]*CONSTANTS[53]))*( STATES[2]*STATES[6] - ( CONSTANTS[5]*ALGEBRAIC[1])/CONSTANTS[51]))/( (1.00000+ALGEBRAIC[1]/CONSTANTS[52]+STATES[6]/CONSTANTS[53])*(1.00000+STATES[2]/CONSTANTS[55]+CONSTANTS[5]/CONSTANTS[54]));
ALGEBRAIC[23] = (CONDVAR[6]>=0.00000 ? ALGEBRAIC[9] : ALGEBRAIC[9]);
ALGEBRAIC[10] = ( (( CONSTANTS[61]*STATES[0])/CONSTANTS[60])*(1.00000 - STATES[3]/( STATES[0]*CONSTANTS[58]))*pow(STATES[0]/CONSTANTS[60]+STATES[3]/CONSTANTS[59], CONSTANTS[57] - 1.00000))/(1.00000+pow(STATES[0]/CONSTANTS[60]+STATES[3]/CONSTANTS[59], CONSTANTS[57]));
ALGEBRAIC[24] = (CONDVAR[7]>=0.00000 ? ALGEBRAIC[10] : ALGEBRAIC[10]);
ALGEBRAIC[11] = (( CONSTANTS[64]*STATES[3])/CONSTANTS[62])/(1.00000+STATES[3]/CONSTANTS[62]+STATES[4]/CONSTANTS[63]);
ALGEBRAIC[25] = (CONDVAR[8]>=0.00000 ? ALGEBRAIC[11] : ALGEBRAIC[11]);
ALGEBRAIC[12] = (( CONSTANTS[66]*STATES[4])/CONSTANTS[65])/(1.00000+STATES[4]/CONSTANTS[65]);
ALGEBRAIC[26] = (CONDVAR[9]>=0.00000 ? ALGEBRAIC[12] : ALGEBRAIC[12]);
ALGEBRAIC[13] = ( (CONSTANTS[72]/( CONSTANTS[71]*CONSTANTS[69]))*( STATES[4]*STATES[6] - ( CONSTANTS[4]*ALGEBRAIC[1])/CONSTANTS[67]))/( (1.00000+STATES[4]/CONSTANTS[71]+CONSTANTS[4]/CONSTANTS[70])*(1.00000+STATES[6]/CONSTANTS[69]+ALGEBRAIC[1]/CONSTANTS[68]));
ALGEBRAIC[27] = (CONDVAR[10]>=0.00000 ? ALGEBRAIC[13] : ALGEBRAIC[13]);
ALGEBRAIC[14] = ( CONSTANTS[75]*pow(STATES[5]/ALGEBRAIC[0], CONSTANTS[73]))/(pow(CONSTANTS[74], CONSTANTS[73])+pow(STATES[5]/ALGEBRAIC[0], CONSTANTS[73]));
ALGEBRAIC[28] = (CONDVAR[11]>=0.00000 ? ALGEBRAIC[14] : ALGEBRAIC[14]);
ALGEBRAIC[15] = (( CONSTANTS[79]*STATES[6]*CONSTANTS[8])/( CONSTANTS[77]*CONSTANTS[78]))/( (1.00000+STATES[6]/CONSTANTS[77]+ALGEBRAIC[1]/CONSTANTS[76])*(1.00000+CONSTANTS[8]/CONSTANTS[78]));
ALGEBRAIC[29] = (CONDVAR[12]>=0.00000 ? ALGEBRAIC[15] : ALGEBRAIC[15]);
}
void
getStateInformation(double* SI)
{
SI[0] = 1.0;
SI[1] = 1.0;
SI[2] = 1.0;
SI[3] = 1.0;
SI[4] = 1.0;
SI[5] = 1.0;
SI[6] = 1.0;
SI[7] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = ALGEBRAIC[3] - 0.00000;
CONDVAR[1] = ALGEBRAIC[4] - 0.00000;
CONDVAR[2] = ALGEBRAIC[5] - 0.00000;
CONDVAR[3] = ALGEBRAIC[6] - 0.00000;
CONDVAR[4] = ALGEBRAIC[7] - 0.00000;
CONDVAR[5] = ALGEBRAIC[8] - 0.00000;
CONDVAR[6] = ALGEBRAIC[9] - 0.00000;
CONDVAR[7] = ALGEBRAIC[10] - 0.00000;
CONDVAR[8] = ALGEBRAIC[11] - 0.00000;
CONDVAR[9] = ALGEBRAIC[12] - 0.00000;
CONDVAR[10] = ALGEBRAIC[13] - 0.00000;
CONDVAR[11] = ALGEBRAIC[14] - 0.00000;
CONDVAR[12] = ALGEBRAIC[15] - 0.00000;
}