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

parallel update

Browse Source
This commit is contained in:
Joachim Schöberl 2016-02-28 05:32:08 +01:00
parent 346e3db54d
commit 9ef2028729
1 changed files with 10 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
libsrc/meshing/refine.cpp
View File

@ -17,6 +17,12 @@ namespace netgen
mesh.ComputeNVertices(); mesh.ComputeNVertices();
mesh.SetNP(mesh.GetNV()); mesh.SetNP(mesh.GetNV());
if (mesh.mlbetweennodes.Size() < mesh.GetNV())
{
mesh.mlbetweennodes.SetSize(mesh.GetNV());
mesh.mlbetweennodes = INDEX_2(PointIndex::BASE-1,PointIndex::BASE-1);
}
INDEX_2_HASHTABLE<PointIndex> between(mesh.GetNP() + 5); INDEX_2_HASHTABLE<PointIndex> between(mesh.GetNP() + 5);
@ -70,7 +76,10 @@ namespace netgen
QuickSort (parents, par_nr); QuickSort (parents, par_nr);
for (int i = 0; i < parents.Size(); i++) for (int i = 0; i < parents.Size(); i++)
between.Set (parents[i], mesh.GetNV()+i+PointIndex::BASE); {
between.Set (parents[i], mesh.GetNV()+i+PointIndex::BASE);
mesh.mlbetweennodes[mesh.GetNV()+i+PointIndex::BASE] = parents[i];
}
mesh.SetNP(mesh.GetNV() + parents.Size()); mesh.SetNP(mesh.GetNV() + parents.Size());
Array<bool, PointIndex::BASE> pointset(mesh.GetNP()); Array<bool, PointIndex::BASE> pointset(mesh.GetNP());