Lines Matching full:0
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()
43 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { 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()
51 const CeedScalar wdetJ = q_data[0][i]; in ElasResidual_Linear()
62 for (CeedInt j = 0; j < 3; j++) { // Component in ElasResidual_Linear()
63 for (CeedInt k = 0; k < 3; k++) { // Derivative in ElasResidual_Linear()
64 grad_u[j][k] = 0; in ElasResidual_Linear()
65 for (CeedInt m = 0; m < 3; m++) grad_u[j][k] += dXdx[m][k] * du[j][m]; 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()
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()
110 for (CeedInt j = 0; j < 3; j++) { // Component in ElasResidual_Linear()
111 for (CeedInt k = 0; k < 3; k++) { // Derivative in ElasResidual_Linear()
112 dvdX[k][j][i] = 0; in ElasResidual_Linear()
113 for (CeedInt m = 0; m < 3; m++) dvdX[k][j][i] += dXdx[k][m] * sigma[j][m] * wdetJ; in ElasResidual_Linear()
126 const CeedScalar(*deltaug)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0], in ElasJacobian_Linear()
132 CeedScalar(*deltadvdX)[3][CEED_Q_VLA] = (CeedScalar(*)[3][CEED_Q_VLA])out[0]; in ElasJacobian_Linear()
140 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { 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()
148 const CeedScalar wdetJ = q_data[0][i]; in ElasJacobian_Linear()
159 for (CeedInt j = 0; j < 3; j++) { // Component in ElasJacobian_Linear()
160 for (CeedInt k = 0; k < 3; k++) { // Derivative in ElasJacobian_Linear()
161 graddeltau[j][k] = 0; in ElasJacobian_Linear()
162 for (CeedInt m = 0; m < 3; m++) graddeltau[j][k] += dXdx[m][k] * deltadu[j][m]; 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()
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()
205 for (CeedInt j = 0; j < 3; j++) { // Component in ElasJacobian_Linear()
206 for (CeedInt k = 0; k < 3; k++) { // Derivative in ElasJacobian_Linear()
207 deltadvdX[k][j][i] = 0; in ElasJacobian_Linear()
208 for (CeedInt m = 0; m < 3; m++) deltadvdX[k][j][i] += dXdx[k][m] * dsigma[j][m] * wdetJ; in ElasJacobian_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()
238 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { 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()
246 const CeedScalar wdetJ = q_data[0][i]; in ElasEnergy_Linear()
257 for (CeedInt j = 0; j < 3; j++) { // Component in ElasEnergy_Linear()
258 for (CeedInt k = 0; k < 3; k++) { // Derivative in ElasEnergy_Linear()
259 grad_u[j][k] = 0; in ElasEnergy_Linear()
260 for (CeedInt m = 0; m < 3; m++) grad_u[j][k] += dXdx[m][k] * du[j][m]; 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 …lambda * strain_vol * strain_vol / 2. + strain_vol * mu + (e[0][1] * e[0][1] + e[0][2] * e[0][2] +… in ElasEnergy_Linear()
288 …const CeedScalar(*u)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[0], (*ug)[3][CEED_Q_VLA] = … in ElasDiagnostic_Linear()
292 CeedScalar(*diagnostic)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in ElasDiagnostic_Linear()
306 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { 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()
324 for (CeedInt j = 0; j < 3; j++) { // Component in ElasDiagnostic_Linear()
325 for (CeedInt k = 0; k < 3; k++) { // Derivative in ElasDiagnostic_Linear()
326 grad_u[j][k] = 0; in ElasDiagnostic_Linear()
327 for (CeedInt m = 0; m < 3; m++) grad_u[j][k] += dXdx[m][k] * du[j][m]; 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()
341 diagnostic[0][i] = u[0][i]; in ElasDiagnostic_Linear()
346 const CeedScalar strain_vol = e[0][0] + e[1][1] + e[2][2]; in ElasDiagnostic_Linear()
351 diagnostic[5][i] = 0.; in ElasDiagnostic_Linear()
352 for (CeedInt j = 0; j < 3; j++) { in ElasDiagnostic_Linear()
353 for (CeedInt m = 0; m < 3; m++) diagnostic[5][i] += e[j][m] * e[m][j]; in ElasDiagnostic_Linear()
359 …lambda * strain_vol * strain_vol / 2. + strain_vol * mu + (e[0][1] * e[0][1] + e[0][2] * e[0][2] +… in ElasDiagnostic_Linear()