Mercurial > repos > PredictPDPS / changeset

--- a/src/PET/PET.jl	Mon Apr 22 09:26:19 2024 +0300
+++ b/src/PET/PET.jl	Mon Apr 22 13:28:58 2024 +0300
@@ -129,7 +129,9 @@
     sino_sparsity = 0.5,            # Percentage of zero parts in the mask
     L = 300.0,
     L_experiment = false,
-    stable_interval = union(Set(1:1000),Set(2500:3000)),
+    stable_interval = Set(0),
+    # stable_interval = union(Set(1:1000),Set(2500:3000)),
+    # stable_interval = union(Set(1000:2000),Set(3000:4000)),
 )
 
 const shepplogan = imgen_shepplogan_radon(p_known₀_pet.sz)

--- a/src/PredictPDPS.jl	Mon Apr 22 09:26:19 2024 +0300
+++ b/src/PredictPDPS.jl	Mon Apr 22 13:28:58 2024 +0300
@@ -135,13 +135,26 @@
     handle_interrupt = true,
     init = :zero,
     plot_movement = false,
-    # stable_interval = Set(),
-    stable_interval = union(Set(1:1500),Set(6000:7500)),
+    stable_interval = Set(0),
 )
 
 const square = imgen_square((200, 300))
 const lighthouse = imgen_shake("lighthouse", (200, 300))
 
+const p_known₀_denoising = (
+    noise_level = 0.5,
+    shake_noise_level = 0.025,
+    shake = 3.0,
+    α = 1.0,
+    ρ̃₀ = 1.0,
+    σ̃₀ = 1.0,
+    δ = 0.9,
+    σ₀ = 0.1,
+    τ₀ = 0.01,
+    # stable_interval = union(Set(1:1500),Set(6000:7500)),
+    # stable_interval = union(Set(2000:4000),Set(8000:10000)),
+)
+
 const p_known₀ = (
     noise_level = 0.5,
     shake_noise_level = 0.05,
@@ -171,17 +184,6 @@
     τ₀ = 0.01,
 )
 
-const p_known₀_denoising = (
-    noise_level = 0.5,
-    shake_noise_level = 0.025,
-    shake = 3.0,
-    α = 1.0,
-    ρ̃₀ = 1.0,
-    σ̃₀ = 1.0,
-    δ = 0.9,
-    σ₀ = 0.1,
-    τ₀ = 0.01,
-)
 
 const experiments_pdps_known = (
     Experiment(Algorithm, DisplacementConstant, lighthouse,