xref: /libCEED/tests/t216-elemrestriction.c (revision 2b730f8b5a9c809740a0b3b302db43a719c636b1)
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];
109560d06aSjeremylt   CeedElemRestriction r, r_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*2b730f8bSJeremy L Thompson   CeedElemRestrictionCreate(ceed, num_elem, 2, 1, comp_stride, num_elem + 1, CEED_MEM_HOST, CEED_USE_POINTER, ind, &r);
19*2b730f8bSJeremy L Thompson   CeedElemRestrictionCreate(ceed, num_elem, 2, 1, comp_stride + 1, num_elem + 1, CEED_MEM_HOST, CEED_USE_POINTER, ind, &r_2);
209560d06aSjeremylt 
219560d06aSjeremylt   CeedElemRestrictionReferenceCopy(r, &r_2);  // This destroys the previous r_2
229560d06aSjeremylt   CeedElemRestrictionDestroy(&r);
239560d06aSjeremylt 
249560d06aSjeremylt   CeedInt comp_stride_2;
259560d06aSjeremylt   CeedElemRestrictionGetCompStride(r_2, &comp_stride_2);
26*2b730f8bSJeremy L Thompson   if (comp_stride_2 != comp_stride) printf("Error copying CeedElemRestriction reference\n");
279560d06aSjeremylt 
289560d06aSjeremylt   CeedElemRestrictionDestroy(&r_2);
299560d06aSjeremylt   CeedDestroy(&ceed);
309560d06aSjeremylt   return 0;
319560d06aSjeremylt }
32