A random matrix is set up, and then result = a**T * 2*b is computed
this is compared with the result of plain scalapack routines and
the maximum norm of "scalapack - elpa" is computed (should be zero on
machine precision)