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_APPLY_31_QF_H
       5              : #define PALACE_LIBCEED_APPLY_31_QF_H
       6              : 
       7            0 : CEED_QFUNCTION(f_apply_31)(void *, CeedInt Q, const CeedScalar *const *in,
       8              :                            CeedScalar *const *out)
       9              : {
      10            0 :   const CeedScalar *__restrict__ qd1 = in[0], *__restrict__ qd2 = in[0] + 9 * Q,
      11            0 :                                  *__restrict__ u1 = in[1], *__restrict__ u2 = in[2];
      12            0 :   CeedScalar *__restrict__ v1 = out[0], *__restrict__ v2 = out[1];
      13              : 
      14            0 :   CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++)
      15              :   {
      16            0 :     const CeedScalar u10 = u1[i + Q * 0];
      17            0 :     const CeedScalar u11 = u1[i + Q * 1];
      18            0 :     const CeedScalar u12 = u1[i + Q * 2];
      19            0 :     v1[i + Q * 0] = qd1[i + Q * 0] * u10 + qd1[i + Q * 3] * u11 + qd1[i + Q * 6] * u12;
      20            0 :     v1[i + Q * 1] = qd1[i + Q * 1] * u10 + qd1[i + Q * 4] * u11 + qd1[i + Q * 7] * u12;
      21            0 :     v1[i + Q * 2] = qd1[i + Q * 2] * u10 + qd1[i + Q * 5] * u11 + qd1[i + Q * 8] * u12;
      22              : 
      23            0 :     v2[i] = qd2[i] * u2[i];
      24              :   }
      25            0 :   return 0;
      26              : }
      27              : 
      28              : #endif  // PALACE_LIBCEED_APPLY_31_QF_H