43   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in ElasResidual_Linear()  local
46         {ug[0][0][i], ug[1][0][i], ug[2][0][i]},  in ElasResidual_Linear()
47         {ug[0][1][i], ug[1][1][i], ug[2][1][i]},  in ElasResidual_Linear()
48         {ug[0][2][i], ug[1][2][i], ug[2][2][i]}  in ElasResidual_Linear()
51     const CeedScalar wdetJ      = q_data[0][i];  in ElasResidual_Linear()
53         {q_data[1][i], q_data[2][i], q_data[3][i]},  in ElasResidual_Linear()
54         {q_data[4][i], q_data[5][i], q_data[6][i]},  in ElasResidual_Linear()
55         {q_data[7][i], q_data[8][i], q_data[9][i]}  in ElasResidual_Linear()
112         dvdX[k][j][i] = 0;  in ElasResidual_Linear()
113         for (CeedInt m = 0; m < 3; m++) dvdX[k][j][i] += dXdx[k][m] * sigma[j][m] * wdetJ;  in ElasResidual_Linear()
140   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in ElasJacobian_Linear()  local
143         {deltaug[0][0][i], deltaug[1][0][i], deltaug[2][0][i]},  in ElasJacobian_Linear()
144         {deltaug[0][1][i], deltaug[1][1][i], deltaug[2][1][i]},  in ElasJacobian_Linear()
145         {deltaug[0][2][i], deltaug[1][2][i], deltaug[2][2][i]}  in ElasJacobian_Linear()
148     const CeedScalar wdetJ      = q_data[0][i];  in ElasJacobian_Linear()
150         {q_data[1][i], q_data[2][i], q_data[3][i]},  in ElasJacobian_Linear()
151         {q_data[4][i], q_data[5][i], q_data[6][i]},  in ElasJacobian_Linear()
152         {q_data[7][i], q_data[8][i], q_data[9][i]}  in ElasJacobian_Linear()
207         deltadvdX[k][j][i] = 0;  in ElasJacobian_Linear()
208         for (CeedInt m = 0; m < 3; m++) deltadvdX[k][j][i] += dXdx[k][m] * dsigma[j][m] * wdetJ;  in ElasJacobian_Linear()
238   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in ElasEnergy_Linear()  local
241         {ug[0][0][i], ug[1][0][i], ug[2][0][i]},  in ElasEnergy_Linear()
242         {ug[0][1][i], ug[1][1][i], ug[2][1][i]},  in ElasEnergy_Linear()
243         {ug[0][2][i], ug[1][2][i], ug[2][2][i]}  in ElasEnergy_Linear()
246     const CeedScalar wdetJ      = q_data[0][i];  in ElasEnergy_Linear()
248         {q_data[1][i], q_data[2][i], q_data[3][i]},  in ElasEnergy_Linear()
249         {q_data[4][i], q_data[5][i], q_data[6][i]},  in ElasEnergy_Linear()
250         {q_data[7][i], q_data[8][i], q_data[9][i]}  in ElasEnergy_Linear()
275     energy[i] =  in ElasEnergy_Linear()
306   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) {  in ElasDiagnostic_Linear()  local
309         {ug[0][0][i], ug[1][0][i], ug[2][0][i]},  in ElasDiagnostic_Linear()
310         {ug[0][1][i], ug[1][1][i], ug[2][1][i]},  in ElasDiagnostic_Linear()
311         {ug[0][2][i], ug[1][2][i], ug[2][2][i]}  in ElasDiagnostic_Linear()
315         {q_data[1][i], q_data[2][i], q_data[3][i]},  in ElasDiagnostic_Linear()
316         {q_data[4][i], q_data[5][i], q_data[6][i]},  in ElasDiagnostic_Linear()
317         {q_data[7][i], q_data[8][i], q_data[9][i]}  in ElasDiagnostic_Linear()
341     diagnostic[0][i] = u[0][i];  in ElasDiagnostic_Linear()
342     diagnostic[1][i] = u[1][i];  in ElasDiagnostic_Linear()
343     diagnostic[2][i] = u[2][i];  in ElasDiagnostic_Linear()
347     diagnostic[3][i]            = -lambda * strain_vol;  in ElasDiagnostic_Linear()
350     diagnostic[4][i] = strain_vol;  in ElasDiagnostic_Linear()
351     diagnostic[5][i] = 0.;  in ElasDiagnostic_Linear()
353       for (CeedInt m = 0; m < 3; m++) diagnostic[5][i] += e[j][m] * e[m][j];  in ElasDiagnostic_Linear()
355     diagnostic[6][i] = 1 + strain_vol;  in ElasDiagnostic_Linear()
358     diagnostic[7][i] =  in ElasDiagnostic_Linear()