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/mytypes.hpp"
#include "Headers/mydebug.hpp"
#include "Headers/LocalRandomMatrix.hpp"
#include "Headers/mersenne_twister.hpp"
#include "Headers/matrix_util.hpp"
#include "Headers/translation.hpp"
#ifdef GPU
    #include "magma_v2.h"
    #include <cuda_runtime_api.h>
#endif
 
std::vector<Real_t*> prefactor;
static inline bool   Initialize(int const argc, char const* const* argv, int const Nargs_common) {
      return true;
}
static inline bool Finalize(int const argc, char const* const* argv) { return true; }
 
static inline Integer_t SectorDimension(TransSector const& Sector) {
    return (Sector.momentum() == 0) ? Sector.dim() - Sector.parityOp().numPairs : Sector.dim();
}
 
template<typename Type>
static inline void generateLocal_h(matrix<Type>& mat, Integer_t const d_mat, Integer_t const seed) {
    return;
}
template<typename Type>
static inline void generateLocal_op(matrix<Type>& mat, Integer_t const d_mat,
                                    Integer_t const seed) {
    return;
}
 
static inline Integer_t constructGlobal_h(matrix<Complex_t>& global, matrix<Complex_t> const& local,
                                          Integer_t const numLoc, TransSector const& Sector,
                                          void* GPUconf = nullptr) {
    RandomMatrix_GUE(global, SectorDimension(Sector));
    return SectorDimension(Sector);
};
 
static inline Integer_t constructGlobal_op(matrix<Complex_t>&       global,
                                           matrix<Complex_t> const& local, Integer_t const numLoc,
                                           TransSector const& Sector, void* GPUconf = nullptr) {
    return constructGlobal_h(global, local, numLoc, Sector);
};
 
#endif