[Octave-bug-tracker] [bug #45542] griddatan unable to perform 4-dimensional interpolation.

[anonymous]

URL:
  <http://savannah.gnu.org/bugs/?45542>

                 Summary: griddatan unable to perform 4-dimensional
interpolation.
                 Project: GNU Octave
            Submitted by: None
            Submitted on: Mon 13 Jul 2015 02:25:03 PM UTC
                Category: Octave Function
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: Crash
                  Status: None
             Assigned to: None
         Originator Name: Juho Iipponen
        Originator Email: address@hidden
             Open/Closed: Open
         Discussion Lock: Any
                 Release: 4.0.0
        Operating System: GNU/Linux

    _______________________________________________________

Details:

The griddatan function seems to be unable to perform 4-dimensional
interpolation (in e.g. (x, y, z, time)-coordinates). While running the
attached script "test_griddatan.m" with test input data "test_griddatan.mat",
the function fails while attempting to fill in the interpolated values for
output (last line of griddatan.m to be executed for linear interpolation). The
error message is:


error: griddatan: A(I) = X: X must have the same size as I
error: called from
    griddatan at line 84 column 15
    test_griddatan at line 2 column 26


To me, this seems to be caused by different shape of vectors
"y(tri(tri_list,:))" and "bary_list", which are multiplied element-wise. In
this example, the former is of size 5x1, while the latter is of size 1x5. If I
transpose "y(tri(tri_list,:))" before multiplication, the last line yields a
result of 3.1866e-14, which is very close to the one given by Matlab
(3.1858e-14). However, I have no experience with multidimensional geometry, so
I'm not sure if taking the transpose is the correct solution. 



    _______________________________________________________

File Attachments:


-------------------------------------------------------
Date: Mon 13 Jul 2015 02:25:03 PM UTC  Name: test_griddatan.m  Size: 106B  
By: None

<http://savannah.gnu.org/bugs/download.php?file_id=34436>
-------------------------------------------------------
Date: Mon 13 Jul 2015 02:25:03 PM UTC  Name: test_griddatan.mat  Size: 1kB  
By: None

<http://savannah.gnu.org/bugs/download.php?file_id=34437>

    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?45542>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/