using interpn for octave in the photon dose calculation significantly…

… boosts performance
HomolkaN · Sep 17, 2021 · 1db68c1 · 1db68c1
1 parent 11759ed
commit 1db68c1
Showing 1 changed file with 10 additions and 6 deletions.
diff --git a/matRad_calcPhotonDose.m b/matRad_calcPhotonDose.m
@@ -202,9 +202,11 @@
             Interp_kernel2 = griddedInterpolant(convMx_X',convMx_Z',convMx2','linear','none');
             Interp_kernel3 = griddedInterpolant(convMx_X',convMx_Z',convMx3','linear','none');
         else
-            Interp_kernel1 = @(x,y)interp2(convMx_X(1,:),convMx_Z(:,1),convMx1,x,y,'linear',NaN);
-            Interp_kernel2 = @(x,y)interp2(convMx_X(1,:),convMx_Z(:,1),convMx2,x,y,'linear',NaN);
-            Interp_kernel3 = @(x,y)interp2(convMx_X(1,:),convMx_Z(:,1),convMx3,x,y,'linear',NaN);
+            %For some reason the use of interpn here is much faster
+            %than using interp2 in Octave
+            Interp_kernel1 = @(x,y)interpn(convMx_X(1,:),convMx_Z(:,1),convMx1,x,y,'linear',NaN);
+            Interp_kernel2 = @(x,y)interpn(convMx_X(1,:),convMx_Z(:,1),convMx2,x,y,'linear',NaN);
+            Interp_kernel3 = @(x,y)interpn(convMx_X(1,:),convMx_Z(:,1),convMx3,x,y,'linear',NaN);
         end
     end
 
@@ -243,9 +245,11 @@
                 Interp_kernel2 = griddedInterpolant(convMx_X',convMx_Z',convMx2','linear','none');
                 Interp_kernel3 = griddedInterpolant(convMx_X',convMx_Z',convMx3','linear','none');
             else
-                Interp_kernel1 = @(x,y)interp2(convMx_X(1,:),convMx_Z(:,1),convMx1,x,y,'linear',NaN);
-                Interp_kernel2 = @(x,y)interp2(convMx_X(1,:),convMx_Z(:,1),convMx2,x,y,'linear',NaN);
-                Interp_kernel3 = @(x,y)interp2(convMx_X(1,:),convMx_Z(:,1),convMx3,x,y,'linear',NaN);
+                %For some reason the use of interpn here is much faster
+                %than using interp2 in Octave
+                Interp_kernel1 = @(x,y)interpn(convMx_X(1,:),convMx_Z(:,1),convMx1,x,y,'linear',NaN);
+                Interp_kernel2 = @(x,y)interpn(convMx_X(1,:),convMx_Z(:,1),convMx2,x,y,'linear',NaN);
+                Interp_kernel3 = @(x,y)interpn(convMx_X(1,:),convMx_Z(:,1),convMx3,x,y,'linear',NaN);
             end
 
         end