From 9c054580c8353fc7f62b09d6fd65aba34f1751e7 Mon Sep 17 00:00:00 2001
From: "Ryan@Home" <goluckyryan@gmail.com>
Date: Sun, 23 Feb 2025 06:05:24 -0500
Subject: [PATCH] the bottleneck is the 9jsymbol

---
 Raphael/DWBA_ZR.py | 2 +-
 Raphael/dwba_zr.py | 4 ++--
 2 files changed, 3 insertions(+), 3 deletions(-)

diff --git a/Raphael/DWBA_ZR.py b/Raphael/DWBA_ZR.py
index df47c03..de7b0c7 100755
--- a/Raphael/DWBA_ZR.py
+++ b/Raphael/DWBA_ZR.py
@@ -33,7 +33,7 @@ angList = []
 xsec = []
 
 start_time = time.time()
-for i in range(0, 181, 5):
+for i in range(0, 181, 10):
   angList.append(i)
   kaka = haha.AngDist(i)
   xsec.append(kaka)
diff --git a/Raphael/dwba_zr.py b/Raphael/dwba_zr.py
index 673859d..7d726a6 100755
--- a/Raphael/dwba_zr.py
+++ b/Raphael/dwba_zr.py
@@ -299,11 +299,11 @@ class DWBA_ZR:
       plt.show(block=False)
       input("Press Enter to continue...")
 
-  def Gamma(self, L1:int, J1, L2:int, J2, m, ma, mb):
+  def Gamma(self, L1:int, J1, L2:int, J2, m:int, ma, mb):
     if  int(L1 + L2 + self.l)%2 != 0: #check if the sum of L1 + L2 + l is even
         return 0
     else:
-        fact0 = wigner_9j(S(2*self.j)/2, S(2*self.l)/2, S(2*self.s)/2, S(2*J1)/2, S(2*L1)/2, S(2*self.spin_a)/2, S(2*J2)/2, S(2*L2)/2, S(2*self.spin_b)/2).evalf()
+        fact0 = wigner_9j(S(2*self.j)/2, self.l, S(2*self.s)/2, S(2*J1)/2, L1, S(2*self.spin_a)/2, J2, S(2*L2)/2, S(2*self.spin_b)/2).evalf()
         if fact0 == 0:
             return 0
         else: