R/load.Hst.ls.2Zs.R

In widals: Weighting by Inverse Distance with Adaptive Least Squares

load.Hst.ls.2Zs <-
function (Z, Z.na, Hst.ls.Z, xwhich, rgr.lags = c(0)) 
{
    tau <- tau
    min.ndx <- max(1, -min(rgr.lags) + 1)
    max.ndx <- min(tau, tau - max(rgr.lags))
    for (i in min.ndx:max.ndx) {
        zi.na <- Z.na[i, ]
        Hst.ls.Z[[i]][!zi.na, 2 * xwhich - 1] <- Z[i + rgr.lags, 
            !zi.na]
        Hst.ls.Z[[i]][zi.na, 2 * xwhich - 1] <- 0
        Hst.ls.Z[[i]][zi.na, 2 * xwhich] <- Z[i + rgr.lags, zi.na]
        Hst.ls.Z[[i]][!zi.na, 2 * xwhich] <- 0
    }
    return(Hst.ls.Z)
}