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

[1D] Added a new log level for saving intermediate solutions

Browse source 
Running Sim1D::solve with loglevel set to 7 or higher will now save the solution
after each attempted steady-state solution, after timestepping, and after
regridding to the file 'debug_sim1d.xml'. This can be useful for diagnosing
convergence issues with some problems.
This commit is contained in:
Ray Speth 2012-12-04 20:28:59 +00:00
parent 2735802eb9
commit 351423bcda
1 changed files with 14 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

14
src/oneD/Sim1D.cpp
View file

@ -301,16 +301,27 @@ void Sim1D::solve(int loglevel, bool refine_grid)
writelog("]"); writelog("]");
writelog(" point grid(s).\n"); writelog(" point grid(s).\n");
} }
if (loglevel > 6) {
save("debug_sim1d.xml", "debug",
"After successful Newton solve");
}
ok = true; ok = true;
soln_number++; soln_number++;
} else { } else {
char buf[100]; char buf[100];
if (loglevel > 0) { if (loglevel > 0) {
writelog(" failure. \n"); writelog(" failure. \n");
if (loglevel > 6) {
save("debug_sim1d.xml", "debug",
"After unsuccessful Newton solve");
}
writelog("Take "+int2str(nsteps)+" timesteps "); writelog("Take "+int2str(nsteps)+" timesteps ");
} }
dt = timeStep(nsteps, dt, DATA_PTR(m_x), DATA_PTR(m_xnew), dt = timeStep(nsteps, dt, DATA_PTR(m_x), DATA_PTR(m_xnew),
loglevel-1); loglevel-1);
if (loglevel > 6) {
save("debug_sim1d.xml", "debug", "After timestepping");
}
if (loglevel == 1) { if (loglevel == 1) {
sprintf(buf, " %10.4g %10.4g \n", dt, sprintf(buf, " %10.4g %10.4g \n", dt,
log10(ssnorm(DATA_PTR(m_x), DATA_PTR(m_xnew)))); log10(ssnorm(DATA_PTR(m_x), DATA_PTR(m_xnew))));
@ -337,6 +348,9 @@ void Sim1D::solve(int loglevel, bool refine_grid)
if (refine_grid) { if (refine_grid) {
new_points = refine(loglevel); new_points = refine(loglevel);
if (new_points && loglevel > 6) {
save("debug_sim1d.xml", "debug", "After regridding");
}
if (new_points < 0) { if (new_points < 0) {
writelog("Maximum number of grid points reached."); writelog("Maximum number of grid points reached.");
new_points = 0; new_points = 0;