xref: /libCEED/tests/t306-basis-f.f90 (revision 144df584ee8fd0c153e6370f8ed4ec2a94d9bc56)
1*144df584Sjeremylt!-----------------------------------------------------------------------
2*144df584Sjeremylt      program test
3*144df584Sjeremylt
4*144df584Sjeremylt      include 'ceedf.h'
5*144df584Sjeremylt
6*144df584Sjeremylt      integer ceed,err
7*144df584Sjeremylt      integer b
8*144df584Sjeremylt      integer p,q
9*144df584Sjeremylt
10*144df584Sjeremylt      character arg*32
11*144df584Sjeremylt
12*144df584Sjeremylt      call getarg(1,arg)
13*144df584Sjeremylt
14*144df584Sjeremylt      call ceedinit(trim(arg)//char(0),ceed,err)
15*144df584Sjeremylt      call ceedbasiscreatetensorh1lagrange(ceed,3,1,4,5,ceed_gauss_lobatto,b,&
16*144df584Sjeremylt     & err)
17*144df584Sjeremylt
18*144df584Sjeremylt      call ceedbasisgetnumnodes(b,p,err)
19*144df584Sjeremylt      call ceedbasisgetnumquadraturepoints(b,q,err)
20*144df584Sjeremylt
21*144df584Sjeremylt
22*144df584Sjeremylt      if (p .NE. 64) then
23*144df584Sjeremylt! LCOV_EXCL_START
24*144df584Sjeremylt        write(*,*) 'Error ',p,' != 64 '
25*144df584Sjeremylt! LCOV_EXCL_STOP
26*144df584Sjeremylt      endif
27*144df584Sjeremylt      if (q .NE. 125) then
28*144df584Sjeremylt! LCOV_EXCL_START
29*144df584Sjeremylt        write(*,*) 'Error ',q,' != 125 '
30*144df584Sjeremylt! LCOV_EXCL_STOP
31*144df584Sjeremylt      endif
32*144df584Sjeremylt
33*144df584Sjeremylt      call ceedbasisdestroy(b,err)
34*144df584Sjeremylt      call ceeddestroy(ceed,err)
35*144df584Sjeremylt
36*144df584Sjeremylt      end
37*144df584Sjeremylt!-----------------------------------------------------------------------
38