1
subroutine fpsuev(idim,tu,nu,tv,nv,c,u,mu,v,mv,f,wu,wv,lu,lv)
3
integer idim,nu,nv,mu,mv
6
real*8 tu(nu),tv(nv),c((nu-4)*(nv-4)*idim),u(mu),v(mv),
7
* f(mu*mv*idim),wu(mu,4),wv(mv,4)
9
integer i,i1,j,j1,k,l,l1,l2,l3,m,nuv,nu4,nv4
13
c ..subroutine references..
23
if(arg.lt.tb) arg = tb
24
if(arg.gt.te) arg = te
25
10 if(arg.lt.tu(l1) .or. l.eq.nu4) go to 20
29
20 call fpbspl(tu,nu,3,arg,l,h)
42
if(arg.lt.tb) arg = tb
43
if(arg.gt.te) arg = te
44
50 if(arg.lt.tv(l1) .or. l.eq.nv4) go to 60
48
60 call fpbspl(tv,nv,3,arg,l,h)
70
sp = sp+c(l2)*h(i1)*wv(j,j1)