Lines Matching refs:grad_u
65 CEED_QFUNCTION_HELPER CeedScalar computeJM1(const CeedScalar grad_u[3][3]) { in computeJM1()
66 return grad_u[0][0] * (grad_u[1][1] * grad_u[2][2] - grad_u[1][2] * grad_u[2][1]) + in computeJM1()
67 grad_u[0][1] * (grad_u[1][2] * grad_u[2][0] - grad_u[1][0] * grad_u[2][2]) + in computeJM1()
68 …grad_u[0][2] * (grad_u[1][0] * grad_u[2][1] - grad_u[2][0] * grad_u[1][1]) + grad_u[0][0] + grad_u… in computeJM1()
69 …grad_u[0][0] * grad_u[1][1] + grad_u[0][0] * grad_u[2][2] + grad_u[1][1] * grad_u[2][2] - grad_u[0… in computeJM1()
70 grad_u[0][2] * grad_u[2][0] - grad_u[1][2] * grad_u[2][1]; in computeJM1()
98 … commonFS(const CeedScalar lambda, const CeedScalar mu, const CeedScalar grad_u[3][3], CeedScalar … in commonFS()
105 E2work[m] = grad_u[indj[m]][indk[m]] + grad_u[indk[m]][indj[m]]; in commonFS()
106 for (CeedInt n = 0; n < 3; n++) E2work[m] += grad_u[n][indj[m]] * grad_u[n][indk[m]]; in commonFS()
114 const CeedScalar Jm1 = computeJM1(grad_u); in commonFS()
154 CeedScalar(*grad_u)[3][CEED_Q_VLA] = (CeedScalar(*)[3][CEED_Q_VLA])out[1]; 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()
208 {grad_u[0][0][i] + 1, grad_u[0][1][i], grad_u[0][2][i] }, in ElasFSResidual_NH()
209 {grad_u[1][0][i], grad_u[1][1][i] + 1, grad_u[1][2][i] }, in ElasFSResidual_NH()
210 {grad_u[2][0][i], grad_u[2][1][i], grad_u[2][2][i] + 1} in ElasFSResidual_NH()
216 {grad_u[0][0][i], grad_u[0][1][i], grad_u[0][2][i]}, in ElasFSResidual_NH()
217 {grad_u[1][0][i], grad_u[1][1][i], grad_u[1][2][i]}, in ElasFSResidual_NH()
218 {grad_u[2][0][i], grad_u[2][1][i], grad_u[2][2][i]} in ElasFSResidual_NH()
258 const CeedScalar(*grad_u)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[2]; in ElasFSJacobian_NH() local
304 {grad_u[0][0][i] + 1, grad_u[0][1][i], grad_u[0][2][i] }, in ElasFSJacobian_NH()
305 {grad_u[1][0][i], grad_u[1][1][i] + 1, grad_u[1][2][i] }, in ElasFSJacobian_NH()
306 {grad_u[2][0][i], grad_u[2][1][i], grad_u[2][2][i] + 1} in ElasFSJacobian_NH()
312 {grad_u[0][0][i], grad_u[0][1][i], grad_u[0][2][i]}, in ElasFSJacobian_NH()
313 {grad_u[1][0][i], grad_u[1][1][i], grad_u[1][2][i]}, in ElasFSJacobian_NH()
314 {grad_u[2][0][i], grad_u[2][1][i], grad_u[2][2][i]} in ElasFSJacobian_NH()
433 CeedScalar grad_u[3][3]; 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()
446 E2work[m] = grad_u[indj[m]][indk[m]] + grad_u[indk[m]][indj[m]]; in ElasFSEnergy_NH()
447 for (CeedInt n = 0; n < 3; n++) E2work[m] += grad_u[n][indj[m]] * grad_u[n][indk[m]]; in ElasFSEnergy_NH()
454 const CeedScalar Jm1 = computeJM1(grad_u); in ElasFSEnergy_NH()
503 CeedScalar grad_u[3][3]; 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()
516 E2work[m] = grad_u[indj[m]][indk[m]] + grad_u[indk[m]][indj[m]]; in ElasFSDiagnostic_NH()
517 for (CeedInt n = 0; n < 3; n++) E2work[m] += grad_u[n][indj[m]] * grad_u[n][indk[m]]; in ElasFSDiagnostic_NH()
531 const CeedScalar Jm1 = computeJM1(grad_u); in ElasFSDiagnostic_NH()