Merge branch 'main' into compathelper/new_version/2022-09-14-02-42-48…

…-613-00127700948

.gitignore

@@ -1 +1,2 @@
 /Manifest.toml
+.vscode

Project.toml

@@ -1,7 +1,7 @@
 name = "Jadex"
 uuid = "302b6ea9-03d1-47d0-b3e6-7bc5ec1c11d5"
 authors = ["Brian Svoboda <[email protected]> and contributors"]
-version = "0.2.1"
+version = "0.3.0"
 
 [deps]
 CSV = "336ed68f-0bac-5ca0-87d4-7b16caf5d00b"
@@ -15,7 +15,9 @@ RecursiveFactorization = "f2c3362d-daeb-58d1-803e-2bc74f2840b4"
 [compat]
 CSV = "0.10"
 DataFrames = "1"
+LoopVectorization = "0.12"
 Interpolations = "0.14"
+RecursiveFactorization = "0.2"
 julia = "1.7"
 
 [extras]

src/gridded.jl

@@ -3,7 +3,6 @@ module GridRunner
 export runseq, rungrid, get_interp, interp_errors
 
 using Base.Threads
-using Distributed
 using DataFrames
 using Interpolations
 

test/runtests.jl

@@ -235,9 +235,13 @@ end
             # LoopVectorization will warn on `BigFloat` that certain
             # optimizations are unavailable.
             @test_logs (:warn,) match_mode=:any solve!(sol_big, rdf)
-            @test maxabsdiff(sol_f64.τl, sol_big.τl) ≈ 0.599_831_102
-            @test maxabsdiff(sol_f64.tex, sol_big.tex) ≈ 0.0222_790_875
-            @test maxabsdiff(sol_f64.xpop, sol_big.xpop) ≈ 0.273_337_065
+            # The last values are numerically unstable and produce
+            # hardware dependent outputs.
+            slice = 1:10
+            ϵ = √eps()
+            @test maxabsdiff(sol_f64.τl[slice], sol_big.τl[slice]) < ϵ
+            @test maxabsdiff(sol_f64.tex[slice], sol_big.tex[slice]) < ϵ
+            @test maxabsdiff(sol_f64.xpop[slice], sol_big.xpop[slice]) < ϵ
         end
 
         @testset "HCO+ grid" begin