anchor.diff:        ik=ik_table(ikb)
         ib=ib_table(ikb)
         idx = idx_b + ikb1
  -
         call hpsi_omp_KB_RT(ik,ztpsi(:,4,idx),ztpsi(:,1,idx))
         call hpsi_omp_KB_RT(ik,ztpsi(:,1,idx),ztpsi(:,2,idx))
         call hpsi_omp_KB_RT(ik,ztpsi(:,2,idx),ztpsi(:,3,idx))
         call hpsi_omp_KB_RT(ik,ztpsi(:,3,idx),ztpsi(:,4,idx))
  +    enddo
  +
  +#if 1
  +    call update_LBLK(zfac,ztpsi(:,:,idx_b),zu(:,:,:), ikb_s,ikb_e)
  +#else
  +    do ikb1 = 0, num_ikb1-1
  +      ikb = ikb0 + ikb1
  +      ik=ik_table(ikb)
  +      ib=ib_table(ikb)
  +      idx = idx_b + ikb1
         call update(zfac,ztpsi(:,:,idx),zu(:,ib,ik))
  +    enddo
  +#endif
   
   #ifdef ARTED_CURRENT_OPTIMIZED
  +    do ikb1 = 0, num_ikb1-1
  +      ikb = ikb0 + ikb1
  +      ik=ik_table(ikb)
  +      ib=ib_table(ikb)
         call current_omp_KB_ST(ib,ik,zu(:,ib,ik))
  -#endif
  -    enddo
  +    enddo
  +#endif
   
   #ifdef ARTED_SC
       loop_count = loop_count + num_ikb1
  
positive.code:     end do
      TIMELOG_END(LOG_HPSI_UPDATE)
    end subroutine
  end subroutine
  
  
negative.code:     end do
      TIMELOG_END(LOG_HPSI_INIT)
    end subroutine
  
    subroutine update(zfac,tpsi,zu)
      use Global_Variables, only: NLx,NLy,NLz
  