xref: /libCEED/backends/blocked/ceed-blocked.c (revision ba4bcd917872bbba2ca5b7fccd7fd5f16d3aa2e4)
1 // Copyright (c) 2017-2018, Lawrence Livermore National Security, LLC.
2 // Produced at the Lawrence Livermore National Laboratory. LLNL-CODE-734707.
3 // All Rights reserved. See files LICENSE and NOTICE for details.
4 //
5 // This file is part of CEED, a collection of benchmarks, miniapps, software
6 // libraries and APIs for efficient high-order finite element and spectral
7 // element discretizations for exascale applications. For more information and
8 // source code availability see http://github.com/ceed.
9 //
10 // The CEED research is supported by the Exascale Computing Project 17-SC-20-SC,
11 // a collaborative effort of two U.S. Department of Energy organizations (Office
12 // of Science and the National Nuclear Security Administration) responsible for
13 // the planning and preparation of a capable exascale ecosystem, including
14 // software, applications, hardware, advanced system engineering and early
15 // testbed platforms, in support of the nation's exascale computing imperative.
16 
17 #include <ceed-impl.h>
18 #include <string.h>
19 #include "ceed-blocked.h"
20 
21 static int CeedDestroy_Blocked(Ceed ceed) {
22   int ierr;
23   Ceed_Blocked *impl = ceed->data;
24   Ceed ceedref = impl->ceedref;
25   ierr = CeedFree(&ceedref); CeedChk(ierr);
26   ierr = CeedFree(&impl); CeedChk(ierr);
27   return 0;
28 }
29 
30 static int CeedInit_Blocked(const char *resource, Ceed ceed) {
31   if (strcmp(resource, "/cpu/self")
32       && strcmp(resource, "/cpu/self/blocked"))
33     return CeedError(ceed, 1, "Blocked backend cannot use resource: %s", resource);
34 
35   int ierr;
36   Ceed_Blocked *impl;
37   Ceed ceedref;
38 
39   // Create refrence CEED that implementation will be dispatched
40   //   through unless overridden
41   ierr = CeedCalloc(1, &impl); CeedChk(ierr);
42   CeedInit("/cpu/self/ref", &ceedref);
43   ceed->data = impl;
44   impl->ceedref = ceedref;
45 
46   ceed->VecCreate = CeedVectorCreate_Blocked;
47   ceed->BasisCreateTensorH1 = CeedBasisCreateTensorH1_Blocked;
48   ceed->BasisCreateH1 = CeedBasisCreateH1_Blocked;
49   ceed->ElemRestrictionCreate = CeedElemRestrictionCreate_Blocked;
50   ceed->ElemRestrictionCreateBlocked = CeedElemRestrictionCreate_Blocked;
51   ceed->QFunctionCreate = CeedQFunctionCreate_Blocked;
52   ceed->OperatorCreate = CeedOperatorCreate_Blocked;
53   ceed->Destroy = CeedDestroy_Blocked;
54 
55   return 0;
56 }
57 
58 __attribute__((constructor))
59 static void Register(void) {
60   CeedRegister("/cpu/self/blocked", CeedInit_Blocked, 10);
61 }
62