xref: /libCEED/tests/t307-basis.c (revision 9560d06a92c065fb7d600a8c20ade8d9a4cda324)
1*9560d06aSjeremylt /// @file
2*9560d06aSjeremylt /// Test creation, copying, and distruction of a H1Lagrange basis
3*9560d06aSjeremylt /// \test Test creation, copying, and distruction of a H1Lagrange basis
4*9560d06aSjeremylt #include <ceed.h>
5*9560d06aSjeremylt 
6*9560d06aSjeremylt int main(int argc, char **argv) {
7*9560d06aSjeremylt   Ceed ceed;
8*9560d06aSjeremylt   CeedBasis b, b_2;
9*9560d06aSjeremylt   CeedInt P_1d = 4;
10*9560d06aSjeremylt 
11*9560d06aSjeremylt   CeedInit(argv[1], &ceed);
12*9560d06aSjeremylt 
13*9560d06aSjeremylt   CeedBasisCreateTensorH1Lagrange(ceed, 1, 1, P_1d, 4, CEED_GAUSS_LOBATTO, &b);
14*9560d06aSjeremylt   CeedBasisCreateTensorH1Lagrange(ceed, 1, 1, P_1d+1, 4, CEED_GAUSS_LOBATTO,
15*9560d06aSjeremylt                                   &b_2);
16*9560d06aSjeremylt 
17*9560d06aSjeremylt   CeedBasisReferenceCopy(b, &b_2); // This destroys the previous b_2
18*9560d06aSjeremylt   CeedBasisDestroy(&b);
19*9560d06aSjeremylt 
20*9560d06aSjeremylt   CeedInt P_1d_2;
21*9560d06aSjeremylt   CeedBasisGetNumNodes1D(b_2, &P_1d_2);
22*9560d06aSjeremylt   if (P_1d != P_1d_2)
23*9560d06aSjeremylt     // LCOV_EXCL_START
24*9560d06aSjeremylt     printf("Error copying CeedBasis reference.");
25*9560d06aSjeremylt   // LCOV_EXCL_STOP
26*9560d06aSjeremylt 
27*9560d06aSjeremylt   CeedBasisDestroy(&b_2);
28*9560d06aSjeremylt   CeedDestroy(&ceed);
29*9560d06aSjeremylt   return 0;
30*9560d06aSjeremylt }
31