Line data    Source code

       1              : // Copyright Amazon.com, Inc. or its affiliates. All Rights Reserved.
       2              : // SPDX-License-Identifier: Apache-2.0
       3              : 
       4              : #ifndef PALACE_LIBCEED_GEOM_21_QF_H
       5              : #define PALACE_LIBCEED_GEOM_21_QF_H
       6              : 
       7              : #include "utils_21_qf.h"
       8              : 
       9        19858 : CEED_QFUNCTION(f_build_geom_factor_21)(void *, CeedInt Q, const CeedScalar *const *in,
      10              :                                        CeedScalar *const *out)
      11              : {
      12        19858 :   const CeedScalar *attr = in[0], *qw = in[1], *J = in[2];
      13        19858 :   CeedScalar *qd_attr = out[0], *qd_wdetJ = out[0] + Q, *qd_adjJt = out[0] + 2 * Q;
      14              : 
      15       463690 :   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++)
      16              :   {
      17              :     CeedScalar J_loc[2], adjJt_loc[2];
      18       443832 :     MatUnpack21(J + i, Q, J_loc);
      19              :     const CeedScalar detJ = AdjJt21<true>(J_loc, adjJt_loc);
      20              : 
      21       443832 :     qd_attr[i] = attr[i];
      22       443832 :     qd_wdetJ[i] = qw[i] * detJ;
      23       443832 :     qd_adjJt[i + Q * 0] = adjJt_loc[0] / detJ;
      24       443832 :     qd_adjJt[i + Q * 1] = adjJt_loc[1] / detJ;
      25              :   }
      26        19858 :   return 0;
      27              : }
      28              : 
      29              : #endif  // PALACE_LIBCEED_GEOM_21_QF_H