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

set occparams and generate triangulation with BRepMesh

Browse Source
This commit is contained in:
Christoph Wintersteiger 2017-10-20 16:39:37 +02:00
parent 8f11b84911
commit a90ee2e437
2 changed files with 16 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
libsrc/occ/occgenmesh.cpp
View File

@ -1137,7 +1137,12 @@ namespace netgen
Handle(Geom_Surface) surf = BRep_Tool::Surface (face); Handle(Geom_Surface) surf = BRep_Tool::Surface (face);
Handle(Poly_Triangulation) triangulation = BRep_Tool::Triangulation (face, loc); Handle(Poly_Triangulation) triangulation = BRep_Tool::Triangulation (face, loc);
if (triangulation.IsNull()) continue; if (triangulation.IsNull())
{
BRepTools::Clean (geom.shape);
BRepMesh_IncrementalMesh (geom.shape, 0.01, true);
triangulation = BRep_Tool::Triangulation (face, loc);
}
BRepAdaptor_Surface sf(face, Standard_True); BRepAdaptor_Surface sf(face, Standard_True);
BRepLProp_SLProps prop(sf, 2, 1e-5); BRepLProp_SLProps prop(sf, 2, 1e-5);
@ -1286,6 +1291,10 @@ namespace netgen
int OCCGenerateMesh (OCCGeometry & geom, shared_ptr<Mesh> & mesh, MeshingParameters & mparam) int OCCGenerateMesh (OCCGeometry & geom, shared_ptr<Mesh> & mesh, MeshingParameters & mparam)
{ {
cout << "MeshingParameters: " << mparam << endl;
cout << "OCCParameters: " << endl;
occparam.Print(cout);
multithread.percent = 0; multithread.percent = 0;
if (mparam.perfstepsstart <= MESHCONST_ANALYSE) if (mparam.perfstepsstart <= MESHCONST_ANALYSE)

6
libsrc/occ/python_occ.cpp
View File

@ -33,6 +33,12 @@ DLL_HEADER void ExportNgOCC(py::module &m)
SetGlobalMesh(mesh); SetGlobalMesh(mesh);
mesh->SetGeometry(geo); mesh->SetGeometry(geo);
ng_geometry = geo; ng_geometry = geo;
occparam.resthcloseedgefac = 2.0;
occparam.resthcloseedgeenable = 1;
occparam.resthminedgelen = 0.2;
occparam.resthminedgelenenable = 1;
try try
{ {
geo->GenerateMesh(mesh,param); geo->GenerateMesh(mesh,param);