20   CeedScalar nu;  // Poisson's ratio
21   CeedScalar E;   // Young's Modulus
28 CEED_QFUNCTION(ElasResidual_Linear)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *…  in ElasResidual_Linear()
30 …const CeedScalar(*ug)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0], (*q_data)[CEED_Q…  in ElasResidual_Linear()
33   CeedScalar(*dvdX)[3][CEED_Q_VLA] = (CeedScalar(*)[3][CEED_Q_VLA])out[0];  in ElasResidual_Linear()
39   const CeedScalar E       = context->E;  in ElasResidual_Linear()
40   const CeedScalar nu      = context->nu;  in ElasResidual_Linear()
45     const CeedScalar du[3][3] = {  in ElasResidual_Linear()
51     const CeedScalar wdetJ      = q_data[0][i];  in ElasResidual_Linear()
52     const CeedScalar dXdx[3][3] = {  in ElasResidual_Linear()
61     CeedScalar grad_u[3][3];  in ElasResidual_Linear()
72     const CeedScalar e[3][3] = {  in ElasResidual_Linear()
98     const CeedScalar ss      = E / ((1 + nu) * (1 - 2 * nu));  in ElasResidual_Linear()
99     const CeedScalar sigma00 = ss * ((1 - nu) * e[0][0] + nu * e[1][1] + nu * e[2][2]),  in ElasResidual_Linear()
103     const CeedScalar sigma[3][3] = {  in ElasResidual_Linear()
124 CEED_QFUNCTION(ElasJacobian_Linear)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *…  in ElasJacobian_Linear()
126   const CeedScalar(*deltaug)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0],  in ElasJacobian_Linear()
127         (*q_data)[CEED_Q_VLA]               = (const CeedScalar(*)[CEED_Q_VLA])in[1];  in ElasJacobian_Linear()
132   CeedScalar(*deltadvdX)[3][CEED_Q_VLA] = (CeedScalar(*)[3][CEED_Q_VLA])out[0];  in ElasJacobian_Linear()
136   const CeedScalar E       = context->E;  in ElasJacobian_Linear()
137   const CeedScalar nu      = context->nu;  in ElasJacobian_Linear()
142     const CeedScalar deltadu[3][3] = {  in ElasJacobian_Linear()
148     const CeedScalar wdetJ      = q_data[0][i];  in ElasJacobian_Linear()
149     const CeedScalar dXdx[3][3] = {  in ElasJacobian_Linear()
158     CeedScalar graddeltau[3][3];  in ElasJacobian_Linear()
168     const CeedScalar de[3][3] = {  in ElasJacobian_Linear()
193     const CeedScalar ss       = E / ((1 + nu) * (1 - 2 * nu));  in ElasJacobian_Linear()
194     const CeedScalar dsigma00 = ss * ((1 - nu) * de[0][0] + nu * de[1][1] + nu * de[2][2]),  in ElasJacobian_Linear()
198     const CeedScalar dsigma[3][3] = {  in ElasJacobian_Linear()
219 CEED_QFUNCTION(ElasEnergy_Linear)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *co…  in ElasEnergy_Linear()
221 …const CeedScalar(*ug)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0], (*q_data)[CEED_Q…  in ElasEnergy_Linear()
224   CeedScalar(*energy) = (CeedScalar(*))out[0];  in ElasEnergy_Linear()
228   const CeedScalar E       = context->E;  in ElasEnergy_Linear()
229   const CeedScalar nu      = context->nu;  in ElasEnergy_Linear()
232   const CeedScalar TwoMu  = E / (1 + nu);  in ElasEnergy_Linear()
233   const CeedScalar mu     = TwoMu / 2;  in ElasEnergy_Linear()
234   const CeedScalar Kbulk  = E / (3 * (1 - 2 * nu));  // Bulk Modulus  in ElasEnergy_Linear()
235   const CeedScalar lambda = (3 * Kbulk - TwoMu) / 3;  in ElasEnergy_Linear()
240     const CeedScalar du[3][3] = {  in ElasEnergy_Linear()
246     const CeedScalar wdetJ      = q_data[0][i];  in ElasEnergy_Linear()
247     const CeedScalar dXdx[3][3] = {  in ElasEnergy_Linear()
256     CeedScalar grad_u[3][3];  in ElasEnergy_Linear()
267     const CeedScalar e[3][3] = {  in ElasEnergy_Linear()
274     const CeedScalar strain_vol = e[0][0] + e[1][1] + e[2][2];  in ElasEnergy_Linear()
286 CEED_QFUNCTION(ElasDiagnostic_Linear)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar…  in ElasDiagnostic_Linear()
288 …const CeedScalar(*u)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0], (*ug)[3][CEED_Q_VLA] = …  in ElasDiagnostic_Linear()
289         (*q_data)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[2];  in ElasDiagnostic_Linear()
292   CeedScalar(*diagnostic)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0];  in ElasDiagnostic_Linear()
296   const CeedScalar E       = context->E;  in ElasDiagnostic_Linear()
297   const CeedScalar nu      = context->nu;  in ElasDiagnostic_Linear()
300   const CeedScalar TwoMu  = E / (1 + nu);  in ElasDiagnostic_Linear()
301   const CeedScalar mu     = TwoMu / 2;  in ElasDiagnostic_Linear()
302   const CeedScalar Kbulk  = E / (3 * (1 - 2 * nu));  // Bulk Modulus  in ElasDiagnostic_Linear()
303   const CeedScalar lambda = (3 * Kbulk - TwoMu) / 3;  in ElasDiagnostic_Linear()
308     const CeedScalar du[3][3] = {  in ElasDiagnostic_Linear()
314     const CeedScalar dXdx[3][3] = {  in ElasDiagnostic_Linear()
323     CeedScalar grad_u[3][3];  in ElasDiagnostic_Linear()
334     const CeedScalar e[3][3] = {  in ElasDiagnostic_Linear()
346     const CeedScalar strain_vol = e[0][0] + e[1][1] + e[2][2];  in ElasDiagnostic_Linear()