anchor.diff:    
     !DEC$ FORCEINLINE
     call multiply_poly(Ix_pol,n_Ix,Iy_pol,n_Iy,d_poly,n_pt_tmp)
  +  if (n_pt_tmp == -1) then
  +    return
  +  endif
     n_pt_out = n_pt_tmp+n_Iz
     do i=0,n_pt_out
       d1(i)=0.d0
  
positive.code:       L_power(p) = ao_power(l,p)
      enddo
      double  precision              :: ERI
      do_p = .True.
      do p = 1, ao_prim_num(i)
        if (.not.do_p) exit
        do_p = .False.
        do_q = .True.
        coef1 = ao_coef_transp(p,i)
        do q = 1, ao_prim_num(j)
          if (.not.do_q) exit
          do_q = .False.
          do_r = .True.
          coef2 = coef1*ao_coef_transp(q,j)
          do r = 1, ao_prim_num(k)
            if (.not.do_r) exit
            do_r = .False.
            coef3 = coef2*ao_coef_transp(r,k)
            do s = 1, ao_prim_num(l)
              coef4 = coef3*ao_coef_transp(s,l)
  
negative.code:   num_k = ao_nucl(k)
    num_l = ao_nucl(l)
    ao_bielec_integral = 0.d0
    
    if (num_i /= num_j .or. num_k /= num_l .or. num_j /= num_k)then
      do p = 1, 3
  