Added a symmetric eigensolve which uses pdgetrf as a subcall for 1D QRs and added
a separate test for this code, as well as a non-separate benchmark. Also fixed corner
value in QR_scala_2D, and added as well as fixed-up profiling.