fix remaining end

MagneticResonanceImaging · Dec 11, 2024 · 3d3d735 · 3d3d735
1 parent f5511ab
commit 3d3d735
Showing 1 changed file with 12 additions and 13 deletions.
diff --git a/src/Reconstruction/IterativeReconstruction.jl b/src/Reconstruction/IterativeReconstruction.jl
@@ -215,22 +215,21 @@ function reconstruction_multiCoil(acqData::AcquisitionData{T}
         E = encodingOps_parallel(acqData, reconSize, senseMapsUnCorr; slice=k, encParams...)
       end
 
-        W = WeightingOp(Complex{T}; weights=weights[j], rep=numChan)
-        kdata = arrayType((multiCoilData(acqData, j, k, rep=l))) .* repeat(weights[j], numChan)
-        if !isnothing(L_inv)
-          kdata = vec(reshape(kdata, :, numChan) * L_inv')
-        end
+      W = WeightingOp(Complex{T}; weights=weights[j], rep=numChan)
+      kdata = arrayType((multiCoilData(acqData, j, k, rep=l))) .* repeat(weights[j], numChan)
+      if !isnothing(L_inv)
+        kdata = vec(reshape(kdata, :, numChan) * L_inv')
+      end
 
-        EFull = ∘(W, E[j])
-        EFullᴴEFull = normalOperator(EFull; normalOpParams(arrayType)...)
-        solv = createLinearSolver(solver, EFull; AHA=EFullᴴEFull, reg=reg, params...)
-        I = solve!(solv, kdata)
+      EFull = ∘(W, E[j])
+      EFullᴴEFull = normalOperator(EFull; normalOpParams(arrayType)...)
+      solv = createLinearSolver(solver, EFull; AHA=EFullᴴEFull, reg=reg, params...)
+      I = solve!(solv, kdata)
 
-        if isCircular( trajectory(acqData, j) )
-          circularShutter!(reshape(I, reconSize), 1.0)
-        end
-        Ireco[:,k,j,l] = Array(I)
+      if isCircular( trajectory(acqData, j) )
+        circularShutter!(reshape(I, reconSize), 1.0)
       end
+      Ireco[:,k,j,l] = Array(I)
     end
   end
   Ireco_ = reshape(Ireco, volumeSize(reconSize, numSl)..., numContr, 1,numRep)