From 77dcd1da7a505aadff3a7784149947dadfa7bef3 Mon Sep 17 00:00:00 2001
From: Nuwan Yapa <yapanuwan@gmail.com>
Date: Tue, 28 Jan 2025 12:34:43 -0500
Subject: [PATCH] Renaming

---
 dirac.jl => nucleons.jl                               | 4 ++--
 test/{Pb208_flds.jl => Pb208_meson_flds.jl}           | 0
 test/{Pb208_basic.jl => Pb208_nucleon_basic.jl}       | 2 +-
 test/{Pb208_density.jl => Pb208_nucleon_dens.jl}      | 2 +-
 test/{Pb208_spectrum.jl => Pb208_nucleon_spectrum.jl} | 2 +-
 test/{Pb208_dens.jl => Pb208_nucleon_wf.jl}           | 4 ++--
 6 files changed, 7 insertions(+), 7 deletions(-)
 rename dirac.jl => nucleons.jl (96%)
 rename test/{Pb208_flds.jl => Pb208_meson_flds.jl} (100%)
 rename test/{Pb208_basic.jl => Pb208_nucleon_basic.jl} (97%)
 rename test/{Pb208_density.jl => Pb208_nucleon_dens.jl} (98%)
 rename test/{Pb208_spectrum.jl => Pb208_nucleon_spectrum.jl} (96%)
 rename test/{Pb208_dens.jl => Pb208_nucleon_wf.jl} (87%)

diff --git a/dirac.jl b/nucleons.jl
similarity index 96%
rename from dirac.jl
rename to nucleons.jl
index a2c0437..f69766b 100644
--- a/dirac.jl
+++ b/nucleons.jl
@@ -27,7 +27,7 @@ end
     divs is the number of mesh divisions so solution would be returned as a 2×(1+divs) matrix,
     shooting method divides the interval into two partitions at r_max/2, ensuring convergence at both r=0 and r=r_max,
     the other parameters are the same from dirac!(...)."
-function solveDiracWf(κ, p, E, Φ0, W0, B0, A0, r_max, divs; shooting=true, normalize=true)
+function solveNucleonWf(κ, p, E, Φ0, W0, B0, A0, r_max, divs; shooting=true, normalize=true)
     saveat = r_max / divs
 
     if shooting
@@ -124,7 +124,7 @@ function calculateNucleonDensity(N, p, Φ0, W0, B0, A0, r_max, divs, E_min=0, E_
     ρr2 = zeros(2, divs + 1) # ρ×r² values
     
     for (κ, E, occ) in zip(κs, Es, occs)
-        wf = solveDiracWf(κ, p, E, Φ0, W0, B0, A0, r_max, divs; shooting=true, normalize=true)
+        wf = solveNucleonWf(κ, p, E, Φ0, W0, B0, A0, r_max, divs; shooting=true, normalize=true)
         wf2 = wf .* wf
         ρr2 += (occ / (4 * pi)) * wf2 # 2j+1 factor is accounted in the occupancy number
     end
diff --git a/test/Pb208_flds.jl b/test/Pb208_meson_flds.jl
similarity index 100%
rename from test/Pb208_flds.jl
rename to test/Pb208_meson_flds.jl
diff --git a/test/Pb208_basic.jl b/test/Pb208_nucleon_basic.jl
similarity index 97%
rename from test/Pb208_basic.jl
rename to test/Pb208_nucleon_basic.jl
index 6ab004a..51e0bab 100644
--- a/test/Pb208_basic.jl
+++ b/test/Pb208_nucleon_basic.jl
@@ -1,5 +1,5 @@
 using DelimitedFiles, Interpolations, Plots
-include("../dirac.jl")
+include("../nucleons.jl")
 
 # test data generated from Hartree.f
 # format: x S(x) V(x) R(x) A(x)
diff --git a/test/Pb208_density.jl b/test/Pb208_nucleon_dens.jl
similarity index 98%
rename from test/Pb208_density.jl
rename to test/Pb208_nucleon_dens.jl
index 5e3df8e..f869a79 100644
--- a/test/Pb208_density.jl
+++ b/test/Pb208_nucleon_dens.jl
@@ -1,5 +1,5 @@
 using DelimitedFiles, Interpolations, Plots
-include("../dirac.jl")
+include("../nucleons.jl")
 
 # test data generated from Hartree.f
 # format: x S(x) V(x) R(x) A(x)
diff --git a/test/Pb208_spectrum.jl b/test/Pb208_nucleon_spectrum.jl
similarity index 96%
rename from test/Pb208_spectrum.jl
rename to test/Pb208_nucleon_spectrum.jl
index 1d96389..282015b 100644
--- a/test/Pb208_spectrum.jl
+++ b/test/Pb208_nucleon_spectrum.jl
@@ -1,5 +1,5 @@
 using DelimitedFiles, Interpolations, Plots
-include("../dirac.jl")
+include("../nucleons.jl")
 
 # test data generated from Hartree.f
 # format: x S(x) V(x) R(x) A(x)
diff --git a/test/Pb208_dens.jl b/test/Pb208_nucleon_wf.jl
similarity index 87%
rename from test/Pb208_dens.jl
rename to test/Pb208_nucleon_wf.jl
index 46da8ab..0a13905 100644
--- a/test/Pb208_dens.jl
+++ b/test/Pb208_nucleon_wf.jl
@@ -1,5 +1,5 @@
 using DelimitedFiles, Interpolations, Plots
-include("../dirac.jl")
+include("../nucleons.jl")
 
 # test data generated from Hartree.f
 # format: x S(x) V(x) R(x) A(x)
@@ -25,7 +25,7 @@ groundE = findEs(κ, p, S_interp, V_interp, R_interp, A_interp, r_max, E_min, E_
 println("ground state E = $groundE")
 
 divs = 400
-wf = solveDiracWf(κ, p, groundE, S_interp, V_interp, R_interp, A_interp, r_max, divs)
+wf = solveNucleonWf(κ, p, groundE, S_interp, V_interp, R_interp, A_interp, r_max, divs)
 rs = range(0, r_max, length=divs+1)
 gs = wf[1, :]
 fs = wf[2, :]