From 54ac3ced340f589b03f79259c02c69d8b72c9ca8 Mon Sep 17 00:00:00 2001
From: Michiel Cottaar <michiel.cottaar@ndcn.ox.ac.uk>
Date: Thu, 1 Feb 2024 15:03:21 +0000
Subject: [PATCH] Get excitation pulse to converge

---
 src/helper_functions.jl | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/helper_functions.jl b/src/helper_functions.jl
index b58fb81..c5d0e48 100644
--- a/src/helper_functions.jl
+++ b/src/helper_functions.jl
@@ -48,7 +48,7 @@ function excitation_pulse(; flip_angle=90, phase=0., frequency=0., shape=:sinc,
         return pulse
     end
 
-    grad = TrapezoidGradient(pulse=pulse, rise_time=:min, slice_thickness=[Inf, Inf, min_slice_thickness], orientation=[0, 0, 1.], rotate=rotate_grad)
+    grad = TrapezoidGradient(pulse=pulse, duration=:min, slice_thickness=[Inf, Inf, min_slice_thickness], orientation=[0, 0, 1.], rotate=rotate_grad)
     if !rephase
         return grad
     end
@@ -57,7 +57,7 @@ function excitation_pulse(; flip_angle=90, phase=0., frequency=0., shape=:sinc,
         TrapezoidGradient(orientation=[0, 0, -1.], rotate=rotate_grad, duration=:min);
         TR=Inf
     )
-    @constraint get_global_model() qvec(grad, 1, nothing)[2] == -qvec(seq[2])[2]
+    @constraint get_global_model() qvec(grad, 1, nothing)[3] == -qvec(seq[2])[3]
     return seq
 end
 
-- 
GitLab