convertpw4qmc.cpp File Reference

Include dependency graph for convertpw4qmc.cpp:

Go to the source code of this file.

Functions
void	convertToVecStrings (int argc, char *argv[], std::vector< std::string > &vec)
bool	file_exists (const std::string &name)
int	main (int argc, char *argv[])

Function Documentation

convertToVecStrings()

void convertToVecStrings	(	int	argc,
		char *	argv[],
		std::vector< std::string > &	vec
	)

Definition at line 24 of file convertpw4qmc.cpp.

Referenced by main().

{
  for (int i = 1; i < argc; i++)
  {
    vec.push_back(argv[i]);
  }
}

file_exists()

bool file_exists ( const std::string &  name )

inline

Definition at line 33 of file convertpw4qmc.cpp.

Referenced by LCAOrbitalBuilder::createSPOSetFromXML(), and main().

{
  ifstream ifile(name.c_str());
  return (bool)ifile;
}

main()

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 39 of file convertpw4qmc.cpp.

References convertToVecStrings(), file_exists(), EshdfFile::writeQboxAtoms(), EshdfFile::writeQboxBoilerPlate(), EshdfFile::writeQboxElectrons(), EshdfFile::writeQboxSupercell(), EshdfFile::writeQEAtoms(), EshdfFile::writeQEBoilerPlate(), EshdfFile::writeQEElectrons(), and EshdfFile::writeQESupercell().

{
#ifdef HAVE_MPI
  int rank = 0, world_size = 1;
  int provided;
  MPI_Init_thread(NULL, NULL, MPI_THREAD_FUNNELED, &provided);
  MPI_Comm_rank(MPI_COMM_WORLD, &rank);
  MPI_Comm_size(MPI_COMM_WORLD, &world_size);
#endif

  vector<string> vecParams;
  convertToVecStrings(argc, argv, vecParams);
  string fname;
  string ofname = "eshdf.h5";
  for (int i = 0; i < vecParams.size(); i++)
  {
    if (vecParams[i] == "--output" || vecParams[i] == "-o")
    {
      ofname = vecParams[i + 1];
      i++;
    }
    else
    {
      fname = vecParams[i];
    }
  }
  if (file_exists(fname) == false)
  {
    cerr << "must specify a valid xml file name as an argument" << endl;
    cerr << fname << " did not work" << endl;
    return (1);
  }

  ifstream is(fname.c_str());

  XmlNode inFile(&is, 0, true);
  string name = inFile.getName();
  string directory("./");
  const int last_slash_index = fname.find_last_of('/');
  if (last_slash_index > -1)
  {
    directory = fname.substr(0, last_slash_index + 1);
  }

  if (name == "qes:espresso")
  {
    // may be good to test different versions of espresso to see
    // if this changes at all
    cout << "xml file comes from quantum espresso" << endl;
    EshdfFile outFile(ofname);
    outFile.writeQEBoilerPlate(inFile);
    outFile.writeQESupercell(inFile);
    outFile.writeQEAtoms(inFile);
    outFile.writeQEElectrons(inFile, directory);
  }
  else if (name == "fpmd:sample")
  {
    cout << "xml file comes from qbox" << endl;
    EshdfFile outFile(ofname);
    outFile.writeQboxBoilerPlate(inFile);
    outFile.writeQboxSupercell(inFile);
    outFile.writeQboxAtoms(inFile);
    outFile.writeQboxElectrons(inFile);
  }
  else
  {
    cerr << "xml format is unrecognized" << endl;
    return (1);
  }

#ifdef HAVE_MPI
  MPI_Finalize();
#endif
  return 0;
}