Skip to content
Snippets Groups Projects
Verified Commit 955ccb5f authored by Michiel Cottaar's avatar Michiel Cottaar
Browse files

Convert InstantGradient1D to InstantGradient3D

parent 7227f0b8
No related branches found
No related tags found
No related merge requests found
Pipeline #23466 passed
......@@ -7,7 +7,7 @@ include("readouts/readouts.jl")
import .AbstractTypes: BaseComponent, GradientWaveform, EventComponent, RFPulseComponent, ReadoutComponent, split_timestep
import .GradientWaveforms: ConstantGradient, ChangingGradient, NoGradient, split_gradient
import .InstantGradients: InstantGradient, InstantGradient3D
import .InstantGradients: InstantGradient, InstantGradient3D, InstantGradient1D
import .Pulses: GenericPulse, InstantPulse, SincPulse, ConstantPulse
import .Readouts: ADC, SingleReadout
......
module Linearise
import StaticArrays: SVector
import ...Components: GradientWaveform, split_timestep, InstantPulse, InstantGradient3D
import ...Variables: amplitude, phase, gradient_strength3, duration
import ...Components: GradientWaveform, split_timestep, InstantPulse, InstantGradient3D, InstantGradient1D
import ...Variables: amplitude, phase, gradient_strength3, duration, qval
import ..Abstract: edge_times, start_time, end_time, ContainerBlock, iter_instant_gradients, iter_instant_pulses
import ..BaseSequences: BaseSequence, Sequence
import ..BuildingBlocks: BaseBuildingBlock
......@@ -181,6 +181,9 @@ function LinearSequence(container::BaseSequence, times::AbstractVector{<:Number}
if !(tstart <= real_time < tfinal)
continue
end
if pulse isa InstantGradient1D
pulse = InstantGradient3D(qval(pulse) .* pulse.orientation, pulse.group)
end
index = findmin(t -> abs(t - real_time), times)[2]
to_store[index] = pulse
end
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment