199       for (CeedInt k = 0; k < 3; k++) {  // Derivative  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()
232       for (CeedInt k = 0; k < 3; k++) {  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()
240       for (CeedInt k = 0; k < 3; k++) {  // Derivative  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()
295       for (CeedInt k = 0; k < 3; k++) {  // Derivative  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()
351       for (CeedInt k = 0; k < 3; k++) Cinv_contract_E += C_inv[j][k] * deltaE[j][k];  in ElasFSJacobian_NH()  local
356       for (CeedInt k = 0; k < 3; k++) {  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()
364       for (CeedInt k = 0; k < 3; k++) {  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()
371 …for (CeedInt k = 0; k < 3; k++) deltaS[j][k] = lambda * Cinv_contract_E * C_inv[j][k] - 2. * (lamb…  in ElasFSJacobian_NH()  local
377       for (CeedInt k = 0; k < 3; k++) {  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()
385       for (CeedInt k = 0; k < 3; k++) {  // Derivative  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()
435       for (int k = 0; k < 3; k++) {  // Derivative  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()
505       for (int k = 0; k < 3; k++) {  // Derivative  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()