514 CEED_QFUNCTION_HELPER void VelocityGradient(const State grad_s[3], CeedScalar grad_velocity[3][3]) {  in VelocityGradient()
515   grad_velocity[0][0] = grad_s[0].Y.velocity[0];  in VelocityGradient()
516   grad_velocity[0][1] = grad_s[1].Y.velocity[0];  in VelocityGradient()
517   grad_velocity[0][2] = grad_s[2].Y.velocity[0];  in VelocityGradient()
518   grad_velocity[1][0] = grad_s[0].Y.velocity[1];  in VelocityGradient()
519   grad_velocity[1][1] = grad_s[1].Y.velocity[1];  in VelocityGradient()
520   grad_velocity[1][2] = grad_s[2].Y.velocity[1];  in VelocityGradient()
521   grad_velocity[2][0] = grad_s[0].Y.velocity[2];  in VelocityGradient()
522   grad_velocity[2][1] = grad_s[1].Y.velocity[2];  in VelocityGradient()
523   grad_velocity[2][2] = grad_s[2].Y.velocity[2];  in VelocityGradient()
526 CEED_QFUNCTION_HELPER void KMStrainRate(const CeedScalar grad_velocity[3][3], CeedScalar strain_rat…  in KMStrainRate()
528   strain_rate[0]          = grad_velocity[0][0];  in KMStrainRate()
529   strain_rate[1]          = grad_velocity[1][1];  in KMStrainRate()
530   strain_rate[2]          = grad_velocity[2][2];  in KMStrainRate()
531   strain_rate[3]          = weight * (grad_velocity[1][2] + grad_velocity[2][1]);  in KMStrainRate()
532   strain_rate[4]          = weight * (grad_velocity[0][2] + grad_velocity[2][0]);  in KMStrainRate()
533   strain_rate[5]          = weight * (grad_velocity[0][1] + grad_velocity[1][0]);  in KMStrainRate()
538   CeedScalar grad_velocity[3][3];  in KMStrainRate_State()  local
539   VelocityGradient(grad_s, grad_velocity);  in KMStrainRate_State()
540   KMStrainRate(grad_velocity, strain_rate);  in KMStrainRate_State()
544 CEED_QFUNCTION_HELPER void RotationRate(const CeedScalar grad_velocity[3][3], CeedScalar rotation_r…  in RotationRate()
548   rotation_rate[1][2] = 0.5 * (grad_velocity[1][2] - grad_velocity[2][1]);  in RotationRate()
549   rotation_rate[0][2] = 0.5 * (grad_velocity[0][2] - grad_velocity[2][0]);  in RotationRate()
550   rotation_rate[0][1] = 0.5 * (grad_velocity[0][1] - grad_velocity[1][0]);  in RotationRate()
579   CeedScalar grad_velocity[3][3];  in Vorticity()  local
580   VelocityGradient(grad_s, grad_velocity);  in Vorticity()
581   Curl3(grad_velocity, vorticity);  in Vorticity()