/*
* Copy constructor.
*/
CoordinateComponent::CoordinateComponent(const CoordinateComponent& source)
: SBase(source)
, mBoundaryMin(NULL)
, mBoundaryMax(NULL)
{
this->mSpatialId=source.mSpatialId;
this->mComponentType=source.mComponentType;
this->mSbmlUnit=source.mSbmlUnit;
this->mIndex=source.mIndex;
this->mIsSetIndex=source.mIsSetIndex;
if (source.isSetBoundaryMin())
this->mBoundaryMin = source.mBoundaryMin->clone();
if (source.isSetBoundaryMax())
this->mBoundaryMax = source.mBoundaryMax->clone();
connectToChild();
}
void readSpatialSBML() {
SBMLDocument *document2 = readSBML("spatial_example2.xml");
Model *model2 = document2->getModel();
Compartment *comp;
SpatialCompartmentPlugin* cplugin;
for (unsigned int i = 0; i < model2->getNumCompartments(); i++) {
comp = model2->getCompartment(i);
cout << "Compartment" << i << ": " << comp->getId() << endl;
cplugin = static_cast<SpatialCompartmentPlugin*>(comp->getPlugin("spatial"));
if (cplugin->getCompartmentMapping()->isSetId()) {
cout << "Comp" << i << " CMSpId: " << cplugin->getCompartmentMapping()->getId() << endl;
cout << "Comp" << i << " CM_DType: " << cplugin->getCompartmentMapping()->getDomainType() << endl;
cout << "Comp" << i << " CM_UnitSz: " << cplugin->getCompartmentMapping()->getUnitSize() << endl;
}
}
Species *sp;
SpatialSpeciesPlugin* srplugin;
for (unsigned int i = 0; i < model2->getNumSpecies(); i++) {
sp = model2->getSpecies(i);
cout << "Species" << i << ": " << sp->getId() << endl;
srplugin = static_cast<SpatialSpeciesPlugin*>(sp->getPlugin("spatial"));
if (srplugin->getIsSpatial()) {
cout << "species" << i << " isSpatial: " << srplugin->getIsSpatial() << endl;
}
}
Parameter *param;
SpatialParameterPlugin* pplugin;
for (unsigned int i = 0; i < model2->getNumParameters(); i++) {
param = model2->getParameter(i);
cout << "Parameter" << i << ": " << param->getId() << endl;
pplugin = static_cast<SpatialParameterPlugin*>(param->getPlugin("spatial"));
if (pplugin->isSetSpatialSymbolReference()) {
cout << "Parameter" << i << " SpRefId: " << pplugin->getSpatialSymbolReference()->getSpatialRef() << endl;
}
if (pplugin->isSetDiffusionCoefficient()) {
cout << "Diff_" << i << " SpeciesVarId: " << pplugin->getDiffusionCoefficient()->getVariable() << endl;
cout << "Diff_" << i << " Type: " << DiffusionKind_toString(pplugin->getDiffusionCoefficient()->getType()) << endl;
for (unsigned int j = 0; j < pplugin->getDiffusionCoefficient()->getNumCoordinateReferences(); ++j)
cout << "Diff_" << i << " SpCoordIndex " << j << " : " << CoordinateKind_toString(pplugin->getDiffusionCoefficient()->getCoordinateReference(j) ->getCoordinate()) << endl;
}
if (pplugin->isSetAdvectionCoefficient()) {
cout << "Adv_" << i << " SpeciesVarId: " << pplugin->getAdvectionCoefficient()->getVariable() << endl;
cout << "Adv_" << i << " SpCoordIndex: " << CoordinateKind_toString(pplugin->getAdvectionCoefficient()->getCoordinate()) << endl;
}
if (pplugin->isSetBoundaryCondition()) {
cout << "BC_" << i << " SpeciesVarId: " << pplugin->getBoundaryCondition()->getVariable() << endl;
cout << "BC_" << i << " SpCoordBoundary: " << pplugin->getBoundaryCondition()->getCoordinateBoundary() << endl;
cout << "BC_" << i << " SpBoundaryType: " << pplugin->getBoundaryCondition()->getType() << endl;
}
}
Reaction *rxn;
SpatialReactionPlugin* rplugin;
for (unsigned int i = 0; i < model2->getNumReactions(); i++) {
rxn = model2->getReaction(i);
cout << "Reaction" << i << ": " << rxn->getId() << endl;
rplugin = static_cast<SpatialReactionPlugin*>(rxn->getPlugin("spatial"));
if (rplugin->getIsLocal()) {
cout << "rxn" << i << " isLocal: " << rplugin->getIsLocal() << endl;
}
}
Rule *rule;
for (unsigned int i = 0; i < model2->getNumRules(); i++) {
rule = model2->getRule(i);
cout << "Rule" << i << ": " << rule->getVariable() << endl;
}
//
// Get a SpatialModelPlugin object plugged in the model object.
//
// The type of the returned value of SBase::getPlugin() function is
// SBasePlugin*, and thus the value needs to be cast for the
// corresponding derived class.
//
SpatialModelPlugin* mplugin2;
mplugin2 = static_cast<SpatialModelPlugin*>(model2->getPlugin("spatial"));
cout << "URI: " << mplugin2->getURI() << endl;
cout << "prefix: " << mplugin2->getPrefix() << endl;
// get a Geometry object via SpatialModelPlugin object.
Geometry* geometry2 = mplugin2->getGeometry();
cout << "Geometry coordSystem: " << geometry2->getCoordinateSystem() << endl;
// get a CoordComponent object via the Geometry object.
CoordinateComponent* coordComp = geometry2->getCoordinateComponent(0);
std::cout << "CoordComponent Id: " << coordComp->getId() << std::endl;
std::cout << "CoordComponent type: " << CoordinateKind_toString( coordComp->getType()) << std::endl;
std::cout << "CoordComponent sbmlUnit: " << coordComp->getUnit() << std::endl;
if (coordComp->isSetBoundaryMin())
{
Boundary* minX = coordComp->getBoundaryMin();
std::cout << "minX name: " << minX->getId() << std::endl;
std::cout << "minX value: " << minX->getValue() << std::endl;
}
if (coordComp->isSetBoundaryMax())
{
Boundary* maxX = coordComp->getBoundaryMax();
std::cout << "maxX name: " << maxX->getId() << std::endl;
std::cout << "maxX value: " << maxX->getValue() << std::endl;
}
// get a DomainType object via the Geometry object.
DomainType* domainType2 = geometry2->getDomainType(0);
std::cout << "DomainType Id: " << domainType2->getId() << std::endl;
std::cout << "DomainType spatialDim: " << domainType2->getSpatialDimension() << std::endl;
// get a Domain object via the Geometry object.
Domain* domain = geometry2->getDomain(0);
std::cout << "Domain1 Id: " << domain->getId() << std::endl;
std::cout << "Domain1 domainType: " << domain->getDomainType() << std::endl;
// get an internal point via the domain object
InteriorPoint* internalPt = domain->getInteriorPoint(0);
std::cout << "InternalPt_1 coord1: " << internalPt->getCoord1() << std::endl;
// get a Domain object via the Geometry object.
domain = geometry2->getDomain(1);
std::cout << "Domain2 Id: " << domain->getId() << std::endl;
std::cout << "Domain2 domainType: " << domain->getDomainType() << std::endl;
// get an internal point via the domain object
internalPt = domain->getInteriorPoint(0);
std::cout << "InternalPt_2 coord1: " << internalPt->getCoord1() << std::endl;
// get an AdjacentDomains object via the Geometry object.
AdjacentDomains* adjDomain = geometry2->getAdjacentDomains(0);
std::cout << "AdjDomain Id: " << adjDomain->getId() << std::endl;
std::cout << "AdjDomain domain1: " << adjDomain->getDomain1() << std::endl;
std::cout << "AdjDomain domain2: " << adjDomain->getDomain2() << std::endl;
// get an AnalyticGeometry object via the Geometry object.
GeometryDefinition* gd;
for (unsigned int i = 0; i < geometry2->getNumGeometryDefinitions(); i++) {
gd = geometry2->getGeometryDefinition(i);
if (gd->isAnalyticGeometry()) {
AnalyticGeometry* analyticalGeom = static_cast<AnalyticGeometry*>(gd);
std::cout << "AnalGeom Id: " << analyticalGeom->getId() << std::endl;
// analVol from analGeom.
AnalyticVolume* av = analyticalGeom->getAnalyticVolume(0);
std::cout << "AnalVol Id: " << av->getId() << std::endl;
std::cout << "AnalVol domainType: " << av->getDomainType() << std::endl;
std::cout << "AnalVol funcType: " << av->getFunctionType() << std::endl;
std::cout << "AnalVol ordinal: " << av->getOrdinal() << std::endl;
const ASTNode* mathNode = av->getMath();
char* mathStr = writeMathMLToString(mathNode);
std::cout << "AnalVol math: " << mathStr << std::endl;
}
if (gd->isSampledFieldGeometry()) {
SampledFieldGeometry* sfGeom = static_cast<SampledFieldGeometry*>(gd);
std::cout << "SampledFieldGeom Id: " << sfGeom->getId() << std::endl;
// sampledField from sfGeom
SampledField* sf = sfGeom->getSampledField();
std::cout << "SampledField Id: " << sf->getId() << std::endl;
std::cout << "SampledField dataType: " << sf->getDataType() << std::endl;
std::cout << "SampledField interpolation: " << sf->getInterpolationType() << std::endl;
std::cout << "SampledField encoding: " << sf->getEncoding() << std::endl;
std::cout << "SampledField numSamples1: " << sf->getNumSamples1() << std::endl;
std::cout << "SampledField numSamples2: " << sf->getNumSamples2() << std::endl;
std::cout << "SampledField numSamples3: " << sf->getNumSamples3() << std::endl;
const ImageData* id = sf->getImageData();
int* samples = new int[id->getSamplesLength()];
id->getSamples(samples);
std::cout << "ImageData samples[0]: " << samples[0] << std::endl;
std::cout << "ImageData dtype: " << id->getDataType() << std::endl;
std::cout << "ImageData samplesLen: " << id->getSamplesLength() << std::endl;
// sampledVolVol from sfGeom.
SampledVolume* sv = sfGeom->getSampledVolume(0);
std::cout << "SampledVol Id: " << sv->getId() << std::endl;
std::cout << "SampledVol domainType: " << sv->getDomainType() << std::endl;
std::cout << "SampledVol sampledVal: " << sv->getSampledValue() << std::endl;
std::cout << "SampledVol min: " << sv->getMinValue() << std::endl;
std::cout << "SampledVol max: " << sv->getMaxValue() << std::endl;
}
}
delete document2;
}