[Octave-bug-tracker] [bug #60818] delaunayn - 2D code path vectorization doesn't match nD algorithm

[Nicholas Jankowski]

Follow-up Comment #48, bug #60818 (project octave):

We could arbitrarily split up the vector equations by simplexes, but in what
way? half? quarter? When we set out I was wondering if there was a way to loop
over the dimensions. With volume calculation being the main task I figured
some aspect of iteration increasing the volume projection would work.  not
something I could see with the LU decomp though.

I'm still not sure why the method does so poorly with low point count, but
other than that, with the tweaked version that removes calls to 'dot' and
'cross' which is fairly quick anyway, the laplace methed is faster with lower
memory up to and including 6D.  So if LU is just used for degree > 6, maybe
it's natural to partition by degree number. 

    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?60818>

_______________________________________________
  Message sent via Savannah
  https://savannah.gnu.org/