1 subroutine mkspl2pb(fun,x,nx,th,nth,fspl,nf3,
2 > ibcxmin,bcxmin,ibcxmax,bcxmax,
3 > wk,inwk,ilinx,ilinth,ier)
22 real fspl(4,4,nf3,nth)
39 write(6,
'('.lt.
' ?mkspl2pb -- array dim error, nf3nx'')')
42 if(inwk.lt.5*max(nx,nth))
then
43 write(6,
'('' ?mkspl2pb -- array dim error, inwk too small'')')
49 fspl(1,1,ix,ith)=fun(x(ix),th(ith))
53 call bpsplinb(x,nx,th,nth,fspl,nx,ibcxmin,bcxmin,ibcxmax,bcxmax,
54 > wk,inwk,ilinx,ilinth,ier)