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 3 entries in the algebraic variable array.
   There are a total of 1 entries in each of the rate and state variable arrays.
   There are a total of 7 entries in the constant variable array.
 */
/*
 * VOI is x in component main (dimensionless).
 * ALGEBRAIC[0] is sin in component sin (dimensionless).
 * STATES[0] is sin in component sin (dimensionless).
 * CONSTANTS[0] is deriv_approx_initial_value in component main (dimensionless).
 * ALGEBRAIC[2] is sin in component sin (dimensionless).
 * CONSTANTS[2] is k2_oPi in component sin (dimensionless).
 * CONSTANTS[3] is k2Pi in component sin (dimensionless).
 * CONSTANTS[4] is kPi_2 in component sin (dimensionless).
 * CONSTANTS[5] is kPi in component sin (dimensionless).
 * CONSTANTS[6] is kPi_32 in component sin (dimensionless).
 * ALGEBRAIC[1] is z in component sin (dimensionless).
 * CONSTANTS[1] is C in component sin (dimensionless).
 * RATES[0] is d/dt sin in component sin (dimensionless).
 * There are a total of 6 condition variables.
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 0;
CONSTANTS[1] = 0.75;
CONSTANTS[2] = 2.00000/ 3.14159265358979;
CONSTANTS[3] =  2.00000* 3.14159265358979;
CONSTANTS[4] =  3.14159265358979/2.00000;
CONSTANTS[5] =  3.14159265358979;
CONSTANTS[6] = ( 3.00000* 3.14159265358979)/2.00000;
STATES[0] = CONSTANTS[0];
RATES[0] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
                 double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[0] - cos(VOI);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] =  sin(VOI);
ALGEBRAIC[1] = (CONDVAR[0]<0.00000 ?  VOI*CONSTANTS[2] - 0.500000 : CONDVAR[1]<0.00000 ?  ( 3.14159265358979 - VOI)*CONSTANTS[2] - 0.500000 : CONDVAR[2]<0.00000 ?  (VOI -  3.14159265358979)*CONSTANTS[2] - 0.500000 :  (CONSTANTS[3] - VOI)*CONSTANTS[2] - 0.500000);
ALGEBRAIC[2] = (CONDVAR[3]<0.00000 ? - ( ALGEBRAIC[1]*ALGEBRAIC[1])+CONSTANTS[1]+ALGEBRAIC[1] : CONDVAR[4]<0.00000 ? - ( ALGEBRAIC[1]*ALGEBRAIC[1])+CONSTANTS[1]+ALGEBRAIC[1] : CONDVAR[5]<0.00000 ? ( ALGEBRAIC[1]*ALGEBRAIC[1] - CONSTANTS[1]) - ALGEBRAIC[1] : ( ALGEBRAIC[1]*ALGEBRAIC[1] - CONSTANTS[1]) - ALGEBRAIC[1]);
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
}
void
getStateInformation(double* SI)
{
SI[0] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = VOI - CONSTANTS[4];
CONDVAR[1] = VOI - CONSTANTS[5];
CONDVAR[2] = VOI - CONSTANTS[6];
CONDVAR[3] = VOI - CONSTANTS[4];
CONDVAR[4] = VOI - CONSTANTS[5];
CONDVAR[5] = VOI - CONSTANTS[6];
}