Update test_NgIteration_HS.jl

IlianPihlajamaa · Oct 3, 2024 · b8744da · b8744da
1 parent d964e95
commit b8744da
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/test/test_NgIteration_HS.jl b/test/test_NgIteration_HS.jl
@@ -57,7 +57,7 @@ kBT = 1.0
 pot = HardSpheres(1.0)
 system = SimpleLiquid(dims, ρ, kBT, pot)
 closure = PercusYevick()
-method = NgIteration(tolerance=10^-10, N_stages=5, M=M, verbose=false, max_iterations=10^3)
+method = NgIteration(tolerance=10^-10, N_stages=5, M=M, verbose=false, max_iterations=10^3, dr=sqrt(π/(M+1))/(2π))
 sol = solve(system, closure, method)
 sol2 = solve(system, closure, Exact(M=M))
 
@@ -78,7 +78,7 @@ kBT = 1.0
 pot = HardSpheres(1.0)
 system = SimpleLiquid(dims, ρ, kBT, pot)
 closure = PercusYevick()
-method = NgIteration(tolerance=10^-10, N_stages=5, M=M, verbose=false, max_iterations=10^3)
+method = NgIteration(tolerance=10^-10, N_stages=5, M=M, verbose=false, max_iterations=10^3, dr=sqrt(π/(M+1))/(2π))
 sol = solve(system, closure, method)
 sol2 = solve(system, closure, Exact(M=M))
 
@@ -92,4 +92,4 @@ atol = 0.1
 
 end 
 
-test2()
+test2()