Skip to content
Snippets Groups Projects
Commit 68711f27 authored by Amaury Johnen's avatar Amaury Johnen
Browse files

replacing tab for space character (continued) + adding anisotropy measure to MElement

parent dfb7dee1
No related branches found
No related tags found
No related merge requests found
......@@ -25,6 +25,7 @@
#include "Numeric.h"
#include "CondNumBasis.h"
#include "Context.h"
#include "qualityMeasuresJacobian.h"
#define SQU(a) ((a)*(a))
......@@ -233,6 +234,11 @@ double MElement::maxDistToStraight() const
return maxdx;
}
double MElement::minAnisotropyMeasure()
{
return jacobianBasedQuality::minAnisotropyMeasure(this);
}
void MElement::scaledJacRange(double &jmin, double &jmax, GEntity *ge) const
{
jmin = jmax = 1.0;
......@@ -667,7 +673,7 @@ void MElement::getSignedJacobian(fullVector<double> &jacobian, int o) const
void MElement::getNodesCoord(fullMatrix<double> &nodesXYZ) const
{
const int numNodes = getNumShapeFunctions();
const int numNodes = getNumVertices();
for (int i = 0; i < numNodes; i++) {
const MVertex *v = getShapeFunctionNode(i);
nodesXYZ(i,0) = v->x();
......
......@@ -216,6 +216,7 @@ class MElement
signedInvCondNumRange(sICNMin, sICNMax);
return sICNMin;
}
double minAnisotropyMeasure();
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);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment