In an attempt to speed up Mathematica's PseudoInverse function, I wonder if helping it zero out small numbers would work.
With cubic scaling for CPU time and quadratic scaling for memory, I don't know that this will help much but I would like to try.
Can the Tolerance parameter be set to do this zeroing of small numbers?
To be more accurate, these are numbers close to zero in the calculation, not the solution.
