few optim in cplx Boys fct

Author: AbdAmmar

src/ao_one_e_ints/one_e_Coul_integrals_cosgtos.irp.f

@@ -111,8 +111,9 @@
   complex*16                   :: accu, P_center(3)
   complex*16                   :: d(0:n_pt_in)
 
-  complex*16                   :: V_n_e_cosgtos
-  complex*16                   :: crint_2
+  complex*16, external         :: V_n_e_cosgtos
+  complex*16, external         :: crint_2
+  complex*16, external         :: crint_sum_2
 
   if ( (A_center(1)/=B_center(1)) .or. (A_center(2)/=B_center(2)) .or. (A_center(3)/=B_center(3)) .or. &
        (A_center(1)/=C_center(1)) .or. (A_center(2)/=C_center(2)) .or. (A_center(3)/=C_center(3)) ) then
@@ -162,22 +163,22 @@
     return
   endif
 
-  call give_cpolynomial_mult_center_one_e( A_center, B_center, alpha, beta &
-                                         , power_A, power_B, C_center, n_pt_in, d, n_pt_out)
+  call give_cpolynomial_mult_center_one_e(A_center, B_center, alpha, beta, &
+                                          power_A, power_B, C_center, n_pt_in, d, n_pt_out)
 
   if(n_pt_out < 0) then
     NAI_pol_mult_cosgtos = (0.d0, 0.d0)
     return
   endif
 
-  accu = (0.d0, 0.d0)
-  do i = 0, n_pt_out, 2
-    accu += crint_2(shiftr(i, 1), const) * d(i)
-
-!    print *, shiftr(i, 1), real(const), real(d(i)), real(crint_2(shiftr(i, 1), const))
-  enddo
+  !accu = (0.d0, 0.d0)
+  !do i = 0, n_pt_out, 2
+  !  accu += crint_2(shiftr(i, 1), const) * d(i)
+  !enddo
+  accu = crint_sum_2(n_pt_out, const, d)
   NAI_pol_mult_cosgtos = accu * coeff
 
+  return
 end
 
 ! ---