1*9ba83ac0SJeremy L Thompson // Copyright (c) 2017-2026, 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 11c0b5abf0SJeremy L Thompson #include <ceed/types.h> 12c9c2c079SJeremy L Thompson 13cb32e2e7SValeria Barra // ----------------------------------------------------------------------------- 142b730f8bSJeremy L Thompson CEED_QFUNCTION(Error)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) { 15d4d45553Srezgarshakeri const CeedScalar *u = in[0], *target = in[1], *q_data = in[2]; 169b072555Sjeremylt CeedScalar *error = out[0]; 17cb32e2e7SValeria Barra for (CeedInt i = 0; i < Q; i++) { 1838f32c05Srezgarshakeri error[i] = (u[i] - target[i]) * (u[i] - target[i]) * q_data[i]; 19cb32e2e7SValeria Barra } 20cb32e2e7SValeria Barra return 0; 21cb32e2e7SValeria Barra } 22cb32e2e7SValeria Barra 23cb32e2e7SValeria Barra // ----------------------------------------------------------------------------- 242b730f8bSJeremy L Thompson CEED_QFUNCTION(Error3)(void *ctx, CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) { 25d4d45553Srezgarshakeri const CeedScalar *u = in[0], *target = in[1], *q_data = in[2]; 269b072555Sjeremylt CeedScalar *error = out[0]; 276c88e6a2Srezgarshakeri for (CeedInt i = 0; i < Q; i++) { 2838f32c05Srezgarshakeri error[i + 0 * Q] = (u[i + 0 * Q] - target[i + 0 * Q]) * (u[i + 0 * Q] - target[i + 0 * Q]) * q_data[i]; 2938f32c05Srezgarshakeri error[i + 1 * Q] = (u[i + 1 * Q] - target[i + 1 * Q]) * (u[i + 1 * Q] - target[i + 1 * Q]) * q_data[i]; 3038f32c05Srezgarshakeri error[i + 2 * Q] = (u[i + 2 * Q] - target[i + 2 * Q]) * (u[i + 2 * Q] - target[i + 2 * Q]) * q_data[i]; 31cb32e2e7SValeria Barra } 32cb32e2e7SValeria Barra return 0; 33cb32e2e7SValeria Barra } 34cb32e2e7SValeria Barra // ----------------------------------------------------------------------------- 35