2       subroutine eval(dimn,x,rslt)  argument
3       integer dimn  local
13       do d=1,dimn
27       integer dimn,d  variable
51       do dimn=1,maxdim
52         qdim=q**dimn
53         xdim=2**dimn
55         do d=0,dimn-1
65         call ceedvectorcreate(ceed,xdim*dimn,x,err)
68         call ceedvectorcreate(ceed,qdim*dimn,xq,err)
74         call ceedbasiscreatetensorh1lagrange(ceed,dimn,dimn,2,q,&
76         call ceedbasiscreatetensorh1lagrange(ceed,dimn,1,q,q,&
83           do d=0,dimn-1
86           call eval(dimn,xxx,uuq(i))
95         call ceedbasiscreatetensorh1lagrange(ceed,dimn,dimn,2,q,ceed_gauss,bxg,&
97         call ceedbasiscreatetensorh1lagrange(ceed,dimn,1,q,q,ceed_gauss,bug,err)
105           do d=0,dimn-1
108           call eval(dimn,xxx,fx)
112           write(*,*) uuq(i+offset2),' not equal to ',fx,dimn,i