Relaxed tests on tensorproduct

Author: theogf

test/kernels/tensorproduct.jl

@@ -32,27 +32,27 @@
         tmp = Matrix{Float64}(undef, 10, 10)
 
         for kernel in (kernel1, kernel2)
-            @test kernelmatrix(kernel, X) == trueX
-            @test kernelmatrix(kernel, X'; obsdim = 1) == trueX
+            @test kernelmatrix(kernel, X) ≈ trueX
+            @test kernelmatrix(kernel, X'; obsdim = 1) ≈ trueX
 
-            @test kernelmatrix(kernel, X, Y) == trueXY
-            @test kernelmatrix(kernel, X', Y'; obsdim = 1) == trueXY
+            @test kernelmatrix(kernel, X, Y) ≈ trueXY
+            @test kernelmatrix(kernel, X', Y'; obsdim = 1) ≈ trueXY
 
             fill!(tmp, 0)
             kernelmatrix!(tmp, kernel, X)
-            @test tmp == trueX
+            @test tmp ≈ trueX
 
             fill!(tmp, 0)
             kernelmatrix!(tmp, kernel, X'; obsdim = 1)
-            @test tmp == trueX
+            @test tmp ≈ trueX
 
             fill!(tmp, 0)
             kernelmatrix!(tmp, kernel, X, Y)
-            @test tmp == trueXY
+            @test tmp ≈ trueXY
 
             fill!(tmp, 0)
             kernelmatrix!(tmp, kernel, X', Y'; obsdim = 1)
-            @test tmp == trueXY
+            @test tmp ≈ trueXY
         end
     end
 
@@ -96,10 +96,10 @@
             tmp = Matrix{Float64}(undef, 10, 10)
 
             @test kernelmatrix(kernel, X) == trueX
-            @test kernelmatrix(kernel, X'; obsdim = 1) == trueX
+            @test kernelmatrix(kernel, X'; obsdim = 1) ≈ trueX
 
-            @test kernelmatrix(kernel, X, Y) == trueXY
-            @test kernelmatrix(kernel, X', Y'; obsdim = 1) == trueXY
+            @test kernelmatrix(kernel, X, Y) ≈ trueXY
+            @test kernelmatrix(kernel, X', Y'; obsdim = 1) ≈ trueXY
 
             fill!(tmp, 0)
             kernelmatrix!(tmp, kernel, X)