LNRF test

Also removed Dilaton-Axion from tests (see issue #52).

test/unit/gradusbase.geometry.jl

@@ -1,4 +1,3 @@
-
 @testset "tetradframe" begin
     all_metrics = (
         # can't do first order yet since no four velocity
@@ -8,7 +7,7 @@
         BoyerLindquistAD(1.0, -0.998),
         JohannsenAD(M = 1.0, a = 0.998, α13 = 1.0),
         JohannsenAD(M = 1.0, a = 0.998, α22 = 1.0),
-        DilatonAxionAD(M = 1.0, a = 0.998, β = 0.2, b = 1.0),
+        # DilatonAxionAD(M = 1.0, a = 0.998, β = 0.2, b = 1.0),
     )
     radii = 5.0:0.8:10.0
     angles = 0.1:0.5:2π
@@ -30,4 +29,19 @@
         # ensure it gives minkowski
         @test isapprox(res, minkowski, atol = 1e-13)
     end
+
+
+    @testset "lnrf" begin
+        for m in all_metrics, r in radii, θ in angles
+            u = @SVector([0.0, r, θ, 0.0])
+
+            # function that we are testing
+            M = GradusBase.lnrframe(m, u)
+
+            m_mat = Gradus.metric(m, u)
+            @tullio res[a, b] := m_mat[i, j] * M[a][i] * M[b][j]
+            # ensure it gives minkowski
+            @test isapprox(res, minkowski, atol = 1e-13)
+        end
+    end
 end