some work on the decomposition. there is three implementations now:
the old one for testing purposes, a faster one in pure sage, and one
that is in compiled code. however, due to the complex types the compiled
version only gives a marginal speed-up.
"compiled" version of matrix decomposition