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

ellipsoid

Browse Source
This commit is contained in:
Joachim Schoeberl 2023-06-08 16:38:34 +03:00
parent ef393472b2
commit ec9d028c60
2 changed files with 32 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
libsrc/occ/python_occ_basic.cpp
View File

@ -151,6 +151,11 @@ DLL_HEADER void ExportNgOCCBasic(py::module &m)
.def(py::init([](gp_Pnt p, gp_Dir N, gp_Dir Vx) { .def(py::init([](gp_Pnt p, gp_Dir N, gp_Dir Vx) {
return gp_Ax3(p,N, Vx); return gp_Ax3(p,N, Vx);
}), py::arg("p")=gp_Pnt(0,0,0), py::arg("n")=gp_Vec(0,0,1), py::arg("h")=gp_Vec(1,0,0)) }), py::arg("p")=gp_Pnt(0,0,0), py::arg("n")=gp_Vec(0,0,1), py::arg("h")=gp_Vec(1,0,0))
.def(py::init([](gp_Ax1 ax1) {
gp_Ax3 ax;
ax.SetAxis(ax1);
return ax;
}), py::arg("axis"))
.def(py::init<gp_Ax2>()) .def(py::init<gp_Ax2>())
.def_property("p", [](gp_Ax3 & ax) { return ax.Location(); }, [](gp_Ax3&ax, gp_Pnt p) { ax.SetLocation(p); }) .def_property("p", [](gp_Ax3 & ax) { return ax.Location(); }, [](gp_Ax3&ax, gp_Pnt p) { ax.SetLocation(p); })
; ;
@ -364,6 +369,9 @@ DLL_HEADER void ExportNgOCCBasic(py::module &m)
py::implicitly_convertible<py::tuple, gp_Vec2d>(); py::implicitly_convertible<py::tuple, gp_Vec2d>();
py::implicitly_convertible<py::tuple, gp_Dir2d>(); py::implicitly_convertible<py::tuple, gp_Dir2d>();
py::implicitly_convertible<gp_Ax1, gp_Ax3>();
py::implicitly_convertible<gp_Ax3, gp_Ax1>();
py::implicitly_convertible<gp_Ax3, gp_Ax2>(); py::implicitly_convertible<gp_Ax3, gp_Ax2>();

24
libsrc/occ/python_occ_shapes.cpp
View File

@ -21,6 +21,7 @@
#include <BRepBuilderAPI_MakeVertex.hxx> #include <BRepBuilderAPI_MakeVertex.hxx>
#include <BRepBuilderAPI_MakeWire.hxx> #include <BRepBuilderAPI_MakeWire.hxx>
#include <BRepBuilderAPI_Transform.hxx> #include <BRepBuilderAPI_Transform.hxx>
#include <BRepBuilderAPI_GTransform.hxx>
#include <BRepExtrema_DistShapeShape.hxx> #include <BRepExtrema_DistShapeShape.hxx>
#include <BRepFilletAPI_MakeChamfer.hxx> #include <BRepFilletAPI_MakeChamfer.hxx>
#include <BRepFilletAPI_MakeFillet.hxx> #include <BRepFilletAPI_MakeFillet.hxx>
@ -1815,10 +1816,33 @@ DLL_HEADER void ExportNgOCCShapes(py::module &m)
BRepPrimAPI_MakeHalfSpace builder(face, refpnt); BRepPrimAPI_MakeHalfSpace builder(face, refpnt);
return builder.Shape(); return builder.Shape();
}, py::arg("p"), py::arg("n"), "Create a half space threw point p normal to n"); }, py::arg("p"), py::arg("n"), "Create a half space threw point p normal to n");
m.def("Sphere", [] (gp_Pnt cc, double r) { m.def("Sphere", [] (gp_Pnt cc, double r) {
return BRepPrimAPI_MakeSphere (cc, r).Solid(); return BRepPrimAPI_MakeSphere (cc, r).Solid();
}, py::arg("c"), py::arg("r"), "create sphere with center 'c' and radius 'r'"); }, py::arg("c"), py::arg("r"), "create sphere with center 'c' and radius 'r'");
m.def("Ellipsoid", [] (gp_Ax3 ax, double r1, double r2, optional<double> hr3) {
auto sp = BRepPrimAPI_MakeSphere (gp_Pnt(0,0,0), 1).Solid();
gp_GTrsf gtrafo;
double r3 = hr3.value_or(r2);
gtrafo.SetVectorialPart({ r2, 0, 0, 0, r3, 0, 0, 0, r1 });
gtrafo.SetTranslationPart( { 0.0, 0.0, 0.0 } );
BRepBuilderAPI_GTransform gbuilder(sp, gtrafo, true);
PropagateProperties(gbuilder, sp, occ2ng(gtrafo));
auto gsp = gbuilder.Shape();
gp_Trsf trafo;
trafo.SetTransformation(ax, gp_Ax3());
BRepBuilderAPI_Transform builder(gsp, trafo, true);
PropagateProperties(builder, gsp, occ2ng(trafo));
return builder.Shape();
}, py::arg("axes"), py::arg("r1"), py::arg("r2"), py::arg("r3")=std::nullopt,
"create ellipsoid with local coordinates given by axes, radi 'r1', 'r2', 'r3'");
m.def("Cylinder", [] (gp_Pnt cpnt, gp_Dir cdir, double r, double h, m.def("Cylinder", [] (gp_Pnt cpnt, gp_Dir cdir, double r, double h,
optional<string> bot, optional<string> top, optional<string> mantle) { optional<string> bot, optional<string> top, optional<string> mantle) {
auto builder = BRepPrimAPI_MakeCylinder (gp_Ax2(cpnt, cdir), r, h); auto builder = BRepPrimAPI_MakeCylinder (gp_Ax2(cpnt, cdir), r, h);