From b0c34f2100761296e95cac7cb9482ba32ab73d23 Mon Sep 17 00:00:00 2001
From: Stanislas Rapacchi <12797255+srapacchi@users.noreply.github.com>
Date: Thu, 4 Sep 2025 14:27:19 +0200
Subject: [PATCH] Fix CalculatekspaceUnfunc with waveforms_and_times

Function gradient_waveforms is deprecated.
---
 matlab/+mr/@Sequence/Sequence.m | 14 ++++++++++++--
 1 file changed, 12 insertions(+), 2 deletions(-)

diff --git a/matlab/+mr/@Sequence/Sequence.m b/matlab/+mr/@Sequence/Sequence.m
index a638aca4..57991f7f 100644
--- a/matlab/+mr/@Sequence/Sequence.m
+++ b/matlab/+mr/@Sequence/Sequence.m
@@ -1481,8 +1481,17 @@ function setBlock(obj, index, varargin)
             end
             
             % now calculate the actual k-space trajectory based on the
-            % gradient waveforms
-            gw=obj.gradient_waveforms();
+            % gradient waveforms (interp on grad raster time)
+            wave_data=obj.waveforms_and_times();
+            tmax=max([wave_data{1}(1,end) wave_data{2}(1,end) wave_data{3}(1,end)]);
+            dt=obj.gradRasterTime; % time raster
+            nt=ceil(tmax/dt);
+            %tmax=nt*dt;
+
+            gw=zeros(3,nt);
+            for i=1:3
+                gw(i,:)=interp1(wave_data{i}(1,:),wave_data{i}(2,:),((1:nt)-0.5)*dt,'linear',0);
+            end
             i_excitation=round(t_excitation/obj.gradRasterTime);
             i_refocusing=round(t_refocusing/obj.gradRasterTime);
 %             ii_next_excitation=min(length(i_excitation),1);
@@ -2807,3 +2816,4 @@ function setExtensionStringAndID(obj, str, id)
         end
     end
 end % classdef
+