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

occ: box, cyl

Browse Source
This commit is contained in:
Joachim Schoeberl 2021-07-24 13:14:21 +02:00
parent d783b71d49
commit 157b246f7f
3 changed files with 28 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
libsrc/occ/occgenmesh.cpp
View File

@ -642,7 +642,8 @@ namespace netgen
Box<3> bb = geom.GetBoundingBox(); Box<3> bb = geom.GetBoundingBox();
int projecttype = PLANESPACE; // int projecttype = PLANESPACE;
int projecttype = PARAMETERSPACE;
static Timer tinit("init"); static Timer tinit("init");
tinit.Start(); tinit.Start();

3
libsrc/occ/occgeom.hpp
View File

@ -290,7 +290,8 @@ namespace netgen
Point<3> & newp, PointGeomInfo & newgi) const override; Point<3> & newp, PointGeomInfo & newgi) const override;
void BuildFMap(); void BuildFMap();
auto GetShape() const { return shape; }
Box<3> GetBoundingBox() const Box<3> GetBoundingBox() const
{ return boundingbox; } { return boundingbox; }

26
libsrc/occ/python_occ.cpp
View File

@ -8,9 +8,15 @@
#include <meshing.hpp> #include <meshing.hpp>
#include <occgeom.hpp> #include <occgeom.hpp>
#include <Standard_Version.hxx> #include <Standard_Version.hxx>
#include <gp_Ax2.hxx>
#include <BRepPrimAPI_MakeSphere.hxx> #include <BRepPrimAPI_MakeSphere.hxx>
#include <BRepPrimAPI_MakeCylinder.hxx>
#include <BRepPrimAPI_MakeBox.hxx>
#include <BOPAlgo_MakerVolume.hxx> #include <BOPAlgo_MakerVolume.hxx>
#include <BOPAlgo_Section.hxx> #include <BOPAlgo_Section.hxx>
#include <BRepAlgoAPI_Cut.hxx>
using namespace netgen; using namespace netgen;
@ -206,6 +212,7 @@ DLL_HEADER void ExportNgOCC(py::module &m)
}, py::arg("mp") = nullptr, }, py::arg("mp") = nullptr,
py::call_guard<py::gil_scoped_release>(), py::call_guard<py::gil_scoped_release>(),
(meshingparameter_description + occparameter_description).c_str()) (meshingparameter_description + occparameter_description).c_str())
.def_property_readonly("shape", [](const OCCGeometry & self) { return self.GetShape(); })
; ;
py::enum_<TopAbs_ShapeEnum>(m, "TopAbs_ShapeEnum", "Enumeration of all supported TopoDS_Shapes") py::enum_<TopAbs_ShapeEnum>(m, "TopAbs_ShapeEnum", "Enumeration of all supported TopoDS_Shapes")
@ -238,14 +245,17 @@ DLL_HEADER void ExportNgOCC(py::module &m)
{ {
// https://dev.opencascade.org/doc/occt-7.3.0/overview/html/occt_user_guides__boolean_operations.html#occt_algorithms_10a // https://dev.opencascade.org/doc/occt-7.3.0/overview/html/occt_user_guides__boolean_operations.html#occt_algorithms_10a
BOPAlgo_MakerVolume aMV; BOPAlgo_MakerVolume aMV;
// BOPAlgo_Section aMV; // BOPAlgo_Section aMV; // only vertices + edges
// BOPAlgo_Builder aBuilder; // BOPAlgo_Builder aMV;
// BRepAlgoAPI_Cut aMV;
TopTools_ListOfShape aLSObjects; TopTools_ListOfShape aLSObjects;
aLSObjects.Append (shape1); aLSObjects.Append (shape1);
aLSObjects.Append (shape2); aLSObjects.Append (shape2);
// aBuilder.SetArguments(aLSObjects); // aBuilder.SetArguments(aLSObjects);
aMV.SetArguments(aLSObjects); aMV.SetArguments(aLSObjects);
// aMV.SetIntersect(true);
aMV.Perform(); aMV.Perform();
// aMV.Build();
return aMV.Shape(); return aMV.Shape();
}); });
; ;
@ -255,6 +265,18 @@ DLL_HEADER void ExportNgOCC(py::module &m)
gp_Pnt cc { py::cast<double> (c[0]), py::cast<double>(c[1]), py::cast<double>(c[2]) }; gp_Pnt cc { py::cast<double> (c[0]), py::cast<double>(c[1]), py::cast<double>(c[2]) };
return BRepPrimAPI_MakeSphere (cc, r).Shape(); return BRepPrimAPI_MakeSphere (cc, r).Shape();
}); });
m.def("Cylinder", [] (py::tuple pnt, py::tuple dir, double r, double h)
{
gp_Pnt cpnt { py::cast<double> (pnt[0]), py::cast<double>(pnt[1]), py::cast<double>(pnt[2]) };
gp_Dir cdir { py::cast<double> (dir[0]), py::cast<double>(dir[1]), py::cast<double>(dir[2]) };
return BRepPrimAPI_MakeCylinder (gp_Ax2(cpnt, cdir), r, h).Shape();
});
m.def("Box", [] (py::tuple p1, py::tuple p2)
{
gp_Pnt cp1 { py::cast<double> (p1[0]), py::cast<double>(p1[1]), py::cast<double>(p1[2]) };
gp_Pnt cp2 { py::cast<double> (p2[0]), py::cast<double>(p2[1]), py::cast<double>(p2[2]) };
return BRepPrimAPI_MakeBox (cp1, cp2).Shape();
});
m.def("LoadOCCGeometry",[] (const string & filename) m.def("LoadOCCGeometry",[] (const string & filename)
{ {