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

mesh.AddElements also for 1D segments

Browse Source
This commit is contained in:
Joachim Schoeberl 2023-05-07 17:45:20 +02:00
parent 8a79f91cb1
commit 9c0eeb5370
1 changed files with 16 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
libsrc/meshing/python_mesh.cpp
View File

@ -1075,6 +1075,21 @@ DLL_HEADER void ExportNetgenMeshing(py::module &m)
// throw std::runtime_error("AddPoints needs buffer of type int"); // throw std::runtime_error("AddPoints needs buffer of type int");
int * ptr = static_cast<int*> (info.ptr); int * ptr = static_cast<int*> (info.ptr);
if (dim == 1)
{
ELEMENT_TYPE type;
int np = info.shape[1];
self.LineSegments().SetAllocSize(self.LineSegments().Size()+info.shape[0]);
for (auto i : Range(info.shape[0]))
{
Segment el;
for (int j = 0; j < np; j++)
el[j] = ptr[j]+PointIndex::BASE-base;
el.si = index;
self.AddSegment(el);
ptr += info.strides[0]/sizeof(int);
}
}
if (dim == 2) if (dim == 2)
{ {
ELEMENT_TYPE type; ELEMENT_TYPE type;
@ -1092,7 +1107,7 @@ DLL_HEADER void ExportNetgenMeshing(py::module &m)
for (auto i : Range(info.shape[0])) for (auto i : Range(info.shape[0]))
{ {
Element2d el(type); Element2d el(type);
for (int j = 0; j < np;j ++) for (int j = 0; j < np; j++)
el[j] = ptr[j]+PointIndex::BASE-base; el[j] = ptr[j]+PointIndex::BASE-base;
el.SetIndex(index); el.SetIndex(index);
self.AddSurfaceElement (el); self.AddSurfaceElement (el);