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

[1D] Exceeding the max number of grid points should be an error

Browse source
This commit is contained in:
Bryan W. Weber 2016-11-10 19:27:27 -05:00 • committed by Bryan W. Weber
parent a82dd5a483
commit 0a363a108d
2 changed files with 2 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
src/oneD/Sim1D.cpp
View file

@ -331,10 +331,6 @@ void Sim1D::solve(int loglevel, bool refine_grid)
saveResidual("debug_sim1d.xml", "residual",
"After regridding");
}
if (new_points < 0) {
writelog("Maximum number of grid points reached.");
new_points = 0;
}
} else {
debuglog("grid refinement disabled.\n", loglevel);
new_points = 0;

5
src/oneD/refine.cpp
View file

@ -47,8 +47,7 @@ int Refiner::analyze(size_t n, const doublereal* z,
const doublereal* x)
{
if (n >= m_npmax) {
writelog("max number of grid points reached ({}).\n", m_npmax);
return -2;
throw CanteraError("Refiner::analyze", "max number of grid points reached ({}).", m_npmax);
}
if (m_domain->nPoints() <= 1) {
@ -66,7 +65,7 @@ int Refiner::analyze(size_t n, const doublereal* z,
// check consistency
if (n != m_domain->nPoints()) {
throw CanteraError("analyze","inconsistent");
throw CanteraError("Refiner::analyze", "inconsistent");
}
// find locations where cell size ratio is too large.