prevent refinement from shifting two-point control points

Cantera · Jan 6, 2024 · a93ae70 · a93ae70
1 parent b63b666
commit a93ae70
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/src/oneD/refine.cpp b/src/oneD/refine.cpp
@@ -184,6 +184,12 @@ int Refiner::analyze(size_t n, const double* z, const double* x)
         if (fflame && fflame->isFree() && z[j] == fflame->m_zfixed) {
             m_keep[j] = 1;
         }
+
+        // Keep the point if it is a control point used for two-point flame control
+        if (fflame && fflame->twoPointControlEnabled() && (z[j] == fflame->m_zLeft || z[j] == fflame->m_zRight)) {
+            m_keep[j] = 1;
+        }
+
     }
 
     // Don't allow pruning to remove multiple adjacent grid points