511 CEED_QFUNCTION_HELPER void VelocityGradient(const State grad_s[3], CeedScalar grad_velocity[3][3]) {  in VelocityGradient()
512   grad_velocity[0][0] = grad_s[0].Y.velocity[0];  in VelocityGradient()
513   grad_velocity[0][1] = grad_s[1].Y.velocity[0];  in VelocityGradient()
514   grad_velocity[0][2] = grad_s[2].Y.velocity[0];  in VelocityGradient()
515   grad_velocity[1][0] = grad_s[0].Y.velocity[1];  in VelocityGradient()
516   grad_velocity[1][1] = grad_s[1].Y.velocity[1];  in VelocityGradient()
517   grad_velocity[1][2] = grad_s[2].Y.velocity[1];  in VelocityGradient()
518   grad_velocity[2][0] = grad_s[0].Y.velocity[2];  in VelocityGradient()
519   grad_velocity[2][1] = grad_s[1].Y.velocity[2];  in VelocityGradient()
520   grad_velocity[2][2] = grad_s[2].Y.velocity[2];  in VelocityGradient()
523 CEED_QFUNCTION_HELPER void KMStrainRate(const CeedScalar grad_velocity[3][3], CeedScalar strain_rat…  in KMStrainRate()
525   strain_rate[0]          = grad_velocity[0][0];  in KMStrainRate()
526   strain_rate[1]          = grad_velocity[1][1];  in KMStrainRate()
527   strain_rate[2]          = grad_velocity[2][2];  in KMStrainRate()
528   strain_rate[3]          = weight * (grad_velocity[1][2] + grad_velocity[2][1]);  in KMStrainRate()
529   strain_rate[4]          = weight * (grad_velocity[0][2] + grad_velocity[2][0]);  in KMStrainRate()
530   strain_rate[5]          = weight * (grad_velocity[0][1] + grad_velocity[1][0]);  in KMStrainRate()
535   CeedScalar grad_velocity[3][3];  in KMStrainRate_State()  local
536   VelocityGradient(grad_s, grad_velocity);  in KMStrainRate_State()
537   KMStrainRate(grad_velocity, strain_rate);  in KMStrainRate_State()
541 CEED_QFUNCTION_HELPER void RotationRate(const CeedScalar grad_velocity[3][3], CeedScalar rotation_r…  in RotationRate()
545   rotation_rate[1][2] = 0.5 * (grad_velocity[1][2] - grad_velocity[2][1]);  in RotationRate()
546   rotation_rate[0][2] = 0.5 * (grad_velocity[0][2] - grad_velocity[2][0]);  in RotationRate()
547   rotation_rate[0][1] = 0.5 * (grad_velocity[0][1] - grad_velocity[1][0]);  in RotationRate()
576   CeedScalar grad_velocity[3][3];  in Vorticity()  local
577   VelocityGradient(grad_s, grad_velocity);  in Vorticity()
578   Curl3(grad_velocity, vorticity);  in Vorticity()