Prevent errors with optional gradient definition

MagneticParticleImaging · Jun 19, 2024 · f168899 · f168899
1 parent a5d1138
commit f168899
Showing 1 changed file with 6 additions and 4 deletions.
diff --git a/src/Utils.jl b/src/Utils.jl
@@ -50,10 +50,12 @@ end
 generateHeaderDict(bSF::MPIFile, b::Vector) = generateHeaderDict(bSF,b[1])
 
 function consistenceCheck(bSF::MPIFile, bMeas::MPIFile)
-  gSF = acqGradient(bSF)[:,:,1,1]
-  gMeas = acqGradient(bMeas)[:,:,1,1]
-  if gSF != gMeas
-    @warn("The gradient strength of the system matrix ($gSF T/m) does not fit to the measurements ($gMeas T/m)!")
+  if !(isnothing(acqGradient(bSF)) || isnothing(acqGradient(bMeas))) # /acquisition/gradient is optional
+    gSF = acqGradient(bSF)[:,:,1,1]
+    gMeas = acqGradient(bMeas)[:,:,1,1]
+    if gSF != gMeas
+      @warn("The gradient strength of the system matrix ($gSF T/m) does not fit to the measurements ($gMeas T/m)!")
+    end
   end
 
   dfSF = dfStrength(bSF)