Closes #9 : added fallback for diag.

diag in srb-matrix_mod is a wrapper around pdiag and sdiag,
parallel and serial solvers, respectively.

Also swapped out zhengst for zhegst.  Because UPLO = 'U', they are
this same and this saves some logic setting up workspaces...