fix operator for multicoilmultiecho

MagneticResonanceImaging · Jun 26, 2023 · 36d32de · 36d32de
1 parent 86dbb0d
commit 36d32de
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/MRIOperators/src/EncodingOp.jl b/MRIOperators/src/EncodingOp.jl
@@ -90,7 +90,8 @@ function encodingOp_multiEcho_parallel(acqData::AcquisitionData{T,D}, shape::NTu
  # fourier operators
  ft = encodingOps_simple(acqData, shape; kargs...)
  S = SensitivityOp(reshape(smaps,:,numChan),numContrasts(acqData))
- return DiagOp(repeat(ft, inner=numChan)...) ∘ S
+ ops2 = [copy(ft[n]) for j=1:numChan,n=eachindex(ft)]
+ return DiagOp(ops2...) ∘ S
 end
 
 ###################################