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