[GiNaC-list] Inverse of complex matrix

[Sheplyakov Alexei]

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.