anchor.diff:    
           C_center(1:3) = nucl_coord(k,1:3)
     
  -        v_k_dump = v_k(k,1:klocmax)
  -        n_k_dump = n_k(k,1:klocmax)
  -        dz_k_dump =  dz_k(k,1:klocmax)
  +        v_k_dump = pseudo_v_k(k,1:pseudo_klocmax)
  +        n_k_dump = pseudo_n_k(k,1:pseudo_klocmax)
  +        dz_k_dump =  pseudo_dz_k(k,1:pseudo_klocmax)
   
  -        c = c + Vloc(klocmax, v_k_dump,n_k_dump, dz_k_dump, &
  +        c = c + Vloc(pseudo_klocmax, v_k_dump,n_k_dump, dz_k_dump, &
                        A_center,power_A,alpha,B_center,power_B,beta,C_center)
     
  +
  +        n_kl_dump = pseudo_n_kl(k,1:pseudo_kmax,0:pseudo_lmax)
  +        v_kl_dump = pseudo_v_kl(k,1:pseudo_kmax,0:pseudo_lmax)
  +        dz_kl_dump = pseudo_dz_kl(k,1:pseudo_kmax,0:pseudo_lmax)
  +
  +        c = c + Vpseudo(pseudo_lmax,pseudo_kmax,v_kl_dump,n_kl_dump,dz_kl_dump,A_center,power_A,alpha,B_center,power_B,beta,C_center)
  +  
         enddo
  -      ao_nucl_elec_integral_pseudo_local(i,j) = ao_nucl_elec_integral_pseudo_local(i,j) + &
  +      ao_nucl_elec_integral_pseudo(i,j) = ao_nucl_elec_integral_pseudo(i,j) + &
                                      ao_coef_normalized_ordered_transp(l,j)*ao_coef_normalized_ordered_transp(m,i)*c
        enddo
        enddo
  
positive.code:   !$OMP DEFAULT (NONE) &
    !$OMP PRIVATE (i,j,k,l,m,alpha,beta,A_center,B_center,C_center,power_A,power_B, &
    !$OMP          num_A,num_B,Z,c,n_pt_in, &
    !$OMP          v_k_dump,n_k_dump, dz_k_dump, &
    !$OMP          wall_0,wall_2,thread_num, output_monoints) & 
    !$OMP SHARED (ao_num,ao_prim_num,ao_expo_ordered_transp,ao_power,ao_nucl,nucl_coord,ao_coef_normalized_ordered_transp, &
    !$OMP         ao_nucl_elec_integral_pseudo_local,nucl_num,nucl_charge, &
    !$OMP         klocmax,v_k,n_k, dz_k, &
    !$OMP         wall_1)
    
    !$OMP DO SCHEDULE (guided)
  
negative.code:  BEGIN_PROVIDER [ double precision, ao_nucl_elec_integral_pseudo, (ao_num_align,ao_num)]
   BEGIN_DOC
  ! interaction nuclear electron
   END_DOC
   implicit none
   double precision  :: alpha, beta, gama, delta
   integer           :: num_A,num_B
   double precision  :: A_center(3),B_center(3),C_center(3)
  