anchor.diff:    
     
     implicit none
  +  include 'Utils/constants.include.F'
     double precision               :: p,q
     integer                        :: a_x,b_x,c_x,d_x,a_y,b_y,c_y,d_y,a_z,b_z,c_z,d_z
     integer                        :: i, n_iter, n_pt, j
  
positive.code:   p01_1 = 0.5d0/q
    p10_2 = 0.5d0 *  q /(p * q + p * p)
    p01_2 = 0.5d0 *  p /(q * q + q * p)
    double precision               :: B10(n_pt), B01(n_pt), B00(n_pt)
    double precision               :: t1(n_pt), t2(n_pt)
    !DIR$ ATTRIBUTES ALIGN : $IRP_ALIGN :: t1, t2, B10, B01, B00
    ix = a_x+b_x
    jx = c_x+d_x
  
negative.code:   BEGIN_DOC
    !  recursive function involved in the bielectronic integral
    END_DOC
    include 'Utils/constants.include.F'
    integer, intent(in)            :: c, n_pt
    double precision, intent(in)   :: B_10(n_pt_max_integrals),B_01(n_pt_max_integrals),B_00(n_pt_max_integrals)
    double precision, intent(out)  :: res(n_pt_max_integrals)
    integer                        :: i
    
    if(c==1)then
  