how about we use same method for invert(A) as for A^^-1 ?
Subject: how about we use same method for invert(A) as for A^^-1 ?
From: Robert Dodier
Date: Mon, 10 Dec 2012 16:36:36 +0000 (UTC)
Hello, invert(A) computes the matrix inverse via the adjoint method,
while A^^-1 uses Gaussian elimination. Both of them can handle symbolic
expressions and CRE in addition to numbers. Is there any reason to
prefer the adjoint method? The adjoint method is much slower.
In any event, it seems messy to have invert(A) and A^^-1 use different
methods ....
There are other methods to compute a matrix inverse -- there are, at
least, tminverse (undocumented core function), invert_by_lu (share), and
probably one in LAPACK.
best,
Robert Dodier