https://github.com/cran/fields
Tip revision: ce722edae3c1b9e1af2985ce3500b11058facf0e authored by Doug Nychka on 24 August 2006, 01:46:17 UTC
version 3.04
version 3.04
Tip revision: ce722ed
Wimage.cov.R
"Wimage.cov" <-
function( Y=NULL,ind=NULL, H.obj, find.row=FALSE ){
# local multipcation by H function
Mult.H<-function( u, H.obj){
IND<- H.obj$ind0
u[IND] <- H.obj$H0%*%c(u[IND])
u*H.obj$H1
}
#
# This is a shortcut option to find a single row of the covariance matrix
#
if( find.row){ Y<-1}
#
# if ind is passed fill complete image with the Y subset
#
if(is.null(ind)){
tmp<- Y}
else{
tmp<- matrix(0,nrow=H.obj$m,ncol=H.obj$n)
tmp[ind]<- Y
}
#
#
hold<- Wtransform.image( tmp, cut.min=H.obj$cut.min, inv=TRUE,transpose=TRUE)
#
hold<- Mult.H( hold, H.obj=H.obj) # Note if H not symmetric this should be
# multiplication by transpose.
hold<- Mult.H( hold, H.obj=H.obj)
Wtransform.image( hold, cut.min=H.obj$cut.min, inv=TRUE)
}