yapanuwan
/
NuclearRMF
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Elaborate docstring

This commit is contained in:
Nuwan Yapa 2024-06-18 12:32:04 -04:00
parent 48f757f907
commit a6b9234f10
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
dirac.jl
View File

@ -22,8 +22,8 @@ end
M is the mass in MeV/c2, M is the mass in MeV/c2,
E in the energy in MeV, E in the energy in MeV,
S(r) & V(r) are functions corresponding to scalar and vector potentials in MeV, S(r) & V(r) are functions corresponding to scalar and vector potentials in MeV,
r_max is the outer boundary, r_max is the outer boundary in fm,
r_min (=r_max/1000) is inside boundary which cannot be 0 due to the centrifugal term." r_min (=r_max/1000) is inside boundary in fm which cannot be 0 due to the centrifugal term."
function boundaryValue(κ, M, E, S, V, r_max, r_min=r_max/1000) function boundaryValue(κ, M, E, S, V, r_max, r_min=r_max/1000)
prob = ODEProblem(dirac!, [0, 1], (r_min, r_max)) prob = ODEProblem(dirac!, [0, 1], (r_min, r_max))
sol = solve(prob, RK4(), p=(κ, M, E, S, V)) sol = solve(prob, RK4(), p=(κ, M, E, S, V))
@ -34,8 +34,8 @@ end
κ is the generalized angular momentum, κ is the generalized angular momentum,
M is the mass in MeV/c2, M is the mass in MeV/c2,
S(r) & V(r) are functions corresponding to scalar and vector potentials in MeV, S(r) & V(r) are functions corresponding to scalar and vector potentials in MeV,
r_max is the outer boundary, r_max is the outer boundary in fm,
r_min (=r_max/1000) is inside boundary which cannot be 0 due to the centrifugal term." r_min (=r_max/1000) is inside boundary in fm which cannot be 0 due to the centrifugal term."
function findEs(κ, M, S, V, r_max, r_min=r_max/1000, E_min=0, E_max=M) function findEs(κ, M, S, V, r_max, r_min=r_max/1000, E_min=0, E_max=M)
f(E) = boundaryValue(κ, M, E, S, V, r_max, r_min) f(E) = boundaryValue(κ, M, E, S, V, r_max, r_min)
return find_zeros(f, (E_min, E_max)) return find_zeros(f, (E_min, E_max))