https://github.com/cran/fields
Tip revision: 292e80962d92c497484ad0735f3f13a4b8d3e0ac authored by Doug Nychka on 29 November 2007, 14:39:12 UTC
version 4.1
version 4.1
Tip revision: 292e809
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)
}