start test design

.github/workflows/ci.yml

@@ -0,0 +1,29 @@
+name: CI
+on:
+  - push
+
+jobs:
+  CI:
+    runs-on: ubuntu-latest
+    steps:
+      - uses: actions/checkout@v4
+      - uses: julia-actions/setup-julia@v1
+        with:
+          version: 1.9
+          arch: x64
+      - name: instantiate environment
+        run: |
+            pwd
+            julia --project=df -e '
+              using Pkg
+              Pkg.add("BenchmarkTools")
+              Pkg.add("Test")
+              Pkg.add(url="https://github.com/JuliaStellarDynamics/OrbitalElements.jl.git")
+              Pkg.add(url="https://github.com/JuliaStellarDynamics/DistributionFunctions.jl.git")
+            ls
+            julia --project=df -e '
+              using Pkg
+              Pkg.develop(PackageSpec(path="/home/runner/work/DistributionFunctions.jl/DistributionFunctions.jl"))'
+            julia --project=df test/runtests.jl
+
+

src/Analytic/Isochrone/isotropic.jl

@@ -22,7 +22,7 @@ function F(EL::Tuple{Float64,Float64}, df::IsotropicIsochrone)
     # rescale dimensionless energy to goes between 0 and 1/2
     mE = E/scaleEnergy
 
-    return (sqrt(mE)*M*(27.0+2.0*mE*(-1.0+4.0*mE)*(33.0+4.0*mE*(-7.0+2.0*mE))+
+    return (sqrt(mE)*df.potential.M*(27.0+2.0*mE*(-1.0+4.0*mE)*(33.0+4.0*mE*(-7.0+2.0*mE))+
     (3.0*(-9.0+4.0*mE*(7.0+4.0*mE))*asin(sqrt(mE)))/sqrt(-((-1.0+mE)*mE))))/
     (128.0*sqrt(2.0)*(-1.0+mE)^(4)*scaleDF*(pi)^(3))
 

test/runtests.jl

@@ -0,0 +1,6 @@
+using DistributionFunctions
+using OrbitalElements
+
+using Test
+
+include("test_spheres.jl")

test/test_spheres.jl

@@ -9,8 +9,16 @@ n1,n2 = -1,2
 
 @testset "distributionfunctions" begin
     @testset "isochrone" begin
-        IsoDF = IsotropicIsochrone()
-        EL = EL_from_ae(1.0,0.5,IsoDF.potential)
-        resonance = Resonance(n1,n2,IsoDF.potential)
+        @testset "isotropic" begin
+            IsoDF = IsotropicIsochrone()
+            EL = EL_from_ae(1.0,0.5,IsoDF.potential)
+            resonance = Resonance(n1,n2,IsoDF.potential)
+            @test F(EL, IsoDF) ≈ 0.023905 atol=1e-6
+        end
+        @testset "osipkovmerritt" begin
+
+        end
+    end
+    @testset "plummer" begin
     end
 end