From 753b8d52d8f64147285120f1dcbc9ba629fb74f6 Mon Sep 17 00:00:00 2001
From: Amaury Johnan <amjohnen@gmail.com>
Date: Tue, 26 Apr 2016 11:52:49 +0000
Subject: [PATCH] testing a special measure improving hex-dominant meshes
---
Geo/MElement.cpp | 17 +++++++++++++++++
Geo/MElement.h | 1 +
2 files changed, 18 insertions(+)
diff --git a/Geo/MElement.cpp b/Geo/MElement.cpp
index bb9ecbe93a..49f377c3c4 100644
--- a/Geo/MElement.cpp
+++ b/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;
diff --git a/Geo/MElement.h b/Geo/MElement.h
index e504498a8b..21ef1d7afe 100644
--- a/Geo/MElement.h
+++ b/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);
--
GitLab