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_HCURL_21_QF_H
       5              : #define PALACE_LIBCEED_HCURL_21_QF_H
       6              : 
       7              : #include "../coeff/coeff_2_qf.h"
       8              : #include "utils_21_qf.h"
       9              : 
      10         3450 : CEED_QFUNCTION(f_apply_hcurl_21)(void *__restrict__ ctx, CeedInt Q,
      11              :                                  const CeedScalar *const *in, CeedScalar *const *out)
      12              : {
      13         3450 :   const CeedScalar *attr = in[0], *wdetJ = in[0] + Q, *adjJt = in[0] + 2 * Q, *u = in[1];
      14         3450 :   CeedScalar *v = out[0];
      15              : 
      16        72594 :   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++)
      17              :   {
      18        69144 :     const CeedScalar u_loc[1] = {u[i + Q * 0]};
      19              :     CeedScalar coeff[4], adjJt_loc[2], v_loc[1];
      20        69144 :     CoeffUnpack2((const CeedIntScalar *)ctx, (CeedInt)attr[i], coeff);
      21        69144 :     MatUnpack21(adjJt + i, Q, adjJt_loc);
      22              :     MultAtBCx21(adjJt_loc, coeff, adjJt_loc, u_loc, v_loc);
      23              : 
      24        69144 :     v[i + Q * 0] = wdetJ[i] * v_loc[0];
      25              :   }
      26         3450 :   return 0;
      27              : }
      28              : 
      29              : #endif  // PALACE_LIBCEED_HCURL_21_QF_H