198     for (CeedInt j = 0; j < 3; j++) {    // Component  in ElasFSResidual_NH()  local
200         grad_u[j][k][i] = 0;  in ElasFSResidual_NH()
201         for (CeedInt m = 0; m < 3; m++) grad_u[j][k][i] += dXdx[m][k] * du[j][m];  in ElasFSResidual_NH()
231     for (CeedInt j = 0; j < 3; j++) {  in ElasFSResidual_NH()  local
233         P[j][k] = 0;  in ElasFSResidual_NH()
234         for (CeedInt m = 0; m < 3; m++) P[j][k] += F[j][m] * S[m][k];  in ElasFSResidual_NH()
239     for (CeedInt j = 0; j < 3; j++) {    // Component  in ElasFSResidual_NH()  local
241         dvdX[k][j][i] = 0;  in ElasFSResidual_NH()
242         for (CeedInt m = 0; m < 3; m++) dvdX[k][j][i] += dXdx[k][m] * P[j][m] * wdetJ;  in ElasFSResidual_NH()
294     for (CeedInt j = 0; j < 3; j++) {    // Component  in ElasFSJacobian_NH()  local
296         graddeltau[j][k] = 0;  in ElasFSJacobian_NH()
297         for (CeedInt m = 0; m < 3; m++) graddeltau[j][k] += dXdx[m][k] * deltadu[j][m];  in ElasFSJacobian_NH()
350     for (CeedInt j = 0; j < 3; j++) {  in ElasFSJacobian_NH()  local
351       for (CeedInt k = 0; k < 3; k++) Cinv_contract_E += C_inv[j][k] * deltaE[j][k];  in ElasFSJacobian_NH()
355     for (CeedInt j = 0; j < 3; j++) {  in ElasFSJacobian_NH()  local
357         deltaECinv[j][k] = 0;  in ElasFSJacobian_NH()
358         for (CeedInt m = 0; m < 3; m++) deltaECinv[j][k] += deltaE[j][m] * C_inv[m][k];  in ElasFSJacobian_NH()
363     for (CeedInt j = 0; j < 3; j++) {  in ElasFSJacobian_NH()  local
365         deltaS[j][k] = 0;  in ElasFSJacobian_NH()
366         for (CeedInt m = 0; m < 3; m++) deltaS[j][k] += C_inv[j][m] * deltaECinv[m][k];  in ElasFSJacobian_NH()
370     for (CeedInt j = 0; j < 3; j++) {  in ElasFSJacobian_NH()  local
371 …CeedInt k = 0; k < 3; k++) deltaS[j][k] = lambda * Cinv_contract_E * C_inv[j][k] - 2. * (lambda * …  in ElasFSJacobian_NH()
376     for (CeedInt j = 0; j < 3; j++) {  in ElasFSJacobian_NH()  local
378         deltaP[j][k] = 0;  in ElasFSJacobian_NH()
379 …for (CeedInt m = 0; m < 3; m++) deltaP[j][k] += graddeltau[j][m] * S[m][k] + F[j][m] * deltaS[m][k…  in ElasFSJacobian_NH()
384     for (CeedInt j = 0; j < 3; j++) {    // Component  in ElasFSJacobian_NH()  local
386         deltadvdX[k][j][i] = 0;  in ElasFSJacobian_NH()
387         for (CeedInt m = 0; m < 3; m++) deltadvdX[k][j][i] += dXdx[k][m] * deltaP[j][m] * wdetJ;  in ElasFSJacobian_NH()
434     for (int j = 0; j < 3; j++) {    // Component  in ElasFSEnergy_NH()  local
436         grad_u[j][k] = 0;  in ElasFSEnergy_NH()
437         for (int m = 0; m < 3; m++) grad_u[j][k] += dXdx[m][k] * du[j][m];  in ElasFSEnergy_NH()
504     for (int j = 0; j < 3; j++) {    // Component  in ElasFSDiagnostic_NH()  local
506         grad_u[j][k] = 0;  in ElasFSDiagnostic_NH()
507         for (int m = 0; m < 3; m++) grad_u[j][k] += dXdx[m][k] * du[j][m];  in ElasFSDiagnostic_NH()
538     for (CeedInt j = 0; j < 3; j++) {  in ElasFSDiagnostic_NH()  local
539       for (CeedInt m = 0; m < 3; m++) diagnostic[5][i] += E2[j][m] * E2[m][j] / 4.;  in ElasFSDiagnostic_NH()