Lines Matching full:0
37 CeedScalar sum = 0; in log1p_series_shifted()
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()
83 A[0][0] * A[2][2] - A[0][2] * A[2][0], /* *NOPAD* */ in computeMatinvSym()
84 A[0][0] * A[1][1] - A[0][1] * A[1][0], /* *NOPAD* */ in computeMatinvSym()
85 A[0][2] * A[1][0] - A[0][0] * A[1][2], /* *NOPAD* */ in computeMatinvSym()
86 A[0][1] * A[1][2] - A[0][2] * A[1][1], /* *NOPAD* */ in computeMatinvSym()
87 A[0][2] * A[2][1] - A[0][1] * A[2][2] /* *NOPAD* */ in computeMatinvSym()
89 for (CeedInt m = 0; m < 6; m++) Ainv[m] = B[m] / (detA); in computeMatinvSym()
102 const CeedInt indj[6] = {0, 1, 2, 1, 0, 0}, indk[6] = {0, 1, 2, 2, 2, 1}; in commonFS()
104 for (CeedInt m = 0; m < 6; 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()
109 {E2work[0], E2work[5], E2work[4]}, in commonFS()
119 {1 + E2[0][0], E2[0][1], E2[0][2] }, in commonFS()
120 {E2[0][1], 1 + E2[1][1], E2[1][2] }, in commonFS()
121 {E2[0][2], E2[1][2], 1 + E2[2][2]} in commonFS()
129 {Cinvwork[0], Cinvwork[5], Cinvwork[4]}, in commonFS()
136 for (CeedInt m = 0; m < 6; m++) { in commonFS()
138 for (CeedInt n = 0; n < 3; n++) Swork[m] += mu * C_inv[indj[m]][n] * E2[n][indk[m]]; in commonFS()
149 …const CeedScalar(*ug)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0], (*q_data)[CEED_Q… in ElasFSResidual_NH()
152 CeedScalar(*dvdX)[3][CEED_Q_VLA] = (CeedScalar(*)[3][CEED_Q_VLA])out[0]; in ElasFSResidual_NH()
180 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ElasFSResidual_NH()
183 {ug[0][0][i], ug[1][0][i], ug[2][0][i]}, in ElasFSResidual_NH()
184 {ug[0][1][i], ug[1][1][i], ug[2][1][i]}, in ElasFSResidual_NH()
185 {ug[0][2][i], ug[1][2][i], ug[2][2][i]} in ElasFSResidual_NH()
188 const CeedScalar wdetJ = q_data[0][i]; in ElasFSResidual_NH()
198 for (CeedInt j = 0; j < 3; j++) { // Component in ElasFSResidual_NH()
199 for (CeedInt k = 0; k < 3; k++) { // Derivative in ElasFSResidual_NH()
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()
224 {Swork[0], Swork[5], Swork[4]}, in ElasFSResidual_NH()
231 for (CeedInt j = 0; j < 3; j++) { in ElasFSResidual_NH()
232 for (CeedInt k = 0; k < 3; k++) { in ElasFSResidual_NH()
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()
240 for (CeedInt k = 0; k < 3; k++) { // Derivative in ElasFSResidual_NH()
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()
255 const CeedScalar(*deltaug)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0], in ElasFSJacobian_NH()
261 CeedScalar(*deltadvdX)[3][CEED_Q_VLA] = (CeedScalar(*)[3][CEED_Q_VLA])out[0]; in ElasFSJacobian_NH()
275 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ElasFSJacobian_NH()
278 {deltaug[0][0][i], deltaug[1][0][i], deltaug[2][0][i]}, in ElasFSJacobian_NH()
279 {deltaug[0][1][i], deltaug[1][1][i], deltaug[2][1][i]}, in ElasFSJacobian_NH()
280 {deltaug[0][2][i], deltaug[1][2][i], deltaug[2][2][i]} in ElasFSJacobian_NH()
283 const CeedScalar wdetJ = q_data[0][i]; in ElasFSJacobian_NH()
294 for (CeedInt j = 0; j < 3; j++) { // Component in ElasFSJacobian_NH()
295 for (CeedInt k = 0; k < 3; k++) { // Derivative in ElasFSJacobian_NH()
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()
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()
319 const CeedInt indj[6] = {0, 1, 2, 1, 0, 0}, indk[6] = {0, 1, 2, 2, 2, 1}; in ElasFSJacobian_NH()
321 for (CeedInt m = 0; m < 6; m++) { in ElasFSJacobian_NH()
322 deltaEwork[m] = 0; in ElasFSJacobian_NH()
323 …for (CeedInt n = 0; n < 3; n++) deltaEwork[m] += (graddeltau[n][indj[m]] * F[n][indk[m]] + F[n][in… in ElasFSJacobian_NH()
326 {deltaEwork[0], deltaEwork[5], deltaEwork[4]}, in ElasFSJacobian_NH()
334 {Cinvwork[0], Cinvwork[5], Cinvwork[4]}, in ElasFSJacobian_NH()
341 {Swork[0], Swork[5], Swork[4]}, in ElasFSJacobian_NH()
349 CeedScalar Cinv_contract_E = 0; in ElasFSJacobian_NH()
350 for (CeedInt j = 0; j < 3; j++) { in ElasFSJacobian_NH()
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()
356 for (CeedInt k = 0; k < 3; k++) { in ElasFSJacobian_NH()
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()
364 for (CeedInt k = 0; k < 3; k++) { in ElasFSJacobian_NH()
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()
371 …for (CeedInt k = 0; k < 3; k++) deltaS[j][k] = lambda * Cinv_contract_E * C_inv[j][k] - 2. * (lamb… in ElasFSJacobian_NH()
376 for (CeedInt j = 0; j < 3; j++) { in ElasFSJacobian_NH()
377 for (CeedInt k = 0; k < 3; k++) { in ElasFSJacobian_NH()
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()
385 for (CeedInt k = 0; k < 3; k++) { // Derivative in ElasFSJacobian_NH()
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()
400 …const CeedScalar(*ug)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0], (*q_data)[CEED_Q… in ElasFSEnergy_NH()
403 CeedScalar(*energy) = (CeedScalar(*))out[0]; in ElasFSEnergy_NH()
415 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ElasFSEnergy_NH()
418 {ug[0][0][i], ug[1][0][i], ug[2][0][i]}, in ElasFSEnergy_NH()
419 {ug[0][1][i], ug[1][1][i], ug[2][1][i]}, in ElasFSEnergy_NH()
420 {ug[0][2][i], ug[1][2][i], ug[2][2][i]} in ElasFSEnergy_NH()
423 const CeedScalar wdetJ = q_data[0][i]; in ElasFSEnergy_NH()
434 for (int j = 0; j < 3; j++) { // Component in ElasFSEnergy_NH()
435 for (int k = 0; k < 3; k++) { // Derivative in ElasFSEnergy_NH()
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()
443 const CeedInt indj[6] = {0, 1, 2, 1, 0, 0}, indk[6] = {0, 1, 2, 2, 2, 1}; in ElasFSEnergy_NH()
445 for (CeedInt m = 0; m < 6; 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()
450 {E2work[0], E2work[5], E2work[4]}, in ElasFSEnergy_NH()
458 …energy[i] = (lambda * logJ * logJ / 2. - mu * logJ + mu * (E2[0][0] + E2[1][1] + E2[2][2]) / 2.) *… in ElasFSEnergy_NH()
470 …const CeedScalar(*u)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0], (*ug)[3][CEED_Q_VLA] = … in ElasFSDiagnostic_NH()
474 CeedScalar(*diagnostic)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in ElasFSDiagnostic_NH()
486 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ElasFSDiagnostic_NH()
489 {ug[0][0][i], ug[1][0][i], ug[2][0][i]}, in ElasFSDiagnostic_NH()
490 {ug[0][1][i], ug[1][1][i], ug[2][1][i]}, in ElasFSDiagnostic_NH()
491 {ug[0][2][i], ug[1][2][i], ug[2][2][i]} in ElasFSDiagnostic_NH()
504 for (int j = 0; j < 3; j++) { // Component in ElasFSDiagnostic_NH()
505 for (int k = 0; k < 3; k++) { // Derivative in ElasFSDiagnostic_NH()
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()
513 const CeedInt indj[6] = {0, 1, 2, 1, 0, 0}, indk[6] = {0, 1, 2, 2, 2, 1}; in ElasFSDiagnostic_NH()
515 for (CeedInt m = 0; m < 6; 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()
520 {E2work[0], E2work[5], E2work[4]}, in ElasFSDiagnostic_NH()
526 diagnostic[0][i] = u[0][i]; in ElasFSDiagnostic_NH()
536 diagnostic[4][i] = (E2[0][0] + E2[1][1] + E2[2][2]) / 2.; in ElasFSDiagnostic_NH()
537 diagnostic[5][i] = 0.; in ElasFSDiagnostic_NH()
538 for (CeedInt j = 0; j < 3; j++) { in ElasFSDiagnostic_NH()
539 for (CeedInt m = 0; m < 3; m++) diagnostic[5][i] += E2[j][m] * E2[m][j] / 4.; in ElasFSDiagnostic_NH()
544 …diagnostic[7][i] = (lambda * logJ * logJ / 2. - mu * logJ + mu * (E2[0][0] + E2[1][1] + E2[2][2]) … in ElasFSDiagnostic_NH()