using isapprox for floating point comparison (#33596)

one way to fix JuliaLang/julia#33580
possibly due to mkl/julia handling the matrix product of the lazy `'` of AA and `copy(AA')` differently

Author: daviehh

test/runtests.jl

@@ -213,7 +213,7 @@ d[5,1:2:4,8] .= 19
 AA = rand(4,2)
 A = @inferred(convert(SharedArray, AA))
 B = @inferred(convert(SharedArray, copy(AA')))
-@test B*A == AA'*AA
+@test B*A ≈ AA'*AA
 
 d=SharedArray{Int64,2}((10,10); init = D->fill!(D.loc_subarr_1d, myid()), pids=[id_me, id_other])
 d2 = map(x->1, d)