convertMatrixFile.cpp File Reference

Functions
int	main (int argc, char **argv)

Function Documentation

main()

int main	(	int	argc,
		char **	argv
	)

                                {
  if(argc < 2) {
    std::cerr << "Usage: 1.This 2.Matrix file (opt1).UPLO" << std::endl;
    std::exit(EX_USAGE);
  }
  std::string   InFilename(argv[1]);
  std::ifstream InFile(InFilename); checkIsFileOpen(InFile, InFilename);
  debug_print("Input: " << InFilename);
 
  matrix<Complex_t<double>> mat;
  auto [dim1, dim2, inUPLO] = readMatrix(InFile, mat); InFile.close();
  if(dim1 != dim2) {
    std::cerr << "Error: dim1(" << dim1 <<") != dim2(" << dim2 << ")" << std::endl;
    std::exit(EX_USAGE);
  }
 
  debug_print(inUPLO);
  for(int i = 0;i < dim1; ++i) mat.at(i,i) = real(mat.at(i,i));
  if(inUPLO=="Upper") {
      #pragma omp parallel for schedule(dynamic, 1)
      for(int i = 0;i < dim1; ++i) for(int j = 0;j < i; ++j) mat.at(i,j) = conj(mat.at(j,i));
  } else if (inUPLO=="Lower") {
      #pragma omp parallel for schedule(dynamic, 1)
      for(int i = 0;i < dim1; ++i) for(int j = i+1;j < dim1; ++j) mat.at(i,j) = conj(mat.at(j,i));
  } else {
      #pragma omp parallel for schedule(dynamic, 1)
      for(Integer_t i = 0;i < dim1; ++i) for(Integer_t j = i+1;j < dim1; ++j) {
        mat.at(i,j) = ( mat.at(i,j) + conj(mat.at(j,i)) )/2.0;
        mat.at(j,i) = conj(mat.at(i,j));
      }
  }
 
  std::string outUPLO(argc>=3 ? argv[2]: "All");
  std::cout << std::scientific; mat.print(std::cout, dim1, dim2, outUPLO);
 
  return 0;
}