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
10 #include "ex-common.h"
12 /// libCEED Q-function for building quadrature data for a diffusion operator
22   switch (build_data->dim + 10 * build_data->space_dim) {  in build_diff()
32         // J: 0 2   q_data: 0 2   adj(J):  J11 -J01  in build_diff()
33         //    1 3           2 1           -J10  J00  in build_diff()
38         const CeedScalar qw  = w[i] / (J00 * J11 - J10 * J01);  in build_diff()
42         q_data[2][i] = -qw * (J00 * J01 + J10 * J11);  in build_diff()
50         CeedScalar A[3][3];  in build_diff()  local
54             // Equivalent code with J as a VLA and no mod operations:  in build_diff()
55             // 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()
56             A[k][j] =  in build_diff()
57 …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()
62 … 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()
69         q_data[0][i] = qw * (A[0][0] * A[0][0] + A[0][1] * A[0][1] + A[0][2] * A[0][2]);  in build_diff()
70         q_data[1][i] = qw * (A[1][0] * A[1][0] + A[1][1] * A[1][1] + A[1][2] * A[1][2]);  in build_diff()
71         q_data[2][i] = qw * (A[2][0] * A[2][0] + A[2][1] * A[2][1] + A[2][2] * A[2][2]);  in build_diff()
72         q_data[3][i] = qw * (A[1][0] * A[2][0] + A[1][1] * A[2][1] + A[1][2] * A[2][2]);  in build_diff()
73         q_data[4][i] = qw * (A[0][0] * A[2][0] + A[0][1] * A[2][1] + A[0][2] * A[2][2]);  in build_diff()
74         q_data[5][i] = qw * (A[0][0] * A[1][0] + A[0][1] * A[1][1] + A[0][2] * A[1][2]);  in build_diff()
81 /// libCEED Q-function for applying a diff operator
88   switch (build_data->dim) {  in apply_diff()