anchor.diff:            !id(4) = j4 ! 1234 are only used once and we know l=4, hence ijk must not be 4
             iax = d*(j3-1+c*(j2-1+b*(j1-1)))
             ibx = d*(id(k)-1+jd(k)*(id(j)-1+jd(j)*(id(i)-1)))
  +!dec$ vector always nontemporal
  +!dir$ simd
             do j4 = 1,d
              ia = j4+iax
              ib = j4+ibx
  
positive.code:           id(3) = j3
            id(4) = j4
            ibx = a*(id(k)-1+jd(k)*(id(j)-1+jd(j)*(id(i)-1)))
            do j1 = 1,a
             ia = j4+d*(j3-1+c*(j2-1+b*(j1-1)))
             ib = j1+ibx
  
negative.code:          enddo
          enddo
         enddo
  !$omp end parallel do
        else if (l.eq.2) then
  !$omp  parallel do private(j1,j2,j3,j4,iax,ibx,ia,ib,id) shared(s,u)
  !$omp& firstprivate(a,b,c,d,i,j,k,l,f,jd)
         do j1 = 1,a
          do j3 = 1,c
           do j4 = 1,d
  