   double relError_S, relError_D, relError_C, relError_Z;
 
   // Machine precision
-  double eps_float = std::numeric_limits<float>::epsilon();
+  double eps_float  = std::numeric_limits<float>::epsilon();
   double eps_double = std::numeric_limits<double>::epsilon();
 
   // Float (S)

// TODO: handle side='L' and diag='U'

#include <cstdlib>
#include <cstring>
#include <string>


  // Check for help command line argument
  for(int i=1; i<argc; ++i) {
    if(strcmp(argv[i],"--help")==0 || strcmp(argv[i],"-h")==0) {

      // Display help page
      std::cout << "\n"
