static PyObject* py_xml_get_XML_File(PyObject *self, PyObject *args) { char* file; int debug; if (!PyArg_ParseTuple(args, "si:xml_get_XML_File", &file, &debug)) return NULL; int n = xml_get_XML_File(file, debug); if (n < 0) return reportError(n); PyObject* pn = Py_BuildValue("i",n); return pn; }
int main(int argc, char** argv) { int ret; int xml_file = xml_get_XML_File("gri30.xml", 0); assert(xml_file > 0); int phase_node = xml_findID(xml_file, "gri30_mix"); assert(phase_node > 0); int thermo = thermo_newFromXML(phase_node); assert(thermo > 0); int nsp = thermo_nSpecies(thermo); assert(nsp == 53); ret = thermo_setTemperature(thermo, 500); assert(ret == 0); ret = thermo_setPressure(thermo, 5 * 101325); assert(ret == 0); ret = thermo_setMoleFractionsByName(thermo, "CH4:1.0, O2:2.0, N2:7.52"); assert(ret == 0); ret = thermo_equilibrate(thermo, "HP", 0, 1e-9, 50000, 1000, 0); assert(ret == 0); double T = thermo_temperature(thermo); assert(T > 2200 && T < 2300); ret = thermo_print(thermo, 1, 0); assert(ret == 0); int kin = kin_newFromXML(phase_node, thermo, 0, 0, 0, 0); assert(kin > 0); size_t nr = kin_nReactions(kin); assert(nr == 325 ); ret = thermo_setTemperature(thermo, T - 200); assert(ret == 0); char buf [1000]; double ropf[325]; printf("\n Reaction Forward ROP\n"); kin_getFwdRatesOfProgress(kin, 325, ropf); size_t n; // declare this here for C89 compatibility for (n = 0; n < nr; n++) { kin_getReactionString(kin, n, 1000, buf); printf("%35s %8.6e\n", buf, ropf[n]); } printf("\n Species Mix diff coeff\n"); int tran = trans_new("Mix", thermo, 0); double dkm[53]; trans_getMixDiffCoeffs(tran, 53, dkm); int k; // declare this here for C89 compatibility for (k = 0; k < nsp; k++) { thermo_getSpeciesName(thermo, k, 1000, buf); printf("%10s %8.6e\n", buf, dkm[k]); } ret = thermo_setTemperature(thermo, 1050); assert(ret == 0); ret = thermo_setPressure(thermo, 5 * 101325); assert(ret == 0); ret = thermo_setMoleFractionsByName(thermo, "CH4:1.0, O2:2.0, N2:7.52"); assert(ret == 0); printf("\ntime Temperature\n"); int reactor = reactor_new(5); int net = reactornet_new(); ret = reactor_setThermoMgr(reactor, thermo); assert(ret == 0); ret = reactor_setKineticsMgr(reactor, kin); assert(ret == 0); ret = reactornet_addreactor(net, reactor); assert(ret == 0); double t = 0.0; while (t < 0.1 && ret == 0) { double T = reactor_temperature(reactor); t = reactornet_time(net); printf("%.2e %.3f\n", t, T); ret = reactornet_advance(net, t + 5e-3); assert(ret == 0); } ct_appdelete(); return 0; }
void xmlmethods(int nlhs, mxArray* plhs[], int nrhs, const mxArray* prhs[]) { int j, m, iok = 0; char* file, *key, *val, *nm; int job = getInt(prhs[1]); int i = getInt(prhs[2]); // Check for proper number of arguments if (!nargs_ok(job,nrhs-1)) { mexErrMsgTxt("Wrong number of inputs."); return; } else if (nlhs > 1) { mexErrMsgTxt("Too many output arguments"); } // options that do not return a value if (job < 20) { switch (job) { case 0: nm = getString(prhs[3]); iok = xml_new(nm); break; case 1: iok = xml_del(i); break; case 2: iok = xml_copy(i); break; case 4: file = getString(prhs[3]); iok = xml_build(i, file); break; case 5: key = getString(prhs[3]); val = getString(prhs[4]); iok = xml_addAttrib(i, key, val); break; case 6: key = getString(prhs[3]); iok = xml_child(i, key); break; case 7: m = getInt(prhs[3]); iok = xml_child_bynumber(i, m); break; case 8: key = getString(prhs[3]); iok = xml_findID(i, key); break; case 9: key = getString(prhs[3]); iok = xml_findByName(i, key); break; case 10: iok = xml_nChildren(i); break; case 11: key = getString(prhs[3]); val = getString(prhs[4]); iok = xml_addChild(i, key, val); break; case 12: key = getString(prhs[3]); j = getInt(prhs[4]); iok = xml_addChildNode(i, j); break; case 13: file = getString(prhs[3]); iok = xml_write(i, file); break; case 14: j = getInt(prhs[3]); iok = xml_removeChild(i, j); break; case 15: file = getString(prhs[3]); iok = xml_get_XML_File(file, 0); break; default: mexErrMsgTxt("unknown job parameter"); } plhs[0] = mxCreateNumericMatrix(1,1,mxDOUBLE_CLASS,mxREAL); double* h = mxGetPr(plhs[0]); *h = double(iok); if (iok < 0) { reportError(); } return; } // options that return strings char* v = (char*)mxCalloc(80, sizeof(char)); switch (job) { case 20: // return an attribute key = getString(prhs[3]); iok = xml_attrib(i, key, v); break; case 21: // return the value of the node iok = xml_value(i, v); break; case 22: iok = xml_tag(i, v); break; default: mexErrMsgTxt("unknown job parameter"); } if (iok < 0) { plhs[0] = mxCreateNumericMatrix(1,1,mxDOUBLE_CLASS,mxREAL); double* h = mxGetPr(plhs[0]); *h = double(iok); if (iok < 0) { reportError(); } } else { plhs[0] = mxCreateString(v); } }