anchor.diff:    PROVIDE ref_bitmask_energy
   
     allocate (shortcut(0:n+1,2), sort_idx(n,2), sorted(Nint,n,2), version(Nint,n,2))
  -  allocate(ut(N_st_8,n))
  +  allocate( ut(N_st_8,n))
   
     v_0 = 0.d0
     s_0 = 0.d0
  
positive.code:   enddo
    !$OMP END DO
  
    !$OMP DO SCHEDULE(dynamic)
    do sh=1,shortcut(0,1)
      do sh2=1,shortcut(0,1)
        if (sh==sh2) cycle
  
negative.code:   enddo
    !$OMP END DO
  
    !$OMP DO SCHEDULE(dynamic)
    do sh=1,shortcut(0,2)
      do i=shortcut(sh,2),shortcut(sh+1,2)-1
        org_i = sort_idx(i,2)
  