further tests

Author: mboberg

src/MPISphericalHarmonics.jl

@@ -47,6 +47,10 @@ function MagneticFieldCoefficients(L::Int)
   end
   return MagneticFieldCoefficients(reshape([SphericalHarmonicCoefficients(L),SphericalHarmonicCoefficients(L),SphericalHarmonicCoefficients(L)],3,1))
 end
+
+# constructor using t-design
+MagneticFieldCoefficients(coeffs::Array{SphericalHarmonicCoefficients,2}, tDesign::SphericalTDesign, ffp=nothing) = MagneticFieldCoefficients(coeffs,ustrip(Unitful.m(tDesign.radius)), ustrip.(Unitful.m.(tDesign.center)), ffp)
+  
 
 # read coefficients from an HDF5-file
 function MagneticFieldCoefficients(path::String)
@@ -69,8 +73,8 @@ function MagneticFieldCoefficients(path::String)
   else
     # convert file of SphericalHarmonicCoefficients into MagneticFieldCoefficients
     # -> set all additional informations to 0 or nothing
-    # use radius = 0.042 as default value
-    return MagneticFieldCoefficients(shcoeffs, 0.042)
+    # use radius = 0.01 as default value
+    return MagneticFieldCoefficients(shcoeffs, 0.01)
   end
 end
 
@@ -163,7 +167,7 @@ function magneticField(coords::AbstractArray{T, 2}, field::Union{AbstractArray{T
 end
 
 #TODO: This should be merged with and moved to MPIFiles
-function loadTDesign(filename::String)
+function loadTDesignCoefficients(filename::String)
   field, radius, N, t, center, correction  = h5open(filename, "r") do file
     field = read(file,"/fields") 		# measured field (size: 3 x #points x #patches)
     radius = read(file,"/positions/tDesign/radius")	# radius of the measured ball
@@ -203,7 +207,7 @@ function SphericalHarmonicsDefinedField(filename::String)
     func = fastfunc.(coeffs_MF.coeffs)
   else
     # load measured field
-    coeffs_MF, expansion, func = loadTDesign(filename)
+    coeffs_MF, expansion, func = loadTDesignCoefficients(filename)
   end
 
   return SphericalHarmonicsDefinedField(func=func)

test/idealGradientField.h5

@@ -14,7 +14,7 @@ using Aqua
 
     ## Calculate some field values
     # load a spherical t-design
-    tDes = loadTDesign(6,26,42u"mm")
+    tDes = loadTDesign(6,26,42.0u"mm")
     # ideal gradient field with FFP
     idealField = IdealFFP([-1,-1,2])
     # get field values
@@ -34,6 +34,88 @@ using Aqua
       @test isapprox(coeffs[j,1].c, ideal[j], atol = 1e-10)
     end
 
+    @testset "MagneticFieldCoefficients" begin
+
+      ## Test Constructor
+      # Errors
+      @test_throws DomainError MagneticFieldCoefficients(-2)
+      @test_throws DimensionMismatch MagneticFieldCoefficients(coeffs[1:2,:])
+
+      # standard constructors
+      coeffsMF = MagneticFieldCoefficients(coeffs)
+      @test coeffsMF.coeffs == coeffs
+      @test coeffsMF.radius == 0.0
+      @test coeffsMF.center == zeros(Float64,3)
+      @test coeffsMF.ffp === nothing
+
+      # constructor with t-design
+      coeffsMF = MagneticFieldCoefficients(coeffs, tDes, zeros(Float64,3,1))
+      @test coeffsMF.radius == 0.042
+      @test coeffsMF.center == zeros(Float64,3)
+      @test coeffsMF.ffp == zeros(Float64,3,1) 
+    
+      # constructor with wrong FFP sizes
+      @test_throws DimensionMismatch MagneticFieldCoefficients(coeffs, tDes, zeros(2,1))
+      @test_throws DimensionMismatch MagneticFieldCoefficients(coeffs, tDes, zeros(3,2))
+    end
+
+    @testset "Load data from file" begin
+      ɛ = eps(Float64)
+
+      ## measurement data (without coefficients)
+      filename = "idealGradientField.h5"
+      func = SphericalHarmonicsDefinedField(filename)
+      @test isapprox(func[0.01,0.01,0.01], [-0.01,-0.01,0.02], atol=ε)
+
+      # get coefficients
+      coeffsMF, = MPISphericalHarmonics.loadTDesignCoefficients(filename)
+      @test isapprox(coeffsMF.radius, 0.042, atol=ε) # radius
+      @test isapprox(coeffsMF.coeffs[1][1,1], -1.0, atol=1e-10) # gradient (x)
+      @test isapprox(coeffsMF.coeffs[2][1,-1], -1.0, atol=1e-10) # gradient (y)
+      @test isapprox(coeffsMF.coeffs[3][1,0], 2.0, atol=1e-10) # gradient (z)
+
+      ## load coefficients from file
+      # write coefficients in file
+      filename2 = "idealGradientFieldwithCoeffs.h5"
+      filename3 = "idealGradientFieldwithCoeffs2.h5"
+      filename4 = "idealGradientFieldwithCoeffs3.h5"
+      cp(filename, filename2)
+      write(filename2, coeffsMF.coeffs)
+      # add radius and center
+      cp(filename2, filename3)
+      h5open(filename3, "cw") do file
+        write(file, "/radius", coeffsMF.radius)
+        write(file, "/center", coeffsMF.center)
+      end
+      # add FFP
+      cp(filename3, filename4)
+      h5open(filename4, "cw") do file
+        write(file, "/ffp", zeros(3,1))
+      end
+
+      # only coefficients (no further informations given)
+      coeffsTest = MagneticFieldCoefficients(filename2)
+      @test isapprox(coeffsTest.radius, 0.01, atol=ε) # radius
+ 
+      # with given radius & center
+      coeffsTest = MagneticFieldCoefficients(filename3)
+      @test isapprox(coeffsTest.radius, 0.042, atol=ε) # radius
+      @test isapprox(coeffsTest.center, zeros(3), atol=ε) # center
+      @test coeffsTest.ffp === nothing # FFP
+      
+      # with given FFP
+      coeffsTest = MagneticFieldCoefficients(filename4)
+      @test coeffsTest.ffp == zeros(3,1) # FFP
+
+      # remove test files
+      rm(filename2)
+      rm(filename3)
+      rm(filename4)
+    end
+
+    @testset "Multiple patches" begin
+      #TODO
+    end
   end
 
 end