dlaed3 seems to be working. TODO test with larger matrices