added output routines for diagnostics of complex solutions in double double
and quad double precision