xref: /libCEED/tests/t307-basis.c (revision 4fee36f0a30516a0b5ad51bf7eb3b32d83efd623)
19560d06aSjeremylt /// @file
29560d06aSjeremylt /// Test creation, copying, and distruction of a H1Lagrange basis
39560d06aSjeremylt /// \test Test creation, copying, and distruction of a H1Lagrange basis
49560d06aSjeremylt #include <ceed.h>
59560d06aSjeremylt 
69560d06aSjeremylt int main(int argc, char **argv) {
79560d06aSjeremylt   Ceed      ceed;
8*4fee36f0SJeremy L Thompson   CeedBasis basis, basis_2;
9*4fee36f0SJeremy L Thompson   CeedInt   p = 4;
109560d06aSjeremylt 
119560d06aSjeremylt   CeedInit(argv[1], &ceed);
129560d06aSjeremylt 
13*4fee36f0SJeremy L Thompson   CeedBasisCreateTensorH1Lagrange(ceed, 1, 1, p, 4, CEED_GAUSS_LOBATTO, &basis);
14*4fee36f0SJeremy L Thompson   CeedBasisCreateTensorH1Lagrange(ceed, 1, 1, p + 1, 4, CEED_GAUSS_LOBATTO, &basis_2);
159560d06aSjeremylt 
16*4fee36f0SJeremy L Thompson   CeedBasisReferenceCopy(basis, &basis_2);  // This destroys the previous basis_2
17*4fee36f0SJeremy L Thompson   CeedBasisDestroy(&basis);
189560d06aSjeremylt 
19*4fee36f0SJeremy L Thompson   CeedInt p_2;
20*4fee36f0SJeremy L Thompson   CeedBasisGetNumNodes1D(basis_2, &p_2);
21*4fee36f0SJeremy L Thompson   if (p != p_2) printf("Error copying CeedBasis reference\n");
229560d06aSjeremylt 
23*4fee36f0SJeremy L Thompson   CeedBasisDestroy(&basis_2);
249560d06aSjeremylt   CeedDestroy(&ceed);
259560d06aSjeremylt   return 0;
269560d06aSjeremylt }
27