svd first uses dgesdd_ and if that fails to converge it then
uses dgesvd_. Same for sgesdd_, cgesdd_, and zgesdd_

Geometry::write(...) now does nothing because
valgrind complains of undefined behavior here.
We don't need Geometry to write anyway

IoSimple::Out::write adds an equal
if T is arithmetic type or enum type

__is_enum is a compiler extension.
When we move to C++11 this will become std::is_enum