1*52bfb9bbSJeremy L Thompson!----------------------------------------------------------------------- 2*52bfb9bbSJeremy L Thompson! 3*52bfb9bbSJeremy L Thompson! Header with common subroutine 4*52bfb9bbSJeremy L Thompson! 5*52bfb9bbSJeremy L Thompson include 't320-basis-f.h' 6*52bfb9bbSJeremy L Thompson!----------------------------------------------------------------------- 7*52bfb9bbSJeremy L Thompson program test 8*52bfb9bbSJeremy L Thompson 9*52bfb9bbSJeremy L Thompson include 'ceedf.h' 10*52bfb9bbSJeremy L Thompson 11*52bfb9bbSJeremy L Thompson integer ceed,err 12*52bfb9bbSJeremy L Thompson integer p,q,d 13*52bfb9bbSJeremy L Thompson parameter(p=6) 14*52bfb9bbSJeremy L Thompson parameter(q=4) 15*52bfb9bbSJeremy L Thompson parameter(d=2) 16*52bfb9bbSJeremy L Thompson 17*52bfb9bbSJeremy L Thompson real*8 qref(d*q) 18*52bfb9bbSJeremy L Thompson real*8 qweight(q) 19*52bfb9bbSJeremy L Thompson real*8 interp(p*q) 20*52bfb9bbSJeremy L Thompson real*8 grad(d*p*q) 21*52bfb9bbSJeremy L Thompson 22*52bfb9bbSJeremy L Thompson integer b 23*52bfb9bbSJeremy L Thompson 24*52bfb9bbSJeremy L Thompson character arg*32 25*52bfb9bbSJeremy L Thompson 26*52bfb9bbSJeremy L Thompson call getarg(1,arg) 27*52bfb9bbSJeremy L Thompson 28*52bfb9bbSJeremy L Thompson call ceedinit(trim(arg)//char(0),ceed,err) 29*52bfb9bbSJeremy L Thompson 30*52bfb9bbSJeremy L Thompson call buildmats(qref,qweight,interp,grad) 31*52bfb9bbSJeremy L Thompson 32*52bfb9bbSJeremy L Thompson call ceedbasiscreateh1(ceed,ceed_triangle,1,p,q,interp,grad,qref,qweight,& 33*52bfb9bbSJeremy L Thompson & b,err) 34*52bfb9bbSJeremy L Thompson call ceedbasisview(b,err) 35*52bfb9bbSJeremy L Thompson 36*52bfb9bbSJeremy L Thompson call ceedbasisdestroy(b,err) 37*52bfb9bbSJeremy L Thompson call ceeddestroy(ceed,err) 38*52bfb9bbSJeremy L Thompson 39*52bfb9bbSJeremy L Thompson end 40*52bfb9bbSJeremy L Thompson!----------------------------------------------------------------------- 41