improve UDV storage to handle larger scales
@Hofmann has the issue that his matrices quickly exceed the capabilities of the double precision storage in D of a UDV decomposition. Let's try whether using 2x64bit for the double resolves that issue. This should not degrade the performance since all operations with D scale linearly.