libSBML C++ API  5.18.0
spec_example1.cpp

Example #1 from the SBML Level 3 Hierarchical Model Composition specification.

/**
* @file spec_example1.cpp
* @brief SBML hierarchical composition example
* @author Lucian Smith
*
* <!--------------------------------------------------------------------------
* This sample program is distributed under a different license than the rest
* of libSBML. This program uses the open-source MIT license, as follows:
*
* Copyright (c) 2013-2018 by the California Institute of Technology
* (California, USA), the European Bioinformatics Institute (EMBL-EBI, UK)
* and the University of Heidelberg (Germany), with support from the National
* Institutes of Health (USA) under grant R01GM070923. All rights reserved.
*
* Permission is hereby granted, free of charge, to any person obtaining a
* copy of this software and associated documentation files (the "Software"),
* to deal in the Software without restriction, including without limitation
* the rights to use, copy, modify, merge, publish, distribute, sublicense,
* and/or sell copies of the Software, and to permit persons to whom the
* Software is furnished to do so, subject to the following conditions:
*
* The above copyright notice and this permission notice shall be included in
* all copies or substantial portions of the Software.
*
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
* IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
* THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
* FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
* DEALINGS IN THE SOFTWARE.
*
* Neither the name of the California Institute of Technology (Caltech), nor
* of the European Bioinformatics Institute (EMBL-EBI), nor of the University
* of Heidelberg, nor the names of any contributors, may be used to endorse
* or promote products derived from this software without specific prior
* written permission.
* ------------------------------------------------------------------------ -->
*/
#include <iostream>
#include <sstream>
#include <sbml/SBMLTypes.h>
#ifdef WIN32
#include <conio.h>
#endif
using namespace std;
int main(int argc,char** argv)
{
int retval = 0;
SBMLNamespaces sbmlns(3,1,"comp",1);
// create the document
SBMLDocument *document = new SBMLDocument(&sbmlns);
//Create our submodel
= static_cast<CompSBMLDocumentPlugin*>(document->getPlugin("comp"));
compdoc->setRequired(true);
mod1->setId("enzyme");
mod1->setName("enzyme");
comp->setSpatialDimensions((unsigned int)3);
comp->setConstant(true);
comp->setId("comp");
comp->setSize(1L);
Species spec(&sbmlns);
spec.setCompartment("comp");
spec.setConstant(false);
spec.setBoundaryCondition(false);
spec.setId("S");
mod1->addSpecies(&spec);
spec.setId("E");
mod1->addSpecies(&spec);
spec.setId("D");
mod1->addSpecies(&spec);
spec.setId("ES");
mod1->addSpecies(&spec);
Reaction rxn(&sbmlns);
rxn.setReversible(true);
rxn.setFast(false);
Reaction rxn2(rxn);
rxn.setId("J0");
rxn2.setId("J1");
SpeciesReference sr(&sbmlns);
sr.setConstant(true);
sr.setSpecies("S");
rxn.addReactant(&sr);
sr.setSpecies("E");
rxn.addReactant(&sr);
rxn2.addProduct(&sr);
sr.setSpecies("ES");
rxn.addProduct(&sr);
rxn2.addReactant(&sr);
sr.setSpecies("D");
rxn2.addProduct(&sr);
mod1->addReaction(&rxn);
mod1->addReaction(&rxn2);
// create the Model
Model* model=document->createModel();
model->setId("aggregate");
// Create a submodel
CompModelPlugin* mplugin = static_cast<CompModelPlugin*>(model->getPlugin("comp"));
Submodel* submod1 = mplugin->createSubmodel();
submod1->setId("submod1");
submod1->setModelRef("enzyme");
Submodel submod2;
submod2.setId("submod2");
submod2.setModelRef("enzyme");
mplugin->addSubmodel(&submod2);
writeSBMLToFile(document,"eg-simple-aggregate.xml");
writeSBMLToFile(document,"enzyme_model.xml");
delete document;
document = readSBMLFromFile("enzyme_model.xml");
if (document == NULL)
{
cout << "Error reading back in file." << endl;
retval = -1;
}
else
{
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, "enzyme_model_rt.xml");
delete document;
}
#ifdef WIN32
if (retval != 0)
{
cout << "(Press any key to exit.)" << endl;
_getch();
}
#endif
return retval;
}