From a8dd587eaaca7ca9912ce37adaf9c4e880ea37a3 Mon Sep 17 00:00:00 2001
From: Van Dung Nguyen <vandung.nguyen@ulg.ac.be>
Date: Fri, 30 Mar 2012 10:15:53 +0000
Subject: [PATCH] add third derivatives in Element
---
Geo/MElement.cpp | 11 +++++++++--
Geo/MElement.h | 2 ++
2 files changed, 11 insertions(+), 2 deletions(-)
diff --git a/Geo/MElement.cpp b/Geo/MElement.cpp
index 4c6dbae508..7dff4ec5f2 100644
--- a/Geo/MElement.cpp
+++ b/Geo/MElement.cpp
@@ -142,6 +142,13 @@ void MElement::getHessShapeFunctions(double u, double v, double w, double s[][3]
else Msg::Error("Function space not implemented for this type of element");
}
+void MElement::getThirdDerivativeShapeFunctions(double u, double v, double w, double s[][3][3][3],
+ int o){
+ const polynomialBasis* fs = getFunctionSpace(o);
+ if(fs) fs->dddf(u, v, w, s);
+ else Msg::Error("Function space not implemented for this type of element");
+};
+
SPoint3 MElement::barycenter()
{
SPoint3 p(0., 0., 0.);
@@ -295,7 +302,7 @@ double MElement::getJacobian(double u, double v, double w, double jac[3][3])
jac[j][1] += ver->y() * gg[j];
jac[j][2] += ver->z() * gg[j];
}
- // printf("GSF (%d,%g %g) = %g %g \n",i,u,v,gg[0],gg[1]);
+ // printf("GSF (%d,%g %g) = %g %g \n",i,u,v,gg[0],gg[1]);
}
return _computeDeterminantAndRegularize(this, jac);
}
@@ -1120,7 +1127,7 @@ MElement *MElement::copy(std::map<int, MVertex*> &vertexMap,
if(vertexMap.count(numV))
vmv.push_back(vertexMap[numV]);
else {
- MVertex *mv = new MVertex(v->x(), v->y(), v->z(), 0, numV);
+ MVertex *mv = new MVertex(v->x(), v->y(), v->z(), 0, numV);
vmv.push_back(mv);
vertexMap[numV] = mv;
}
diff --git a/Geo/MElement.h b/Geo/MElement.h
index 4d3cc54a50..f4dc4dfde8 100644
--- a/Geo/MElement.h
+++ b/Geo/MElement.h
@@ -235,6 +235,8 @@ class MElement
int order=-1);
virtual void getHessShapeFunctions(double u, double v, double w, double s[][3][3],
int order=-1);
+ virtual void getThirdDerivativeShapeFunctions(double u, double v, double w, double s[][3][3][3],
+ int order=-1);
// return the Jacobian of the element evaluated at point (u,v,w) in
// parametric coordinates
double getJacobian(const fullMatrix<double> &gsf, double jac[3][3]);
--
GitLab