END_TEST
START_TEST (test_SBMLDocument_setLevelAndVersion_UnitsError)
{
SBMLDocument_t *d = SBMLDocument_create();
SBMLDocument_setLevelAndVersionNonStrict(d, 2, 4);
Model_t *m1 = SBMLDocument_createModel(d);
Compartment_t *c = Model_createCompartment(m1);
Compartment_setId(c, "c");
Parameter_t *p = Model_createParameter(m1);
Parameter_setId(p, "p");
Parameter_setUnits(p, "mole");
Rule_t * r = Model_createAssignmentRule(m1);
Rule_setVariable(r, "c");
Rule_setFormula(r, "p*p");
fail_unless(SBMLDocument_setLevelAndVersionNonStrict(d,2,2) == 1);
fail_unless(SBMLDocument_setLevelAndVersionNonStrict(d,2,3) == 1);
fail_unless(SBMLDocument_setLevelAndVersionNonStrict(d,1,2) == 1);
fail_unless(SBMLDocument_setLevelAndVersionNonStrict(d,1,1) == 0);
SBMLDocument_free(d);
}
END_TEST
START_TEST (test_SBMLConvertStrict_convertL1ParamRule)
{
SBMLDocument_t *d = SBMLDocument_createWithLevelAndVersion(1, 2);
Model_t * m = SBMLDocument_createModel(d);
/* create a compartment */
Compartment_t * c = Model_createCompartment(m);
Compartment_setId(c, "c");
/* create a parameter */
Parameter_t * p = Model_createParameter(m);
Parameter_setId(p, "p");
Parameter_t * p1 = Model_createParameter(m);
Parameter_setId(p1, "p1");
/* create a math element */
ASTNode_t *math = SBML_parseFormula("p");
/* create an assignment rule */
Rule_t *ar = Model_createAssignmentRule(m);
Rule_setVariable(ar, "p1");
Rule_setMath(ar, math);
Rule_setUnits(ar, "mole");
fail_unless( SBMLDocument_setLevelAndVersionStrict(d, 2, 1) == 1 );
fail_unless( SBMLDocument_getLevel (d) == 2, NULL );
fail_unless( SBMLDocument_getVersion(d) == 1, NULL );
Rule_t * r1 = Model_getRule(SBMLDocument_getModel(d), 0);
fail_unless (Rule_getUnits(r1) == NULL );
SBMLDocument_free(d);
}
END_TEST
START_TEST (test_SBMLConvert_convertToL3_parameter)
{
SBMLDocument_t *d = SBMLDocument_createWithLevelAndVersion(2, 2);
Model_t *m = SBMLDocument_createModel(d);
const char *sid = "C";
Parameter_t *p = Model_createParameter(m);
Parameter_t *p1;
Parameter_setId ( p, sid );
fail_unless( SBMLDocument_setLevelAndVersionNonStrict(d, 3, 1) == 1, NULL);
p1 = Model_getParameter(m, 0);
fail_unless(Parameter_hasRequiredAttributes(p1) == 1);
SBMLDocument_free(d);
}
END_TEST
START_TEST (test_SBMLConvert_convertFromL3_conversionFactor)
{
SBMLDocument_t *d = SBMLDocument_createWithLevelAndVersion(3, 1);
Model_t *m = SBMLDocument_createModel(d);
const char *sid = "P";
Model_setConversionFactor(m, sid);
Parameter_t *c = Model_createParameter(m);
Parameter_setId ( c, sid );
Parameter_setConstant( c, 1);
fail_unless(SBMLDocument_setLevelAndVersion(d, 1, 1) == 0);
fail_unless(SBMLDocument_setLevelAndVersion(d, 1, 2) == 0);
fail_unless(SBMLDocument_setLevelAndVersion(d, 2, 1) == 0);
fail_unless(SBMLDocument_setLevelAndVersion(d, 2, 2) == 0);
fail_unless(SBMLDocument_setLevelAndVersion(d, 2, 3) == 0);
fail_unless(SBMLDocument_setLevelAndVersion(d, 2, 4) == 0);
fail_unless(SBMLDocument_setLevelAndVersion(d, 3, 1) == 1);
}
END_TEST
START_TEST (test_create_l1v1_units)
{
Model_t *m;
Compartment_t *c;
KineticLaw_t *kl;
Parameter_t *p;
Reaction_t *r;
Species_t *s;
SpeciesReference_t *sr;
Unit_t *u;
UnitDefinition_t *ud;
/**
* <sbml level="1" version="1">
*/
m = Model_create(2, 4);
/**
* <unitDefinition name="substance">
* <listOfUnits>
* <unit kind="mole" scale="-3"/>
* </listOfUnits>
* </unitDefinition>
*/
ud = Model_createUnitDefinition(m);
UnitDefinition_setName(ud, "substance");
u = Model_createUnit(m);
Unit_setKind(u,UNIT_KIND_MOLE);
Unit_setScale(u, -3);
/**
* <unitDefinition name="mls">
* <listOfUnits>
* <unit kind="mole" scale="-3"/>
* <unit kind="liter" exponent="-1"/>
* <unit kind="second" exponent="-1"/>
* </listOfUnits>
* </unitDefinition>
*/
ud = Model_createUnitDefinition(m);
UnitDefinition_setName(ud, "mls");
u = Model_createUnit(m);
Unit_setKind(u, UNIT_KIND_MOLE);
Unit_setScale(u, -3);
u = Model_createUnit(m);
Unit_setKind(u, UNIT_KIND_LITER);
Unit_setExponent(u, -1);
u = Model_createUnit(m);
Unit_setKind(u, UNIT_KIND_SECOND);
Unit_setExponent(u, -1);
/**
* <listOfCompartments>
* <compartment name="cell"/>
* </listOfCompartments>
*/
c = Model_createCompartment(m);
Compartment_setName(c, "cell");
/**
* <listOfSpecies>
* <specie name="x0" compartment="cell" initialAmount="1"/>
* <specie name="x1" compartment="cell" initialAmount="1"/>
* <specie name="s1" compartment="cell" initialAmount="1"/>
* <specie name="s2" compartment="cell" initialAmount="1"/>
* </listOfSpecies>
*/
s = Model_createSpecies(m);
Species_setName(s, "x0");
Species_setCompartment(s, "cell");
Species_setInitialAmount(s, 1);
s = Model_createSpecies(m);
Species_setName(s, "x1");
Species_setCompartment(s, "cell");
Species_setInitialAmount(s, 1);
s = Model_createSpecies(m);
Species_setName(s, "s1");
Species_setCompartment(s, "cell");
Species_setInitialAmount(s, 1);
s = Model_createSpecies(m);
Species_setName(s, "s2");
Species_setCompartment(s, "cell");
Species_setInitialAmount(s, 1);
/**
* <listOfParameters>
* <parameter name="vm" value="2" units="mls"/>
* <parameter name="km" value="2"/>
* </listOfParameters>
*/
p = Model_createParameter(m);
Parameter_setName (p, "vm");
Parameter_setUnits(p, "mls");
Parameter_setValue(p, 2);
p = Model_createParameter(m);
Parameter_setName (p, "km");
Parameter_setValue(p, 2);
/**
* <reaction name="v1">
* <listOfReactants>
* <specieReference specie="x0"/>
* </listOfReactants>
* <listOfProducts>
* <specieReference specie="s1"/>
* </listOfProducts>
* <kineticLaw formula="(vm * s1)/(km + s1)"/>
* </reaction>
*/
r = Model_createReaction(m);
Reaction_setName(r, "v1");
sr = Model_createReactant(m);
SpeciesReference_setSpecies(sr, "x0");
sr = Model_createProduct(m);
SpeciesReference_setSpecies(sr, "s1");
kl = Model_createKineticLaw(m);
KineticLaw_setFormula(kl, "(vm * s1)/(km + s1)");
/**
* <reaction name="v2">
* <listOfReactants>
* <specieReference specie="s1"/>
* </listOfReactants>
* <listOfProducts>
* <specieReference specie="s2"/>
* </listOfProducts>
* <kineticLaw formula="(vm * s2)/(km + s2)"/>
* </reaction>
*/
r = Model_createReaction(m);
Reaction_setName(r, "v2");
sr = Model_createReactant(m);
SpeciesReference_setSpecies(sr, "s1");
sr = Model_createProduct(m);
SpeciesReference_setSpecies(sr, "s2");
kl = Model_createKineticLaw(m);
KineticLaw_setFormula(kl, "(vm * s2)/(km + s2)");
/**
* <reaction name="v3">
* <listOfReactants>
* <specieReference specie="s2"/>
* </listOfReactants>
* <listOfProducts>
* <specieReference specie="x1"/>
* </listOfProducts>
* <kineticLaw formula="(vm * s1)/(km + s1)"/>
* </reaction>
*/
r = Model_createReaction(m);
Reaction_setName(r, "v3");
sr = Model_createReactant(m);
SpeciesReference_setSpecies(sr, "s2");
sr = Model_createProduct(m);
SpeciesReference_setSpecies(sr, "x1");
kl = Model_createKineticLaw(m);
KineticLaw_setFormula(kl, "(vm * s1)/(km + s1)");
Model_free(m);
}
