Lines Matching full:0
14 // in[0] is Jacobians with shape [dim, dim, Q] in build_mass_diff()
17 CeedScalar(*q_data)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[0]; in build_mass_diff()
24 const CeedScalar(*J)[1][CEED_Q_VLA] = (const CeedScalar(*)[1][CEED_Q_VLA])in[0]; in build_mass_diff()
26 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in build_mass_diff()
28 q_data[0][i] = w[i] * J[0][0][i]; in build_mass_diff()
31 q_data[1][i] = w[i] / J[0][0][i]; in build_mass_diff()
35 const CeedScalar(*J)[2][CEED_Q_VLA] = (const CeedScalar(*)[2][CEED_Q_VLA])in[0]; in build_mass_diff()
37 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in build_mass_diff()
38 // J: 0 2 q_data: 0 2 adj(J): J22 -J12 in build_mass_diff()
40 const CeedScalar J00 = J[0][0][i]; in build_mass_diff()
41 const CeedScalar J10 = J[0][1][i]; in build_mass_diff()
42 const CeedScalar J01 = J[1][0][i]; in build_mass_diff()
47 q_data[0][i] = w[i] * (J00 * J11 - J10 * J01); in build_mass_diff()
56 const CeedScalar(*J)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0]; in build_mass_diff()
58 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in build_mass_diff()
61 for (CeedInt j = 0; j < 3; j++) { in build_mass_diff()
62 for (CeedInt k = 0; k < 3; k++) { in build_mass_diff()
69 … 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()
72 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()
79 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()
80 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()
81 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()
82 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()
83 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()
84 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()
94 // in[0], out[0] solution values with shape [1, 1, Q] in apply_mass_diff()
101 const CeedScalar *u = in[0], *ug = in[1]; in apply_mass_diff()
102 CeedScalar *v = out[0], *vg = out[1]; in apply_mass_diff()
104 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in apply_mass_diff()
106 v[i] = q_data[0][i] * u[i]; in apply_mass_diff()
113 const CeedScalar *u = in[0]; in apply_mass_diff()
115 CeedScalar *v = out[0]; in apply_mass_diff()
118 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in apply_mass_diff()
120 v[i] = q_data[0][i] * u[i]; in apply_mass_diff()
133 for (int j = 0; j < 2; j++) { in apply_mass_diff()
134 vg[j][i] = (ug[0][i] * dXdxdXdx_T[0][j] + ug[1][i] * dXdxdXdx_T[1][j]); in apply_mass_diff()
139 const CeedScalar *u = in[0]; in apply_mass_diff()
141 CeedScalar *v = out[0]; in apply_mass_diff()
144 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in apply_mass_diff()
146 v[i] = q_data[0][i] * u[i]; in apply_mass_diff()
161 for (int j = 0; j < 3; j++) { in apply_mass_diff()
162 …vg[j][i] = (ug[0][i] * dXdxdXdx_T[0][j] + ug[1][i] * dXdxdXdx_T[1][j] + ug[2][i] * dXdxdXdx_T[2][j… in apply_mass_diff()