swh:1:snp:dfbb8ae2fb3632e8a8608d675a49ab1f110b7c6d
Tip revision: 56c6d241a6642cc8bd7ee1b4b209bf9888daa74c authored by Doug Nychka on 20 October 2008, 00:00:00 UTC
version 5.01
version 5.01
Tip revision: 56c6d24
cvrf.f
c** finds value of h minimizing the generalized cross-validation
subroutine cvrf
* (h,nobs,x,y,wt,sy,trace,diag,din,dout,cv,ierr)
implicit double precision (a-h,o-z)
REAL*8 h,trace,cv
REAL*8 x(nobs),y(nobs),wt(nobs)
REAL*8 sy(nobs),diag(nobs),dumm1(1),dumm2(1)
real*8 din(10), dout(10)
integer ngrid, ierr, job(3),ideriv
data job/3,0,0/
data ideriv,ngrid/0,1/
call rcss(h,nobs,x,y,wt,sy,trace,diag,cv,
+ ngrid,dumm1,dumm2,job,ideriv,din,dout,ierr)
nit= int( dout(1))
trace=dout(3)
c
return
end