From 6aa90d9ee793c0ba981eb2ed1522a2a3470d97d4 Mon Sep 17 00:00:00 2001
From: Michiel Cottaar <michiel.cottaar@ndcn.ox.ac.uk>
Date: Mon, 26 Feb 2024 13:23:59 +0000
Subject: [PATCH] Ensure times are positive
---
src/components/pulses/sinc_pulses.jl | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/components/pulses/sinc_pulses.jl b/src/components/pulses/sinc_pulses.jl
index db5e47f..53ffb9d 100644
--- a/src/components/pulses/sinc_pulses.jl
+++ b/src/components/pulses/sinc_pulses.jl
@@ -88,7 +88,7 @@ effective_time(pulse::SincPulse) = N_left(pulse) * lobe_duration(pulse)
function make_generic(block::SincPulse)
normed_times = -N_left(block):0.1:N_right(block) + 1e-5
- times = (normed_times .+ N_left(block)) .* lobe_duration(block)
+ times = max.(0., (normed_times .+ N_left(block))) .* lobe_duration(block)
amplitudes = amplitude(block) .* (normalised_function.(normed_times; apodise=block.apodise))
phases = [frequency(block) .* lobe_duration(block)] .* normed_times .* 360
return GenericPulse(times, amplitudes, phases, effective_time(block))
--
GitLab