1*5aed82e4SJeremy L Thompson // Copyright (c) 2017-2024, 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. 35754ecacSJeremy L Thompson // 43d8e8822SJeremy L Thompson // SPDX-License-Identifier: BSD-2-Clause 55754ecacSJeremy L Thompson // 63d8e8822SJeremy L Thompson // This file is part of CEED: http://github.com/ceed 75754ecacSJeremy L Thompson 85754ecacSJeremy L Thompson /// @file 95754ecacSJeremy L Thompson /// Linear elasticity manufactured solution true solution for solid mechanics example using PETSc 105754ecacSJeremy L Thompson 115754ecacSJeremy L Thompson #ifndef MANUFACTURED_TRUE_H 125754ecacSJeremy L Thompson #define MANUFACTURED_TRUE_H 135754ecacSJeremy L Thompson 14c9c2c079SJeremy L Thompson #include <ceed.h> 155754ecacSJeremy L Thompson #include <math.h> 165754ecacSJeremy L Thompson 175754ecacSJeremy L Thompson // ----------------------------------------------------------------------------- 185754ecacSJeremy L Thompson // True solution for linear elasticity manufactured solution 195754ecacSJeremy L Thompson // ----------------------------------------------------------------------------- 202b730f8bSJeremy L Thompson CEED_QFUNCTION(MMSTrueSoln)(void *ctx, const CeedInt Q, const CeedScalar *const *in, CeedScalar *const *out) { 215754ecacSJeremy L Thompson // Inputs 225754ecacSJeremy L Thompson const CeedScalar *coords = in[0]; 235754ecacSJeremy L Thompson 245754ecacSJeremy L Thompson // Outputs 255754ecacSJeremy L Thompson CeedScalar *true_soln = out[0]; 265754ecacSJeremy L Thompson 275754ecacSJeremy L Thompson // Quadrature Point Loop 282b730f8bSJeremy L Thompson CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { 295754ecacSJeremy L Thompson // Setup 305754ecacSJeremy L Thompson CeedScalar x = coords[i + 0 * Q], y = coords[i + 1 * Q], z = coords[i + 2 * Q]; 315754ecacSJeremy L Thompson 325754ecacSJeremy L Thompson // True solution 335754ecacSJeremy L Thompson // -- Component 1 345754ecacSJeremy L Thompson true_soln[i + 0 * Q] = exp(2 * x) * sin(3 * y) * cos(4 * z) / 1e8; 355754ecacSJeremy L Thompson 365754ecacSJeremy L Thompson // -- Component 2 375754ecacSJeremy L Thompson true_soln[i + 1 * Q] = exp(3 * y) * sin(4 * z) * cos(2 * x) / 1e8; 385754ecacSJeremy L Thompson 395754ecacSJeremy L Thompson // -- Component 3 405754ecacSJeremy L Thompson true_soln[i + 2 * Q] = exp(4 * z) * sin(2 * x) * cos(3 * y) / 1e8; 415754ecacSJeremy L Thompson 425754ecacSJeremy L Thompson } // End of Quadrature Point Loop 435754ecacSJeremy L Thompson 445754ecacSJeremy L Thompson return 0; 455754ecacSJeremy L Thompson } 465754ecacSJeremy L Thompson // ----------------------------------------------------------------------------- 475754ecacSJeremy L Thompson 485754ecacSJeremy L Thompson #endif // End MANUFACTURED_TRUE_H 49