__global__ void vreman_nut_compute(const T *__restrict__ a11, const T *__restrict__ a12, const T *__restrict__ a13, const T *__restrict__ a21, const T *__restrict__ a22, const T *__restrict__ a23, const T *__restrict__ a31, const T *__restrict__ a32, const T *__restrict__ a33, const T *__restrict__ delta, T *__restrict__ nut, const T *__restrict__ mult, const T c, const T eps, const int n)