generateRM.hpp

Go to the documentation of this file.

// Generate a random matrix used by all programs in this directory
#ifndef GENERATE_RM
#define GENERATE_RM
 
#include "Headers/LocalRandomMatrix.hpp"
#include "partWithMax.hpp"
#include "Headers/mersenne_twister.hpp"
#include <iostream>
#include <sysexits.h>
#define _USE_MATH_DEFINES
#include <cmath>
#ifndef M_PI
  #define M_PI 3.14159265358979323846
#endif
 
int const N_ARGS        = 6;
std::string const USAGE_sp = "";
 
static inline bool Initialize(int argc, char** argv) { return true; }
static inline bool Finalize(int argc, char** argv) { return true; }
 
static inline void generateLocal_h(matrix<Complex_t>& mat, Integer_t const d_lmat, Integer_t const seed){
  if(seed >= 0) {
    init_genrand(seed);
    std::cout << "# init_genrand(" << seed << ")\n";
  }
  // RandomMatrix_GUE(mat, d_lmat);
};
 
static inline void generateLocal_op(matrix<Complex_t>& mat, Integer_t const d_lmat, Integer_t const seed){
  if(seed >= 0) init_genrand(seed);
  // RandomMatrix_GUE(mat, d_lmat);
};
 
static inline void constructGlobal_h(matrix<Complex_t>& global, matrix<Complex_t> const& local, Integer_t const m, Integer_t const n, Integer_t const dim_loc){
  if(dim_loc != 2) {
    std::cerr << "Error: constructGlobal_h (dim_loc=" << dim_loc << ")" << std::endl;
    std::exit(EX_USAGE);
  }
  RandomMatrix_mBody_lLocal(global, m, n, dim_loc);
};
 
static inline void constructGlobal_op(matrix<Complex_t>& global, matrix<Complex_t> const& local, Integer_t const m, Integer_t const n, Integer_t const dim_loc){
  constructGlobal_h(global, local, m, n, dim_loc);
};
 
#endif