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

Layer parameter for RestrictH

Browse Source
This commit is contained in:
Matthias Hochsteger 2025-03-10 10:05:11 +01:00
parent 7b13db740d
commit 3b79dbc8ff
4 changed files with 12 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
libsrc/csg/genmesh.cpp
View File

@ -718,7 +718,7 @@ namespace netgen
mesh -> LoadLocalMeshSize (mparam.meshsizefilename);
for (auto mspnt : mparam.meshsize_points)
mesh -> RestrictLocalH (mspnt.pnt, mspnt.h);
mesh -> RestrictLocalH (mspnt.pnt, mspnt.h, mspnt.layer);
}
spoints.SetSize(0);

2
libsrc/meshing/basegeom.cpp
View File

@ -474,7 +474,7 @@ namespace netgen
}
for(const auto& mspnt : mparam.meshsize_points)
mesh.RestrictLocalH(mspnt.pnt, mspnt.h);
mesh.RestrictLocalH(mspnt.pnt, mspnt.h, mspnt.layer);
mesh.LoadLocalMeshSize(mparam.meshsizefilename);
}

2
libsrc/meshing/meshtype.hpp
View File

@ -1635,7 +1635,7 @@ namespace netgen
Point<3> pnt;
double h;
int layer = 1;
MeshSizePoint (Point<3> _pnt, double _h) : pnt(_pnt), h(_h) { ; }
MeshSizePoint (Point<3> pnt_, double h_, int layer_ = 1) : pnt(pnt_), h(h_), layer(layer_) { ; }
MeshSizePoint () = default;
MeshSizePoint (const MeshSizePoint &) = default;
MeshSizePoint (MeshSizePoint &&) = default;

18
libsrc/meshing/python_mesh.cpp
View File

@ -1683,25 +1683,25 @@ py::arg("point_tolerance") = -1.)
return mp;
}), py::arg("mp")=nullptr, meshingparameter_description.c_str())
.def("__str__", &ToString<MP>)
.def("RestrictH", [](MP & mp, double x, double y, double z, double h)
.def("RestrictH", [](MP & mp, double x, double y, double z, double h, int layer)
{
mp.meshsize_points.Append ( MeshingParameters::MeshSizePoint(Point<3> (x,y,z), h));
}, py::arg("x"), py::arg("y"), py::arg("z"), py::arg("h")
mp.meshsize_points.Append ( MeshingParameters::MeshSizePoint(Point<3> (x,y,z), h, layer));
}, py::arg("x"), py::arg("y"), py::arg("z"), py::arg("h"), py::arg("layer")=1
)
.def("RestrictH", [](MP & mp, const Point<3>& p, double h)
.def("RestrictH", [](MP & mp, const Point<3>& p, double h, int layer)
{
mp.meshsize_points.Append ({p, h});
}, py::arg("p"), py::arg("h"))
mp.meshsize_points.Append ({p, h, layer});
}, py::arg("p"), py::arg("h"), py::arg("layer")=1)
.def("RestrictHLine", [](MP& mp, const Point<3>& p1, const Point<3>& p2,
double maxh)
double maxh, int layer)
{
int steps = int(Dist(p1, p2) / maxh) + 2;
auto v = p2 - p1;
for (int i = 0; i <= steps; i++)
{
mp.meshsize_points.Append({p1 + double(i)/steps * v, maxh});
mp.meshsize_points.Append({p1 + double(i)/steps * v, maxh, layer});
}
}, py::arg("p1"), py::arg("p2"), py::arg("maxh"))
}, py::arg("p1"), py::arg("p2"), py::arg("maxh"), py::arg("layer")=1)
;
m.def("SetTestoutFile", FunctionPointer ([] (const string & filename)