Displaying 3 results from an estimated 3 matches for "hinv3".
Did you mean:
hinv
2009 Jun 17
1
Inverting a square matrix using solve() with LAPACK=TRUE (PR#13762)
...nvert a matrix using "LAPACK" or "LINPACK"
if (nrow(A) != ncol(A)) stop("Matrix muxt be square")
qrA <- qr(A, LAPACK=LAPACK, tol=tol)
if (LAPACK) {
apply(diag(1, ncol(A)), 2, function(x) solve(qrA, x))
} else solve(qrA)
}
hinv3 <- solve.lapack(h5)
all.equal(hinv1, hinv3) # Now, they are equal
2009 Jun 18
1
Inverting a square... (PR#13762)
...INPACK"
> if (nrow(A) !=3D ncol(A)) stop("Matrix muxt be square")
> qrA <- qr(A, LAPACK=3DLAPACK, tol=3Dtol)
> if (LAPACK) {
> apply(diag(1, ncol(A)), 2, function(x) solve(qrA, x))=20
> } else solve(qrA)
> }
>=20
> hinv3 <- solve.lapack(h5)
> all.equal(hinv1, hinv3) # Now, they are equal
>=20
> ______________________________________________
> R-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
--=20
O__ ---- Peter Dalgaard =C3=98ster Farimagsga...
2009 Jun 18
0
Inverting a square matrix using solve() with LAPACK=TRUE (PR#13765)
...or "LINPACK"
> if (nrow(A) != ncol(A)) stop("Matrix muxt be square")
> qrA <- qr(A, LAPACK=LAPACK, tol=tol)
> if (LAPACK) {
> apply(diag(1, ncol(A)), 2, function(x) solve(qrA, x))
> } else solve(qrA)
> }
>
> hinv3 <- solve.lapack(h5)
> all.equal(hinv1, hinv3) # Now, they are equal
>
> ______________________________________________
> R-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
--
O__ ---- Peter Dalgaard ?ster Farimagsgade 5, En...