anchor.diff:      do j = 1, k
           call xmm2(tm1(1,1,j), dy, tm2(1,1,j), xargs)
       enddo
  -    CALL xmm3(tm2, dy, c(1,1,1,i), xargs)
  +    CALL xmm3(tm2, dz, c(1,1,1,i), xargs)
     END DO
     !$OMP MASTER
     duration = libxsmm_timer_duration(start, libxsmm_timer_tick())
  
positive.code:     enddo
    END DO 
    dx = 1.
    dy = transpose(dx)
  
    WRITE(*, "(6(A,I0),A,I0)") "m=", m, " n=", n, " k=", k, " mm=", mm, " nn=", nn, " kk=", kk, " size=", UBOUND(a, 4) 
  
  
negative.code:     do j = 1, k
          call mxmf2(tm1(:,:,j), mm, dy, n, tm2(:,:,j), nn, xargs)
      enddo
      call mxmf2(tm2, mm*nn, dz, k, c(:,:,:,i), kk, xargs)
    END DO
    !$OMP MASTER
    duration = libxsmm_timer_duration(start, libxsmm_timer_tick())
  