Merge pull request #511 from pvillacorta/fix-spin-reset

Adjust reset flag interpolation to trigger at jump
JuliaHealth · Nov 22, 2024 · c6d2183 · c6d2183
2 parents ad6f5d2 + 835ffc8
commit c6d2183
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/KomaMRICore/src/simulation/Flow.jl b/KomaMRICore/src/simulation/Flow.jl
@@ -82,11 +82,11 @@ function replace_view(replace_by, idx)
 end
 
 function get_mask(spin_reset, t::Real)
-   itp  = KomaMRIBase.interpolate(spin_reset, KomaMRIBase.Gridded(KomaMRIBase.Constant{KomaMRIBase.Previous}()), Val(size(spin_reset, 1)), t)
+   itp  = KomaMRIBase.interpolate(spin_reset, KomaMRIBase.Gridded(KomaMRIBase.Constant{KomaMRIBase.Next}()), Val(size(spin_reset, 1)), t)
    return KomaMRIBase.resample(itp, t)
 end
 function get_mask(spin_reset, t::AbstractArray)
-   itp  = KomaMRIBase.interpolate(spin_reset, KomaMRIBase.Gridded(KomaMRIBase.Constant{KomaMRIBase.Previous}()), Val(size(spin_reset, 1)), t)
+   itp  = KomaMRIBase.interpolate(spin_reset, KomaMRIBase.Gridded(KomaMRIBase.Constant{KomaMRIBase.Next}()), Val(size(spin_reset, 1)), t)
    mask = KomaMRIBase.resample(itp, t)
    mask .= (cumsum(mask; dims=2) .== 1)
    return mask