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) { 19d4d45553Srezgarshakeri 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*6c88e6a2Srezgarshakeri error[i] = (u[i] - target[i])*q_data[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) { 30d4d45553Srezgarshakeri const CeedScalar *u = in[0], *target = in[1], *q_data = in[2]; 319b072555Sjeremylt CeedScalar *error = out[0]; 32*6c88e6a2Srezgarshakeri for (CeedInt i=0; i<Q; i++) { 33*6c88e6a2Srezgarshakeri error[i+0*Q] = (u[i+0*Q] - target[i+0*Q])*q_data[i]*q_data[i]; 34*6c88e6a2Srezgarshakeri error[i+1*Q] = (u[i+1*Q] - target[i+1*Q])*q_data[i]*q_data[i]; 35*6c88e6a2Srezgarshakeri error[i+2*Q] = (u[i+2*Q] - target[i+2*Q])*q_data[i]*q_data[i]; 36cb32e2e7SValeria Barra } 37cb32e2e7SValeria Barra return 0; 38cb32e2e7SValeria Barra } 39cb32e2e7SValeria Barra // ----------------------------------------------------------------------------- 40f6b55d2cSvaleriabarra 41f6b55d2cSvaleriabarra #endif // common_h 42