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

[Nicholas Jankowski]

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

so all that said:

biggest memory jump seems to happen inside lu.  no real way to work around
that. 

laplace does scale worse than lu at high dimension and high number of points,
but seems laplace still edges it out to about 5D (ignoring hte small 10pt
case), after which laplace scaling gets pretty poor. 

just based on timing, I'd recommend using laplace expansion up to 5D, then
switch to the LU approach.  unfortunately that will run into memory issues for
>5D. I don't really know any way around that for those running larger, higher
dimension triangulations. Personally, I don't know anyone running such high
dimension triangulations. 

Attached is the m file which would do that. Haven't worked up a patch yet, as
this should still get some BISTs added to it. 



(file #51724)
    _______________________________________________________

Additional Item Attachment:

File name: delaunayn.m                    Size:8 KB
    <https://file.savannah.gnu.org/file/delaunayn.m?file_id=51724>



    _______________________________________________________

Reply to this item at:

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

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