13d8e8822SJeremy L Thompson // Copyright (c) 2017-2022, Lawrence Livermore National Security, LLC and other CEED contributors. 23d8e8822SJeremy L Thompson // All Rights Reserved. See the top-level LICENSE and NOTICE files for details. 3cb32e2e7SValeria Barra // 43d8e8822SJeremy L Thompson // SPDX-License-Identifier: BSD-2-Clause 5cb32e2e7SValeria Barra // 63d8e8822SJeremy 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 14c9c2c079SJeremy L Thompson #include <ceed.h> 15c9c2c079SJeremy L Thompson 16cb32e2e7SValeria Barra // ----------------------------------------------------------------------------- 17cb32e2e7SValeria Barra CEED_QFUNCTION(Error)(void *ctx, CeedInt Q, 18cb32e2e7SValeria Barra const CeedScalar *const *in, CeedScalar *const *out) { 19*d4d45553Srezgarshakeri const CeedScalar *u = in[0], *target = in[1], *q_data = in[2]; 209b072555Sjeremylt CeedScalar *error = out[0]; 21cb32e2e7SValeria Barra for (CeedInt i=0; i<Q; i++) { 22*d4d45553Srezgarshakeri error[i] = (u[i] - target[i])*q_data[i]; 23cb32e2e7SValeria Barra } 24cb32e2e7SValeria Barra return 0; 25cb32e2e7SValeria Barra } 26cb32e2e7SValeria Barra 27cb32e2e7SValeria Barra // ----------------------------------------------------------------------------- 28cb32e2e7SValeria Barra CEED_QFUNCTION(Error3)(void *ctx, CeedInt Q, 29cb32e2e7SValeria Barra const CeedScalar *const *in, CeedScalar *const *out) { 30*d4d45553Srezgarshakeri const CeedScalar *u = in[0], *target = in[1], *q_data = in[2]; 319b072555Sjeremylt CeedScalar *error = out[0]; 32cb32e2e7SValeria Barra for (CeedInt i=0; i<3*Q; i++) { 33*d4d45553Srezgarshakeri error[i] = (u[i] - target[i])*q_data[i]; 34cb32e2e7SValeria Barra } 35cb32e2e7SValeria Barra return 0; 36cb32e2e7SValeria Barra } 37cb32e2e7SValeria Barra // ----------------------------------------------------------------------------- 38f6b55d2cSvaleriabarra 39f6b55d2cSvaleriabarra #endif // common_h 40