2       subroutine eval(dimn,x,rslt)  argument
3       integer dimn  local
8       if (dimn>1) then
11       if (dimn>2) then
24       integer dimn,d  variable
56       do dimn=1,maxdim
57         qdimn=q**dimn
58         pdimn=p**dimn
59         xdim=2**dimn
60         dimxqdimn=dimn*qdimn
64         do d=0,dimn-1
74         call ceedvectorcreate(ceed,xdim*dimn,x,err)
77         call ceedvectorcreate(ceed,pdimn*dimn,xq,err)
80         call ceedvectorcreate(ceed,qdimn*dimn,uq,err)
82         call ceedvectorcreate(ceed,qdimn*dimn,ones,err)
87         call ceedbasiscreatetensorh1lagrange(ceed,dimn,dimn,2,p,&
93           do d=0,dimn-1
96           call eval(dimn,xxx,uu(i))
102         call ceedbasiscreatetensorh1lagrange(ceed,dimn,1,p,q,ceed_gauss,bug,err)
123           write(*,'(A,I1,A,F12.6,A,F12.6)')'[',dimn,'] Error: ',sum1,  ' != ',&