Displaying 1 result from an estimated 1 matches for "2ecall".
Did you mean:
recall
2009 Jun 18
1
Inverting a square... (PR#13762)
...adhan
> Version: 2.8.1
> OS: Windows
> Submission from: (NULL) (162.129.251.19)
>=20
>=20
> Inverting a matrix with solve(), but using LAPACK=3DTRUE, gives erroneo=
us
> results:
Thanks, but there seems to be a much easier fix.
Inside coef.qr, we have
coef[qr$pivot, ] <-
=2ECall("qr_coef_real", qr, y, PACKAGE =3D "base")[seq_len(p)]
which should be [seq_len(p),]
(otherwise, in the matrix case, the RHS will recycle only the 1st p
elements, i.e., the 1st column).
>=20
> Here is an example:
>=20
> hilbert <- function(n) { i <- 1:n...