From 1311db53f7198758d30e72d257a94c68dfec2203 Mon Sep 17 00:00:00 2001
From: Nuwan Yapa <yapanuwan@gmail.com>
Date: Wed, 16 Apr 2025 16:48:03 -0400
Subject: [PATCH] Missed boundary conditions during merge

---
 nucleons.jl | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/nucleons.jl b/nucleons.jl
index b3fde98..095c212 100644
--- a/nucleons.jl
+++ b/nucleons.jl
@@ -68,12 +68,12 @@ function solveNucleonWf(κ, p::Bool, E, s::system; normalize=true, algo=Vern9())
     right_r = s.r_mesh.r[mid_idx:end]
 
     # left partition
-    prob = ODEProblem(dirac!, [0, 1], (0, r_mid))
+    prob = ODEProblem(dirac!, init_BC(), (0, r_mid))
     sol = solve(prob, algo, p=(κ, E, f1, f2), saveat=left_r)
     wf_left = hcat(sol.u...)
 
     # right partition
-    prob = ODEProblem(dirac!, [0, 1], (s.r_mesh.r_max, r_mid))
+    prob = ODEProblem(dirac!, asymp_BC(κ, p, E, s.r_mesh.r_max), (s.r_mesh.r_max, r_mid))
     sol = solve(prob, algo, p=(κ, E, f1, f2), saveat=right_r)
     wf_right = reverse(hcat(sol.u...); dims=2)