From fcc1d97e7b03f5a51f70b93c72eabe00d77d7cd8 Mon Sep 17 00:00:00 2001
From: Nuwan Yapa <yapanuwan@gmail.com>
Date: Wed, 3 Jul 2024 15:16:14 -0400
Subject: [PATCH] Minor change

---
 calculations/3body_Berggren_B2R_EC.jl | 16 +++++++++-------
 calculations/3body_Berggren_R2R_EC.jl |  8 ++++----
 2 files changed, 13 insertions(+), 11 deletions(-)

diff --git a/calculations/3body_Berggren_B2R_EC.jl b/calculations/3body_Berggren_B2R_EC.jl
index aa38b17..4b55290 100644
--- a/calculations/3body_Berggren_B2R_EC.jl
+++ b/calculations/3body_Berggren_B2R_EC.jl
@@ -6,9 +6,10 @@ include("../berggren.jl")
 
 println("No of threads = ", Threads.nthreads())
 
-atol = 10^-5
-maxevals = 10^5
-R_cutoff = 16
+training_c = [2.0, 1.9, 1.8]
+extrapolating_c = 0.0 : 0.2 : 1.2
+
+training_ref = -0.72763 # complete list not needed because identification is simple
 
 exact_ref = reverse([4.076662025307587-0.012709842443350328im,
     3.613318119833891-0.007335804709990623im,
@@ -23,6 +24,10 @@ m = 1.0
 Va_of_r(r) = 2 * exp(-(r-3)^2 / (1.5)^2)
 Vb_of_r(r) = -exp(-(r/3)^2)
 
+atol = 10^-5
+maxevals = 10^5
+R_cutoff = 16
+
 # due to Jacobi coordinates
 μ1 = m * 1/2
 μ2 = m * 2/3
@@ -82,10 +87,7 @@ display(test_evals)
 Es = n1s = l1s = n2s = l2s = mask1 = mask2 = T1 = T2 = V1_cache = V_relative_cache = V1 = V_relative = U = V2 = nothing
 GC.gc()
 
-current_E = -0.72763
-
-training_c = [2.0, 1.9, 1.8]
-extrapolating_c = 0.0 : 0.2 : 1.2
+current_E = training_ref
 
 exact = ComplexF64[]
 training = ComplexF64[]
diff --git a/calculations/3body_Berggren_R2R_EC.jl b/calculations/3body_Berggren_R2R_EC.jl
index 77894e1..6fa5cc8 100644
--- a/calculations/3body_Berggren_R2R_EC.jl
+++ b/calculations/3body_Berggren_R2R_EC.jl
@@ -6,10 +6,6 @@ include("../berggren.jl")
 
 println("No of threads = ", Threads.nthreads())
 
-atol = 10^-5
-maxevals = 10^5
-R_cutoff = 16
-
 training_c = [1.1, 0.9, 0.7, 0.5]
 extrapolating_c = 0.0 : 0.2 : 1.2
 
@@ -31,6 +27,10 @@ m = 1.0
 Va_of_r(r) = 2 * exp(-(r-3)^2 / (1.5)^2)
 Vb_of_r(r) = -exp(-(r/3)^2)
 
+atol = 10^-5
+maxevals = 10^5
+R_cutoff = 16
+
 # due to Jacobi coordinates
 μ1 = m * 1/2
 μ2 = m * 2/3