END_TEST
START_TEST (test_UnitDefinition_isVariantOfVolume_2)
{
Unit_t *dim = Unit_create(2, 4);
Unit_setKind( dim , UnitKind_forName("dimensionless") );
Unit_t *u = UnitDefinition_createUnit(UD);
fail_unless( !UnitDefinition_isVariantOfVolume(UD) );
Unit_setKind(u, UNIT_KIND_METRE);
Unit_setExponent(u, 3);
fail_unless( UnitDefinition_isVariantOfVolume(UD) );
Unit_setScale(u, 100);
fail_unless( UnitDefinition_isVariantOfVolume(UD) );
Unit_setMultiplier(u, 5);
fail_unless( UnitDefinition_isVariantOfVolume(UD) );
Unit_setOffset(u, -5);
fail_unless( UnitDefinition_isVariantOfVolume(UD) );
Unit_setExponent(u, 2);
fail_unless( !UnitDefinition_isVariantOfVolume(UD) );
Unit_setExponent(u, 3);
UnitDefinition_addUnit( UD, dim );
fail_unless( UnitDefinition_isVariantOfVolume(UD) );
}
END_TEST
START_TEST (test_UnitDefinition_isVariantOfTime)
{
Unit_t *dim = Unit_create(2, 4);
Unit_setKind( dim , UnitKind_forName("dimensionless") );
Unit_t *u = UnitDefinition_createUnit(UD);
fail_unless( !UnitDefinition_isVariantOfTime(UD) );
Unit_setKind(u, UNIT_KIND_SECOND);
Unit_setExponent(u, 1);
fail_unless( UnitDefinition_isVariantOfTime(UD) );
Unit_setScale(u, -10);
fail_unless( UnitDefinition_isVariantOfTime(UD) );
Unit_setMultiplier(u, 10);
fail_unless( UnitDefinition_isVariantOfTime(UD) );
Unit_setOffset(u, 30);
fail_unless( UnitDefinition_isVariantOfTime(UD) );
Unit_setExponent(u, 2);
fail_unless( !UnitDefinition_isVariantOfTime(UD) );
Unit_setExponent(u, 1);
UnitDefinition_addUnit( UD, dim );
fail_unless( UnitDefinition_isVariantOfTime(UD) );
Unit_free(dim);
}
END_TEST
START_TEST (test_Unit_setExponent1)
{
int i = Unit_setExponent(U, 2);
fail_unless( i == LIBSBML_OPERATION_SUCCESS );
fail_unless( Unit_getExponent(U) == 2 );
}
END_TEST
START_TEST (test_UnitDefinition_printUnits)
{
UnitDefinition_t *ud = UnitDefinition_create(2, 4);
UnitDefinition_setId(ud, "mmls");
Unit_t *perTime = UnitDefinition_createUnit(ud);
Unit_setKind( perTime , UnitKind_forName("second") );
Unit_setExponent( perTime, -1);
char * ud_str = UnitDefinition_printUnits(ud, 0);
fail_unless(!strcmp(ud_str,
"second (exponent = -1, multiplier = 1, scale = 0)"));
char * ud_str1 = UnitDefinition_printUnits(ud, 1);
fail_unless(!strcmp(ud_str1, "(1 second)^-1"));
UnitDefinition_t *ud1 = UnitDefinition_create(2, 4);
UnitDefinition_setId(ud1, "mmls");
Unit_t *u = UnitDefinition_createUnit(ud1);
Unit_setKind(u, UNIT_KIND_KILOGRAM);
Unit_setExponent(u, 1);
Unit_setScale(u, 2);
Unit_setMultiplier(u, 3.0);
char * ud_str2 = UnitDefinition_printUnits(ud1, 0);
fail_unless(!strcmp(ud_str2,
"kilogram (exponent = 1, multiplier = 3, scale = 2)"));
char * ud_str3 = UnitDefinition_printUnits(ud1, 1);
fail_unless(!strcmp(ud_str3, "(300 kilogram)^1"));
safe_free(ud_str);
safe_free(ud_str1);
safe_free(ud_str2);
safe_free(ud_str3);
UnitDefinition_free(ud);
UnitDefinition_free(ud1);
}
END_TEST
START_TEST (test_UnitDefinition_getUnit)
{
Unit_t *mole = Unit_create(2, 4);
Unit_t *litre = Unit_create(2, 4);
Unit_t *second = Unit_create(2, 4);
Unit_setKind( mole , UnitKind_forName("mole") );
Unit_setKind( litre , UnitKind_forName("litre") );
Unit_setKind( second, UnitKind_forName("second") );
Unit_setScale (mole , -3);
Unit_setExponent(litre , -1);
Unit_setExponent(second, -1);
UnitDefinition_addUnit( UD, mole );
UnitDefinition_addUnit( UD, litre );
UnitDefinition_addUnit( UD, second );
Unit_free(mole);
Unit_free(litre);
Unit_free(second);
fail_unless( UnitDefinition_getNumUnits(UD) == 3 );
mole = UnitDefinition_getUnit(UD, 0);
litre = UnitDefinition_getUnit(UD, 1);
second = UnitDefinition_getUnit(UD, 2);
fail_unless( Unit_getKind(mole) == UNIT_KIND_MOLE );
fail_unless( Unit_getKind(litre) == UNIT_KIND_LITRE );
fail_unless( Unit_getKind(second) == UNIT_KIND_SECOND );
fail_unless( Unit_getScale(mole) == -3 );
fail_unless( Unit_getExponent(litre) == -1 );
fail_unless( Unit_getExponent(second) == -1 );
}
END_TEST
START_TEST (test_Unit_set_get_unset)
{
Unit_t *u = Unit_create(2, 4);
// defaults
fail_unless( Unit_getKind (u) == UNIT_KIND_INVALID );
fail_unless( Unit_getExponent (u) == 1 );
fail_unless( Unit_isSetExponent(u) );
fail_unless( Unit_getScale (u) == 0 );
fail_unless( Unit_isSetScale(u) );
fail_unless( Unit_getMultiplier(u) == 1.0 );
fail_unless( Unit_isSetMultiplier(u) );
fail_unless( !Unit_isSetKind(u) );
Unit_setKind(u, UNIT_KIND_WATT);
fail_unless( Unit_getKind (u) == UNIT_KIND_WATT );
Unit_setExponent(u, 3);
fail_unless( Unit_getExponent (u) == 3 );
Unit_setScale(u, 4);
fail_unless( Unit_getScale (u) == 4 );
Unit_setMultiplier(u, 3.2);
fail_unless( Unit_getMultiplier(u) == 3.2 );
fail_unless( Unit_unsetKind(u) == LIBSBML_OPERATION_SUCCESS);
fail_unless( Unit_unsetExponent(u) == LIBSBML_UNEXPECTED_ATTRIBUTE);
fail_unless( Unit_unsetScale(u) == LIBSBML_UNEXPECTED_ATTRIBUTE);
fail_unless( Unit_unsetMultiplier(u) == LIBSBML_UNEXPECTED_ATTRIBUTE);
fail_unless( Unit_getKind (u) == UNIT_KIND_INVALID );
fail_unless( Unit_getExponent (u) == 1 );
fail_unless( Unit_isSetExponent(u) );
fail_unless( Unit_getScale (u) == 0 );
fail_unless( Unit_isSetScale(u) );
fail_unless( Unit_getMultiplier(u) == 1.0 );
fail_unless( Unit_isSetMultiplier(u) );
fail_unless( !Unit_isSetKind(u) );
Unit_free(u);
}
END_TEST
START_TEST (test_Unit_unset)
{
Unit_t *u = Unit_create(2, 4);
Unit_setKind(u, UNIT_KIND_WATT);
Unit_setExponent(u, 3);
Unit_setScale(u, 4);
Unit_setMultiplier(u, 3.2);
fail_unless( Unit_getKind (u) == UNIT_KIND_WATT );
fail_unless( Unit_getExponent (u) == 3 );
fail_unless( Unit_getScale (u) == 4 );
fail_unless( Unit_getMultiplier(u) == 3.2 );
fail_unless (Unit_isSetKind(u) == 1);
fail_unless (Unit_isSetExponent(u) == 1);
fail_unless (Unit_isSetScale(u) == 1);
fail_unless (Unit_isSetMultiplier(u) == 1);
fail_unless (Unit_isSetOffset(u) == 0);
Unit_unsetKind(u);
Unit_unsetExponent(u);
Unit_unsetScale(u);
Unit_unsetMultiplier(u);
fail_unless (Unit_isSetKind(u) == 0);
fail_unless (Unit_isSetExponent(u) == 1);
fail_unless (Unit_isSetScale(u) == 1);
fail_unless (Unit_isSetMultiplier(u) == 1);
fail_unless (Unit_isSetOffset(u) == 0);
fail_unless( Unit_getKind (u) == UNIT_KIND_INVALID );
fail_unless( Unit_getExponent (u) == 1 );
fail_unless( Unit_getScale (u) == 0 );
fail_unless( Unit_getMultiplier(u) == 1.0 );
Unit_free(u);
}
END_TEST
START_TEST (test_UnitDefinition_isVariantOfSubstancePerTime_3)
{
UnitDefinition_t *ud =
UnitDefinition_create(2, 2);
Unit_t *dim = Unit_create(2, 2);
Unit_setKind( dim , UnitKind_forName("dimensionless") );
Unit_t *perTime = UnitDefinition_createUnit(ud);
Unit_setKind( perTime , UnitKind_forName("second") );
Unit_setExponent( perTime, -1);
Unit_t *u = UnitDefinition_createUnit(ud);
fail_unless( !UnitDefinition_isVariantOfSubstancePerTime(ud) );
Unit_setKind(u, UNIT_KIND_GRAM);
Unit_setExponent(u, 1);
fail_unless( UnitDefinition_isVariantOfSubstancePerTime(ud) );
Unit_setScale(u, -1);
Unit_setScale(perTime, -1);
fail_unless( UnitDefinition_isVariantOfSubstancePerTime(ud) );
Unit_setMultiplier(u, 2);
fail_unless( UnitDefinition_isVariantOfSubstancePerTime(ud) );
Unit_setOffset(u, 3);
fail_unless( UnitDefinition_isVariantOfSubstancePerTime(ud) );
Unit_setExponent(u, -3);
fail_unless( !UnitDefinition_isVariantOfSubstancePerTime(ud) );
Unit_setExponent(u, 1);
Unit_setExponent(perTime, -3);
fail_unless( !UnitDefinition_isVariantOfSubstancePerTime(ud) );
Unit_setExponent(perTime, -1);
UnitDefinition_addUnit( ud, dim );
fail_unless( UnitDefinition_isVariantOfSubstancePerTime(ud) );
UnitDefinition_free(ud);
Unit_free(dim);
}
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);
}
