testing a special measure improving hex-dominant meshes

Geo/MElement.cpp

@@ -243,6 +243,23 @@ double MElement::minAnisotropyMeasure()
 #endif
 }
 
+double MElement::specialQuality()
+{
+#if defined(HAVE_MESH)
+  double minJ, maxJ;
+  jacobianBasedQuality::minMaxJacobianDeterminant(this, minJ, maxJ);
+  if (minJ < 0 && maxJ >= 0) { // accept -inf as a answer
+    return minJ/maxJ;
+  }
+  if (minJ < 0 && maxJ < 0) {
+    return -std::numeric_limits<double>::infinity();
+  }
+  return jacobianBasedQuality::minAnisotropyMeasure(this);
+#else
+  return 0.;
+#endif
+}
+
 void MElement::scaledJacRange(double &jmin, double &jmax, GEntity *ge) const
 {
   jmin = jmax = 1.0;

Geo/MElement.h

@@ -217,6 +217,7 @@ class MElement
     return sICNMin;
   }
   double minAnisotropyMeasure();
+  double specialQuality();
   virtual double angleShapeMeasure() { return 1.0; }
   virtual void scaledJacRange(double &jmin, double &jmax, GEntity *ge = 0) const;
   virtual void idealJacRange(double &jmin, double &jmax, GEntity *ge = 0);