ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2025-01-27 13:20:34 +05:00
Code Issues Packages Projects Releases Wiki Activity

don't invert

Browse Source
This commit is contained in:
Joachim Schoeberl 2011-04-19 16:04:07 +00:00
parent 162293dde7
commit b1bd297775
1 changed files with 15 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
libsrc/meshing/meshclass.cpp
View File

@ -863,13 +863,21 @@ namespace netgen
PrintMessage (3, n, " surface elements"); PrintMessage (3, n, " surface elements");
for (i = 1; i <= n; i++) for (i = 1; i <= n; i++)
{ {
int j;
int surfnr, bcp, domin, domout, nep, faceind = 0; int surfnr, bcp, domin, domout, nep, faceind = 0;
infile >> surfnr >> bcp >> domin >> domout; infile >> surfnr >> bcp >> domin >> domout;
surfnr--; surfnr--;
for (j = 1; j <= facedecoding.Size(); j++) bool invert_el = false;
/*
if (domin == 0)
{
invert_el = true;
Swap (domin, domout);
}
*/
for (int j = 1; j <= facedecoding.Size(); j++)
if (GetFaceDescriptor(j).SurfNr() == surfnr && if (GetFaceDescriptor(j).SurfNr() == surfnr &&
GetFaceDescriptor(j).BCProperty() == bcp && GetFaceDescriptor(j).BCProperty() == bcp &&
GetFaceDescriptor(j).DomainIn() == domin && GetFaceDescriptor(j).DomainIn() == domin &&
@ -888,19 +896,21 @@ namespace netgen
Element2d tri(nep); Element2d tri(nep);
tri.SetIndex(faceind); tri.SetIndex(faceind);
for (j = 1; j <= nep; j++) for (int j = 1; j <= nep; j++)
infile >> tri.PNum(j); infile >> tri.PNum(j);
if (strcmp (str, "surfaceelementsgi") == 0) if (strcmp (str, "surfaceelementsgi") == 0)
for (j = 1; j <= nep; j++) for (int j = 1; j <= nep; j++)
infile >> tri.GeomInfoPi(j).trignum; infile >> tri.GeomInfoPi(j).trignum;
if (strcmp (str, "surfaceelementsuv") == 0) if (strcmp (str, "surfaceelementsuv") == 0)
for (j = 1; j <= nep; j++) for (int j = 1; j <= nep; j++)
infile >> tri.GeomInfoPi(j).u >> tri.GeomInfoPi(j).v; infile >> tri.GeomInfoPi(j).u >> tri.GeomInfoPi(j).v;
if (invertsurf) if (invertsurf)
tri.Invert(); tri.Invert();
if (invert_el)
tri.Invert();
AddSurfaceElement (tri); AddSurfaceElement (tri);
} }