35   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in ComputeSGS_DDAnisotropicTrainingDataNodal()  local
36     const CeedScalar qi[5]                 = {q[0][i], q[1][i], q[2][i], q[3][i], q[4][i]};  in ComputeSGS_DDAnisotropicTrainingDataNodal()
38         {grad_velo[0][0][i], grad_velo[0][1][i], grad_velo[0][2][i]},  in ComputeSGS_DDAnisotropicTrainingDataNodal()
39         {grad_velo[1][0][i], grad_velo[1][1][i], grad_velo[1][2][i]},  in ComputeSGS_DDAnisotropicTrainingDataNodal()
40         {grad_velo[2][0][i], grad_velo[2][1][i], grad_velo[2][2][i]}  in ComputeSGS_DDAnisotropicTrainingDataNodal()
42 …ar km_A_ij[6] = {A_ij_delta[0][i], A_ij_delta[1][i], A_ij_delta[2][i], A_ij_delta[3][i], A_ij_delt…  in ComputeSGS_DDAnisotropicTrainingDataNodal()
43     const CeedScalar delta      = A_ij_delta[6][i];  in ComputeSGS_DDAnisotropicTrainingDataNodal()
50     for (int j = 0; j < 6; j++) v[j][i] = inv_multiplicity[i] * inputs[j];  in ComputeSGS_DDAnisotropicTrainingDataNodal()
52 …v[0 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_XX][i] - Square(s.Y.velocity[0])) * inv_mult…  in ComputeSGS_DDAnisotropicTrainingDataNodal()
53 …v[1 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_YY][i] - Square(s.Y.velocity[1])) * inv_mult…  in ComputeSGS_DDAnisotropicTrainingDataNodal()
54 …v[2 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_ZZ][i] - Square(s.Y.velocity[2])) * inv_mult…  in ComputeSGS_DDAnisotropicTrainingDataNodal()
55 …v[3 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_YZ][i] - s.Y.velocity[1] * s.Y.velocity[2]) …  in ComputeSGS_DDAnisotropicTrainingDataNodal()
56 …v[4 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_XZ][i] - s.Y.velocity[0] * s.Y.velocity[2]) …  in ComputeSGS_DDAnisotropicTrainingDataNodal()
57 …v[5 + 6][i] = (velo_prod[DIFF_FILTER_VELOCITY_SQUARED_XY][i] - s.Y.velocity[0] * s.Y.velocity[1]) …  in ComputeSGS_DDAnisotropicTrainingDataNodal()