# HG changeset patch
# User Neil Dizon <neil.dizon@helsinki.fi>
# Date 1713692764 -10800
# Node ID c286925c0f35fcb536ed060dcd1418cb30b26193
# Parent  dc6339756e35cdf75a13dff1e99fd3f7243ec326
Change rand to randn in PET

diff -r dc6339756e35 -r c286925c0f35 src/PET/ImGenerate.jl
--- a/src/PET/ImGenerate.jl	Sat Apr 20 13:54:22 2024 +0300
+++ b/src/PET/ImGenerate.jl	Sun Apr 21 12:46:04 2024 +0300
@@ -108,7 +108,7 @@
 
 
 function rotatebytheta(params)
-    r = params.rotation_factor*(2*rand(rng)-1)
+    r = params.rotation_factor*randn(rng)
     return r
 end
 
@@ -272,8 +272,8 @@
 end
 
 
-function imgen_shepplogan_radon(origsize,sz)
-    im = convert(Array{Float64},TestImages.shepp_logan(origsize, highContrast=true))
+function imgen_shepplogan_radon(sz)
+    im = convert(Array{Float64},TestImages.shepp_logan(sz[1], highContrast=true))
     dynrange = maximum(im)
     return ImGen(curry(generate_radon, im, sz), sz, 1, dynrange, "shepplogan$(sz[1])x$(sz[2])")
 end
diff -r dc6339756e35 -r c286925c0f35 src/PET/PET.jl
--- a/src/PET/PET.jl	Sat Apr 20 13:54:22 2024 +0300
+++ b/src/PET/PET.jl	Sun Apr 21 12:46:04 2024 +0300
@@ -101,8 +101,8 @@
     noise_level = 0.5,
     shake_noise_level = 0.1,
     shake = 1.0, 
-    rotation_factor = 0.15,                  
-    rotation_noise_level = 0.0025,          
+    rotation_factor = 0.075,                  
+    rotation_noise_level = 0.0075,         
     α = 1.0,
     ρ̃₀ = 1.0,
     σ̃₀ = 1.0,
@@ -120,7 +120,7 @@
     L_experiment = false,
 )
 
-const petscan = imgen_shepplogan_radon(p_known₀_pet.origsize, p_known₀_pet.sz)
+const petscan = imgen_shepplogan_radon(p_known₀_pet.sz)
 
 const pet_experiments_pdps_known = (
     Experiment(AlgorithmDualScaling, DisplacementConstant, petscan,