src/ao_two_e_erf_ints/two_e_integrals_erf.irp.f

@@ -16,8 +16,6 @@ double precision function ao_two_e_integral_erf(i,j,k,l)
   integer                        :: iorder_p(3), iorder_q(3)
   double precision               :: ao_two_e_integral_schwartz_accel_erf
 
-  provide mu_erf
-
    if (ao_prim_num(i) * ao_prim_num(j) * ao_prim_num(k) * ao_prim_num(l) > 1024 ) then
      ao_two_e_integral_erf = ao_two_e_integral_schwartz_accel_erf(i,j,k,l)
      return