anchor.diff:             id(4) = j4
              !id(5) = j5
              id(6) = j6
  +           iax = -1+e*(j4-1+d*(j3-1+c*(j2-1+b*(j1-1))))
  +           ibx = e*(id(m)-1+jd(m)*(id(l)-1+jd(l)*(id(k)-1
  +     &           +jd(k)*(id(j)-1+jd(j)*(id(i)-1)))))
              do j5 = 1,e
  -            ia = j6+f*(j5-1+e*(j4-1+d*(j3-1+c*(j2-1+b*(j1-1)))))
  -            ib = j5+e*(id(m)-1+jd(m)*(id(l)-1+jd(l)
  -     1         *(id(k)-1+jd(k)*(id(j)-1+jd(j)*(id(i)-1)))))
  +            ia = j6+f*(j5+iax)
  +            ib = j5+ibx
               s(ib) = u(ia) * factor
              enddo
             enddo
  
positive.code:            id(4) = j4
             id(5) = j5
             !id(6) = j6 ! 123456 are only used once and we know n=6, hence ijklm must not be 6
             ia12345 = f*(j5-1+e*(j4-1+d*(j3-1+c*(j2-1+b*(j1-1)))))
             ib12345 = f*(id(m)-1+jdm*(id(l)-1+jdl*(id(k)-1+jdk
       &                *(id(j)-1+jdj*(id(i)-1)))))
             do j6 = 1,f
              ia = j6+ia12345
              ib = j6+ib12345
              s(ib) = u(ia) * factor
             enddo
            enddo
  
negative.code:            id(4) = j4
             id(5) = j5
             id(6) = j6
             ibx = a*(id(m)-1+jd(m)*(id(l)-1+jd(l)*(id(k)-1
       &           +jd(k)*(id(j)-1+jd(j)*(id(i)-1)))))
             do j1 = 1,a
              ia = j6+f*(j5-1+e*(j4-1+d*(j3-1+c*(j2-1+b*(j1-1)))))
              ib = j1+ibx
              s(ib) = u(ia) * factor
             enddo
            enddo
  