Normalization fix
This commit is contained in:
parent
990d5a6aeb
commit
2ee76f2c1b
|
|
@ -53,7 +53,7 @@ function solveNucleonWf(κ, p, E, Φ0, W0, B0, A0, r_max, divs; shooting=true, n
|
||||||
|
|
||||||
if normalize
|
if normalize
|
||||||
norm = sum(wf .* wf) * r_max / divs # integration by Reimann sum
|
norm = sum(wf .* wf) * r_max / divs # integration by Reimann sum
|
||||||
wf = wf ./ sqrt(norm)
|
wf = wf ./ sqrt(norm * 2) # WHY FACTOR OF 2?
|
||||||
end
|
end
|
||||||
|
|
||||||
return wf
|
return wf
|
||||||
|
|
|
||||||
Loading…
Reference in New Issue