Lines Matching refs:CeedScalar
24 CeedScalar mu_1;
25 CeedScalar mu_2;
26 CeedScalar lambda;
40 CEED_QFUNCTION_HELPER CeedScalar log1p_series_shifted(CeedScalar x) { in log1p_series_shifted()
41 const CeedScalar left = sqrt(2.) / 2 - 1, right = sqrt(2.) - 1; in log1p_series_shifted()
42 CeedScalar sum = 0; in log1p_series_shifted()
52 CeedScalar y = x / (2. + x); in log1p_series_shifted()
53 const CeedScalar y2 = y * y; in log1p_series_shifted()
70 CEED_QFUNCTION_HELPER CeedScalar computeJM1(const CeedScalar grad_u[3][3]) { in computeJM1()
84 CEED_QFUNCTION_HELPER int computeMatinvSym(const CeedScalar A[3][3], const CeedScalar detA, CeedSca… in computeMatinvSym()
86 CeedScalar B[6] = { in computeMatinvSym()
102 …UNCTION_HELPER int commonFSMR(const CeedScalar mu_1, const CeedScalar mu_2, const CeedScalar lambd… in commonFSMR()
103 … CeedScalar Swork[6], CeedScalar Cwork[6], CeedScalar Cinvwork[6], CeedScalar *logJ) { in commonFSMR()
107 CeedScalar E2work[6]; in commonFSMR()
112 CeedScalar E2[3][3] = { in commonFSMR()
120 const CeedScalar C[3][3] = { in commonFSMR()
134 const CeedScalar I_1 = C[0][0] + C[1][1] + C[2][2]; in commonFSMR()
136 const CeedScalar Jm1 = computeJM1(grad_u); in commonFSMR()
139 const CeedScalar detC = (Jm1 + 1.) * (Jm1 + 1.); in commonFSMR()
158 CEED_QFUNCTION(ElasFSResidual_MR)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *co… in ElasFSResidual_MR()
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()
165 CeedScalar(*grad_u)[3][CEED_Q_VLA] = (CeedScalar(*)[3][CEED_Q_VLA])out[1]; in ElasFSResidual_MR()
169 const CeedScalar mu_1 = context->mu_1; in ElasFSResidual_MR()
170 const CeedScalar mu_2 = context->mu_2; in ElasFSResidual_MR()
171 const CeedScalar lambda = context->lambda; in ElasFSResidual_MR()
184 const CeedScalar du[3][3] = { in ElasFSResidual_MR()
190 const CeedScalar wdetJ = q_data[0][i]; in ElasFSResidual_MR()
191 const CeedScalar dXdx[3][3] = { in ElasFSResidual_MR()
209 const CeedScalar F[3][3] = { in ElasFSResidual_MR()
215 const CeedScalar tempgradu[3][3] = { in ElasFSResidual_MR()
222 CeedScalar Swork[6], Cwork[6], Cinvwork[6], logJ; in ElasFSResidual_MR()
226 const CeedScalar S[3][3] = { in ElasFSResidual_MR()
233 CeedScalar P[3][3]; in ElasFSResidual_MR()
256 CEED_QFUNCTION(ElasFSJacobian_MR)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *co… in ElasFSJacobian_MR()
258 const CeedScalar(*deltaug)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0], in ElasFSJacobian_MR()
259 (*q_data)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1]; in ElasFSJacobian_MR()
261 const CeedScalar(*grad_u)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[2]; in ElasFSJacobian_MR()
264 CeedScalar(*deltadvdX)[3][CEED_Q_VLA] = (CeedScalar(*)[3][CEED_Q_VLA])out[0]; in ElasFSJacobian_MR()
268 const CeedScalar mu_1 = context->mu_1; in ElasFSJacobian_MR()
269 const CeedScalar mu_2 = context->mu_2; in ElasFSJacobian_MR()
270 const CeedScalar lambda = context->lambda; in ElasFSJacobian_MR()
275 const CeedScalar deltadu[3][3] = { in ElasFSJacobian_MR()
281 const CeedScalar wdetJ = q_data[0][i]; in ElasFSJacobian_MR()
282 const CeedScalar dXdx[3][3] = { in ElasFSJacobian_MR()
292 CeedScalar graddeltau[3][3]; in ElasFSJacobian_MR()
302 const CeedScalar F[3][3] = { in ElasFSJacobian_MR()
308 const CeedScalar tempgradu[3][3] = { in ElasFSJacobian_MR()
315 CeedScalar Swork[6], Cwork[6], Cinvwork[6], logJ; in ElasFSJacobian_MR()
320 CeedScalar dEwork[6]; in ElasFSJacobian_MR()
325 CeedScalar dE[3][3] = { in ElasFSJacobian_MR()
332 const CeedScalar C[3][3] = { in ElasFSJacobian_MR()
337 const CeedScalar C_inv[3][3] = { in ElasFSJacobian_MR()
343 CeedScalar Cinv_contract_dE = 0; in ElasFSJacobian_MR()
349 CeedScalar C_contract_dE = 0; in ElasFSJacobian_MR()
355 CeedScalar dE_Cinv[3][3]; in ElasFSJacobian_MR()
364 CeedScalar Cinv_dE_Cinv[3][3]; in ElasFSJacobian_MR()
377 CeedScalar tr_dE = dE[0][0] + dE[1][1] + dE[2][2]; in ElasFSJacobian_MR()
378 CeedScalar dSwork[6]; in ElasFSJacobian_MR()
384 CeedScalar dS[3][3] = { in ElasFSJacobian_MR()
390 const CeedScalar S[3][3] = { in ElasFSJacobian_MR()
396 CeedScalar dP[3][3]; in ElasFSJacobian_MR()
419 CEED_QFUNCTION(ElasFSEnergy_MR)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *cons… in ElasFSEnergy_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()
428 const CeedScalar mu_1 = context->mu_1; in ElasFSEnergy_MR()
429 const CeedScalar mu_2 = context->mu_2; in ElasFSEnergy_MR()
430 const CeedScalar lambda = context->lambda; in ElasFSEnergy_MR()
435 const CeedScalar du[3][3] = { in ElasFSEnergy_MR()
441 const CeedScalar wdetJ = q_data[0][i]; in ElasFSEnergy_MR()
442 const CeedScalar dXdx[3][3] = { in ElasFSEnergy_MR()
450 CeedScalar grad_u[3][3]; in ElasFSEnergy_MR()
461 CeedScalar E2work[6]; in ElasFSEnergy_MR()
466 CeedScalar E2[3][3] = { in ElasFSEnergy_MR()
474 const CeedScalar C[3][3] = { in ElasFSEnergy_MR()
480 CeedScalar CC[3][3]; in ElasFSEnergy_MR()
488 const CeedScalar Jm1 = computeJM1(grad_u); 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()
496 const CeedScalar I_2 = 0.5 * (I_1 * I_1 - tr_CC); in ElasFSEnergy_MR()
498 const CeedScalar logJ = log1p_series_shifted(Jm1); in ElasFSEnergy_MR()
510 CEED_QFUNCTION(ElasFSDiagnostic_MR)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *… in ElasFSDiagnostic_MR()
512 …const CeedScalar(*u)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0], (*ug)[3][CEED_Q_VLA] = … in ElasFSDiagnostic_MR()
513 (*q_data)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[2]; in ElasFSDiagnostic_MR()
516 CeedScalar(*diagnostic)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in ElasFSDiagnostic_MR()
520 const CeedScalar mu_1 = context->mu_1; in ElasFSDiagnostic_MR()
521 const CeedScalar mu_2 = context->mu_2; in ElasFSDiagnostic_MR()
522 const CeedScalar lambda = context->lambda; in ElasFSDiagnostic_MR()
527 const CeedScalar du[3][3] = { in ElasFSDiagnostic_MR()
533 const CeedScalar dXdx[3][3] = { in ElasFSDiagnostic_MR()
542 CeedScalar grad_u[3][3]; in ElasFSDiagnostic_MR()
553 CeedScalar E2work[6]; in ElasFSDiagnostic_MR()
558 CeedScalar E2[3][3] = { in ElasFSDiagnostic_MR()
570 const CeedScalar Jm1 = computeJM1(grad_u); in ElasFSDiagnostic_MR()
571 const CeedScalar logJ = log1p_series_shifted(Jm1); in ElasFSDiagnostic_MR()
584 const CeedScalar C[3][3] = { in ElasFSDiagnostic_MR()
590 CeedScalar CC[3][3]; 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()
605 const CeedScalar I_2 = 0.5 * (pow(I_1, 2) - tr_CC); in ElasFSDiagnostic_MR()