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

Fixed 1D grid refinement based on curvature requirement

Browse source 
The variable "dz" was being used before being calculated, causing
the algorithm to never insert points based on the curvature criterion.
This commit is contained in:
Ray Speth 2012-04-10 03:13:12 +00:00
parent d8db89e708
commit eb337eac8d
1 changed files with 5 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
src/oneD/refine.cpp
View file

@ -81,12 +81,16 @@ int Refiner::analyze(size_t n, const doublereal* z,
* find locations where cell size ratio is too large. * find locations where cell size ratio is too large.
*/ */
size_t j; size_t j;
vector_fp dz(n-1, 0.0);
string name; string name;
doublereal vmin, vmax, smin, smax, aa, ss; doublereal vmin, vmax, smin, smax, aa, ss;
doublereal dmax, r; doublereal dmax, r;
vector_fp v(n), s(n-1); vector_fp v(n), s(n-1);
vector_fp dz(n-1);
for (j = 0; j < n-1; j++) {
dz[j] = z[j+1] - z[j];
}
for (size_t i = 0; i < m_nv; i++) { for (size_t i = 0; i < m_nv; i++) {
if (m_active[i]) { if (m_active[i]) {
name = m_domain->componentName(i); name = m_domain->componentName(i);
@ -178,9 +182,7 @@ int Refiner::analyze(size_t n, const doublereal* z,
} }
} }
dz[0] = z[1] - z[0];
for (j = 1; j < n-1; j++) { for (j = 1; j < n-1; j++) {
dz[j] = z[j+1] - z[j];
if (dz[j] > m_ratio*dz[j-1]) { if (dz[j] > m_ratio*dz[j-1]) {
m_loc[j] = 1; m_loc[j] = 1;
m_c["point "+int2str(j)] = 1; m_c["point "+int2str(j)] = 1;