Displaying 2 results from an estimated 2 matches for "qr_qy_real".
2012 Dec 03
1
qr.qy and qr.qty give an error message when y is integer and LAPACK=TRUE
...torage.mode of y is not set to "double" (in contrast to when Linpack QR has been used).
I assume that the error issued when LAPACK=TRUE is not intended.
Berend
Suggested code change in qr.qy
Replace
if(!is.null(a) && is.logical(a) && a)
return(.Call("qr_qy_real", qr, as.matrix(y), 0, PACKAGE = "base"))
with
if(!is.null(a) && is.logical(a) && a) {
storage.mode(y) <- "double"
return(.Call("qr_qy_real", qr, as.matrix(y), 0, PACKAGE = "base"))
}
and a similar change in...
2004 Mar 01
0
se.contrast ....too hard??? .... Too easy????? .....too trivial???? ...... Too boring.....too????????
...as.matrix(y)
if (!is.complex(y))
y[] <- as.complex(y)
return(.Call("qr_qy_cmplx", qr, y, 1, PACKAGE = "base"))
}
a <- attr(qr, "useLAPACK")
if (!is.null(a) && is.logical(a) && a)
return(.Call("qr_qy_real", qr, as.matrix(y), 1, PACKAGE =
"base"))
n <- nrow(qr$qr)
k <- as.integer(qr$rank)
ny <- NCOL(y)
if (NROW(y) != n)
stop("qr and y must have the same number of rows")
storage.mode(y) <- "double"
.Fortran("dqrqty&qu...