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

Merge branch 'delaunay_tree_with_double' into 'master'

Browse Source 
Use new delaunay tree with double

See merge request jschoeberl/netgen!277
This commit is contained in:
Joachim Schöberl 2019-10-15 08:10:19 +00:00
commit 6b99cf3c21
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
libsrc/meshing/delaunay.cpp
View File

@ -3,7 +3,7 @@
namespace netgen
{
template<int dim, typename T=INDEX>
template<int dim, typename T=INDEX, typename TSCAL=double>
class DelaunayTree
{
public:
@ -14,7 +14,7 @@ namespace netgen
struct Leaf
{
Point<2*dim, float> p[N];
Point<2*dim, TSCAL> p[N];
T index[N];
int n_elements;
int nr;
@ -208,7 +208,7 @@ namespace netgen
{
// add two new nodes and one new leaf
int n_elements = leaf->n_elements;
ArrayMem<float, N> coords(n_elements);
ArrayMem<TSCAL, N> coords(n_elements);
ArrayMem<int, N> order(n_elements);
// separate points in two halves, first sort all coordinates in direction dir
@ -282,8 +282,8 @@ namespace netgen
}
};
typedef BoxTree<3> DTREE;
// typedef DelaunayTree<3> DTREE;
// typedef BoxTree<3> DTREE;
typedef DelaunayTree<3> DTREE;
static const int deltetfaces[][3] =
{ { 1, 2, 3 },