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 // ----------------------------------------------------------------------------- 17*2b730f8bSJeremy L Thompson CEED_QFUNCTION(Error)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) { 18d4d45553Srezgarshakeri const CeedScalar *u = in[0], *target = in[1], *q_data = in[2]; 199b072555Sjeremylt CeedScalar *error = out[0]; 20cb32e2e7SValeria Barra for (CeedInt i = 0; i < Q; i++) { 2138f32c05Srezgarshakeri error[i] = (u[i] - target[i]) * (u[i] - target[i]) * q_data[i]; 22cb32e2e7SValeria Barra } 23cb32e2e7SValeria Barra return 0; 24cb32e2e7SValeria Barra } 25cb32e2e7SValeria Barra 26cb32e2e7SValeria Barra // ----------------------------------------------------------------------------- 27*2b730f8bSJeremy L Thompson CEED_QFUNCTION(Error3)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) { 28d4d45553Srezgarshakeri const CeedScalar *u = in[0], *target = in[1], *q_data = in[2]; 299b072555Sjeremylt CeedScalar *error = out[0]; 306c88e6a2Srezgarshakeri for (CeedInt i = 0; i < Q; i++) { 3138f32c05Srezgarshakeri error[i + 0 * Q] = (u[i + 0 * Q] - target[i + 0 * Q]) * (u[i + 0 * Q] - target[i + 0 * Q]) * q_data[i]; 3238f32c05Srezgarshakeri error[i + 1 * Q] = (u[i + 1 * Q] - target[i + 1 * Q]) * (u[i + 1 * Q] - target[i + 1 * Q]) * q_data[i]; 3338f32c05Srezgarshakeri error[i + 2 * Q] = (u[i + 2 * Q] - target[i + 2 * Q]) * (u[i + 2 * Q] - target[i + 2 * Q]) * q_data[i]; 34cb32e2e7SValeria Barra } 35cb32e2e7SValeria Barra return 0; 36cb32e2e7SValeria Barra } 37cb32e2e7SValeria Barra // ----------------------------------------------------------------------------- 38f6b55d2cSvaleriabarra 39f6b55d2cSvaleriabarra #endif // common_h 40