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

disable parallel vertex2element table, needs sorting

Browse Source
This commit is contained in:
Joachim Schöberl 2019-09-09 10:55:27 +02:00
parent 5ffab34d13
commit cbd129bf53
1 changed files with 8 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
libsrc/meshing/topology.cpp
View File

@ -394,14 +394,14 @@ namespace netgen
}); });
vert2element = TABLE<ElementIndex,PointIndex::BASE> (cnt); vert2element = TABLE<ElementIndex,PointIndex::BASE> (cnt);
/*
for (ElementIndex ei = 0; ei < ne; ei++) for (ElementIndex ei = 0; ei < ne; ei++)
{ {
const Element & el = (*mesh)[ei]; const Element & el = (*mesh)[ei];
for (int j = 0; j < el.GetNV(); j++) for (int j = 0; j < el.GetNV(); j++)
vert2element.AddSave (el[j], ei); vert2element.AddSave (el[j], ei);
} }
*/
/*
ParallelForRange ParallelForRange
(tm, ne, (tm, ne,
[&] (size_t begin, size_t end) [&] (size_t begin, size_t end)
@ -413,6 +413,8 @@ namespace netgen
vert2element.ParallelAdd (el[j], ei); vert2element.ParallelAdd (el[j], ei);
} }
}); });
requires sorting !!!!
*/
cnt = 0; cnt = 0;
/* /*
@ -438,14 +440,14 @@ namespace netgen
vert2surfelement = TABLE<SurfaceElementIndex,PointIndex::BASE> (cnt); vert2surfelement = TABLE<SurfaceElementIndex,PointIndex::BASE> (cnt);
/*
for (SurfaceElementIndex sei = 0; sei < nse; sei++) for (SurfaceElementIndex sei = 0; sei < nse; sei++)
{ {
const Element2d & el = (*mesh)[sei]; const Element2d & el = (*mesh)[sei];
for (int j = 0; j < el.GetNV(); j++) for (int j = 0; j < el.GetNV(); j++)
vert2surfelement.AddSave (el[j], sei); vert2surfelement.AddSave (el[j], sei);
} }
*/ /*
ParallelForRange ParallelForRange
(tm, nse, (tm, nse,
[&] (size_t begin, size_t end) [&] (size_t begin, size_t end)
@ -457,7 +459,8 @@ namespace netgen
vert2surfelement.ParallelAdd (el[j], sei); vert2surfelement.ParallelAdd (el[j], sei);
} }
}); });
requires sorting !!!
*/
cnt = 0; cnt = 0;
for (SegmentIndex si = 0; si < nseg; si++) for (SegmentIndex si = 0; si < nseg; si++)