research/cantera
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Added a solver variable.

Browse source
This commit is contained in:
Harry Moffat 2007-11-11 23:27:06 +00:00
parent 9b84ab2f55
commit ab100d892e
1 changed files with 6 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

11
test_problems/ChemEquil_gri_matrix/gri_matrix.cpp
View file

@ -20,6 +20,7 @@ int main(int argc, char **argv) {
try { try {
int retnSub; int retnSub;
double T = 500.; double T = 500.;
int solver = -1;
IdealGasMix g("gri30.xml", "gri30_mix"); IdealGasMix g("gri30.xml", "gri30_mix");
double pres = OneAtm; double pres = OneAtm;
@ -39,7 +40,7 @@ int main(int argc, char **argv) {
Xmol[iN2] = 0.4; Xmol[iN2] = 0.4;
g.setState_TPX(T, pres, DATA_PTR(Xmol)); g.setState_TPX(T, pres, DATA_PTR(Xmol));
try { try {
retnSub = equilibrate(g, "TP", -1); retnSub = equilibrate(g, "TP", solver);
cout << g; cout << g;
if (retnSub != 1) { if (retnSub != 1) {
cerr << "ERROR: ChemEquil equilibration step failed at " cerr << "ERROR: ChemEquil equilibration step failed at "
@ -75,7 +76,7 @@ int main(int argc, char **argv) {
Xmol[iN2] = 0.4; Xmol[iN2] = 0.4;
g.setState_TPX(T, pres, DATA_PTR(Xmol)); g.setState_TPX(T, pres, DATA_PTR(Xmol));
try { try {
retnSub = equilibrate(g, "TP", -1); retnSub = equilibrate(g, "TP", solver);
cout << g; cout << g;
if (retnSub != 1) { if (retnSub != 1) {
cerr << "ERROR: ChemEquil equilibration step failed at " cerr << "ERROR: ChemEquil equilibration step failed at "
@ -113,7 +114,7 @@ int main(int argc, char **argv) {
pres = OneAtm; pres = OneAtm;
g.setState_TPX(T, pres, DATA_PTR(Xmol)); g.setState_TPX(T, pres, DATA_PTR(Xmol));
try { try {
retnSub = equilibrate(g, "TP", -1); retnSub = equilibrate(g, "TP", solver);
cout << g; cout << g;
if (retnSub != 1) { if (retnSub != 1) {
cerr << "ERROR: ChemEquil equilibration step failed at " cerr << "ERROR: ChemEquil equilibration step failed at "
@ -151,7 +152,7 @@ int main(int argc, char **argv) {
pres = 1.0; pres = 1.0;
g.setState_TPX(T, pres, DATA_PTR(Xmol)); g.setState_TPX(T, pres, DATA_PTR(Xmol));
try { try {
retnSub = equilibrate(g, "TP", -1); retnSub = equilibrate(g, "TP", solver);
cout << g; cout << g;
if (retnSub != 1) { if (retnSub != 1) {
cerr << "ERROR: ChemEquil equilibration step failed at " cerr << "ERROR: ChemEquil equilibration step failed at "
@ -197,7 +198,7 @@ int main(int argc, char **argv) {
retnSub = 0; retnSub = 0;
try { try {
retnSub = equilibrate(g, "TP", -1); retnSub = equilibrate(g, "TP", solver);
if (retnSub != 1) { if (retnSub != 1) {
cerr << "ERROR: ChemEquil equilibration step failed at " cerr << "ERROR: ChemEquil equilibration step failed at "
<< " T = " << T << " T = " << T