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

[Nicholas Jankowski]

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

A good question. In the extreme each volume’s eqn  diagonal could be
independently calculated, that’s what the for loop did. And I’m pretty
sure  the blocks would independently diagonalize (since you can use matrix row
operations instead of lu decomp to get the diagonal, and row ops to get row
echelon form are block independent), but not sure how to partition that in any
non-arbitrary way.

I’m half convinced this whole problem stems from not making use of what’s
already going on in the initial Delaunay calculation, and we’re doubling
effort for an after the fact check.

    _______________________________________________________

Reply to this item at:

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

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