private void loadFromXml(String fileName)
throws ParserConfigurationException, SAXException, IOException, ParseException {
System.out.println("NeuralNetwork : loading network topology from file " + fileName);
DocumentBuilderFactory factory = DocumentBuilderFactory.newInstance();
DocumentBuilder parser = factory.newDocumentBuilder();
Document doc = parser.parse(fileName);
Node nodeNeuralNetwork = doc.getDocumentElement();
if (!nodeNeuralNetwork.getNodeName().equals("neuralNetwork"))
throw new ParseException(
"[Error] NN-Load: Parse error in XML file, neural network couldn't be loaded.", 0);
// nodeNeuralNetwork ok
// indexNeuralNetworkContent -> indexStructureContent -> indexLayerContent -> indexNeuronContent
// -> indexNeuralInputContent
NodeList nodeNeuralNetworkContent = nodeNeuralNetwork.getChildNodes();
for (int innc = 0; innc < nodeNeuralNetworkContent.getLength(); innc++) {
Node nodeStructure = nodeNeuralNetworkContent.item(innc);
if (nodeStructure.getNodeName().equals("structure")) { // for structure element
NodeList nodeStructureContent = nodeStructure.getChildNodes();
for (int isc = 0; isc < nodeStructureContent.getLength(); isc++) {
Node nodeLayer = nodeStructureContent.item(isc);
if (nodeLayer.getNodeName().equals("layer")) { // for layer element
NeuralLayer neuralLayer = new NeuralLayer(this);
this.listLayers.add(neuralLayer);
NodeList nodeLayerContent = nodeLayer.getChildNodes();
for (int ilc = 0; ilc < nodeLayerContent.getLength(); ilc++) {
Node nodeNeuron = nodeLayerContent.item(ilc);
if (nodeNeuron.getNodeName().equals("neuron")) { // for neuron in layer
Neuron neuron =
new Neuron(
Double.parseDouble(((Element) nodeNeuron).getAttribute("threshold")),
neuralLayer);
neuralLayer.listNeurons.add(neuron);
NodeList nodeNeuronContent = nodeNeuron.getChildNodes();
for (int inc = 0; inc < nodeNeuronContent.getLength(); inc++) {
Node nodeNeuralInput = nodeNeuronContent.item(inc);
// if (nodeNeuralInput==null) System.out.print("-"); else System.out.print("*");
if (nodeNeuralInput.getNodeName().equals("input")) {
// System.out.println("neuron at
// STR:"+innc+" LAY:"+isc+" NEU:"+ilc+" INP:"+inc);
NeuralInput neuralInput =
new NeuralInput(
Double.parseDouble(((Element) nodeNeuralInput).getAttribute("weight")),
neuron);
neuron.listInputs.add(neuralInput);
}
}
}
}
}
}
}
}
}
/**
* Advanced users only; use loadXML() in PApplet.
*
* <p>Added extra code to handle \u2028 (Unicode NLF), which is sometimes inserted by web browsers
* (Safari?) and not distinguishable from a "real" LF (or CRLF) in some text editors (i.e.
* TextEdit on OS X). Only doing this for XML (and not all Reader objects) because LFs are
* essential. https://github.com/processing/processing/issues/2100
*
* @nowebref
*/
public XML(final Reader reader, String options)
throws IOException, ParserConfigurationException, SAXException {
DocumentBuilderFactory factory = DocumentBuilderFactory.newInstance();
// Prevent 503 errors from www.w3.org
try {
factory.setAttribute("http://apache.org/xml/features/nonvalidating/load-external-dtd", false);
} catch (IllegalArgumentException e) {
// ignore this; Android doesn't like it
}
// without a validating DTD, this doesn't do anything since it doesn't know what is ignorable
// factory.setIgnoringElementContentWhitespace(true);
factory.setExpandEntityReferences(false);
// factory.setExpandEntityReferences(true);
// factory.setCoalescing(true);
// builderFactory.setFeature("http://apache.org/xml/features/disallow-doctype-decl", true);
DocumentBuilder builder = factory.newDocumentBuilder();
// builder.setEntityResolver()
// SAXParserFactory spf = SAXParserFactory.newInstance();
// spf.setFeature("http://apache.org/xml/features/disallow-doctype-decl", true);
// SAXParser p = spf.newSAXParser();
// builder = DocumentBuilderFactory.newDocumentBuilder();
// builder = new SAXBuilder();
// builder.setValidation(validating);
Document document =
builder.parse(
new InputSource(
new Reader() {
@Override
public int read(char[] cbuf, int off, int len) throws IOException {
int count = reader.read(cbuf, off, len);
for (int i = 0; i < count; i++) {
if (cbuf[off + i] == '\u2028') {
cbuf[off + i] = '\n';
}
}
return count;
}
@Override
public void close() throws IOException {
reader.close();
}
}));
node = document.getDocumentElement();
}
public void read(Document doc) {
try {
Element root = doc.getDocumentElement();
NodeList nList = doc.getElementsByTagName("user");
for (int i = 0; i < nList.getLength(); i++) {
Node node = nList.item(i);
if (node.getNodeType() == Node.ELEMENT_NODE) {
Element element = (Element) node;
}
}
// Element e1 = (Element) nList.item(0);
} catch (Exception e) {
}
}
/**
* Unlike the loadXML() method in PApplet, this version works with files that are not in UTF-8
* format.
*
* @nowebref
*/
public XML(InputStream input, String options)
throws IOException, ParserConfigurationException, SAXException {
// this(PApplet.createReader(input), options); // won't handle non-UTF8
DocumentBuilderFactory factory = DocumentBuilderFactory.newInstance();
try {
// Prevent 503 errors from www.w3.org
factory.setAttribute("http://apache.org/xml/features/nonvalidating/load-external-dtd", false);
} catch (IllegalArgumentException e) {
// ignore this; Android doesn't like it
}
factory.setExpandEntityReferences(false);
DocumentBuilder builder = factory.newDocumentBuilder();
Document document = builder.parse(new InputSource(input));
node = document.getDocumentElement();
}
/**
* Method that reads a XML-file, and returns a Model that contains the information
*
* @param file
* @return
* @return
*/
public static Model readXML(String file) {
// initialize table to be filled with content of XML file
Model t = new Model();
try {
// Create file to be parsed by document parser
File xmlfile = new File(file);
// create parser
DocumentBuilder parser = DocumentBuilderFactory.newInstance().newDocumentBuilder();
// parse the file
Document parsedfile = parser.parse(xmlfile);
// normalize the parsed file (make it more user-friendly)
parsedfile.getDocumentElement().normalize();
NodeList cells = parsedfile.getElementsByTagName("CELL");
for (int i = 0; i < cells.getLength(); i++) {
// Get cell at list index i
Node currentcell = cells.item(i);
// read the elements "location" attributes row/column
if (Node.ELEMENT_NODE == currentcell.getNodeType()) {
Element cellinfo = (Element) currentcell;
// get the row number from node attribute
int row = Integer.parseInt(cellinfo.getAttribute("row")) - 1;
// get the column number from the node attribute
int col = Integer.parseInt(cellinfo.getAttribute("column")) - 1;
// get content from node
String content = cellinfo.getTextContent();
if (content != null) {
content = content.replace("\n", "");
}
// Make the content an Integer (if it is a number), easier
// for
// using it later on
// put content in table, with row/column inserted as x/y
t.setContent(row, col, (String) content);
}
}
} catch (ParserConfigurationException e) {
System.out.println("Fileparser could not be made");
} catch (IOException f) {
System.out.println("File could not be parsed, did you enter the correct file name?");
} catch (SAXException g) {
System.out.println("Something went wrong in parsing the file");
}
return t;
}
/**
* Unmarshall a Chromosome instance from a given XML Document representation. Its genes will be
* unmarshalled from the gene sub-elements.
*
* @param a_activeConfiguration the current active Configuration object that is to be used during
* construction of the Chromosome instances
* @param a_xmlDocument the XML Document representation of the Chromosome
* @return a new Chromosome instance setup with the data from the XML Document representation
* @throws ImproperXMLException if the given Document is improperly structured or missing data
* @throws InvalidConfigurationException if the given Configuration is in an inconsistent state
* @throws UnsupportedRepresentationException if the actively configured Gene implementation does
* not support the string representation of the alleles used in the given XML document
* @throws GeneCreationException if there is a problem creating or populating a Gene instance
* @author Neil Rotstan
* @since 1.0
*/
public static Chromosome getChromosomeFromDocument(
Configuration a_activeConfiguration, Document a_xmlDocument)
throws ImproperXMLException, InvalidConfigurationException,
UnsupportedRepresentationException, GeneCreationException {
// Extract the root element, which should be a chromosome element.
// After verifying that the root element is not null and that it
// in fact is a chromosome element, then convert it into a Chromosome
// instance.
// ------------------------------------------------------------------
Element rootElement = a_xmlDocument.getDocumentElement();
if (rootElement == null || !(rootElement.getTagName().equals(CHROMOSOME_TAG))) {
throw new ImproperXMLException(
"Unable to build Chromosome instance from XML Document: "
+ "'chromosome' element must be at root of Document.");
}
return getChromosomeFromElement(a_activeConfiguration, rootElement);
}
/** used for cut and paste. */
public void addObjectFromClipboard(String a_value) throws CircularIncludeException {
Reader reader = new StringReader(a_value);
Document document = null;
try {
document = UJAXP.getDocument(reader);
} catch (Exception e) {
e.printStackTrace();
return;
} // try-catch
Element root = document.getDocumentElement();
if (!root.getNodeName().equals("clipboard")) {
return;
} // if
Node child;
for (child = root.getFirstChild(); child != null; child = child.getNextSibling()) {
if (!(child instanceof Element)) {
continue;
} // if
Element element = (Element) child;
IGlyphFactory factory = GlyphFactory.getFactory();
if (XModule.isMatch(element)) {
EModuleInvoke module = (EModuleInvoke) factory.createXModule(element);
addModule(module);
continue;
} // if
if (XContour.isMatch(element)) {
EContour contour = (EContour) factory.createXContour(element);
addContour(contour);
continue;
} // if
if (XInclude.isMatch(element)) {
EIncludeInvoke include = (EIncludeInvoke) factory.createXInclude(element);
addInclude(include);
continue;
} // if
} // while
}
private void call(String script)
throws ParserConfigurationException, TransformerConfigurationException {
File callerScript = this.currentScript;
Document doc = null;
try {
this.currentScript = new File(script);
DocumentBuilder db = DocumentBuilderFactory.newInstance().newDocumentBuilder();
doc = db.parse(this.currentScript);
} catch (Exception ex) {
this.currentScript = callerScript;
Utils.onError(new Error.FileParse(script));
return;
}
NodeList operations = doc.getDocumentElement().getChildNodes();
for (int i = 0; i < operations.getLength(); i++) {
Node operation = operations.item(i);
if (operation.getNodeType() != Node.ELEMENT_NODE) continue;
call(operation);
}
this.currentScript = callerScript;
}
/** @param args the command line arguments */
public static void main(String[] args) {
// TODO code application logic here
Model model;
String insertstmt;
String insertmodel = "",
insertspecies = "",
insertcompartment = "",
insertfunction = "",
insertunitdef = "",
insertunits = "",
insertreaction = "",
insertreactant = "",
insertproduct = "";
String insertmodifier = "",
insertklaw = "",
insertrules = "",
insertconstraint = "",
insertdelay = "",
inserttrigger = "",
insertevent = "",
inserteventassign = "",
insertparameter = "";
String insertstatement = "";
String server, user, password, dbname, filepath;
String Filedata = "";
String cwd = System.getProperty("user.dir");
if (args.length == 0) {
server = "localhost";
user = "******";
password = "******";
dbname = "sbmldb2";
/**
* Path to extract the SBML files from database, where cwd is
* "github\db2sbml\dbtosbml_standalone_Project\dbtosbml" so add a folder in this directory and
* mention folder name instead of extractedbm folder
*/
filepath = cwd + "\\extractedbm\\";
} else {
server = args[0];
user = args[1];
password = args[2];
dbname = args[3];
filepath = args[4];
}
try {
Filedata = readFileAsString(cwd + "\\sbmldbschema.sql");
} catch (Exception e) {
e.printStackTrace();
}
Mysqlconn sql = new Mysqlconn(server, user, password, dbname);
// String modelids.getId() = "MorrisonAllegra" ;
ASTNode math = null;
int level = 0, version = 0;
ArrayList<modellist> modelidlist = sql.getmodels();
insertstatement =
"LOCK TABLES `model` WRITE,`species` WRITE,`compartment` WRITE,`functiondefinition` WRITE,";
insertstatement =
insertstatement
+ "`listofunitdefinitions` WRITE,`listofunits` WRITE,`reaction` WRITE,`simplespeciesreference` WRITE,";
insertstatement =
insertstatement
+ "`modifierspeciesreference` WRITE,`kineticlaw` WRITE,`parameter` WRITE,`sbmlconstraint` WRITE,";
insertstatement =
insertstatement
+ "`event` WRITE,`sbmltrigger` WRITE,`delay` WRITE,`eventassignment` WRITE,`rules` WRITE"
+ ";";
for (modellist modelids : modelidlist) {
ArrayList<modellist> modellevel = sql.getmodeldetails(modelids.getId());
for (modellist modellv : modellevel) {
level = modellv.getlevel();
version = modellv.getversion();
}
SBMLDocument doc = new SBMLDocument(level, version);
ArrayList<modellist> modellists = sql.getmodeldetails(modelids.getId());
if (!modellists.isEmpty())
insertmodel =
insertmodel
+ "\nInsert Into model (id, name,SBML_level,version,notes,annotation) Values";
for (modellist models : modellists) {
insertmodel =
insertmodel
+ "(\'"
+ models.getId()
+ "\',\'"
+ models.getName()
+ "\',"
+ models.getlevel()
+ ","
+ models.getversion()
+ ",\'"
+ models.getnotes()
+ "\',\'"
+ models.getannotation().toString()
+ "\'),";
model = doc.createModel(models.getId());
model.setName(models.getName());
// System.out.println("model : " + models.getId());
// model.setNotes(models.getnotes()); // there is some null exception is command line run
// but run perfectly from netbeans so ommented out
if (!models.getannotation().equals("")) {
Annotation annot = new Annotation(models.getannotation().toString());
model.setAnnotation(annot);
}
doc.setModel(model);
}
if (!modellists.isEmpty()) {
insertmodel = insertmodel.substring(0, insertmodel.length() - 1);
insertmodel = insertmodel + ';';
}
// insertmodel = insertmodel + "\nUNLOCK TABLES;";
// System.out.println(insertmodel);
ArrayList<SpeciesList> specieslist = sql.getspecies(modelids.getId());
if (!specieslist.isEmpty())
insertspecies =
insertspecies
+ "\nInsert Into species (id, name, compartment, initialAmount, initialConcentration,substanceUnits,hasOnlySubstanceUnits,boundaryCondition,constant,conversionFactor,model_id,annotation) Values";
for (SpeciesList species : specieslist) {
insertspecies =
insertspecies
+ "(\'"
+ species.getId()
+ "\',\'"
+ species.getName()
+ "\',\'"
+ species.getcompartment()
+ "\',"
+ species.getia()
+ ","
+ species.getic()
+ ",\'"
+ species.getsu()
+ "\',"
+ species.gethosu()
+ ","
+ species.getbc()
+ ","
+ species.getconstant()
+ ","
+ species.getcf()
+ ",\'"
+ modelids.getId()
+ "\',\'"
+ species.getannotation()
+ "\'),";
Species sp = doc.getModel().createSpecies(species.getId());
sp.setName(species.getName());
sp.setCompartment(species.getcompartment());
sp.setConstant(species.getconstant());
sp.setInitialAmount(species.getia());
sp.setInitialConcentration(species.getic());
sp.setHasOnlySubstanceUnits(species.gethosu());
if (doc.getModel().getLevel() == 3) sp.setConversionFactor(species.getcf());
sp.setBoundaryCondition(species.getbc());
sp.setSubstanceUnits(species.getsu());
if (!species.getannotation().equals("")) {
Annotation annot = new Annotation(species.getannotation().toString());
sp.setAnnotation(annot);
}
// doc.getModel().addSpecies(sp) ;
}
if (!specieslist.isEmpty()) {
insertspecies = insertspecies.substring(0, insertspecies.length() - 1);
insertspecies = insertspecies + ';';
}
ArrayList<CompartmentList> complist = sql.getcompartments(modelids.getId());
if (!complist.isEmpty())
insertcompartment =
insertcompartment
+ "\nInsert Into compartment (id, name,constant,model_id,spacialDimensions,size,units) Values";
for (CompartmentList comp : complist) {
insertcompartment =
insertcompartment
+ "(\'"
+ comp.getId()
+ "\',\'"
+ comp.getName()
+ "\',"
+ comp.getconstant()
+ ",\'"
+ modelids.getId()
+ "\',"
+ comp.getspatialdimensions()
+ ","
+ comp.getsize()
+ ","
+ comp.getunits()
+ "\'),";
Compartment c = doc.getModel().createCompartment(comp.getId());
c.setName(comp.getName());
c.setConstant(comp.getconstant());
c.setSize(comp.getsize());
c.setSpatialDimensions(comp.getspatialdimensions());
if (comp.getspatialdimensions() != 0) c.setUnits(comp.getunits());
// doc.getModel().addSpecies(sp) ;
}
if (!complist.isEmpty()) {
insertcompartment = insertcompartment.substring(0, insertcompartment.length() - 1);
insertcompartment = insertcompartment + ';';
}
ArrayList<functionList> funclist = sql.getfunctions(modelids.getId());
if (!funclist.isEmpty())
insertfunction =
insertfunction + "\nInsert Into functiondefinition (id, xmlns,model_id) Values";
for (functionList func : funclist) {
insertfunction =
insertfunction
+ "(\'"
+ func.getId()
+ "\',\'"
+ func.getxmlns()
+ "\',\'"
+ modelids.getId()
+ "\'),";
FunctionDefinition fd = doc.getModel().createFunctionDefinition(func.getId());
try {
math = ASTNode.parseFormula(func.getxmlns());
fd.setMath(math);
} catch (Exception e) {
e.printStackTrace();
}
}
if (!funclist.isEmpty()) {
insertfunction = insertfunction.substring(0, insertfunction.length() - 1);
insertfunction = insertfunction + ';';
}
ArrayList<unitList> unitlist = sql.getunitlist(modelids.getId());
if (!unitlist.isEmpty())
insertunitdef =
insertunitdef + "\nInsert Into listofunitdefinitions (id,name,model_id) Values";
for (unitList units : unitlist) {
insertunitdef =
insertunitdef
+ "(\'"
+ units.getId()
+ "\',\'"
+ units.getName()
+ "\',\'"
+ modelids.getId()
+ "\'),";
UnitDefinition ud = doc.getModel().createUnitDefinition(units.getId());
ud.setName(units.getName());
ArrayList<unitList> unitdeflist = sql.getunitdef(units.getId());
if (!unitdeflist.isEmpty())
insertunits =
insertunits
+ "\nInsert Into listofunits (listofunitdefinitions_id,kind, scale,exponent,multiplier) Values";
for (unitList unitdef : unitdeflist) {
insertunits =
insertunits
+ "(\'"
+ units.getId()
+ "\',\'"
+ unitdef.getkind()
+ "\',"
+ unitdef.getscale()
+ ","
+ unitdef.getexponent()
+ ","
+ unitdef.getmultiplier()
+ "),";
Unit u = ud.createUnit(Unit.Kind.valueOf(unitdef.getkind()));
u.setScale(unitdef.getscale());
u.setExponent(unitdef.getexponent());
u.setMultiplier(unitdef.getmultiplier());
}
// doc.getModel().addSpecies(sp) ;
if (!unitdeflist.isEmpty()) {
insertunits = insertunits.substring(0, insertunits.length() - 1);
insertunits = insertunits + ';';
}
}
if (!unitlist.isEmpty()) {
insertunitdef = insertunitdef.substring(0, insertunitdef.length() - 1);
insertunitdef = insertunitdef + ';';
}
ArrayList<reactionList> reactionlist = sql.getreactons(modelids.getId());
if (!reactionlist.isEmpty())
insertreaction =
insertreaction
+ "\nInsert Into reaction (id,name, reversible,fast,model_id,compartment,annotation) Values";
for (reactionList reaction : reactionlist) {
insertreaction =
insertreaction
+ "(\'"
+ reaction.getId()
+ "\',\'"
+ reaction.getName()
+ "\',"
+ reaction.getreversible()
+ ","
+ reaction.getfast()
+ ",\'"
+ modelids.getId()
+ "\',\'"
+ reaction.getcompartment()
+ "\',\'"
+ reaction.getannotation()
+ "\'),";
Reaction rn = doc.getModel().createReaction(reaction.getId());
rn.setName(reaction.getName());
if (doc.getModel().getLevel() == 3) rn.setCompartment(reaction.getcompartment());
rn.setFast(reaction.getfast());
rn.setReversible(reaction.getreversible());
if (!reaction.getannotation().equals("")) {
Annotation annot = new Annotation(reaction.getannotation().toString());
rn.setAnnotation(annot);
}
ArrayList<reactionList> reactantlist = sql.getreactants(reaction.getId());
if (!reactantlist.isEmpty())
insertreactant =
insertreactant
+ "\nInsert Into simplespeciesreference (reaction_id,species, sboTerm,stoichiometry,speciestype,constant) Values";
for (reactionList reactant : reactantlist) {
insertreactant =
insertreactant
+ "(\'"
+ reaction.getId()
+ "\',\'"
+ reactant.getspecies()
+ "\',\'"
+ reactant.getsboTerm()
+ "\',"
+ reactant.getstoichometry()
+ ","
+ reactant.getconstant()
+ ",\'reactants\'),";
SpeciesReference rt = new SpeciesReference();
rt.setName(reactant.getspecies());
rt.setSpecies(reactant.getspecies());
// rt.setSBOTerm(reactant.getsboTerm());
rt.setStoichiometry(reactant.getstoichometry());
// rt.setConstant(reactant.getconstant());
rn.addReactant(rt);
}
if (!reactantlist.isEmpty()) {
insertreactant = insertreactant.substring(0, insertreactant.length() - 1);
insertreactant = insertreactant + ';';
}
ArrayList<reactionList> productlist = sql.getproducts(reaction.getId());
if (!productlist.isEmpty())
insertproduct =
insertproduct
+ "\nInsert Into simplespeciesreference (reaction_id,species, sboTerm,stoichiometry,constant,speciestype) Values";
for (reactionList product : productlist) {
insertproduct =
insertproduct
+ "(\'"
+ reaction.getId()
+ "\',\'"
+ product.getspecies()
+ "\',\'"
+ product.getsboTerm()
+ "\',"
+ product.getstoichometry()
+ ","
+ product.getconstant()
+ ",\'products\'),";
SpeciesReference pr = new SpeciesReference();
pr.setName(product.getspecies());
pr.setSpecies(product.getspecies());
// pr.setSBOTerm(product.getsboTerm());
pr.setStoichiometry(product.getstoichometry());
// pr.setConstant(product.getconstant());
rn.addProduct(pr);
}
if (!productlist.isEmpty()) {
insertproduct = insertproduct.substring(0, insertproduct.length() - 1);
insertproduct = insertproduct + ';';
}
ArrayList<reactionList> modifierlist = sql.getmodifiers(reaction.getId());
if (!modifierlist.isEmpty())
insertmodifier =
insertmodifier
+ "\nInsert Into modifierspeciesreference (reaction_id,species, sboTerm,speciestype) Values";
for (reactionList modifier : modifierlist) {
insertmodifier =
insertmodifier
+ "(\'"
+ reaction.getId()
+ "\',\'"
+ modifier.getspecies()
+ "\',\'"
+ modifier.getsboTerm()
+ "\',\'modifiers\'),";
ModifierSpeciesReference m = new ModifierSpeciesReference();
m.setName(modifier.getspecies());
m.setSpecies(modifier.getspecies());
// m.setSBOTerm(modifier.getsboTerm());
rn.addModifier(m);
}
if (!modifierlist.isEmpty()) {
insertmodifier = insertmodifier.substring(0, insertmodifier.length() - 1);
insertmodifier = insertmodifier + ';';
}
ArrayList<reactionList> klawlist = sql.getkineticlaws(reaction.getId());
if (!klawlist.isEmpty())
insertklaw =
insertklaw + "\nInsert Into kineticlaw (reaction_id,kid, math,annotation) Values";
for (reactionList klaw : klawlist) {
insertklaw =
insertklaw
+ "(\'"
+ reaction.getId()
+ "\',\'"
+ klaw.getId()
+ "\',\'"
+ klaw.getmath()
+ "\',\'"
+ klaw.getannotation()
+ "\'),";
KineticLaw kl = rn.createKineticLaw();
try {
math = ASTNode.parseFormula(klaw.getmath());
kl.setMath(math);
if (!klaw.getannotation().equals("")) {
Annotation annot = new Annotation(klaw.getannotation().toString());
kl.setAnnotation(annot);
}
} catch (Exception e) {
e.printStackTrace();
}
}
if (!klawlist.isEmpty()) {
insertklaw = insertklaw.substring(0, insertklaw.length() - 1);
insertklaw = insertklaw + ';';
}
}
if (!reactionlist.isEmpty()) {
insertreaction = insertreaction.substring(0, insertreaction.length() - 1);
insertreaction = insertreaction + ';';
}
ArrayList<parameterList> paralist = sql.getparameters(modelids.getId());
if (!paralist.isEmpty())
insertparameter =
insertparameter
+ "\nInsert Into parameter (id,name,value,units,constant,model_id) Values";
for (parameterList para : paralist) {
insertparameter =
insertparameter
+ "(\'"
+ para.getId()
+ "\',\'"
+ para.getName()
+ "\',"
+ para.getvalue()
+ ","
+ para.getunits()
+ ","
+ para.getconstant()
+ ",\'"
+ modelids.getId()
+ "\'),";
Parameter par = doc.getModel().createParameter(para.getId());
par.setName(para.getId());
par.setConstant(para.getconstant());
par.setUnits(para.getunits());
par.setValue(para.getvalue());
}
if (!paralist.isEmpty()) {
insertparameter = insertparameter.substring(0, insertparameter.length() - 1);
insertparameter = insertparameter + ';';
}
ArrayList<constraintList> conslist = sql.getconstraints(modelids.getId());
if (!conslist.isEmpty())
insertconstraint =
insertconstraint + "\nInsert Into sbmlconstraint (math,message,model_id) Values";
for (constraintList constraint : conslist) {
insertconstraint =
insertconstraint
+ "(\'"
+ constraint.getmath()
+ "\',\'"
+ constraint.getmessage()
+ "\',\'"
+ modelids.getId()
+ "\'),";
Constraint cons = doc.getModel().createConstraint();
try {
math = ASTNode.parseFormula(constraint.getmath());
cons.setMath(math);
cons.setMessage(constraint.getmessage());
} catch (Exception e) {
e.printStackTrace();
}
}
if (!conslist.isEmpty()) {
insertconstraint = insertconstraint.substring(0, insertconstraint.length() - 1);
insertconstraint = insertconstraint + ';';
}
ArrayList<eventsList> eventlist = sql.getevents(modelids.getId());
if (!eventlist.isEmpty())
insertevent =
insertevent + "\nInsert Into event (id,name,UseValuesFromTriggerTime,model_id) Values";
for (eventsList events : eventlist) {
insertevent =
insertevent
+ "(\'"
+ events.getId()
+ "\',\'"
+ events.getName()
+ "\',"
+ events.getuservalues()
+ ",\'"
+ modelids.getId()
+ "\'),";
Event ev = doc.getModel().createEvent(events.getId());
ev.setName(events.getName());
// ev.setUseValuesFromTriggerTime(events.getuservalues());
ArrayList<eventsList> triggerlist = sql.gettriggers(events.getId());
if (!triggerlist.isEmpty())
inserttrigger =
inserttrigger
+ "\nInsert Into sbmltrigger (event_id,initialvalue,persisent,math) Values";
for (eventsList triggers : triggerlist) {
Trigger tr = doc.getModel().createTrigger();
try {
math = ASTNode.parseFormula(triggers.getmath());
tr.setMath(math);
tr.setInitialValue(triggers.getinitialval());
tr.setPersistent(triggers.getpersistent());
} catch (Exception e) {
e.printStackTrace();
}
}
if (!triggerlist.isEmpty()) {
inserttrigger = inserttrigger.substring(0, insertmodel.length() - 1);
inserttrigger = inserttrigger + ';';
}
ArrayList<eventsList> delaylist = sql.getdelays(events.getId());
if (!delaylist.isEmpty())
insertdelay = insertdelay + "\nInsert Into delay (event_id,math) Values";
for (eventsList delays : delaylist) {
Delay d = doc.getModel().createDelay();
try {
math = ASTNode.parseFormula(delays.getmath());
d.setMath(math);
} catch (Exception e) {
e.printStackTrace();
}
}
if (!delaylist.isEmpty()) {
insertdelay = insertdelay.substring(0, insertdelay.length() - 1);
insertdelay = insertdelay + ';';
}
ArrayList<eventsList> evasslist = sql.geteventassignments(events.getId());
if (!evasslist.isEmpty())
inserteventassign =
inserteventassign + "\nInsert Into eventassignment (event_id,variable,math) Values";
for (eventsList evassign : evasslist) {
EventAssignment ea = doc.getModel().createEventAssignment();
try {
math = ASTNode.parseFormula(evassign.getmath());
ea.setMath(math);
} catch (Exception e) {
e.printStackTrace();
}
}
if (!evasslist.isEmpty()) {
inserteventassign = inserteventassign.substring(0, inserteventassign.length() - 1);
inserteventassign = inserteventassign + ';';
}
}
if (!eventlist.isEmpty()) {
insertevent = insertevent.substring(0, insertevent.length() - 1);
insertevent = insertevent + ';';
}
ArrayList<ruleslist> rulelist = sql.getrules(modelids.getId());
if (!rulelist.isEmpty())
insertrules = insertrules + "\nInsert Into rules (id,math,ruletype,model_id) Values";
for (ruleslist rules : rulelist) {
insertrules =
insertrules
+ "(\'"
+ rules.getId()
+ "\',\'"
+ rules.getmath()
+ "\',\'"
+ rules.getruletype()
+ "\',\'"
+ modelids.getId()
+ "\'),";
if (rules.getruletype().equals("assignmentrule")) {
Rule r = doc.getModel().createAssignmentRule();
r.setMetaId(rules.getId());
try {
math = ASTNode.parseFormula(rules.getmath());
r.setMath(math);
} catch (Exception e) {
e.printStackTrace();
}
}
}
if (!rulelist.isEmpty()) {
insertrules = insertrules.substring(0, insertrules.length() - 1);
insertrules = insertrules + ';';
}
SBMLWriter writer = new SBMLWriter();
try {
String Path = filepath + modelids.getId() + ".xml";
writer.write(doc, Path);
DocumentBuilderFactory documentBuilderFactory = DocumentBuilderFactory.newInstance();
DocumentBuilder documentBuilder = documentBuilderFactory.newDocumentBuilder();
Document document = documentBuilder.parse(Path);
Element root = document.getDocumentElement();
Element newdataset = document.createElement("dataset");
root.appendChild(newdataset);
ArrayList<dataset> datasetlist = sql.getdataset(modelids.getId());
for (dataset ds : datasetlist) {
// System.out.println(ds.getexpcond());
Element name = document.createElement("experimentalcondition");
name.setAttribute("bioelement", ds.getbioel());
name.setAttribute("name", ds.getName());
name.setAttribute("descr", ds.getdescr());
name.setAttribute("expcond", ds.getexpcond());
name.setAttribute("value", String.valueOf(ds.getvalue()));
name.setAttribute("type", ds.gettype());
name.setAttribute("uri", ds.geturi());
newdataset.appendChild(name);
}
root.appendChild(newdataset);
DOMSource source = new DOMSource(document);
TransformerFactory transformerFactory = TransformerFactory.newInstance();
Transformer transformer = transformerFactory.newTransformer();
transformer.setOutputProperty(OutputKeys.INDENT, "yes");
transformer.setOutputProperty(OutputKeys.METHOD, "xml");
transformer.setOutputProperty("{http://xml.apache.org/xslt}indent-amount", "4");
StreamResult result = new StreamResult(filepath + modelids.getId() + "d.xml");
transformer.transform(source, result);
System.out.println(
"Files : "
+ modelids.getId()
+ ".xml and "
+ modelids.getId()
+ "d.xml have been generated successfully !!!");
} catch (Exception e) {
e.printStackTrace();
}
insertstatement =
insertstatement
+ "\n\n"
+ insertmodel
+ "\n"
+ insertspecies
+ "\n"
+ insertcompartment
+ "\n"
+ insertfunction;
insertstatement =
insertstatement
+ "\n"
+ insertparameter
+ "\n"
+ insertreaction
+ "\n"
+ insertreactant
+ "\n"
+ insertproduct;
insertstatement =
insertstatement
+ "\n"
+ insertmodifier
+ "\n"
+ insertklaw
+ "\n"
+ insertunitdef
+ "\n"
+ insertunits;
insertstatement =
insertstatement
+ "\n"
+ insertrules
+ "\n"
+ insertconstraint
+ "\n"
+ insertevent
+ "\n"
+ inserttrigger
+ "\n"
+ insertdelay
+ "\n"
+ inserteventassign;
insertcompartment = "";
insertmodel = "";
insertspecies = "";
// System.out.println("document : " + doc);
}
insertstatement = insertstatement + "\nUNLOCK TABLES;";
Filedata = Filedata + "\n\n\n" + insertstatement;
try {
wrtireStringToFile(Filedata, filepath + "sbmldb.sql");
} catch (IOException e) {
e.printStackTrace();
}
// System.out.println(insertstatement);
}
// ## operation readReactorOutputFile(ReactionModel)
public SystemSnapshot readReactorOutputFile(ReactionModel p_reactionModel) {
// #[ operation readReactorOutputFile(ReactionModel)
try {
// open output file and build the DOM tree
String dir = System.getProperty("RMG.workingDirectory");
String filename = "chemkin/reactorOutput.xml";
File inputFile = new File(filename);
DocumentBuilderFactory factory = DocumentBuilderFactory.newInstance();
factory.setValidating(true); // validate the document with the DTD
factory.setIgnoringElementContentWhitespace(true); // ignore whitespace
DocumentBuilder builder = factory.newDocumentBuilder();
Document doc = builder.parse(inputFile);
// get root element and its children
Element root = doc.getDocumentElement();
NodeList rootchildren = root.getChildNodes();
// header is rootchildren.item(0)
// get return message and check for successful run
Element returnmessageElement = (Element) rootchildren.item(1);
Text returnmessageText = (Text) returnmessageElement.getFirstChild();
String returnmessage = returnmessageText.toString();
returnmessage = returnmessage.trim();
if (!returnmessage.contains("SUCCESSFULLY COMPLETED RUN.")) {
System.out.println("External reactor model failed!");
System.out.println("Reactor model error message: " + returnmessage);
System.exit(0);
}
// get outputvalues element and its children
Element outputvaluesElement = (Element) rootchildren.item(2);
NodeList children = outputvaluesElement.getChildNodes();
// get time
Element timeElement = (Element) children.item(0);
Text timeText = (Text) timeElement.getFirstChild();
double time = Double.parseDouble(timeText.getData());
String timeUnits = timeElement.getAttribute("units");
// get systemstate element and its children
Element systemstateElement = (Element) children.item(1);
NodeList states = systemstateElement.getChildNodes();
// get temperature and its units
Element temperatureElement = (Element) states.item(0);
String tempUnits = temperatureElement.getAttribute("units");
Text temperatureText = (Text) temperatureElement.getFirstChild();
double temp = Double.parseDouble(temperatureText.getData());
Temperature T = new Temperature(temp, tempUnits);
// get pressure and its units
Element pressureElement = (Element) states.item(1);
String presUnits = pressureElement.getAttribute("units");
Text pressureText = (Text) pressureElement.getFirstChild();
double pres = Double.parseDouble(pressureText.getData());
Pressure P = new Pressure(pres, presUnits);
// get species amounts (e.g. concentrations)
ArrayList speciesIDs = new ArrayList();
ArrayList amounts = new ArrayList();
ArrayList fluxes = new ArrayList();
String amountUnits = null;
String fluxUnits = null;
// loop thru all the species
// begin at i=2, since T and P take already the first two position of states
int nSpe = (states.getLength() - 2) / 2;
int index = 0;
LinkedHashMap inertGas = new LinkedHashMap();
for (int i = 2; i < nSpe + 2; i++) {
// get amount element and the units
Element amountElement = (Element) states.item(i);
amountUnits = amountElement.getAttribute("units");
Element fluxElement = (Element) states.item(i + nSpe);
fluxUnits = fluxElement.getAttribute("units");
// get speciesid and store in an array list
String thisSpeciesID = amountElement.getAttribute("speciesid");
// get amount (e.g. concentraion) and store in an array list
Text amountText = (Text) amountElement.getFirstChild();
double thisAmount = Double.parseDouble(amountText.getData());
if (thisAmount < 0) {
double aTol = ReactionModelGenerator.getAtol();
// if (Math.abs(thisAmount) < aTol) thisAmount = 0;
// else throw new NegativeConcentrationException("Negative concentration in
// reactorOutput.xml: " + thisSpeciesID);
if (thisAmount < -100.0 * aTol)
throw new NegativeConcentrationException(
"Species "
+ thisSpeciesID
+ " has negative concentration: "
+ String.valueOf(thisAmount));
}
// get amount (e.g. concentraion) and store in an array list
Text fluxText = (Text) fluxElement.getFirstChild();
double thisFlux = Double.parseDouble(fluxText.getData());
if (thisSpeciesID.compareToIgnoreCase("N2") == 0
|| thisSpeciesID.compareToIgnoreCase("Ne") == 0
|| thisSpeciesID.compareToIgnoreCase("Ar") == 0) {
inertGas.put(thisSpeciesID, new Double(thisAmount));
} else {
speciesIDs.add(index, thisSpeciesID);
amounts.add(index, new Double(thisAmount));
fluxes.add(index, new Double(thisFlux));
index++;
}
}
// print results for debugging purposes
/**
* System.out.println(returnmessage); System.out.println("Temp = " + temp + " " + tempUnits);
* System.out.println("Pres = " + pres + " " + presUnits); for (int i = 0; i < amounts.size();
* i++) { System.out.println(speciesIDs.get(i) + " " + amounts.get(i) + " " + amountUnits); }
*/
ReactionTime rt = new ReactionTime(time, timeUnits);
LinkedHashMap speStatus = generateSpeciesStatus(p_reactionModel, speciesIDs, amounts, fluxes);
SystemSnapshot ss = new SystemSnapshot(rt, speStatus, T, P);
ss.inertGas = inertGas;
return ss;
} catch (Exception e) {
System.out.println("Error reading reactor model output: " + e.getMessage());
System.exit(0);
return null;
}
// #]
}
