xref: /libCEED/tests/t216-elemrestriction.c (revision 9560d06a92c065fb7d600a8c20ade8d9a4cda324)
1*9560d06aSjeremylt /// @file
2*9560d06aSjeremylt /// Test creation, copying, and destruction of an element restriction
3*9560d06aSjeremylt /// \test Test creation, copying, and destruction of an element restriction
4*9560d06aSjeremylt #include <ceed.h>
5*9560d06aSjeremylt 
6*9560d06aSjeremylt int main(int argc, char **argv) {
7*9560d06aSjeremylt   Ceed ceed;
8*9560d06aSjeremylt   CeedInt num_elem = 3, comp_stride = 1;
9*9560d06aSjeremylt   CeedInt ind[2*num_elem];
10*9560d06aSjeremylt   CeedElemRestriction r, r_2;
11*9560d06aSjeremylt 
12*9560d06aSjeremylt   CeedInit(argv[1], &ceed);
13*9560d06aSjeremylt 
14*9560d06aSjeremylt   for (CeedInt i=0; i<num_elem; i++) {
15*9560d06aSjeremylt     ind[2*i+0] = i;
16*9560d06aSjeremylt     ind[2*i+1] = i+1;
17*9560d06aSjeremylt   }
18*9560d06aSjeremylt   CeedElemRestrictionCreate(ceed, num_elem, 2, 1, comp_stride, num_elem+1,
19*9560d06aSjeremylt                             CEED_MEM_HOST, CEED_USE_POINTER, ind, &r);
20*9560d06aSjeremylt   CeedElemRestrictionCreate(ceed, num_elem, 2, 1, comp_stride+1, num_elem+1,
21*9560d06aSjeremylt                             CEED_MEM_HOST, CEED_USE_POINTER, ind, &r_2);
22*9560d06aSjeremylt 
23*9560d06aSjeremylt   CeedElemRestrictionReferenceCopy(r, &r_2); // This destroys the previous r_2
24*9560d06aSjeremylt   CeedElemRestrictionDestroy(&r);
25*9560d06aSjeremylt 
26*9560d06aSjeremylt   CeedInt comp_stride_2;
27*9560d06aSjeremylt   CeedElemRestrictionGetCompStride(r_2, &comp_stride_2);
28*9560d06aSjeremylt   if (comp_stride_2 != comp_stride)
29*9560d06aSjeremylt     // LCOV_EXCL_START
30*9560d06aSjeremylt     printf("Error copying CeedElemRestriction reference.");
31*9560d06aSjeremylt   // LCOV_EXCL_STOP
32*9560d06aSjeremylt 
33*9560d06aSjeremylt   CeedElemRestrictionDestroy(&r_2);
34*9560d06aSjeremylt   CeedDestroy(&ceed);
35*9560d06aSjeremylt   return 0;
36*9560d06aSjeremylt }
37