Lines Matching full:1
20 // in[1] is quadrature weights, size (Q) in build_mass_diff()
21 const CeedScalar *w = in[1]; in build_mass_diff()
28 const CeedScalar(*J)[1][CEED_Q_VLA] = (const CeedScalar(*)[1][CEED_Q_VLA])in[0]; in build_mass_diff()
35 q_data[1][i] = w[i] / J[0][0][i]; in build_mass_diff()
43 // 1 3 2 1 -J10 J00 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()
47 const CeedScalar J11 = J[1][1][i]; in build_mass_diff()
54 q_data[1][i] = qw * (J01 * J01 + J11 * J11); in build_mass_diff()
69 // A[k][j] = J[j+1][k+1]*J[j+2][k+2] - J[j+1][k+2]*J[j+2][k+1] in build_mass_diff()
71 …J[(k + 1) % 3][(j + 1) % 3][i] * J[(k + 2) % 3][(j + 2) % 3][i] - J[(k + 2) % 3][(j + 1) % 3][i] *… 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()
83 // 1 6 5 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()
102 // in[1], out[1] solution values with shape [1, 1, Q] in apply_mass_diff()
103 // in[1], out[1] solution gradients with shape [dim, 1, Q] in apply_mass_diff()
108 case 1: { in apply_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()
117 vg[i] = q_data[1][i] * ug[i]; in apply_mass_diff()
122 const CeedScalar(*ug)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1]; in apply_mass_diff()
124 CeedScalar(*vg)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[1]; in apply_mass_diff()
133 // 1 3 in apply_mass_diff()
136 {q_data[1][i], q_data[3][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()
146 const CeedScalar(*ug)[CEED_Q_VLA] = (const CeedScalar(*)[CEED_Q_VLA])in[1]; in apply_mass_diff()
148 CeedScalar(*vg)[CEED_Q_VLA] = (CeedScalar(*)[CEED_Q_VLA])out[1]; in apply_mass_diff()
157 // 1 6 5 in apply_mass_diff()
161 {q_data[1][i], q_data[6][i], q_data[5][i]}, in apply_mass_diff()
167 …t j = 0; j < 3; j++) vg[j][i] = (ug[0][i] * dXdxdXdx_T[0][j] + ug[1][i] * dXdxdXdx_T[1][j] + ug[2]… in apply_mass_diff()