__global__ void sigma_nut_compute(T *__restrict__ g11, T *__restrict__ g12, T *__restrict__ g13, T *__restrict__ g21, T *__restrict__ g22, T *__restrict__ g23, T *__restrict__ g31, T *__restrict__ g32, T *__restrict__ g33, T *__restrict__ delta, T *__restrict__ nut, T *__restrict__ mult, const T c, const T eps, const int n)