Lines Matching full:0
19 // in[0] is Jacobians with shape [dim, dim, Q] in build_mass_diff()
22 CeedScalar(*q_data)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in build_mass_diff()
28 const CeedScalar(*J)[1][CEED_Q_VLA] = (const CeedScalar(*)[1][CEED_Q_VLA])in[0]; in build_mass_diff()
30 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in build_mass_diff()
32 q_data[0][i] = w[i] * J[0][0][i]; in build_mass_diff()
35 q_data[1][i] = w[i] / J[0][0][i]; in build_mass_diff()
39 const CeedScalar(*J)[2][CEED_Q_VLA] = (const CeedScalar(*)[2][CEED_Q_VLA])in[0]; in build_mass_diff()
41 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in build_mass_diff()
42 // J: 0 2 q_data: 0 2 adj(J): J22 -J12 in build_mass_diff()
44 const CeedScalar J00 = J[0][0][i]; in build_mass_diff()
45 const CeedScalar J10 = J[0][1][i]; in build_mass_diff()
46 const CeedScalar J01 = J[1][0][i]; in build_mass_diff()
51 q_data[0][i] = w[i] * (J00 * J11 - J10 * J01); in build_mass_diff()
60 const CeedScalar(*J)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0]; in build_mass_diff()
62 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in build_mass_diff()
66 for (CeedInt j = 0; j < 3; j++) { in build_mass_diff()
67 for (CeedInt k = 0; k < 3; k++) { in build_mass_diff()
76 … const CeedScalar qw = w[i] / (J[0][0][i] * A[0][0] + J[0][1][i] * A[0][1] + J[0][2][i] * A[0][2]); in build_mass_diff()
79 q_data[0][i] = w[i] * (J[0][0][i] * A[0][0] + J[0][1][i] * A[0][1] + J[0][2][i] * A[0][2]); in build_mass_diff()
86 q_data[1][i] = qw * (A[0][0] * A[0][0] + A[0][1] * A[0][1] + A[0][2] * A[0][2]); in build_mass_diff()
87 q_data[2][i] = qw * (A[1][0] * A[1][0] + A[1][1] * A[1][1] + A[1][2] * A[1][2]); in build_mass_diff()
88 q_data[3][i] = qw * (A[2][0] * A[2][0] + A[2][1] * A[2][1] + A[2][2] * A[2][2]); in build_mass_diff()
89 q_data[4][i] = qw * (A[1][0] * A[2][0] + A[1][1] * A[2][1] + A[1][2] * A[2][2]); in build_mass_diff()
90 q_data[5][i] = qw * (A[0][0] * A[2][0] + A[0][1] * A[2][1] + A[0][2] * A[2][2]); in build_mass_diff()
91 q_data[6][i] = qw * (A[0][0] * A[1][0] + A[0][1] * A[1][1] + A[0][2] * A[1][2]); in build_mass_diff()
109 const CeedScalar *u = in[0], *ug = in[1]; in apply_mass_diff()
110 CeedScalar *v = out[0], *vg = out[1]; in apply_mass_diff()
112 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in apply_mass_diff()
114 v[i] = q_data[0][i] * u[i]; in apply_mass_diff()
121 const CeedScalar *u = in[0]; in apply_mass_diff()
123 CeedScalar *v = out[0]; in apply_mass_diff()
126 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in apply_mass_diff()
128 v[i] = q_data[0][i] * u[i]; in apply_mass_diff()
141 …for (int j = 0; j < 2; j++) vg[j][i] = (ug[0][i] * dXdxdXdx_T[0][j] + ug[1][i] * dXdxdXdx_T[1][j]); in apply_mass_diff()
145 const CeedScalar *u = in[0]; in apply_mass_diff()
147 CeedScalar *v = out[0]; in apply_mass_diff()
150 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in apply_mass_diff()
152 v[i] = q_data[0][i] * u[i]; in apply_mass_diff()
167 …for (int j = 0; j < 3; j++) vg[j][i] = (ug[0][i] * dXdxdXdx_T[0][j] + ug[1][i] * dXdxdXdx_T[1][j] … in apply_mass_diff()