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:     do j = 1, k
          call mxmf2(tm1(:,:,j), mm, dy, n, tm2(:,:,j), nn, xargs)
      enddo
      call mxmf2(tm2, mm*nn, dy, k, c(:,:,:,i), kk, xargs)
    END DO
    !$OMP MASTER
    duration = libxsmm_timer_duration(start, libxsmm_timer_tick())
  
negative.code:     do j = 1, k
          call libxsmm_mm(mm, nn, n, tm1(:,:,j), dy, tm2(:,:,j), xargs)
      enddo
      call libxsmm_mm(mm*nn, kk, k, reshape(tm2, (/mm*nn,k/)), dz, c(:,:,1,i), xargs)
    END DO
    !$OMP MASTER
    duration = libxsmm_timer_duration(start, libxsmm_timer_tick())
  