Skip to content
Snippets Groups Projects
Commit e8a34116 authored by Christophe Geuzaine's avatar Christophe Geuzaine
Browse files

pp

parent 28f704ce
No related branches found
No related tags found
No related merge requests found
...@@ -168,12 +168,13 @@ void MElement::getHessShapeFunctions(double u, double v, double w, double s[][3] ...@@ -168,12 +168,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"); 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], void MElement::getThirdDerivativeShapeFunctions(double u, double v, double w,
int o){ double s[][3][3][3], int o)
{
const polynomialBasis* fs = getFunctionSpace(o); const polynomialBasis* fs = getFunctionSpace(o);
if(fs) fs->dddf(u, v, w, s); if(fs) fs->dddf(u, v, w, s);
else Msg::Error("Function space not implemented for this type of element"); else Msg::Error("Function space not implemented for this type of element");
}; }
SPoint3 MElement::barycenter_infty () SPoint3 MElement::barycenter_infty ()
{ {
...@@ -487,6 +488,15 @@ void MElement::xyz2uvw(double xyz[3], double uvw[3]) ...@@ -487,6 +488,15 @@ void MElement::xyz2uvw(double xyz[3], double uvw[3])
} }
} }
void MElement::xyzTouvw(fullMatrix<double> *xu)
{
double _xyz[3] = {(*xu)(0,0),(*xu)(0,1),(*xu)(0,2)}, _uvw[3];
xyz2uvw(_xyz, _uvw);
(*xu)(1,0) = _uvw[0];
(*xu)(1,1) = _uvw[1];
(*xu)(1,2) = _uvw[2];
}
void MElement::movePointFromParentSpaceToElementSpace(double &u, double &v, double &w) void MElement::movePointFromParentSpaceToElementSpace(double &u, double &v, double &w)
{ {
if(!getParent()) return; if(!getParent()) return;
...@@ -971,6 +981,7 @@ void MElement::writeSTL(FILE *fp, bool binary, double scalingFactor) ...@@ -971,6 +981,7 @@ void MElement::writeSTL(FILE *fp, bool binary, double scalingFactor)
} }
} }
} }
void MElement::writePLY2(FILE *fp) void MElement::writePLY2(FILE *fp)
{ {
setVolumePositive(); setVolumePositive();
...@@ -1453,12 +1464,3 @@ MElement *MElementFactory::create(int num, int type, const std::vector<int> &dat ...@@ -1453,12 +1464,3 @@ MElement *MElementFactory::create(int num, int type, const std::vector<int> &dat
return element; return element;
} }
void MElement::xyzTouvw(fullMatrix<double> *xu)
{
double _xyz[3] = {(*xu)(0,0),(*xu)(0,1),(*xu)(0,2)}, _uvw[3];
xyz2uvw(_xyz, _uvw);
(*xu)(1,0) = _uvw[0];
(*xu)(1,1) = _uvw[1];
(*xu)(1,2) = _uvw[2];
}
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment