1 // Copyright (c) 2017-2026, Lawrence Livermore National Security, LLC and other CEED contributors.
2 // All Rights Reserved. See the top-level LICENSE and NOTICE files for details.
4 // SPDX-License-Identifier: BSD-2-Clause
15 /// libCEED Q-function for building quadrature data for a diffusion operator
24   // At every quadrature point, compute w/det(J).adj(J).adj(J)^T and store  in build_diff()
26   switch (build_data->dim + 10 * build_data->space_dim) {  in build_diff()
28       const CeedScalar(*J)[1][CEED_Q_VLA] = (const CeedScalar(*)[1][CEED_Q_VLA])in[0];  in build_diff()  local
30 …CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { q_data[0][i] = w[i] / J[0][0][i]; }  // End of Qu…  in build_diff()
33       const CeedScalar(*J)[2][CEED_Q_VLA] = (const CeedScalar(*)[2][CEED_Q_VLA])in[0];  in build_diff()  local
36         // J: 0 2   q_data: 0 2   adj(J):  J11 -J01  in build_diff()
37         //    1 3           2 1           -J10  J00  in build_diff()
38         const CeedScalar J00 = J[0][0][i];  in build_diff()
39         const CeedScalar J10 = J[0][1][i];  in build_diff()
40         const CeedScalar J01 = J[1][0][i];  in build_diff()
41         const CeedScalar J11 = J[1][1][i];  in build_diff()
42         const CeedScalar qw  = w[i] / (J00 * J11 - J10 * J01);  in build_diff()
46         q_data[2][i] = -qw * (J00 * J01 + J10 * J11);  in build_diff()
50       const CeedScalar(*J)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[0];  in build_diff()  local
56         for (CeedInt j = 0; j < 3; j++) {  in build_diff()  local
58             // Equivalent code with J as a VLA and no mod operations:  in build_diff()
59             // 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_diff()
60             A[k][j] =  in build_diff()
61 …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_diff()
65         // Compute quadrature weight / det(J)  in build_diff()
66 … 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_diff()
85 /// libCEED Q-function for applying a diff operator
93   switch (build_data->dim) {  in apply_diff()
114         // j = direction of vg  in apply_diff()
115 …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_diff()  local
134         // j = direction of vg  in apply_diff()
135 …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_diff()  local