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
16 // -----------------------------------------------------------------------------
37 // -----------------------------------------------------------------------------
40   const CeedScalar(*J)[3][CEED_Q_VLA] = (const CeedScalar(*)[3][CEED_Q_VLA])in[1];  in SetupDiffGeo()  local
50     for (CeedInt j = 0; j < dim; j++) {  in SetupDiffGeo()  local
53         // A[k][j] = J[k+1][j+1]*J[k+2][j+2] - J[k+1][j+2]*J[k+2][j+1]  in SetupDiffGeo()
54         A[k][j] = J[(k + 1) % dim][(j + 1) % dim][i] * J[(k + 2) % dim][(j + 2) % dim][i] -  in SetupDiffGeo()
55                   J[(k + 1) % dim][(j + 2) % dim][i] * J[(k + 2) % dim][(j + 1) % dim][i];  in SetupDiffGeo()
58     const CeedScalar detJ = J[0][0][i] * A[0][0] + J[0][1][i] * A[0][1] + J[0][2][i] * A[0][2];  in SetupDiffGeo()
73 // -----------------------------------------------------------------------------
75 // -----------------------------------------------------------------------------
95 // -----------------------------------------------------------------------------
99 //   ug      - Input vector gradient at quadrature points
100 //   q_data  - Geometric factors
103 //   vg     - Output vector (test functions) gradient at quadrature points
104 // -----------------------------------------------------------------------------
120     for (int j = 0; j < 3; j++) {  // j = direction of vg  in Diff()  local
121 …  vg[i + j * Q] = (du[0] * dXdxdXdx_T[0][j] + du[1] * dXdxdXdx_T[1][j] + du[2] * dXdxdXdx_T[2][j]);  in Diff()
126 // -----------------------------------------------------------------------------