ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2025-01-11 21:50:34 +05:00
Code Issues Packages Projects Releases Wiki Activity

fix identify periodic

Browse Source
This commit is contained in:
Christopher Lackner 2020-03-17 15:32:42 +01:00
parent b8d313f056
commit ff60ca3f55
2 changed files with 7 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
libsrc/meshing/meshclass.cpp
View File

@ -6118,12 +6118,12 @@ namespace netgen
const Transformation<3> &mapping) const Transformation<3> &mapping)
{ {
auto nr = ident->GetMaxNr() + 1; auto nr = ident->GetMaxNr() + 1;
ident->SetType(nr, Identifications::PERIODIC);
double lami[4]; double lami[4];
GetElementOfPoint({0,0,0}, lami, true);
set<int> identified_points; set<int> identified_points;
Point3d pmin, pmax; Point3d pmin, pmax;
GetBox(pmin, pmax); GetBox(pmin, pmax);
auto eps = 1e-10 * (pmax-pmin).Length(); auto eps = 1e-8 * (pmax-pmin).Length();
for(const auto& se : surfelements) for(const auto& se : surfelements)
{ {
if(GetBCName(se.index-1) != s1) if(GetBCName(se.index-1) != s1)
@ -6131,12 +6131,11 @@ namespace netgen
for(const auto& pi : se.PNums()) for(const auto& pi : se.PNums())
{ {
// cout << "pi = " << pi << endl;
if(identified_points.find(pi) != identified_points.end()) if(identified_points.find(pi) != identified_points.end())
continue; continue;
auto pt = (*this)[pi]; auto pt = (*this)[pi];
auto mapped_pt = mapping(pt); auto mapped_pt = mapping(pt);
auto other_nr = GetElementOfPoint(mapped_pt, lami); auto other_nr = GetElementOfPoint(mapped_pt, lami, true);
int index = -1; int index = -1;
auto other_el = VolumeElement(other_nr); auto other_el = VolumeElement(other_nr);
for(auto i : Range(4)) for(auto i : Range(4))
@ -6150,7 +6149,6 @@ namespace netgen
auto other_pi = other_el.PNums()[index]; auto other_pi = other_el.PNums()[index];
identified_points.insert(pi); identified_points.insert(pi);
ident->Add(pi, other_pi, nr); ident->Add(pi, other_pi, nr);
// cout << "other pi = " << other_pi << endl;
} }
} }
return nr; return nr;

4
libsrc/meshing/python_mesh.cpp
View File

@ -856,6 +856,10 @@ DLL_HEADER void ExportNetgenMeshing(py::module &m)
py::arg("identnr"), py::arg("identnr"),
py::arg("type")) py::arg("type"))
.def("IdentifyPeriodicBoundaries", &Mesh::IdentifyPeriodicBoundaries) .def("IdentifyPeriodicBoundaries", &Mesh::IdentifyPeriodicBoundaries)
.def("GetNrIdentifications", [](Mesh& self)
{
return self.GetIdentifications().GetMaxNr();
})
.def ("CalcLocalH", &Mesh::CalcLocalH) .def ("CalcLocalH", &Mesh::CalcLocalH)
.def ("SetMaxHDomain", [] (Mesh& self, py::list maxhlist) .def ("SetMaxHDomain", [] (Mesh& self, py::list maxhlist)
{ {