Lines Matching full:1
30 …CeedScalar(*)[3][CEED_Q_VLA])in[0], (*q_data)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1]; in ElasResidual_Linear()
35 // (*grad_u)[3][CEED_Q_VLA] = (CeedScalar(*)[3][CEED_Q_VLA])out[1]; in ElasResidual_Linear()
46 {ug[0][0][i], ug[1][0][i], ug[2][0][i]}, in ElasResidual_Linear()
47 {ug[0][1][i], ug[1][1][i], ug[2][1][i]}, in ElasResidual_Linear()
48 {ug[0][2][i], ug[1][2][i], ug[2][2][i]} in ElasResidual_Linear()
53 {q_data[1][i], q_data[2][i], q_data[3][i]}, in ElasResidual_Linear()
59 // dXdx = (dx/dX)^(-1) in ElasResidual_Linear()
70 // e = 1/2 (grad u + (grad u)^T) in ElasResidual_Linear()
73 …{(grad_u[0][0] + grad_u[0][0]) / 2., (grad_u[0][1] + grad_u[1][0]) / 2., (grad_u[0][2] + grad_u[2]… in ElasResidual_Linear()
74 …{(grad_u[1][0] + grad_u[0][1]) / 2., (grad_u[1][1] + grad_u[1][1]) / 2., (grad_u[1][2] + grad_u[2]… in ElasResidual_Linear()
75 …{(grad_u[2][0] + grad_u[0][2]) / 2., (grad_u[2][1] + grad_u[1][2]) / 2., (grad_u[2][2] + grad_u[2]… in ElasResidual_Linear()
90 // [1-nu nu nu ] in ElasResidual_Linear()
91 // [ nu 1-nu nu ] in ElasResidual_Linear()
92 // [ nu nu 1-nu ] in ElasResidual_Linear()
93 // S = E/((1+nu)*(1-2*nu)) [ (1-2*nu)/2 ] in ElasResidual_Linear()
94 // [ (1-2*nu)/2 ] in ElasResidual_Linear()
95 // [ (1-2*nu)/2 ] 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()
100 sigma11 = ss * (nu * e[0][0] + (1 - nu) * e[1][1] + nu * e[2][2]), in ElasResidual_Linear()
101 …sigma22 = ss * (nu * e[0][0] + nu * e[1][1] + (1 - nu) * e[2][2]), sigma12 = ss * (1 - 2 * nu) * e… in ElasResidual_Linear()
102 … sigma02 = ss * (1 - 2 * nu) * e[0][2] * 0.5, sigma01 = ss * (1 - 2 * nu) * e[0][1] * 0.5; in ElasResidual_Linear()
127 (*q_data)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1]; in ElasJacobian_Linear()
143 {deltaug[0][0][i], deltaug[1][0][i], deltaug[2][0][i]}, in ElasJacobian_Linear()
144 {deltaug[0][1][i], deltaug[1][1][i], deltaug[2][1][i]}, in ElasJacobian_Linear()
145 {deltaug[0][2][i], deltaug[1][2][i], deltaug[2][2][i]} in ElasJacobian_Linear()
150 {q_data[1][i], q_data[2][i], q_data[3][i]}, in ElasJacobian_Linear()
156 // dXdx = (dx/dX)^(-1) in ElasJacobian_Linear()
167 // e = 1/2 (grad u + (grad u)^T) in ElasJacobian_Linear()
169 …{(graddeltau[0][0] + graddeltau[0][0]) / 2., (graddeltau[0][1] + graddeltau[1][0]) / 2., (graddelt… in ElasJacobian_Linear()
170 …{(graddeltau[1][0] + graddeltau[0][1]) / 2., (graddeltau[1][1] + graddeltau[1][1]) / 2., (graddelt… in ElasJacobian_Linear()
171 …{(graddeltau[2][0] + graddeltau[0][2]) / 2., (graddeltau[2][1] + graddeltau[1][2]) / 2., (graddelt… in ElasJacobian_Linear()
185 // [1-nu nu nu ] in ElasJacobian_Linear()
186 // [ nu 1-nu nu ] in ElasJacobian_Linear()
187 // [ nu nu 1-nu ] in ElasJacobian_Linear()
188 // S = E/((1+nu)*(1-2*nu)) [ (1-2*nu)/2 ] in ElasJacobian_Linear()
189 // [ (1-2*nu)/2 ] in ElasJacobian_Linear()
190 // [ (1-2*nu)/2 ] 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()
195 dsigma11 = ss * (nu * de[0][0] + (1 - nu) * de[1][1] + nu * de[2][2]), in ElasJacobian_Linear()
196 …dsigma22 = ss * (nu * de[0][0] + nu * de[1][1] + (1 - nu) * de[2][2]), dsigma12 = ss * (1 - 2 * nu… in ElasJacobian_Linear()
197 … dsigma02 = ss * (1 - 2 * nu) * de[0][2] / 2, dsigma01 = ss * (1 - 2 * nu) * de[0][1] / 2; in ElasJacobian_Linear()
221 …CeedScalar(*)[3][CEED_Q_VLA])in[0], (*q_data)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1]; in ElasEnergy_Linear()
232 const CeedScalar TwoMu = E / (1 + nu); in ElasEnergy_Linear()
234 const CeedScalar Kbulk = E / (3 * (1 - 2 * nu)); // Bulk Modulus in ElasEnergy_Linear()
241 {ug[0][0][i], ug[1][0][i], ug[2][0][i]}, in ElasEnergy_Linear()
242 {ug[0][1][i], ug[1][1][i], ug[2][1][i]}, in ElasEnergy_Linear()
243 {ug[0][2][i], ug[1][2][i], ug[2][2][i]} in ElasEnergy_Linear()
248 {q_data[1][i], q_data[2][i], q_data[3][i]}, in ElasEnergy_Linear()
254 // dXdx = (dx/dX)^(-1) in ElasEnergy_Linear()
265 // e = 1/2 (grad u + (grad u)^T) in ElasEnergy_Linear()
268 …{(grad_u[0][0] + grad_u[0][0]) / 2., (grad_u[0][1] + grad_u[1][0]) / 2., (grad_u[0][2] + grad_u[2]… in ElasEnergy_Linear()
269 …{(grad_u[1][0] + grad_u[0][1]) / 2., (grad_u[1][1] + grad_u[1][1]) / 2., (grad_u[1][2] + grad_u[2]… in ElasEnergy_Linear()
270 …{(grad_u[2][0] + grad_u[0][2]) / 2., (grad_u[2][1] + grad_u[1][2]) / 2., (grad_u[2][2] + grad_u[2]… in ElasEnergy_Linear()
274 const CeedScalar strain_vol = e[0][0] + e[1][1] + e[2][2]; in ElasEnergy_Linear()
276 …ain_vol * strain_vol / 2. + strain_vol * mu + (e[0][1] * e[0][1] + e[0][2] * e[0][2] + e[1][2] * e… in ElasEnergy_Linear()
288 … CeedScalar(*)[CEED_Q_VLA])in[0], (*ug)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[1], in ElasDiagnostic_Linear()
300 const CeedScalar TwoMu = E / (1 + nu); in ElasDiagnostic_Linear()
302 const CeedScalar Kbulk = E / (3 * (1 - 2 * nu)); // Bulk Modulus in ElasDiagnostic_Linear()
309 {ug[0][0][i], ug[1][0][i], ug[2][0][i]}, in ElasDiagnostic_Linear()
310 {ug[0][1][i], ug[1][1][i], ug[2][1][i]}, in ElasDiagnostic_Linear()
311 {ug[0][2][i], ug[1][2][i], ug[2][2][i]} in ElasDiagnostic_Linear()
315 {q_data[1][i], q_data[2][i], q_data[3][i]}, in ElasDiagnostic_Linear()
321 // dXdx = (dx/dX)^(-1) in ElasDiagnostic_Linear()
332 // e = 1/2 (grad u + (grad u)^T) in ElasDiagnostic_Linear()
335 …{(grad_u[0][0] + grad_u[0][0]) / 2., (grad_u[0][1] + grad_u[1][0]) / 2., (grad_u[0][2] + grad_u[2]… in ElasDiagnostic_Linear()
336 …{(grad_u[1][0] + grad_u[0][1]) / 2., (grad_u[1][1] + grad_u[1][1]) / 2., (grad_u[1][2] + grad_u[2]… in ElasDiagnostic_Linear()
337 …{(grad_u[2][0] + grad_u[0][2]) / 2., (grad_u[2][1] + grad_u[1][2]) / 2., (grad_u[2][2] + grad_u[2]… in ElasDiagnostic_Linear()
342 diagnostic[1][i] = u[1][i]; in ElasDiagnostic_Linear()
346 const CeedScalar strain_vol = e[0][0] + e[1][1] + e[2][2]; in ElasDiagnostic_Linear()
355 diagnostic[6][i] = 1 + strain_vol; in ElasDiagnostic_Linear()
359 …ain_vol * strain_vol / 2. + strain_vol * mu + (e[0][1] * e[0][1] + e[0][2] * e[0][2] + e[1][2] * e… in ElasDiagnostic_Linear()