xref: /libCEED/tests/t216-elemrestriction.c (revision 4fee36f0a30516a0b5ad51bf7eb3b32d83efd623)
19560d06aSjeremylt /// @file
29560d06aSjeremylt /// Test creation, copying, and destruction of an element restriction
39560d06aSjeremylt /// \test Test creation, copying, and destruction of an element restriction
49560d06aSjeremylt #include <ceed.h>
59560d06aSjeremylt 
69560d06aSjeremylt int main(int argc, char **argv) {
79560d06aSjeremylt   Ceed                ceed;
89560d06aSjeremylt   CeedInt             num_elem = 3, comp_stride = 1;
99560d06aSjeremylt   CeedInt             ind[2 * num_elem];
10*4fee36f0SJeremy L Thompson   CeedElemRestriction elem_restriction, elem_restriction_2;
119560d06aSjeremylt 
129560d06aSjeremylt   CeedInit(argv[1], &ceed);
139560d06aSjeremylt 
149560d06aSjeremylt   for (CeedInt i = 0; i < num_elem; i++) {
159560d06aSjeremylt     ind[2 * i + 0] = i;
169560d06aSjeremylt     ind[2 * i + 1] = i + 1;
179560d06aSjeremylt   }
18*4fee36f0SJeremy L Thompson   CeedElemRestrictionCreate(ceed, num_elem, 2, 1, comp_stride, num_elem + 1, CEED_MEM_HOST, CEED_USE_POINTER, ind, &elem_restriction);
19*4fee36f0SJeremy L Thompson   CeedElemRestrictionCreate(ceed, num_elem, 2, 1, comp_stride + 1, num_elem + 1, CEED_MEM_HOST, CEED_USE_POINTER, ind, &elem_restriction_2);
209560d06aSjeremylt 
21*4fee36f0SJeremy L Thompson   CeedElemRestrictionReferenceCopy(elem_restriction, &elem_restriction_2);  // This destroys the previous r_2
22*4fee36f0SJeremy L Thompson   CeedElemRestrictionDestroy(&elem_restriction);
239560d06aSjeremylt 
249560d06aSjeremylt   CeedInt comp_stride_2;
25*4fee36f0SJeremy L Thompson   CeedElemRestrictionGetCompStride(elem_restriction_2, &comp_stride_2);
262b730f8bSJeremy L Thompson   if (comp_stride_2 != comp_stride) printf("Error copying CeedElemRestriction reference\n");
279560d06aSjeremylt 
28*4fee36f0SJeremy L Thompson   CeedElemRestrictionDestroy(&elem_restriction_2);
299560d06aSjeremylt   CeedDestroy(&ceed);
309560d06aSjeremylt   return 0;
319560d06aSjeremylt }
32