xref: /libCEED/rust/libceed-sys/c-src/gallery/mass/ceed-massapply.c (revision 3d8e882215d238700cdceb37404f76ca7fa24eaa)
1*3d8e8822SJeremy L Thompson // Copyright (c) 2017-2022, Lawrence Livermore National Security, LLC and other CEED contributors.
2*3d8e8822SJeremy L Thompson // All Rights Reserved. See the top-level LICENSE and NOTICE files for details.
36a6224a1SJeremy L Thompson //
4*3d8e8822SJeremy L Thompson // SPDX-License-Identifier: BSD-2-Clause
56a6224a1SJeremy L Thompson //
6*3d8e8822SJeremy L Thompson // This file is part of CEED:  http://github.com/ceed
76a6224a1SJeremy L Thompson 
86a6224a1SJeremy L Thompson #include <ceed/ceed.h>
96a6224a1SJeremy L Thompson #include <ceed/backend.h>
106a6224a1SJeremy L Thompson #include <string.h>
116a6224a1SJeremy L Thompson #include "ceed-massapply.h"
126a6224a1SJeremy L Thompson 
136a6224a1SJeremy L Thompson /**
146a6224a1SJeremy L Thompson   @brief Set fields for Ceed QFunction for applying the mass matrix
156a6224a1SJeremy L Thompson **/
166a6224a1SJeremy L Thompson static int CeedQFunctionInit_MassApply(Ceed ceed, const char *requested,
176a6224a1SJeremy L Thompson                                        CeedQFunction qf) {
186a6224a1SJeremy L Thompson   int ierr;
196a6224a1SJeremy L Thompson 
206a6224a1SJeremy L Thompson   // Check QFunction name
216a6224a1SJeremy L Thompson   const char *name = "MassApply";
226a6224a1SJeremy L Thompson   if (strcmp(name, requested))
236a6224a1SJeremy L Thompson     // LCOV_EXCL_START
246a6224a1SJeremy L Thompson     return CeedError(ceed, CEED_ERROR_UNSUPPORTED,
256a6224a1SJeremy L Thompson                      "QFunction '%s' does not match requested name: %s",
266a6224a1SJeremy L Thompson                      name, requested);
276a6224a1SJeremy L Thompson   // LCOV_EXCL_STOP
286a6224a1SJeremy L Thompson 
296a6224a1SJeremy L Thompson   // Add QFunction fields
306a6224a1SJeremy L Thompson   ierr = CeedQFunctionAddInput(qf, "u", 1, CEED_EVAL_INTERP); CeedChk(ierr);
316a6224a1SJeremy L Thompson   ierr = CeedQFunctionAddInput(qf, "qdata", 1, CEED_EVAL_NONE); CeedChk(ierr);
326a6224a1SJeremy L Thompson   ierr = CeedQFunctionAddOutput(qf, "v", 1, CEED_EVAL_INTERP); CeedChk(ierr);
336a6224a1SJeremy L Thompson 
346a6224a1SJeremy L Thompson   return CEED_ERROR_SUCCESS;
356a6224a1SJeremy L Thompson }
366a6224a1SJeremy L Thompson 
376a6224a1SJeremy L Thompson /**
386a6224a1SJeremy L Thompson   @brief Register Ceed QFunction for applying the mass matrix
396a6224a1SJeremy L Thompson **/
406a6224a1SJeremy L Thompson CEED_INTERN int CeedQFunctionRegister_MassApply(void) {
416a6224a1SJeremy L Thompson   return CeedQFunctionRegister("MassApply", MassApply_loc, 1, MassApply,
426a6224a1SJeremy L Thompson                                CeedQFunctionInit_MassApply);
436a6224a1SJeremy L Thompson }
44