qfunctions/bps/bp3.h

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 // -----------------------------------------------------------------------------
32 // Stored: w B^T B detJ = w A^T A / detJ
37 // -----------------------------------------------------------------------------
49     CeedScalar A[3][3];  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()
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()
62     qd[i + Q * 1]       = qw * (A[0][0] * A[0][0] + A[0][1] * A[0][1] + A[0][2] * A[0][2]);  in SetupDiffGeo()
63     qd[i + Q * 2]       = qw * (A[0][0] * A[1][0] + A[0][1] * A[1][1] + A[0][2] * A[1][2]);  in SetupDiffGeo()
64     qd[i + Q * 3]       = qw * (A[0][0] * A[2][0] + A[0][1] * A[2][1] + A[0][2] * A[2][2]);  in SetupDiffGeo()
65     qd[i + Q * 4]       = qw * (A[1][0] * A[1][0] + A[1][1] * A[1][1] + A[1][2] * A[1][2]);  in SetupDiffGeo()
66     qd[i + Q * 5]       = qw * (A[1][0] * A[2][0] + A[1][1] * A[2][1] + A[1][2] * A[2][2]);  in SetupDiffGeo()
67     qd[i + Q * 6]       = qw * (A[2][0] * A[2][0] + A[2][1] * A[2][1] + A[2][2] * A[2][2]);  in SetupDiffGeo()
73 // -----------------------------------------------------------------------------
75 // -----------------------------------------------------------------------------
95 // -----------------------------------------------------------------------------
96 // This QFunction applies the diffusion operator for a scalar field.
99 //   ug      - Input vector gradient at quadrature points
100 //   q_data  - Geometric factors
103 //   vg     - Output vector (test functions) gradient at quadrature points
104 // -----------------------------------------------------------------------------
126 // -----------------------------------------------------------------------------