ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-11-11 16:49:16 +05:00
Code Issues Packages Projects Releases Wiki Activity

Merge branch 'fixmultielementtrafo' into 'master'

Browse Source 
Fixmultielementtrafo

See merge request !49
This commit is contained in:
Matthias Hochsteger 2017-06-01 17:15:38 +02:00
commit dbee019adf
3 changed files with 25 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
libsrc/csg/python_csg.cpp
View File

@ -441,6 +441,9 @@ DLL_HEADER void ExportCSG(py::module &m)
} }
spsol->AddSurfaces(self); spsol->AddSurfaces(self);
int tlonr = self.SetTopLevelObject(spsol->GetSolid(), surf.get()); int tlonr = self.SetTopLevelObject(spsol->GetSolid(), surf.get());
self.GetTopLevelObject(tlonr) -> SetBCProp(surf->GetBase()->GetBCProperty());
self.GetTopLevelObject(tlonr) -> SetBCName(surf->GetBase()->GetBCName());
self.GetTopLevelObject(tlonr) -> SetMaxH(surf->GetBase()->GetMaxH());
for(auto p : surf->GetPoints()) for(auto p : surf->GetPoints())
self.AddUserPoint(p); self.AddUserPoint(p);
self.AddSplineSurface(surf); self.AddSplineSurface(surf);

2
libsrc/interface/nginterface_v2.cpp
View File

@ -757,7 +757,7 @@ namespace netgen
double hdxdxi[4][3]; double hdxdxi[4][3];
for (int j = 0; j<4;j++) for (int j = 0; j<4;j++)
hxi[j][0] = ((double*)&(xi[0]))[j]; hxi[j][0] = ((double*)&(xi[0]))[j];
MultiElementTransformation<1,3> (elnr, 4, &hxi[0][0], 1, &hx[0][0], 2, &hdxdxi[0][0],4); MultiElementTransformation<1,3> (elnr, 4, &hxi[0][0], 1, &hx[0][0], 3, &hdxdxi[0][0],3);
for(int j=0; j<4; j++) for(int j=0; j<4; j++)
for(int k=0; k<3; k++) for(int k=0; k<3; k++)
((double*)&(x[k]))[j] = hx[j][k]; ((double*)&(x[k]))[j] = hx[j][k];

29
py_tutorials/merge.py
View File

@ -3,28 +3,35 @@ from netgen.csg import *
from ngsolve import ngsglobals from ngsolve import ngsglobals
ngsglobals.msg_level = 2 ngsglobals.msg_level = 2
# generate brick and mesh it
geo1 = CSGeometry() geo1 = CSGeometry()
geo1.Add (OrthoBrick( Pnt(0,0,0), Pnt(1,1,1) )) geo1.Add (OrthoBrick( Pnt(0,0,0), Pnt(1,1,1) ))
m1 = geo1.GenerateMesh (maxh=0.1) m1 = geo1.GenerateMesh (maxh=0.1)
m1.Refine() m1.Refine()
m1.Refine()
# generate sphere and mesh it
geo2 = CSGeometry() geo2 = CSGeometry()
geo2.Add (Sphere (Pnt(0.5,0.5,0.5), 0.1)) geo2.Add (Sphere (Pnt(0.5,0.5,0.5), 0.1))
m2 = geo2.GenerateMesh (maxh=0.05) m2 = geo2.GenerateMesh (maxh=0.05)
m2.Refine() m2.Refine()
m2.Refine() m2.Refine()
print ("********************") print ("***************************")
print ("** start merging **") print ("** merging suface meshes **")
print ("********************") print ("***************************")
# create an empty mesh
mesh = Mesh() mesh = Mesh()
fd_outside = mesh.Add (FaceDescriptor(surfnr=1,domin=1))
fd_inside = mesh.Add (FaceDescriptor(surfnr=2,domin=2,domout=1))
# a face-descriptor stores properties associated with a set of surface elements
# bc .. boundary condition marker,
# domin/domout .. domain-number in front/back of surface elements (0 = void),
# surfnr .. number of the surface described by the face-descriptor
fd_outside = mesh.Add (FaceDescriptor(bc=1,domin=1,surfnr=1))
fd_inside = mesh.Add (FaceDescriptor(bc=2,domin=2,domout=1,surfnr=2))
# copy all boundary points from first mesh to new mesh.
# pmap1 maps point-numbers from old to new mesh
pmap1 = { } pmap1 = { }
for e in m1.Elements2D(): for e in m1.Elements2D():
@ -32,11 +39,17 @@ for e in m1.Elements2D():
if (v not in pmap1): if (v not in pmap1):
pmap1[v] = mesh.Add (m1[v]) pmap1[v] = mesh.Add (m1[v])
# copy surface elements from first mesh to new mesh
# we have to map point-numbers:
for e in m1.Elements2D(): for e in m1.Elements2D():
mesh.Add (Element2D (fd_outside, [pmap1[v] for v in e.vertices])) mesh.Add (Element2D (fd_outside, [pmap1[v] for v in e.vertices]))
# same for the second mesh:
pmap2 = { } pmap2 = { }
for e in m2.Elements2D(): for e in m2.Elements2D():
for v in e.vertices: for v in e.vertices: