[GiNaC-list] Inverse of complex matrix

[Jonathan Brandmeyer]

Sometimes, inverting a matrix with complex values does not produce
something that, when multiplied by the original, returns the identity.
I have attached a program that demonstrates this for one case, which was
found by altering checks/check_matricies.cpp to use complex values from
sparse_tree().

Is this a known problem?

Thanks,
-Jonathan Brandmeyer

P.S: GiNaC 1.3.0, G++ 3.3.4+ (Debian), CLN 1.1.8, Debian Sid on i386
-------------- next part --------------
A non-text attachment was scrubbed...
Name: complex_inverse_test.cpp
Type: text/x-c++src
Size: 503 bytes
Desc: not available
Url : http://www.cebix.net/pipermail/ginac-list/attachments/20041105/a4315628/complex_inverse_test.bin