1*3d8e8822SJeremy L Thompson // Copyright (c) 2017-2022, Lawrence Livermore National Security, LLC and other CEED contributors. 2*3d8e8822SJeremy L Thompson // All Rights Reserved. See the top-level LICENSE and NOTICE files for details. 3cb32e2e7SValeria Barra // 4*3d8e8822SJeremy L Thompson // SPDX-License-Identifier: BSD-2-Clause 5cb32e2e7SValeria Barra // 6*3d8e8822SJeremy L Thompson // This file is part of CEED: http://github.com/ceed 7cb32e2e7SValeria Barra 8cb32e2e7SValeria Barra /// @file 9cb32e2e7SValeria Barra /// libCEED QFunctions for BP examples using PETSc 10cb32e2e7SValeria Barra 11f6b55d2cSvaleriabarra #ifndef common_h 12f6b55d2cSvaleriabarra #define common_h 13f6b55d2cSvaleriabarra 14cb32e2e7SValeria Barra // ----------------------------------------------------------------------------- 15cb32e2e7SValeria Barra CEED_QFUNCTION(Error)(void *ctx, CeedInt Q, 16cb32e2e7SValeria Barra const CeedScalar *const *in, CeedScalar *const *out) { 17cb32e2e7SValeria Barra const CeedScalar *u = in[0], *target = in[1]; 189b072555Sjeremylt CeedScalar *error = out[0]; 19cb32e2e7SValeria Barra for (CeedInt i=0; i<Q; i++) { 209b072555Sjeremylt error[i] = u[i] - target[i]; 21cb32e2e7SValeria Barra } 22cb32e2e7SValeria Barra return 0; 23cb32e2e7SValeria Barra } 24cb32e2e7SValeria Barra 25cb32e2e7SValeria Barra // ----------------------------------------------------------------------------- 26cb32e2e7SValeria Barra CEED_QFUNCTION(Error3)(void *ctx, CeedInt Q, 27cb32e2e7SValeria Barra const CeedScalar *const *in, CeedScalar *const *out) { 28cb32e2e7SValeria Barra const CeedScalar *u = in[0], *target = in[1]; 299b072555Sjeremylt CeedScalar *error = out[0]; 30cb32e2e7SValeria Barra for (CeedInt i=0; i<3*Q; i++) { 319b072555Sjeremylt error[i] = u[i] - target[i]; 32cb32e2e7SValeria Barra } 33cb32e2e7SValeria Barra return 0; 34cb32e2e7SValeria Barra } 35cb32e2e7SValeria Barra // ----------------------------------------------------------------------------- 36f6b55d2cSvaleriabarra 37f6b55d2cSvaleriabarra #endif // common_h 38