changeset 27239:bee80e27dcc5

griddatan.m: Fix 4-D interpolation (bug #45542, bug #56515). * griddatan.m: Only use bary_list in calculations if it is not empty.
author Juho Iipponen <juho.iipponen@helsinki.fi>
date Thu, 11 Jul 2019 20:16:58 -0700
parents 177be3c01238
children 336c640c481b
files scripts/geometry/griddatan.m
diffstat 1 files changed, 3 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/scripts/geometry/griddatan.m	Thu Jul 11 19:55:53 2019 -0700
+++ b/scripts/geometry/griddatan.m	Thu Jul 11 20:16:58 2019 -0700
@@ -81,7 +81,9 @@
     yt = y(tri(tri_list,:));
 
     ## Use barycentric coordinate of point to calculate yi
-    yi(valid) = sum (y(tri(tri_list,:)) .* bary_list, 2);
+    if (! isempty (bary_list))
+      yi(valid) = sum (y(tri(tri_list,:))' .* bary_list, 2);
+    endif
 
   else
     error ("griddatan: unknown interpolation METHOD");