Lines Matching refs:CeedScalar
20 CeedScalar nu; // Poisson's ratio
21 CeedScalar E; // Young's Modulus
35 CEED_QFUNCTION_HELPER CeedScalar log1p_series_shifted(CeedScalar x) { in log1p_series_shifted()
36 const CeedScalar left = sqrt(2.) / 2 - 1, right = sqrt(2.) - 1; in log1p_series_shifted()
37 CeedScalar sum = 0; in log1p_series_shifted()
47 CeedScalar y = x / (2. + x); in log1p_series_shifted()
48 const CeedScalar y2 = y * y; in log1p_series_shifted()
65 CEED_QFUNCTION_HELPER CeedScalar computeJM1(const CeedScalar grad_u[3][3]) { in computeJM1()
79 CEED_QFUNCTION_HELPER int computeMatinvSym(const CeedScalar A[3][3], const CeedScalar detA, CeedSca… in computeMatinvSym()
81 CeedScalar B[6] = { in computeMatinvSym()
98 …QFUNCTION_HELPER int commonFS(const CeedScalar lambda, const CeedScalar mu, const CeedScalar grad_… in commonFS()
99 CeedScalar Cinvwork[6], CeedScalar *logJ) { in commonFS()
103 CeedScalar E2work[6]; in commonFS()
108 CeedScalar E2[3][3] = { in commonFS()
114 const CeedScalar Jm1 = computeJM1(grad_u); in commonFS()
118 const CeedScalar C[3][3] = { in commonFS()
125 const CeedScalar detC = (Jm1 + 1.) * (Jm1 + 1.); in commonFS()
128 const CeedScalar C_inv[3][3] = { in commonFS()
147 CEED_QFUNCTION(ElasFSResidual_NH)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *co… in ElasFSResidual_NH()
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()
154 CeedScalar(*grad_u)[3][CEED_Q_VLA] = (CeedScalar(*)[3][CEED_Q_VLA])out[1]; in ElasFSResidual_NH()
158 const CeedScalar E = context->E; in ElasFSResidual_NH()
159 const CeedScalar nu = context->nu; in ElasFSResidual_NH()
160 const CeedScalar TwoMu = E / (1 + nu); in ElasFSResidual_NH()
161 const CeedScalar mu = TwoMu / 2; in ElasFSResidual_NH()
162 const CeedScalar Kbulk = E / (3 * (1 - 2 * nu)); // Bulk Modulus in ElasFSResidual_NH()
163 const CeedScalar lambda = (3 * Kbulk - TwoMu) / 3; in ElasFSResidual_NH()
182 const CeedScalar du[3][3] = { in ElasFSResidual_NH()
188 const CeedScalar wdetJ = q_data[0][i]; in ElasFSResidual_NH()
189 const CeedScalar dXdx[3][3] = { in ElasFSResidual_NH()
207 const CeedScalar F[3][3] = { in ElasFSResidual_NH()
214 CeedScalar Swork[6], Cinvwork[6], logJ; in ElasFSResidual_NH()
215 const CeedScalar tempgradu[3][3] = { in ElasFSResidual_NH()
223 const CeedScalar S[3][3] = { in ElasFSResidual_NH()
230 CeedScalar P[3][3]; in ElasFSResidual_NH()
253 CEED_QFUNCTION(ElasFSJacobian_NH)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *co… in ElasFSJacobian_NH()
255 const CeedScalar(*deltaug)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0], in ElasFSJacobian_NH()
256 (*q_data)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1]; in ElasFSJacobian_NH()
258 const CeedScalar(*grad_u)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[2]; in ElasFSJacobian_NH()
261 CeedScalar(*deltadvdX)[3][CEED_Q_VLA] = (CeedScalar(*)[3][CEED_Q_VLA])out[0]; in ElasFSJacobian_NH()
265 const CeedScalar E = context->E; in ElasFSJacobian_NH()
266 const CeedScalar nu = context->nu; in ElasFSJacobian_NH()
269 const CeedScalar TwoMu = E / (1 + nu); in ElasFSJacobian_NH()
270 const CeedScalar mu = TwoMu / 2; in ElasFSJacobian_NH()
271 const CeedScalar Kbulk = E / (3 * (1 - 2 * nu)); // Bulk Modulus in ElasFSJacobian_NH()
272 const CeedScalar lambda = (3 * Kbulk - TwoMu) / 3; in ElasFSJacobian_NH()
277 const CeedScalar deltadu[3][3] = { in ElasFSJacobian_NH()
283 const CeedScalar wdetJ = q_data[0][i]; in ElasFSJacobian_NH()
284 const CeedScalar dXdx[3][3] = { in ElasFSJacobian_NH()
293 CeedScalar graddeltau[3][3]; in ElasFSJacobian_NH()
303 const CeedScalar F[3][3] = { in ElasFSJacobian_NH()
310 CeedScalar Swork[6], Cinvwork[6], logJ; in ElasFSJacobian_NH()
311 const CeedScalar tempgradu[3][3] = { in ElasFSJacobian_NH()
320 CeedScalar deltaEwork[6]; in ElasFSJacobian_NH()
325 CeedScalar deltaE[3][3] = { in ElasFSJacobian_NH()
333 const CeedScalar C_inv[3][3] = { in ElasFSJacobian_NH()
340 const CeedScalar S[3][3] = { in ElasFSJacobian_NH()
349 CeedScalar Cinv_contract_E = 0; in ElasFSJacobian_NH()
354 CeedScalar deltaECinv[3][3]; in ElasFSJacobian_NH()
362 CeedScalar deltaS[3][3]; in ElasFSJacobian_NH()
375 CeedScalar deltaP[3][3]; in ElasFSJacobian_NH()
398 CEED_QFUNCTION(ElasFSEnergy_NH)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *cons… in ElasFSEnergy_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()
407 const CeedScalar E = context->E; in ElasFSEnergy_NH()
408 const CeedScalar nu = context->nu; in ElasFSEnergy_NH()
409 const CeedScalar TwoMu = E / (1 + nu); in ElasFSEnergy_NH()
410 const CeedScalar mu = TwoMu / 2; in ElasFSEnergy_NH()
411 const CeedScalar Kbulk = E / (3 * (1 - 2 * nu)); // Bulk Modulus in ElasFSEnergy_NH()
412 const CeedScalar lambda = (3 * Kbulk - TwoMu) / 3; in ElasFSEnergy_NH()
417 const CeedScalar du[3][3] = { in ElasFSEnergy_NH()
423 const CeedScalar wdetJ = q_data[0][i]; in ElasFSEnergy_NH()
424 const CeedScalar dXdx[3][3] = { in ElasFSEnergy_NH()
433 CeedScalar grad_u[3][3]; in ElasFSEnergy_NH()
444 CeedScalar E2work[6]; in ElasFSEnergy_NH()
449 CeedScalar E2[3][3] = { in ElasFSEnergy_NH()
454 const CeedScalar Jm1 = computeJM1(grad_u); in ElasFSEnergy_NH()
455 const CeedScalar logJ = log1p_series_shifted(Jm1); in ElasFSEnergy_NH()
468 CEED_QFUNCTION(ElasFSDiagnostic_NH)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *… in ElasFSDiagnostic_NH()
470 …const CeedScalar(*u)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0], (*ug)[3][CEED_Q_VLA] = … in ElasFSDiagnostic_NH()
471 (*q_data)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[2]; in ElasFSDiagnostic_NH()
474 CeedScalar(*diagnostic)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in ElasFSDiagnostic_NH()
478 const CeedScalar E = context->E; in ElasFSDiagnostic_NH()
479 const CeedScalar nu = context->nu; in ElasFSDiagnostic_NH()
480 const CeedScalar TwoMu = E / (1 + nu); in ElasFSDiagnostic_NH()
481 const CeedScalar mu = TwoMu / 2; in ElasFSDiagnostic_NH()
482 const CeedScalar Kbulk = E / (3 * (1 - 2 * nu)); // Bulk Modulus in ElasFSDiagnostic_NH()
483 const CeedScalar lambda = (3 * Kbulk - TwoMu) / 3; in ElasFSDiagnostic_NH()
488 const CeedScalar du[3][3] = { in ElasFSDiagnostic_NH()
494 const CeedScalar dXdx[3][3] = { in ElasFSDiagnostic_NH()
503 CeedScalar grad_u[3][3]; in ElasFSDiagnostic_NH()
514 CeedScalar E2work[6]; in ElasFSDiagnostic_NH()
519 CeedScalar E2[3][3] = { in ElasFSDiagnostic_NH()
531 const CeedScalar Jm1 = computeJM1(grad_u); in ElasFSDiagnostic_NH()
532 const CeedScalar logJ = log1p_series_shifted(Jm1); in ElasFSDiagnostic_NH()