Lines Matching full:0
42 CeedScalar sum = 0; in log1p_series_shifted()
71 return grad_u[0][0] * (grad_u[1][1] * grad_u[2][2] - grad_u[1][2] * grad_u[2][1]) + in computeJM1()
72 grad_u[0][1] * (grad_u[1][2] * grad_u[2][0] - grad_u[1][0] * grad_u[2][2]) + in computeJM1()
73 …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()
74 …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()
75 grad_u[0][2] * grad_u[2][0] - grad_u[1][2] * grad_u[2][1]; in computeJM1()
88 A[0][0] * A[2][2] - A[0][2] * A[2][0], /* *NOPAD* */ in computeMatinvSym()
89 A[0][0] * A[1][1] - A[0][1] * A[1][0], /* *NOPAD* */ in computeMatinvSym()
90 A[0][2] * A[1][0] - A[0][0] * A[1][2], /* *NOPAD* */ in computeMatinvSym()
91 A[0][1] * A[1][2] - A[0][2] * A[1][1], /* *NOPAD* */ in computeMatinvSym()
92 A[0][2] * A[2][1] - A[0][1] * A[2][2] /* *NOPAD* */ in computeMatinvSym()
94 for (CeedInt m = 0; m < 6; m++) Ainv[m] = B[m] / (detA); in computeMatinvSym()
106 const CeedInt indj[6] = {0, 1, 2, 1, 0, 0}, indk[6] = {0, 1, 2, 2, 2, 1}; in commonFSMR()
108 for (CeedInt m = 0; m < 6; m++) { in commonFSMR()
110 for (CeedInt n = 0; n < 3; n++) E2work[m] += grad_u[n][indj[m]] * grad_u[n][indk[m]]; in commonFSMR()
113 {E2work[0], E2work[5], E2work[4]}, in commonFSMR()
121 {1 + E2[0][0], E2[0][1], E2[0][2] }, in commonFSMR()
122 {E2[0][1], 1 + E2[1][1], E2[1][2] }, in commonFSMR()
123 {E2[0][2], E2[1][2], 1 + E2[2][2]} in commonFSMR()
126 Cwork[0] = C[0][0]; in commonFSMR()
130 Cwork[4] = C[0][2]; in commonFSMR()
131 Cwork[5] = C[0][1]; in commonFSMR()
134 const CeedScalar I_1 = C[0][0] + C[1][1] + C[2][2]; in commonFSMR()
144 // *1 for indices 0-2 for I_3 in commonFSMR()
147 for (CeedInt i = 0; i < 6; i++) { in commonFSMR()
160 …const CeedScalar(*ug)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0], (*q_data)[CEED_Q… in ElasFSResidual_MR()
163 CeedScalar(*dvdX)[3][CEED_Q_VLA] = (CeedScalar(*)[3][CEED_Q_VLA])out[0]; in ElasFSResidual_MR()
182 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ElasFSResidual_MR()
185 {ug[0][0][i], ug[1][0][i], ug[2][0][i]}, in ElasFSResidual_MR()
186 {ug[0][1][i], ug[1][1][i], ug[2][1][i]}, in ElasFSResidual_MR()
187 {ug[0][2][i], ug[1][2][i], ug[2][2][i]} in ElasFSResidual_MR()
190 const CeedScalar wdetJ = q_data[0][i]; in ElasFSResidual_MR()
200 for (CeedInt j = 0; j < 3; j++) { // Component in ElasFSResidual_MR()
201 for (CeedInt k = 0; k < 3; k++) { // Derivative in ElasFSResidual_MR()
202 grad_u[j][k][i] = 0; in ElasFSResidual_MR()
203 for (CeedInt m = 0; m < 3; m++) grad_u[j][k][i] += dXdx[m][k] * du[j][m]; in ElasFSResidual_MR()
210 {grad_u[0][0][i] + 1, grad_u[0][1][i], grad_u[0][2][i] }, in ElasFSResidual_MR()
211 {grad_u[1][0][i], grad_u[1][1][i] + 1, grad_u[1][2][i] }, in ElasFSResidual_MR()
212 {grad_u[2][0][i], grad_u[2][1][i], grad_u[2][2][i] + 1} in ElasFSResidual_MR()
216 {grad_u[0][0][i], grad_u[0][1][i], grad_u[0][2][i]}, in ElasFSResidual_MR()
217 {grad_u[1][0][i], grad_u[1][1][i], grad_u[1][2][i]}, in ElasFSResidual_MR()
218 {grad_u[2][0][i], grad_u[2][1][i], grad_u[2][2][i]} in ElasFSResidual_MR()
227 {Swork[0], Swork[5], Swork[4]}, in ElasFSResidual_MR()
234 for (CeedInt j = 0; j < 3; j++) { in ElasFSResidual_MR()
235 for (CeedInt k = 0; k < 3; k++) { in ElasFSResidual_MR()
236 P[j][k] = 0; in ElasFSResidual_MR()
237 for (CeedInt m = 0; m < 3; m++) P[j][k] += F[j][m] * S[m][k]; in ElasFSResidual_MR()
242 for (CeedInt j = 0; j < 3; j++) { // Component in ElasFSResidual_MR()
243 for (CeedInt k = 0; k < 3; k++) { // Derivative in ElasFSResidual_MR()
244 dvdX[k][j][i] = 0; in ElasFSResidual_MR()
245 for (CeedInt m = 0; m < 3; m++) dvdX[k][j][i] += dXdx[k][m] * P[j][m] * wdetJ; in ElasFSResidual_MR()
258 const CeedScalar(*deltaug)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0], in ElasFSJacobian_MR()
264 CeedScalar(*deltadvdX)[3][CEED_Q_VLA] = (CeedScalar(*)[3][CEED_Q_VLA])out[0]; in ElasFSJacobian_MR()
273 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ElasFSJacobian_MR()
276 {deltaug[0][0][i], deltaug[1][0][i], deltaug[2][0][i]}, in ElasFSJacobian_MR()
277 {deltaug[0][1][i], deltaug[1][1][i], deltaug[2][1][i]}, in ElasFSJacobian_MR()
278 {deltaug[0][2][i], deltaug[1][2][i], deltaug[2][2][i]} in ElasFSJacobian_MR()
281 const CeedScalar wdetJ = q_data[0][i]; in ElasFSJacobian_MR()
293 for (CeedInt j = 0; j < 3; j++) { // Component in ElasFSJacobian_MR()
294 for (CeedInt k = 0; k < 3; k++) { // Derivative in ElasFSJacobian_MR()
295 graddeltau[j][k] = 0; in ElasFSJacobian_MR()
296 for (CeedInt m = 0; m < 3; m++) graddeltau[j][k] += dXdx[m][k] * deltadu[j][m]; in ElasFSJacobian_MR()
303 {grad_u[0][0][i] + 1, grad_u[0][1][i], grad_u[0][2][i] }, in ElasFSJacobian_MR()
304 {grad_u[1][0][i], grad_u[1][1][i] + 1, grad_u[1][2][i] }, in ElasFSJacobian_MR()
305 {grad_u[2][0][i], grad_u[2][1][i], grad_u[2][2][i] + 1} in ElasFSJacobian_MR()
309 {grad_u[0][0][i], grad_u[0][1][i], grad_u[0][2][i]}, in ElasFSJacobian_MR()
310 {grad_u[1][0][i], grad_u[1][1][i], grad_u[1][2][i]}, in ElasFSJacobian_MR()
311 {grad_u[2][0][i], grad_u[2][1][i], grad_u[2][2][i]} in ElasFSJacobian_MR()
319 const CeedInt indj[6] = {0, 1, 2, 1, 0, 0}, indk[6] = {0, 1, 2, 2, 2, 1}; in ElasFSJacobian_MR()
321 for (CeedInt m = 0; m < 6; m++) { in ElasFSJacobian_MR()
322 dEwork[m] = 0; in ElasFSJacobian_MR()
323 …for (CeedInt n = 0; n < 3; n++) dEwork[m] += (graddeltau[n][indj[m]] * F[n][indk[m]] + F[n][indj[m… in ElasFSJacobian_MR()
326 {dEwork[0], dEwork[5], dEwork[4]}, in ElasFSJacobian_MR()
333 {Cwork[0], Cwork[5], Cwork[4]}, in ElasFSJacobian_MR()
338 {Cinvwork[0], Cinvwork[5], Cinvwork[4]}, in ElasFSJacobian_MR()
343 CeedScalar Cinv_contract_dE = 0; in ElasFSJacobian_MR()
344 for (CeedInt j = 0; j < 3; j++) { in ElasFSJacobian_MR()
345 for (CeedInt k = 0; k < 3; k++) Cinv_contract_dE += C_inv[j][k] * dE[j][k]; in ElasFSJacobian_MR()
349 CeedScalar C_contract_dE = 0; in ElasFSJacobian_MR()
350 for (CeedInt j = 0; j < 3; j++) { in ElasFSJacobian_MR()
351 for (CeedInt k = 0; k < 3; k++) C_contract_dE += C[j][k] * dE[j][k]; in ElasFSJacobian_MR()
356 for (CeedInt j = 0; j < 3; j++) { in ElasFSJacobian_MR()
357 for (CeedInt k = 0; k < 3; k++) { in ElasFSJacobian_MR()
358 dE_Cinv[j][k] = 0; in ElasFSJacobian_MR()
359 for (CeedInt m = 0; m < 3; m++) dE_Cinv[j][k] += dE[j][m] * C_inv[m][k]; in ElasFSJacobian_MR()
366 for (CeedInt j = 0; j < 3; j++) { in ElasFSJacobian_MR()
367 for (CeedInt k = 0; k < 3; k++) { in ElasFSJacobian_MR()
368 Cinv_dE_Cinv[j][k] = 0; in ElasFSJacobian_MR()
369 for (CeedInt m = 0; m < 3; m++) Cinv_dE_Cinv[j][k] += C_inv[j][m] * dE_Cinv[m][k]; in ElasFSJacobian_MR()
377 CeedScalar tr_dE = dE[0][0] + dE[1][1] + dE[2][2]; in ElasFSJacobian_MR()
379 for (CeedInt i = 0; i < 6; i++) { in ElasFSJacobian_MR()
385 {dSwork[0], dSwork[5], dSwork[4]}, in ElasFSJacobian_MR()
391 {Swork[0], Swork[5], Swork[4]}, in ElasFSJacobian_MR()
397 for (CeedInt j = 0; j < 3; j++) { in ElasFSJacobian_MR()
398 for (CeedInt k = 0; k < 3; k++) { in ElasFSJacobian_MR()
399 dP[j][k] = 0; in ElasFSJacobian_MR()
400 for (CeedInt m = 0; m < 3; m++) dP[j][k] += graddeltau[j][m] * S[m][k] + F[j][m] * dS[m][k]; in ElasFSJacobian_MR()
405 for (CeedInt j = 0; j < 3; j++) { // Component in ElasFSJacobian_MR()
406 for (CeedInt k = 0; k < 3; k++) { // Derivative in ElasFSJacobian_MR()
407 deltadvdX[k][j][i] = 0; in ElasFSJacobian_MR()
408 for (CeedInt m = 0; m < 3; m++) deltadvdX[k][j][i] += dXdx[k][m] * dP[j][m] * wdetJ; in ElasFSJacobian_MR()
421 …const CeedScalar(*ug)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0], (*q_data)[CEED_Q… in ElasFSEnergy_MR()
424 CeedScalar(*energy) = (CeedScalar(*))out[0]; in ElasFSEnergy_MR()
433 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ElasFSEnergy_MR()
436 {ug[0][0][i], ug[1][0][i], ug[2][0][i]}, in ElasFSEnergy_MR()
437 {ug[0][1][i], ug[1][1][i], ug[2][1][i]}, in ElasFSEnergy_MR()
438 {ug[0][2][i], ug[1][2][i], ug[2][2][i]} in ElasFSEnergy_MR()
441 const CeedScalar wdetJ = q_data[0][i]; in ElasFSEnergy_MR()
451 for (int j = 0; j < 3; j++) { // Component in ElasFSEnergy_MR()
452 for (int k = 0; k < 3; k++) { // Derivative in ElasFSEnergy_MR()
453 grad_u[j][k] = 0; in ElasFSEnergy_MR()
454 for (int m = 0; m < 3; m++) grad_u[j][k] += dXdx[m][k] * du[j][m]; in ElasFSEnergy_MR()
460 const CeedInt indj[6] = {0, 1, 2, 1, 0, 0}, indk[6] = {0, 1, 2, 2, 2, 1}; in ElasFSEnergy_MR()
462 for (CeedInt m = 0; m < 6; m++) { in ElasFSEnergy_MR()
464 for (CeedInt n = 0; n < 3; n++) E2work[m] += grad_u[n][indj[m]] * grad_u[n][indk[m]]; in ElasFSEnergy_MR()
467 {E2work[0], E2work[5], E2work[4]}, in ElasFSEnergy_MR()
475 {1 + E2[0][0], E2[0][1], E2[0][2] }, in ElasFSEnergy_MR()
476 {E2[0][1], 1 + E2[1][1], E2[1][2] }, in ElasFSEnergy_MR()
477 {E2[0][2], E2[1][2], 1 + E2[2][2]} in ElasFSEnergy_MR()
481 for (CeedInt j = 0; j < 3; j++) { in ElasFSEnergy_MR()
482 for (CeedInt k = 0; k < 3; k++) { in ElasFSEnergy_MR()
483 CC[j][k] = 0; in ElasFSEnergy_MR()
484 for (CeedInt m = 0; m < 3; m++) CC[j][k] += C[j][m] * C[m][k]; in ElasFSEnergy_MR()
492 const CeedScalar I_1 = C[0][0] + C[1][1] + C[2][2]; in ElasFSEnergy_MR()
494 const CeedScalar tr_CC = CC[0][0] + CC[1][1] + CC[2][2]; in ElasFSEnergy_MR()
512 …const CeedScalar(*u)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0], (*ug)[3][CEED_Q_VLA] = … in ElasFSDiagnostic_MR()
516 CeedScalar(*diagnostic)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in ElasFSDiagnostic_MR()
525 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in ElasFSDiagnostic_MR()
528 {ug[0][0][i], ug[1][0][i], ug[2][0][i]}, in ElasFSDiagnostic_MR()
529 {ug[0][1][i], ug[1][1][i], ug[2][1][i]}, in ElasFSDiagnostic_MR()
530 {ug[0][2][i], ug[1][2][i], ug[2][2][i]} in ElasFSDiagnostic_MR()
543 for (int j = 0; j < 3; j++) { // Component in ElasFSDiagnostic_MR()
544 for (int k = 0; k < 3; k++) { // Derivative in ElasFSDiagnostic_MR()
545 grad_u[j][k] = 0; in ElasFSDiagnostic_MR()
546 for (int m = 0; m < 3; m++) grad_u[j][k] += dXdx[m][k] * du[j][m]; in ElasFSDiagnostic_MR()
552 const CeedInt indj[6] = {0, 1, 2, 1, 0, 0}, indk[6] = {0, 1, 2, 2, 2, 1}; in ElasFSDiagnostic_MR()
554 for (CeedInt m = 0; m < 6; m++) { in ElasFSDiagnostic_MR()
556 for (CeedInt n = 0; n < 3; n++) E2work[m] += grad_u[n][indj[m]] * grad_u[n][indk[m]]; in ElasFSDiagnostic_MR()
559 {E2work[0], E2work[5], E2work[4]}, in ElasFSDiagnostic_MR()
565 diagnostic[0][i] = u[0][i]; in ElasFSDiagnostic_MR()
575 diagnostic[4][i] = (E2[0][0] + E2[1][1] + E2[2][2]) / 2.; in ElasFSDiagnostic_MR()
576 diagnostic[5][i] = 0.; in ElasFSDiagnostic_MR()
577 for (CeedInt j = 0; j < 3; j++) { in ElasFSDiagnostic_MR()
578 for (CeedInt m = 0; m < 3; m++) diagnostic[5][i] += E2[j][m] * E2[m][j] / 4.; in ElasFSDiagnostic_MR()
585 {1 + E2[0][0], E2[0][1], E2[0][2] }, in ElasFSDiagnostic_MR()
586 {E2[0][1], 1 + E2[1][1], E2[1][2] }, in ElasFSDiagnostic_MR()
587 {E2[0][2], E2[1][2], 1 + E2[2][2]} in ElasFSDiagnostic_MR()
591 for (CeedInt j = 0; j < 3; j++) { in ElasFSDiagnostic_MR()
592 for (CeedInt k = 0; k < 3; k++) { in ElasFSDiagnostic_MR()
593 CC[j][k] = 0; in ElasFSDiagnostic_MR()
594 for (CeedInt m = 0; m < 3; m++) CC[j][k] += C[j][m] * C[m][k]; in ElasFSDiagnostic_MR()
601 const CeedScalar I_1 = C[0][0] + C[1][1] + C[2][2]; in ElasFSDiagnostic_MR()
603 const CeedScalar tr_CC = CC[0][0] + CC[1][1] + CC[2][2]; in ElasFSDiagnostic_MR()