# HG changeset patch
# User Neil Dizon <neil.dizon@helsinki.fi>
# Date 1713610462 -10800
# Node ID dc6339756e35cdf75a13dff1e99fd3f7243ec326
# Parent  73f5e9bd63cfe5997d22cf95e5f51a2824f8fd4e
L update fix

diff -r 73f5e9bd63cf -r dc6339756e35 src/PET/AlgorithmDualScaling.jl
--- a/src/PET/AlgorithmDualScaling.jl	Sat Apr 20 12:42:22 2024 +0300
+++ b/src/PET/AlgorithmDualScaling.jl	Sat Apr 20 13:54:22 2024 +0300
@@ -153,11 +153,11 @@
             petgrad!(oldpetgradx, oldx, b, c, S)
             if norm₂(x-oldx)>1e-12
                 L = max(0.9*norm₂(petgradx - oldpetgradx)/norm₂(x-oldx),L)
+                println("Step length parameters: L=$(L)")
+                τ = τ₀/L
+                σ = σ₀*(1-τ₀)/(R_K²*τ)
             end  
-            println("Step length parameters: L=$(L)")
-            τ = τ₀/L
-            σ = σ₀*(1-τ₀)/(R_K²*τ)
-        end     
+        end       
        
         ################################
         # Give function value if needed
diff -r 73f5e9bd63cf -r dc6339756e35 src/PET/AlgorithmGreedy.jl
--- a/src/PET/AlgorithmGreedy.jl	Sat Apr 20 12:42:22 2024 +0300
+++ b/src/PET/AlgorithmGreedy.jl	Sat Apr 20 13:54:22 2024 +0300
@@ -153,11 +153,11 @@
             petgrad!(oldpetgradx, oldx, b, c, S)
             if norm₂(x-oldx)>1e-12
                 L = max(0.9*norm₂(petgradx - oldpetgradx)/norm₂(x-oldx),L)
+                println("Step length parameters: L=$(L)")
+                τ = τ₀/L
+                σ = σ₀*(1-τ₀)/(R_K²*τ)
             end  
-            println("Step length parameters: L=$(L)")
-            τ = τ₀/L
-            σ = σ₀*(1-τ₀)/(R_K²*τ)
-        end     
+        end       
        
         ################################
         # Give function value if needed
diff -r 73f5e9bd63cf -r dc6339756e35 src/PET/AlgorithmNoPrediction.jl
--- a/src/PET/AlgorithmNoPrediction.jl	Sat Apr 20 12:42:22 2024 +0300
+++ b/src/PET/AlgorithmNoPrediction.jl	Sat Apr 20 13:54:22 2024 +0300
@@ -153,11 +153,11 @@
             petgrad!(oldpetgradx, oldx, b, c, S)
             if norm₂(x-oldx)>1e-12
                 L = max(0.9*norm₂(petgradx - oldpetgradx)/norm₂(x-oldx),L)
+                println("Step length parameters: L=$(L)")
+                τ = τ₀/L
+                σ = σ₀*(1-τ₀)/(R_K²*τ)
             end  
-            println("Step length parameters: L=$(L)")
-            τ = τ₀/L
-            σ = σ₀*(1-τ₀)/(R_K²*τ)
-        end     
+        end      
        
         ################################
         # Give function value if needed
diff -r 73f5e9bd63cf -r dc6339756e35 src/PET/AlgorithmPrimalOnly.jl
--- a/src/PET/AlgorithmPrimalOnly.jl	Sat Apr 20 12:42:22 2024 +0300
+++ b/src/PET/AlgorithmPrimalOnly.jl	Sat Apr 20 13:54:22 2024 +0300
@@ -153,11 +153,11 @@
             petgrad!(oldpetgradx, oldx, b, c, S)
             if norm₂(x-oldx)>1e-12
                 L = max(0.9*norm₂(petgradx - oldpetgradx)/norm₂(x-oldx),L)
+                println("Step length parameters: L=$(L)")
+                τ = τ₀/L
+                σ = σ₀*(1-τ₀)/(R_K²*τ)
             end  
-            println("Step length parameters: L=$(L)")
-            τ = τ₀/L
-            σ = σ₀*(1-τ₀)/(R_K²*τ)
-        end     
+        end       
        
         ################################
         # Give function value if needed
diff -r 73f5e9bd63cf -r dc6339756e35 src/PET/AlgorithmProximal.jl
--- a/src/PET/AlgorithmProximal.jl	Sat Apr 20 12:42:22 2024 +0300
+++ b/src/PET/AlgorithmProximal.jl	Sat Apr 20 13:54:22 2024 +0300
@@ -182,10 +182,11 @@
             petgrad!(oldpetgradx, oldx, b, c, S)
             if norm₂(x-oldx)>1e-12
                 L = max(0.9*norm₂(petgradx - oldpetgradx)/norm₂(x-oldx),L)
+                println("Step length parameters: L=$(L)")
+                τ = τ₀/L
+                σ = σ₀*(1-τ₀)/(R_K²*τ)
             end  
-            #println("Step length parameters: L=$(L)")
-            τ, σ, σ̃, ρ̃ = step_lengths(params, γ, R_K², L)
-        end     
+        end       
        
         ################################
         # Give function value if needed
diff -r 73f5e9bd63cf -r dc6339756e35 src/PET/AlgorithmRotation.jl
--- a/src/PET/AlgorithmRotation.jl	Sat Apr 20 12:42:22 2024 +0300
+++ b/src/PET/AlgorithmRotation.jl	Sat Apr 20 13:54:22 2024 +0300
@@ -153,10 +153,10 @@
             petgrad!(oldpetgradx, oldx, b, c, S)
             if norm₂(x-oldx)>1e-12
                 L = max(0.9*norm₂(petgradx - oldpetgradx)/norm₂(x-oldx),L)
+                println("Step length parameters: L=$(L)")
+                τ = τ₀/L
+                σ = σ₀*(1-τ₀)/(R_K²*τ)
             end  
-            println("Step length parameters: L=$(L)")
-            τ = τ₀/L
-            σ = σ₀*(1-τ₀)/(R_K²*τ)
         end     
        
         ################################