Made MElement::scaledJacRange const

Geo/MElement.cpp

@@ -123,7 +123,7 @@ double MElement::metricShapeMeasure2()
   return MetricBasis::boundMinR(this);
 }
 
-double MElement::maxDistToStraight()
+double MElement::maxDistToStraight() const
 {
   const nodalBasis *lagBasis = getFunctionSpace();
   const fullMatrix<double> &uvw = lagBasis->points;
@@ -146,7 +146,7 @@ double MElement::maxDistToStraight()
   return maxdx;
 }
 
-void MElement::scaledJacRange(double &jmin, double &jmax, GEntity *ge)
+void MElement::scaledJacRange(double &jmin, double &jmax, GEntity *ge) const
 {
   jmin = jmax = 1.0;
 #if defined(HAVE_MESH)
@@ -161,7 +161,7 @@ void MElement::scaledJacRange(double &jmin, double &jmax, GEntity *ge)
     SVector3 geoNorm(0.,0.,0.);
     // ... correct Jacobian sign with geometrical normal
     for (int i=0; i<jac->getNumPrimMapNodes(); i++) {
-      MVertex *vert = getVertex(i);
+      const MVertex *vert = getVertex(i);
       if (vert->onWhat() == ge) {
         double u, v;
         vert->getParameter(0,u);

Geo/MElement.h

@@ -176,7 +176,7 @@ class MElement
   virtual double minEdge();
 
   // Max. distance between curved and straight element among all high-order nodes
-  double maxDistToStraight();
+  double maxDistToStraight() const;
 
   // get the quality measures
   double skewness();
@@ -190,7 +190,7 @@ class MElement
     return jmin;
   }
   virtual double angleShapeMeasure() { return 1.0; }
-  virtual void scaledJacRange(double &jmin, double &jmax, GEntity *ge = 0);
+  virtual void scaledJacRange(double &jmin, double &jmax, GEntity *ge = 0) const;
   virtual double metricShapeMeasure();
   virtual double metricShapeMeasure2();