ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-12-24 21:10:33 +05:00
Code Issues Packages Projects Releases Wiki Activity

set correct pointtype

Browse Source
This commit is contained in:
Joachim Schoeberl 2012-07-30 05:27:37 +00:00
parent 14362ac2f0
commit b1bf267ef3
1 changed files with 9 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
libsrc/meshing/secondorder.cpp
View File

@ -47,7 +47,8 @@ namespace netgen
el.epgeominfo[0], el.epgeominfo[1], el.epgeominfo[0], el.epgeominfo[1],
pb, ngi); pb, ngi);
el[2] = mesh.AddPoint (pb); el[2] = mesh.AddPoint (pb, mesh.Point(el[0]).GetLayer(),
EDGEPOINT);
between.Set (i2, el[2]); between.Set (i2, el[2]);
} }
} }
@ -133,7 +134,8 @@ namespace netgen
el.GeomInfoPi (betw[j][1]+1), el.GeomInfoPi (betw[j][1]+1),
pb, newgi); pb, newgi);
newel[onp+j] = mesh.AddPoint (pb); newel[onp+j] = mesh.AddPoint (pb, mesh.Point(pi1).GetLayer(),
SURFACEPOINT);
between.Set (i2, newel[onp+j]); between.Set (i2, newel[onp+j]);
} }
} }
@ -151,7 +153,7 @@ namespace netgen
for (int i = 1; i <= ne; i++) for (int i = 1; i <= ne; i++)
{ {
const Element & el = mesh.VolumeElement(i); const Element & el = mesh.VolumeElement(i);
int onp(0); int onp = 0;
Element newel; Element newel;
newel.SetIndex (el.GetIndex()); newel.SetIndex (el.GetIndex());
@ -213,7 +215,10 @@ namespace netgen
{ {
newel.PNum(onp+1+j) = mesh.AddPoint newel.PNum(onp+1+j) = mesh.AddPoint
(Center (mesh.Point(i2.I1()), (Center (mesh.Point(i2.I1()),
mesh.Point(i2.I2()))); mesh.Point(i2.I2())),
mesh.Point(i2.I1()).GetLayer(),
INNERPOINT);
between.Set (i2, newel.PNum(onp+1+j)); between.Set (i2, newel.PNum(onp+1+j));
} }
} }