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

Patch for second order UNV tetra node order by Takuya OSHIMA.

(Not sure if this is right. But the new ordering seems to match Code_Aster and SYNOISE?)
parent 48823811
No related branches found
No related tags found
No related merge requests found
...@@ -1225,8 +1225,8 @@ int GModel::readUNV(const std::string &name) ...@@ -1225,8 +1225,8 @@ int GModel::readUNV(const std::string &name)
case 118: case 118:
elements[3][elementary].push_back elements[3][elementary].push_back
(new MTetrahedron10(vertices[0], vertices[2], vertices[4], vertices[9], (new MTetrahedron10(vertices[0], vertices[2], vertices[4], vertices[9],
vertices[1], vertices[3], vertices[5], vertices[8], vertices[1], vertices[3], vertices[5], vertices[6],
vertices[6], vertices[7], num)); vertices[8], vertices[7], num));
dim = 3; dim = 3;
break; break;
case 104: case 115: case 104: case 115:
......
...@@ -205,7 +205,7 @@ class MTetrahedron10 : public MTetrahedron { ...@@ -205,7 +205,7 @@ class MTetrahedron10 : public MTetrahedron {
virtual MVertex *getVertex(int num){ return num < 4 ? _v[num] : _vs[num - 4]; } virtual MVertex *getVertex(int num){ return num < 4 ? _v[num] : _vs[num - 4]; }
virtual MVertex *getVertexUNV(int num) virtual MVertex *getVertexUNV(int num)
{ {
static const int map[10] = {0, 4, 1, 5, 2, 6, 8, 9, 7, 3}; static const int map[10] = {0, 4, 1, 5, 2, 6, 7, 9, 8, 3};
return getVertex(map[num]); return getVertex(map[num]);
} }
virtual MVertex *getVertexBDF(int num) virtual MVertex *getVertexBDF(int num)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment