5 Nov
2004
5 Nov
'04
9:59 p.m.
On Thu, Nov 04, 2004 at 06:53:56PM -0500, Jonathan Brandmeyer wrote:
Sometimes, inverting a matrix with complex values does not produce something that, when multiplied by the original, returns the identity.
[huge snip] Actually, the output of your program *is* the unit matrix. The issue is that GiNaC can not factorize polynomials properly, that's why normal() produces a huge expression instead of simply 1. Best regards, Alexei.