public void TestWithParameters()
{
// Create a model
MuCell.Model.SBML.Model model = new MuCell.Model.SBML.Model();
MuCell.Model.SBML.Species species1 = new MuCell.Model.SBML.Species();
MuCell.Model.SBML.Species species2 = new MuCell.Model.SBML.Species();
MuCell.Model.SBML.Parameter parameter = new MuCell.Model.SBML.Parameter();
parameter.ID = "k_1";
parameter.Value = 0.07d;
model.listOfParameters = new List<MuCell.Model.SBML.Parameter>();
model.listOfParameters.Add(parameter);
model.listOfSpecies = new List<MuCell.Model.SBML.Species>();
model.listOfSpecies.Add(species1);
model.listOfSpecies.Add(species2);
// Set some values for species1
species1.ID = "s1";
species1.InitialAmount = 1.0d;
// Set some values for species2
species2.ID = "s2";
species2.InitialAmount = 12.0d;
model.AddId("s1", species1);
model.AddId("s2", species2);
model.AddId("k_1", parameter);
// Set up the reaction
MuCell.Model.SBML.Reaction reaction1 = new MuCell.Model.SBML.Reaction("reaction1");
// Set up the kinetic law
reaction1.KineticLaw = new MuCell.Model.SBML.KineticLaw(model);
reaction1.KineticLaw.Formula = "s1*s2*k_1";
// set up the species reference for the reactants
MuCell.Model.SBML.SpeciesReference ref1 = new MuCell.Model.SBML.SpeciesReference(species1, 1);
// set up the species references for the products
MuCell.Model.SBML.SpeciesReference ref2 = new MuCell.Model.SBML.SpeciesReference(species2, 3);
// Add the references
reaction1.Reactants.Add(ref1);
reaction1.Products.Add(ref2);
// set up the cell definition
MuCell.Model.CellDefinition celldef1 = new MuCell.Model.CellDefinition("cell1");
celldef1.addSBMLModel(model);
// set up a cell
MuCell.Model.CellInstance cell1 = celldef1.createCell();
cell1.localSpeciesDelta.Add("s1", 0.0d);
cell1.localSpeciesDelta.Add("s2", 0.0d);
cell1.localSpeciesVariables.Add("s1", 1.0d);
cell1.localSpeciesVariables.Add("s2", 12.0d);
// set up the statesnaphost
MuCell.Model.StateSnapshot state1 = new MuCell.Model.StateSnapshot();
state1.Cells.Add(cell1);
// get the eval function
MuCell.Model.EffectReactionEvaluationFunction fun1 = reaction1.ToEffectReactionEvaluationFunction();
// Evaluate
fun1(state1, cell1);
cell1.setSpeciesAmount("s1", cell1.UnitTestGetSpeciesAmountsAndDeltas("s1"));
cell1.setSpeciesAmount("s2", cell1.UnitTestGetSpeciesAmountsAndDeltas("s2"));
cell1.UnitTestClearDeltas();
AssertDouble.AreEqual(1.0d-0.84d, cell1.getLocalSimulationSpeciesAmount("s1"));
AssertDouble.AreEqual(12.0d+0.84d*3.0d, cell1.getLocalSimulationSpeciesAmount("s2"));
System.Console.WriteLine("about to eval");
// Evaluate again
fun1(state1, cell1);
cell1.setSpeciesAmount("s1", cell1.UnitTestGetSpeciesAmountsAndDeltas("s1"));
cell1.setSpeciesAmount("s2", cell1.UnitTestGetSpeciesAmountsAndDeltas("s2"));
cell1.UnitTestClearDeltas();
AssertDouble.AreEqual(0.16d-0.162624d, cell1.getLocalSimulationSpeciesAmount("s1"));
AssertDouble.AreEqual(14.52d+0.162624d*3.0d, cell1.getLocalSimulationSpeciesAmount("s2"));
// Check the param is the same
Assert.AreEqual(0.07m, parameter.Value);
}
public void testInterrogateModelForIDs()
{
MuCell.Model.SBML.Model model = new MuCell.Model.SBML.Model();
model.listOfReactions = new List<MuCell.Model.SBML.Reaction>();
model.listOfSpecies = new List<MuCell.Model.SBML.Species>();
MuCell.Model.SBML.Species species1 = new MuCell.Model.SBML.Species();
MuCell.Model.SBML.Species species2 = new MuCell.Model.SBML.Species();
MuCell.Model.SBML.Species species3 = new MuCell.Model.SBML.Species();
species1.ID = "species1";
species2.ID = "species2";
species3.ID = "species 3";
model.AddId("species1", species1);
MuCell.Model.SBML.Reaction r1 = new MuCell.Model.SBML.Reaction();
MuCell.Model.SBML.SpeciesReference sp1r = new MuCell.Model.SBML.SpeciesReference();
sp1r.species = species1;
MuCell.Model.SBML.SpeciesReference sp2r = new MuCell.Model.SBML.SpeciesReference();
sp2r.species = species2;
MuCell.Model.SBML.SpeciesReference sp3r = new MuCell.Model.SBML.SpeciesReference();
sp3r.species = species3;
r1.Reactants.Add(sp1r);
r1.Reactants.Add(sp2r);
r1.Products.Add(sp3r);
MuCell.Model.SBML.KineticLaw kl = new MuCell.Model.SBML.KineticLaw(model);
r1.KineticLaw = kl;
MuCell.Model.SBML.Parameter k1 = new MuCell.Model.SBML.Parameter();
k1.ID = "k1";
r1.Parameters.Add(k1);
r1.Formula = "species1*species_3*k1";
model.listOfSpecies.Add(species1);
model.listOfReactions.Add(r1);
// Assert the unknown entities
Assert.AreEqual(2, r1.KineticLaw.UnknownEntitiesFromFormula().Count);
Assert.AreEqual("species_3", r1.KineticLaw.UnknownEntitiesFromFormula().ToArray()[0].ID);
Assert.AreEqual("k1", r1.KineticLaw.UnknownEntitiesFromFormula().ToArray()[1].ID);
// Assert missing items in model
Assert.IsTrue(model.idExists("species1"));
Assert.IsFalse(model.idExists("species2"));
Assert.IsFalse(model.idExists("species 3"));
Assert.IsFalse(model.idExists("k1"));
model.InterrogateModelForMissingIDs();
// redo formula
r1.Formula = "species1*species_3*k1";
// assert no unknown entities
Assert.AreEqual(0, r1.KineticLaw.UnknownEntitiesFromFormula().Count);
// Assert everything is in the model
Assert.IsTrue(model.idExists("species1"));
Assert.IsTrue(model.idExists("species2"));
Assert.IsTrue(model.idExists("species 3"));
Assert.IsTrue(model.idExists("k1"));
}
