/** @cond doxygenLibsbmlInternal */
int
Replacing::replaceWithAndMaybeDelete(SBase* replacement, bool deleteme, ASTNode* conversionFactor)
{
//All of the following function calls that could result in errors set their own error messages.
SBase* replaced = getReferencedElement();
if (replaced==NULL) return LIBSBML_INVALID_OBJECT;
//Rename things
int ret = updateIDs(replaced, replacement);
if (ret != LIBSBML_OPERATION_SUCCESS) return ret;
//Perform any conversions on references in the submodel.
ret = performConversions(replacement, &conversionFactor);
if (ret != LIBSBML_OPERATION_SUCCESS) return ret;
//Finally, recurse down if there are things the replaced element itself replaced (or used to be replaced by)
CompSBasePlugin* replacedplug = static_cast<CompSBasePlugin*>(replaced->getPlugin(getPrefix()));
if (replacedplug==NULL) {
//assert(false); //Not sure when this situation would come up, so I would like to see an example.
//(sadly, I cannot set an assert in production code.)
return LIBSBML_OPERATION_SUCCESS; //I guess? LS DEBUG
}
for (unsigned int re=0; re<replacedplug->getNumReplacedElements(); re++) {
ret = replacedplug->getReplacedElement(re)->replaceWithAndMaybeDelete(replacement, true, conversionFactor);
if (ret != LIBSBML_OPERATION_SUCCESS) return ret;
}
if (replacedplug->isSetReplacedBy()) {
ret = replacedplug->getReplacedBy()->replaceWithAndMaybeDelete(replacement, deleteme, conversionFactor);
if (ret != LIBSBML_OPERATION_SUCCESS) return ret;
}
return LIBSBML_OPERATION_SUCCESS;
}
int SBaseRef::collectDeletions(set<SBase*>* removed, set<SBase*>* toremove)
{
SBase* todelete = getReferencedElement();
if (todelete==NULL) {
return LIBSBML_INVALID_OBJECT;
}
if (removed) {
if (removed->find(todelete) != removed->end()) {
//Already deleted or replaced.
return LIBSBML_OPERATION_SUCCESS;
}
}
if (toremove) {
toremove->insert(todelete);
}
CompSBasePlugin* todplug = static_cast<CompSBasePlugin*>(todelete->getPlugin(getPrefix()));
if (todplug != NULL) {
for (unsigned long re=0; re<todplug->getNumReplacedElements(); re++) {
todplug->getReplacedElement(re)->collectDeletions(removed, toremove);
}
if (todplug->isSetReplacedBy()) {
todplug->getReplacedBy()->collectDeletions(removed, toremove);
}
}
return LIBSBML_OPERATION_SUCCESS;
}
CompSBasePlugin::CompSBasePlugin(const CompSBasePlugin& orig)
: SBasePlugin(orig)
, mListOfReplacedElements(NULL)
, mReplacedBy(NULL)
{
if (orig.isSetReplacedBy()) {
mReplacedBy = orig.mReplacedBy->clone();
mReplacedBy->connectToParent(getParentSBMLObject());
}
if (orig.getNumReplacedElements() > 0) {
createListOfReplacedElements();
for (unsigned int re=0; re<orig.getNumReplacedElements(); re++) {
mListOfReplacedElements->append(orig.getReplacedElement(re));
}
}
connectToChild();
}
Model* CompModelPlugin::flattenModel() const
{
//First make a copy of our parent (the model to be flattened):
const Model* parent = static_cast<const Model*>(getParentSBMLObject());
if (parent==NULL) {
return NULL;
}
//doc needs to be non-const so that the error messages can be updated. Otherwise, nothing changes.
SBMLDocument* doc = const_cast<SBMLDocument*>(getSBMLDocument());
if (doc==NULL) {
return NULL;
}
//Set the original document so that it can find the model definitions
//and external model definitions while we flatten.
Model* flat = parent->clone();
flat->setSBMLDocument(doc);
CompModelPlugin* flatplug =
static_cast<CompModelPlugin*>(flat->getPlugin(getPrefix()));
// Now instantiate its submodels and
// follow all renaming/deletion/replacement rules.
vector<const Model*> submods;
int success = flatplug->instantiateSubmodels();
if (success != LIBSBML_OPERATION_SUCCESS) {
//instantiateSubmodels sets its own error messages.
delete flat;
return NULL;
}
//Now start the aggregation process.
//This goes from the bottom up, calling 'appendFrom' iteratively
//(from the plugin).
for (unsigned int sm=0; sm<flatplug->getNumSubmodels(); sm++)
{
Model* submodel = flatplug->getSubmodel(sm)->getInstantiation();
if (submodel==NULL) {
//getInstantiation should be calling a cached value by now, but if not, it will set its own error messages.
delete flat;
return NULL;
}
CompModelPlugin* submodplug =
static_cast<CompModelPlugin*>(submodel->getPlugin(getPrefix()));
if (submodplug != NULL) {
//Strip the ports from the submodel, as we no longer need them.
while (submodplug->getNumPorts() > 0)
{
delete submodplug->removePort(0);
}
}
success = flat->appendFrom(submodel);
if (success != LIBSBML_OPERATION_SUCCESS) {
string error = "Unable to flatten model in CompModelPlugin::flattenModel: appending elements from the submodel '" + submodel->getId() + "' to the elements of the parent model failed.";
doc->getErrorLog()->logPackageError("comp", CompModelFlatteningFailed, getPackageVersion(), getLevel(), getVersion(), error);
delete flat;
return NULL;
}
}
// Now we clear the saved referenced elements in the local Port objects,
// but point them to the new object if necessary.
flatplug->resetPorts();
// Next, strip the package info from 'flat'.
// We're going to remove everything but the Ports:
flatplug->mListOfSubmodels.clear();
flatplug->clearReplacedElements();
flatplug->unsetReplacedBy();
List* allelements = flat->getAllElements();
vector<SBase*> nonReplacedElements;
for (unsigned int el=0; el<allelements->getSize(); el++)
{
SBase* element = static_cast<SBase*>(allelements->get(el));
int type = element->getTypeCode();
if (!(type==SBML_COMP_REPLACEDBY ||
type==SBML_COMP_REPLACEDELEMENT ||
type==SBML_COMP_SBASEREF))
{
nonReplacedElements.push_back(element);
}
}
// delete the list
delete allelements;
for (unsigned int el=0; el<nonReplacedElements.size(); el++)
{
SBase* element = nonReplacedElements[el];
CompSBasePlugin* elplug =
static_cast<CompSBasePlugin*>(element->getPlugin(getPrefix()));
if (elplug != NULL)
{
elplug->clearReplacedElements();
elplug->unsetReplacedBy();
}
}
//Finally, unset the document again.
flat->setSBMLDocument(NULL);
return flat;
}
int main(int argc,char** argv)
{
int retval = 0;
SBMLNamespaces sbmlns(3,1,"comp",1);
// create the document
SBMLDocument *document = new SBMLDocument(&sbmlns);
//Define the external model definition
CompSBMLDocumentPlugin* compdoc
= static_cast<CompSBMLDocumentPlugin*>(document->getPlugin("comp"));
compdoc->setRequired(true);
ExternalModelDefinition* extmod = compdoc->createExternalModelDefinition();
extmod->setId("ExtMod1");
extmod->setSource("enzyme_model.xml");
extmod->setModelRef("enzyme");
//Define the 'simple' model
ModelDefinition* mod1 = compdoc->createModelDefinition();
mod1->setId("simple");
Compartment* comp=mod1->createCompartment();
comp->setSpatialDimensions((unsigned int)3);
comp->setConstant(true);
comp->setId("comp");
comp->setSize(1L);
// We have to construct it this way because we get the comp
// plugin from it later.
Species spec(&sbmlns);
spec.setCompartment("comp");
spec.setHasOnlySubstanceUnits(false);
spec.setConstant(false);
spec.setBoundaryCondition(false);
spec.setId("S");
spec.setInitialConcentration(5);
mod1->addSpecies(&spec);
spec.setId("D");
spec.setInitialConcentration(10);
mod1->addSpecies(&spec);
Reaction rxn(&sbmlns);
rxn.setReversible(true);
rxn.setFast(false);
rxn.setId("J0");
SpeciesReference sr(&sbmlns);
sr.setConstant(true);
sr.setStoichiometry(1);
sr.setSpecies("S");
rxn.addReactant(&sr);
sr.setSpecies("D");
rxn.addProduct(&sr);
mod1->addReaction(&rxn);
CompModelPlugin* mod1plug
= static_cast<CompModelPlugin*>(mod1->getPlugin("comp"));
Port port;
port.setId("S_port");
port.setIdRef("S");
mod1plug->addPort(&port);
Port* port2 = mod1plug->createPort();
port2->setId("D_port");
port2->setIdRef("D");
port.setId("comp_port");
port.setIdRef("comp");
mod1plug->addPort(&port);
port.setId("J0_port");
port.setIdRef("J0");
mod1plug->addPort(&port);
// create the Model
Model* model=document->createModel();
model->setId("complexified");
// Set the submodels
CompModelPlugin* mplugin
= static_cast<CompModelPlugin*>(model->getPlugin("comp"));
Submodel* submod1 = mplugin->createSubmodel();
submod1->setId("A");
submod1->setModelRef("ExtMod1");
Submodel* submod2 = mplugin->createSubmodel();
submod2->setId("B");
submod2->setModelRef("simple");
Deletion* del = submod2->createDeletion();
del->setPortRef("J0_port");
// Synchronize the compartments
Compartment* mcomp=model->createCompartment();
mcomp->setSpatialDimensions((unsigned int)3);
mcomp->setConstant(true);
mcomp->setId("comp");
mcomp->setSize(1L);
CompSBasePlugin* compartplug
= static_cast<CompSBasePlugin*>(mcomp->getPlugin("comp"));
ReplacedElement re;
re.setIdRef("comp");
re.setSubmodelRef("A");
compartplug->addReplacedElement(&re);
re.setSubmodelRef("B");
re.unsetIdRef();
re.setPortRef("comp_port");
compartplug->addReplacedElement(&re);
//Synchronize the species
spec.setId("S");
spec.unsetInitialConcentration();
CompSBasePlugin* specplug
= static_cast<CompSBasePlugin*>(spec.getPlugin("comp"));
ReplacedElement* sre = specplug->createReplacedElement();
sre->setSubmodelRef("A");
sre->setIdRef("S");
ReplacedBy* srb = specplug->createReplacedBy();
srb->setSubmodelRef("B");
srb->setPortRef("S_port");
model->addSpecies(&spec);
spec.setId("D");
sre->setIdRef("D");
srb->setPortRef("D_port");
model->addSpecies(&spec);
writeSBMLToFile(document,"eg-ports.xml");
writeSBMLToFile(document,"spec_example3.xml");
delete document;
document = readSBMLFromFile("spec_example3.xml");
if (document == NULL)
{
cout << "Error reading back in file." << endl;
retval = -1;
}
else
{
document->setConsistencyChecks(LIBSBML_CAT_UNITS_CONSISTENCY, false);
document->checkConsistency();
if (document->getErrorLog()->getNumFailsWithSeverity(2) > 0
|| document->getErrorLog()->getNumFailsWithSeverity(3) > 0)
{
stringstream errorstream;
document->printErrors(errorstream);
cout << "Errors encoutered when round-tripping SBML file: \n"
<< errorstream.str() << endl;
retval = -1;
}
writeSBMLToFile(document, "spec_example3_rt.xml");
delete document;
}
#ifdef WIN32
if (retval != 0)
{
cout << "(Press any key to exit.)" << endl;
_getch();
}
#endif
return retval;
}
int main(int argc,char** argv)
{
int retval = 0;
SBMLNamespaces sbmlns(3,1,"comp",1);
// create the document
SBMLDocument *document = new SBMLDocument(&sbmlns);
//Define the external model definitions
CompSBMLDocumentPlugin* compdoc
= static_cast<CompSBMLDocumentPlugin*>(document->getPlugin("comp"));
compdoc->setRequired(true);
ExternalModelDefinition* extmod = compdoc->createExternalModelDefinition();
extmod->setId("ExtMod1");
extmod->setSource("enzyme_model.xml");
extmod->setModelRef("enzyme");
// create the main Model
Model* model=document->createModel();
// Set the submodels
CompModelPlugin* mplugin = static_cast<CompModelPlugin*>(model->getPlugin("comp"));
Submodel* submod1 = mplugin->createSubmodel();
submod1->setId("A");
submod1->setModelRef("ExtMod1");
Submodel* submod2 = mplugin->createSubmodel();
submod2->setId("B");
submod2->setModelRef("ExtMod1");
// create a replacement compartment
Compartment* comp=model->createCompartment();
comp->setSpatialDimensions((unsigned int)3);
comp->setConstant(true);
comp->setId("comp");
comp->setSize(1L);
//Tell the model that this compartment replaces both of the inside ones.
CompSBasePlugin* compartplug = static_cast<CompSBasePlugin*>(comp->getPlugin("comp"));
ReplacedElement re;
re.setIdRef("comp");
re.setSubmodelRef("A");
compartplug->addReplacedElement(&re);
re.setSubmodelRef("B");
compartplug->addReplacedElement(&re);
// create a replacement species
Species* spec = model->createSpecies();
spec->setCompartment("comp");
spec->setHasOnlySubstanceUnits(false);
spec->setConstant(false);
spec->setBoundaryCondition(false);
spec->setId("S");
//Tell the model that this species replaces both of the inside ones.
CompSBasePlugin* spp = static_cast<CompSBasePlugin*>(spec->getPlugin("comp"));
re.setIdRef("S");
re.setSubmodelRef("A");
spp->addReplacedElement(&re);
re.setSubmodelRef("B");
spp->addReplacedElement(&re);
writeSBMLToFile(document,"eg-import-external.xml");
writeSBMLToFile(document,"spec_example2.xml");
delete document;
document = readSBMLFromFile("spec_example2.xml");
if (document == NULL)
{
cout << "Error reading back in file." << endl;
retval = -1;
}
else
{
document->setConsistencyChecks(LIBSBML_CAT_UNITS_CONSISTENCY, false);
document->checkConsistency();
if (document->getErrorLog()->getNumFailsWithSeverity(2) > 0
|| document->getErrorLog()->getNumFailsWithSeverity(3) > 0)
{
stringstream errorstream;
document->printErrors(errorstream);
cout << "Errors encoutered when round-tripping SBML file: \n"
<< errorstream.str() << endl;
retval = -1;
}
writeSBMLToFile(document, "spec_example2_rt.xml");
delete document;
}
#ifdef WIN32
if (retval != 0)
{
cout << "(Press any key to exit.)" << endl;
_getch();
}
#endif
return retval;
}
int ReplacedElement::performReplacementAndCollect(set<SBase*>* removed, set<SBase*>* toremove)
{
SBMLDocument* doc = getSBMLDocument();
if (isSetDeletion()) {
//Deletions don't need to be replaced.
return LIBSBML_OPERATION_SUCCESS;
}
//Find the various objects and plugin objects we need for this to work.
SBase* lore = getParentSBMLObject();
ListOf* lorelist = static_cast<ListOf*>(lore);
if (lore == NULL) {
if (doc) {
string error = "Cannot carry out replacement in ReplacedElement::performReplacement: no parent <listOfReplacedElements> could be found for the given replacement element.";
doc->getErrorLog()->logPackageError("comp", CompModelFlatteningFailed, getPackageVersion(), getLevel(), getVersion(), error, getLine(), getColumn());
}
return LIBSBML_INVALID_OBJECT;
}
if (lore->getTypeCode() != SBML_LIST_OF || lorelist->getItemTypeCode() != SBML_COMP_REPLACEDELEMENT) {
if (doc) {
string error = "Cannot carry out replacement in ReplacedElement::performReplacement: no parent <listOfReplacedElements> could be found for the given replacement element.";
doc->getErrorLog()->logPackageError("comp", CompModelFlatteningFailed, getPackageVersion(), getLevel(), getVersion(), error, getLine(), getColumn());
}
return LIBSBML_INVALID_OBJECT;
}
SBase* parent = lore->getParentSBMLObject();
if (parent==NULL) {
if (doc) {
string error = "Cannot carry out replacement in ReplacedElement::performReplacement: no parent could be found for the parent <listOfReplacedElements> object.";
doc->getErrorLog()->logPackageError("comp", CompModelFlatteningFailed, getPackageVersion(), getLevel(), getVersion(), error, getLine(), getColumn());
}
return LIBSBML_INVALID_OBJECT;
}
SBase* ref = getReferencedElement();
if (ref==NULL) {
//getReferencedElement sets its own error messages.
return LIBSBML_INVALID_OBJECT;
}
if (removed && removed->find(ref)!=removed->end()) {
//Already deleted: can't get the deleted element's ID to
if (doc) {
string error = "Cannot carry out replacement in ReplacedElement::performReplacement: a <" + parent->getElementName() + ">";
if (parent->isSetId()) {
error += " with the ID '" + parent->getId() + "'";
}
error += " has a child <replacedElement> that points to something that has already been deleted, probably because its parent was deleted.";
doc->getErrorLog()->logPackageError("comp", CompDeletedReplacement, getPackageVersion(), getLevel(), getVersion(), error, getLine(), getColumn());
}
return LIBSBML_INVALID_OBJECT;
}
//Update the IDs.
int ret = updateIDs(ref, parent);
if (ret != LIBSBML_OPERATION_SUCCESS) {
return ret;
}
//Perform any conversions on references in the submodel.
ASTNode* blank = NULL;
ret = performConversions(parent, &blank);
if (ret != LIBSBML_OPERATION_SUCCESS)
{
if (blank != NULL)
{
delete blank;
}
return ret;
}
CompSBasePlugin* refplug = static_cast<CompSBasePlugin*>(ref->getPlugin(getPrefix()));
if (refplug != NULL) {
//Now recurse down the 'replace*' tree, renaming IDs and deleting things as we go.
for (unsigned int re=0; re<refplug->getNumReplacedElements(); re++) {
refplug->getReplacedElement(re)->replaceWithAndMaybeDelete(parent, true, blank);
if (toremove) {
toremove->insert(refplug->getReplacedElement(re)->getReferencedElement());
}
}
if (refplug->isSetReplacedBy()) {
//Even if the subelement used to be replaced by something further down, it is now being replaced by the parent. It just can't catch a break, it seems.
refplug->getReplacedBy()->replaceWithAndMaybeDelete(parent, true, blank);
if (toremove) {
toremove->insert(refplug->getReplacedBy()->getReferencedElement());
}
}
}
if (toremove) {
toremove->insert(ref);
}
if (blank != NULL)
{
delete blank;
}
return LIBSBML_OPERATION_SUCCESS;
}
LIBSBML_CPP_NAMESPACE_USE
int main(int argc,char** argv)
{
DynPkgNamespaces sbmlns;
sbmlns.addPackageNamespace("comp", 1, "comp");
// create the document
SBMLDocument *document = new SBMLDocument(&sbmlns);
document->setPackageRequired("dyn", true);
document->setPackageRequired("comp", true);
// create the Model
Model* model=document->createModel();
model->setId("grid2x2");
// create the Compartment
Compartment* compartment = model->createCompartment();
compartment->setId("Loc1");
compartment->setConstant(false);
compartment->setSize(1);
compartment->setSpatialDimensions(2.0);
DynCompartmentPlugin* cplugin =
static_cast<DynCompartmentPlugin*>(compartment->getPlugin("dyn"));
SpatialComponent* component = cplugin->createSpatialComponent();
component->setSpatialIndex(DYN_SPATIALKIND_CARTESIANX);
component->setVariable("q1_X");
component = cplugin->createSpatialComponent();
component->setSpatialIndex(DYN_SPATIALKIND_CARTESIANY);
component->setVariable("q1_Y");
CompSBasePlugin* compPlugin =
static_cast<CompSBasePlugin*>(compartment->getPlugin("comp"));
ReplacedElement* relement = compPlugin->createReplacedElement();
relement->setIdRef("C");
relement->setSubmodelRef("GRID_1_1_cell");
compartment = model->createCompartment();
compartment->setId("Loc2");
compartment->setConstant(false);
compartment->setSize(1);
compartment->setSpatialDimensions(2.0);
cplugin =
static_cast<DynCompartmentPlugin*>(compartment->getPlugin("dyn"));
component = cplugin->createSpatialComponent();
component->setSpatialIndex(DYN_SPATIALKIND_CARTESIANX);
component->setVariable("q2_X");
component = cplugin->createSpatialComponent();
component->setSpatialIndex(DYN_SPATIALKIND_CARTESIANY);
component->setVariable("q2_Y");
compPlugin =
static_cast<CompSBasePlugin*>(compartment->getPlugin("comp"));
relement = compPlugin->createReplacedElement();
relement->setIdRef("C");
relement->setSubmodelRef("GRID_1_2_cell");
compartment = model->createCompartment();
compartment->setId("Loc3");
compartment->setConstant(false);
compartment->setSize(1);
compartment->setSpatialDimensions(2.0);
cplugin =
static_cast<DynCompartmentPlugin*>(compartment->getPlugin("dyn"));
component = cplugin->createSpatialComponent();
component->setSpatialIndex(DYN_SPATIALKIND_CARTESIANX);
component->setVariable("q3_X");
component = cplugin->createSpatialComponent();
component->setSpatialIndex(DYN_SPATIALKIND_CARTESIANY);
component->setVariable("q3_Y");
compPlugin =
static_cast<CompSBasePlugin*>(compartment->getPlugin("comp"));
relement = compPlugin->createReplacedElement();
relement->setIdRef("C");
relement->setSubmodelRef("GRID_2_1_cell");
compartment = model->createCompartment();
compartment->setId("Loc4");
compartment->setConstant(false);
compartment->setSize(1);
compartment->setSpatialDimensions(2.0);
cplugin =
static_cast<DynCompartmentPlugin*>(compartment->getPlugin("dyn"));
component = cplugin->createSpatialComponent();
component->setSpatialIndex(DYN_SPATIALKIND_CARTESIANX);
component->setVariable("q4_X");
component = cplugin->createSpatialComponent();
component->setSpatialIndex(DYN_SPATIALKIND_CARTESIANY);
component->setVariable("q4_Y");
compPlugin =
static_cast<CompSBasePlugin*>(compartment->getPlugin("comp"));
relement = compPlugin->createReplacedElement();
relement->setIdRef("C");
relement->setSubmodelRef("GRID_2_2_cell");
// create Parameters
Parameter* param = model->createParameter();
param->initDefaults();
param->setId("q1_X");
param->setValue(1);
param = model->createParameter();
param->initDefaults();
param->setId("q1_Y");
param->setValue(1);
param = model->createParameter();
param->initDefaults();
param->setId("q2_X");
param->setValue(2);
param = model->createParameter();
param->initDefaults();
param->setId("q2_Y");
param->setValue(1);
param = model->createParameter();
param->initDefaults();
param->setId("q3_X");
param->setValue(1);
param = model->createParameter();
param->initDefaults();
param->setId("q3_Y");
param->setValue(2);
param = model->createParameter();
param->initDefaults();
param->setId("q4_X");
param->setValue(2);
param = model->createParameter();
param->initDefaults();
param->setId("q4_Y");
param->setValue(2);
// create SubModels
CompModelPlugin* mplugin =
static_cast<CompModelPlugin*>(model->getPlugin("comp"));
Submodel* submodel = mplugin->createSubmodel();
submodel->setId("GRID_1_1_cell");
submodel->setModelRef("Cell");
submodel = mplugin->createSubmodel();
submodel->setId("GRID_1_2_cell");
submodel->setModelRef("Cell");
submodel = mplugin->createSubmodel();
submodel->setId("GRID_2_1_cell");
submodel->setModelRef("Cell");
submodel = mplugin->createSubmodel();
submodel->setId("GRID_2_2_cell");
submodel->setModelRef("Cell");
// create the ModelDefinition
CompSBMLDocumentPlugin* dplugin =
static_cast<CompSBMLDocumentPlugin*>(document->getPlugin("comp"));
ModelDefinition* mdef = dplugin->createModelDefinition();
mdef->setId("Cell");
compartment = mdef->createCompartment();
compartment->initDefaults();
compartment->setId("C");
compartment->setSpatialDimensions(2.0);
compartment->setSize(1.0);
Species* species = mdef->createSpecies();
species->setId("R");
species->setCompartment("C");
species->setHasOnlySubstanceUnits(false);
species->setBoundaryCondition(false);
species->setConstant(false);
species = mdef->createSpecies();
species->setId("S");
species->setCompartment("C");
species->setHasOnlySubstanceUnits(false);
species->setBoundaryCondition(false);
species->setConstant(false);
Reaction* reaction = mdef->createReaction();
reaction->setId("Degradation_R");
reaction->setReversible(false);
reaction->setFast(false);
reaction->setCompartment("C");
SpeciesReference* reactant = reaction->createReactant();
reactant->setSpecies("R");
reactant->setStoichiometry(1);
reactant->setConstant(true);
reaction = mdef->createReaction();
reaction->setId("Degradation_S");
reaction->setReversible(false);
reaction->setFast(false);
reaction->setCompartment("C");
reactant = reaction->createReactant();
reactant->setSpecies("S");
reactant->setStoichiometry(1);
reactant->setConstant(true);
document->checkConsistency();
Event* event = mdef->createEvent();
event->setId("event0");
event->setUseValuesFromTriggerTime(false);
DynEventPlugin* eplugin =
static_cast<DynEventPlugin*>(event->getPlugin("dyn"));
eplugin->setApplyToAll(true);
eplugin->setCboTerm("http://cbo.biocomplexity.indiana.edu/svn/cbo/trunk/CBO_1_0.owl#CellDivision");
Trigger* trigger = event->createTrigger();
trigger->setInitialValue(false);
trigger->setPersistent(false);
trigger->setMath(SBML_parseFormula("true"));
if (document->getNumErrors(LIBSBML_SEV_ERROR) > 0)
document->printErrors();
writeSBML(document,"dyn_example2.xml");
delete document;
}
END_TEST
START_TEST (test_comp_sbase)
{
SBMLNamespaces sbmlns(3,1,"comp",1);
Parameter param(&sbmlns);
CompSBasePlugin* pplugin = static_cast<CompSBasePlugin*>(param.getPlugin("comp"));
fail_unless(pplugin->getNumReplacedElements()==0);
ReplacedElement* re = pplugin->createReplacedElement();
re->setMetaId("re1");
fail_unless(pplugin->getNumReplacedElements()==1);
fail_unless(pplugin->addReplacedElement(NULL)==LIBSBML_INVALID_OBJECT);
ReplacedElement re2(3, 1);
re2.setMetaId("re2");
fail_unless(pplugin->addReplacedElement(&re2)==LIBSBML_INVALID_OBJECT);
re2.setDeletion("ID1");
fail_unless(pplugin->addReplacedElement(&re2)==LIBSBML_INVALID_OBJECT);
re2.setSubmodelRef("mod1");
fail_unless(pplugin->addReplacedElement(&re2)==LIBSBML_OPERATION_SUCCESS);
ReplacedElement* reref = pplugin->getReplacedElement(1);
fail_unless(reref != NULL);
fail_unless(reref->getMetaId()=="re2");
reref = pplugin->getReplacedElement(0);
fail_unless(reref != NULL);
fail_unless(reref->getMetaId()=="re1");
re->setDeletion("ID1");
fail_unless(reref->getDeletion()=="ID1");
fail_unless(pplugin->removeReplacedElement(3)==NULL);
fail_unless(pplugin->removeReplacedElement(0)==re);
fail_unless(pplugin->getReplacedElement(1)==NULL);
fail_unless(pplugin->isSetReplacedBy()==false);
ReplacedBy* rb = pplugin->createReplacedBy();
fail_unless(rb != NULL);
fail_unless(pplugin->isSetReplacedBy()==true);
fail_unless(pplugin->setReplacedBy(NULL)==LIBSBML_OPERATION_SUCCESS);
fail_unless(pplugin->isSetReplacedBy()==false);
ReplacedBy rb2(3,1);
fail_unless(pplugin->setReplacedBy(&rb2)==LIBSBML_INVALID_OBJECT);
rb2.setIdRef("ID1");
fail_unless(pplugin->setReplacedBy(&rb2)==LIBSBML_INVALID_OBJECT);
rb2.setSubmodelRef("mod1");
fail_unless(pplugin->setReplacedBy(&rb2)==LIBSBML_OPERATION_SUCCESS);
fail_unless(pplugin->unsetReplacedBy()==LIBSBML_OPERATION_SUCCESS);
}