test fixed

CliMA · Jun 2, 2023 · 4ea2d81 · 4ea2d81
1 parent e1754b2
commit 4ea2d81
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 4 deletions.
diff --git a/test/EnsembleKalmanProcess/runtests.jl b/test/EnsembleKalmanProcess/runtests.jl
@@ -156,7 +156,14 @@ end
         for (scheduler, N_iter) in zip(schedulers, N_iters)
             println("Scheduler: ", nameof(typeof(scheduler)))
             if !(nameof(typeof(process)) == Symbol(Unscented))
-                ekpobj = EKP.EnsembleKalmanProcess(initial_ensemble, y_obs, Γy, process, rng = copy(rng), scheduler = scheduler)
+                ekpobj = EKP.EnsembleKalmanProcess(
+                    initial_ensemble,
+                    y_obs,
+                    Γy,
+                    process,
+                    rng = copy(rng),
+                    scheduler = scheduler,
+                )
             else #no initial ensemble for UKI
                 ekpobj = EKP.EnsembleKalmanProcess(y_obs, Γy, process, rng = copy(rng), scheduler = scheduler)
             end

diff --git a/test/SparseInversion/runtests.jl b/test/SparseInversion/runtests.jl
@@ -191,7 +191,6 @@ include("../EnsembleKalmanProcess/inverse_problem.jl")
             if !isnothing(terminated)
                 break
             end
-            println(i, " ", ekiobj.Δt[end])
 
         end
         push!(init_means, vec(mean(get_u_prior(ekiobj), dims = 2)))
@@ -201,8 +200,6 @@ include("../EnsembleKalmanProcess/inverse_problem.jl")
     for i in 1:length(final_means)
         u_star = transform_constrained_to_unconstrained(prior, ϕ_star)
         inv_sqrt_Γy = sqrt(inv(Γy))
-        println(norm(inv_sqrt_Γy * (y_obs .- G(transform_unconstrained_to_constrained(prior, final_means[i])))) )
-        println(norm(inv_sqrt_Γy * (y_obs .- G(transform_unconstrained_to_constrained(prior, init_means[i])))) )
         #        @test norm(u_star - final_means[i]) < norm(u_star - init_means[i])
         @test norm(inv_sqrt_Γy * (y_obs .- G(transform_unconstrained_to_constrained(prior, final_means[i])))) <
               norm(inv_sqrt_Γy * (y_obs .- G(transform_unconstrained_to_constrained(prior, init_means[i]))))