Fix: Importer - treat gradient delays correctly in combination with RF

CHANGELOG.md

@@ -1,3 +1,5 @@
+- 1.0.9
+  - Fix pulseq importer: correct handling of gradients in RF block that are not slice selection gradients
 - 1.0.8
   - `util.pulseq_plot` defaults to a stacked six-row figure (`stacked=True`) on pypulseq 1.5, 1.4 and pre-1.4
   - stacked plots use a doubled figure height so the six rows stay readable

Cargo.toml

@@ -1,6 +1,6 @@
 [package]
 name = "mrzero_core"
-version = "1.0.8"
+version = "1.0.9"
 edition = "2021"
 
 # See more keys and their definitions at https://doc.rust-lang.org/cargo/reference/manifest.html

pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "maturin"
 
 [project]
 name = "mrzerocore"
-version = "1.0.8"
+version = "1.0.9"
 description = "Core functionality of MRzero"
 authors = [
     {name = "Jonathan Endres", email = "jonathan.endres@uk-erlangen.de"},

python/MRzeroCore/sequence.py

@@ -944,14 +944,23 @@ def integrate_axis(axis: str, t0: float, t1: float) -> float:
             # - .sample(t): return the amplitude at one time point
             # - .integrate(t0, t1): return the gradient moment over the given period
             if blocks[b_idx].gx is not None:
-                rf_hasgrad = True
-                rf_grad[0] = blocks[b_idx].gx.amp
+                if blocks[b_idx].gx.delay > rf.shape_duration + rf.delay:
+                    pass # no graident present during rf 
+                else:
+                    rf_hasgrad = True
+                    rf_grad[0] = blocks[b_idx].gx.amp
             if blocks[b_idx].gy is not None:
-                rf_hasgrad = True
-                rf_grad[1] = blocks[b_idx].gy.amp
+                if blocks[b_idx].gy.delay > rf.shape_duration + rf.delay:
+                    pass # no graident present during rf 
+                else:
+                    rf_hasgrad = True
+                    rf_grad[1] = blocks[b_idx].gy.amp
             if blocks[b_idx].gz is not None:
-                rf_hasgrad = True
-                rf_grad[2] = blocks[b_idx].gz.amp
+                if blocks[b_idx].gz.delay > rf.shape_duration + rf.delay:
+                    pass # no graident present during rf 
+                else:
+                    rf_hasgrad = True
+                    rf_grad[2] = blocks[b_idx].gz.amp
 
             # Shim handling: pulseq-rs always returns a list, with [(1.0, 0.0)]
             # meaning "no shim" - in that case fall back to the default.