assertion r>2M

svretina · May 29, 2024 · 6921b9b · 6921b9b
1 parent 3d6c667
commit 6921b9b
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/src/ODE.jl b/src/ODE.jl
@@ -78,6 +78,7 @@ using ..InitialData
                 end
             else
                 ri = InitialData.r(rstar[i], M) # transform rstar to r for potential
+                @assert ri > 2M
                 dtΠ[i] = drΨ + InitialData.Vpot(l, ri, M) * Φ[i]
                 dtΨ[i] = drΠ
             end