Test simulation waveform event types

Author: cncastillo

KomaMRICore/test/runtests.jl

@@ -247,6 +247,45 @@ end
     end
 end
 
+@testitem "Bloch waveform event types" tags=[:core, :nomotion] begin
+    include("initialize_backend.jl")
+
+    sys = Scanner()
+    obj = Phantom(x=[0.0], ρ=[1.0], T1=[1.0], T2=[1.0])
+    adc = ADC(4, 1e-3)
+    sim_params = Dict{String, Any}(
+        "gpu" => USE_GPU,
+        "return_type" => "mat",
+        "Δt" => 1e-4,
+        "Δt_rf" => 1e-4,
+    )
+
+    grad_events = (
+        Grad(1e-3, 1e-3),
+        Grad([0.0, 1e-3, 0.0], 1e-3),
+        Grad([0.0, 1e-3, 0.0], [0.5e-3, 0.5e-3]),
+    )
+    rf_events = (
+        RF(1e-7, 1e-3),
+        RF([0.0, 1e-7, 0.0], 1e-3),
+        RF([0.0, 1e-7, 0.0], [0.5e-3, 0.5e-3]),
+    )
+
+    for grad in grad_events
+        seq = Sequence([grad;;], [RF(0.0, 0.0);;], [adc])
+        sig = simulate(obj, seq, sys; sim_params, verbose=false)
+        @test size(sig) == (adc.N, 1, 1)
+        @test all(isfinite, sig)
+    end
+
+    for rf in rf_events
+        seq = Sequence([Grad(0.0, 0.0);;], [rf;;], [adc])
+        sig = simulate(obj, seq, sys; sim_params, verbose=false)
+        @test size(sig) == (adc.N, 1, 1)
+        @test all(isfinite, sig)
+    end
+end
+
 @testitem "Bloch_RF_accuracy" tags=[:important, :core, :nomotion] begin
     using OrdinaryDiffEqTsit5
     include("initialize_backend.jl")