anchor.diff:      
     else
       
  +
  +    integer, save :: icount(4) = (/0,0,0,0/)
       !DIR$ LOOP COUNT (1000)
       outer: do i=1,sze
         degree_x2 = 0
  
positive.code:           popcnt(xor( key1(2,2,i), key2(2,2))) +                     &
            popcnt(xor( key1(3,1,i), key2(3,1))) +                     &
            popcnt(xor( key1(3,2,i), key2(3,2)))
        if (degree_x2 > 4) then
          cycle
        else if(degree_x2 /= 0)then
          idx(l) = i
          l = l+1
        endif
  
negative.code:           popcnt(xor( key1(2,1,i), key2(2,1))) +                     &
            popcnt(xor( key1(1,2,i), key2(1,2))) +                     &
            popcnt(xor( key1(2,2,i), key2(2,2)))
        if (degree_x2 <= 4) then
          idx(l) = i
          l = l+1
        endif
  