Changes in src/moleculelist.cpp [35b698:274d45]

File:

• src/moleculelist.cpp (modified) (12 diffs)

src/moleculelist.cpp

@@ -12 +12 @@
 #include "atom.hpp"
 #include "bond.hpp"
-#include "bondgraph.hpp"
 #include "boundary.hpp"
 #include "config.hpp"
@@ -380 +379 @@
  * bonded to the same atom, then we add for this pair a correction term constructed from a Morse
  * potential function fit to QM calculations with respecting to the interatomic hydrogen distance.
- * \param &path path to file
- */
-bool MoleculeListClass::AddHydrogenCorrection(std::string &path)
+ * \param *out output stream for debugging
+ * \param *path path to file
+ */
+bool MoleculeListClass::AddHydrogenCorrection(char *path)
 {
   bond *Binder = NULL;
@@ -400 +400 @@
   // 0a. find dimension of matrices with constants
   line = path;
+  line.append("/");
+  line += FRAGMENTPREFIX;
   line += "1";
   line += FITCONSTANTSUFFIX;
   input.open(line.c_str());
-  if (input.fail()) {
+  if (input == NULL) {
     DoLog(1) && (Log() << Verbose(1) << endl << "Unable to open " << line << ", is the directory correct?" << endl);
     return false;
@@ -567 +569 @@
 
 /** Store force indices, i.e. the connection between the nuclear index in the total molecule config and the respective atom in fragment config.
- * \param &path path to file
+ * \param *out output stream for debugging
+ * \param *path path to file
  * \param *SortIndex Index to map from the BFS labeling to the sequence how of Ion_Type in the config
  * \return true - file written successfully, false - writing failed
  */
-bool MoleculeListClass::StoreForcesFile(std::string &path, int *SortIndex)
+bool MoleculeListClass::StoreForcesFile(char *path,
+    int *SortIndex)
 {
   bool status = true;
-  string filename(path);
-  filename += FORCESFILE;
-  ofstream ForcesFile(filename.c_str());
+  ofstream ForcesFile;
+  stringstream line;
   periodentafel *periode=World::getInstance().getPeriode();
 
   // open file for the force factors
   DoLog(1) && (Log() << Verbose(1) << "Saving  force factors ... ");
-  if (!ForcesFile.fail()) {
+  line << path << "/" << FRAGMENTPREFIX << FORCESFILE;
+  ForcesFile.open(line.str().c_str(), ios::out);
+  if (ForcesFile != NULL) {
     //Log() << Verbose(1) << "Final AtomicForcesList: ";
     //output << prefix << "Forces" << endl;
@@ -606 +611 @@
   } else {
     status = false;
-    DoLog(1) && (Log() << Verbose(1) << "failed to open file " << filename << "." << endl);
+    DoLog(1) && (Log() << Verbose(1) << "failed to open file " << line.str() << "." << endl);
   }
   ForcesFile.close();
@@ -615 +620 @@
 /** Writes a config file for each molecule in the given \a **FragmentList.
  * \param *out output stream for debugging
- * \param &prefix path and prefix to the fragment config files
+ * \param *configuration standard configuration to attach atoms in fragment molecule to.
  * \param *SortIndex Index to map from the BFS labeling to the sequence how of Ion_Type in the config
  * \return true - success (each file was written), false - something went wrong.
  */
-bool MoleculeListClass::OutputConfigForListOfFragments(std::string &prefix, int *SortIndex)
+bool MoleculeListClass::OutputConfigForListOfFragments(config *configuration, int *SortIndex)
 {
   ofstream outputFragment;
-  std::string FragmentName;
+  char FragmentName[MAXSTRINGSIZE];
   char PathBackup[MAXSTRINGSIZE];
   bool result = true;
@@ -640 +645 @@
   for (MoleculeList::iterator ListRunner = ListOfMolecules.begin(); ListRunner != ListOfMolecules.end(); ListRunner++) {
     // save default path as it is changed for each fragment
-    path = World::getInstance().getConfig()->GetDefaultPath();
+    path = configuration->GetDefaultPath();
     if (path != NULL)
       strcpy(PathBackup, path);
@@ -653 +658 @@
     // output xyz file
     FragmentNumber = FixedDigitNumber(ListOfMolecules.size(), FragmentCounter++);
-    FragmentName = prefix + FragmentNumber + ".conf.xyz";
-    outputFragment.open(FragmentName.c_str(), ios::out);
+    sprintf(FragmentName, "%s/%s%s.conf.xyz", configuration->configpath, FRAGMENTPREFIX, FragmentNumber);
+    outputFragment.open(FragmentName, ios::out);
     DoLog(2) && (Log() << Verbose(2) << "Saving bond fragment No. " << FragmentNumber << "/" << FragmentCounter - 1 << " as XYZ ...");
     if ((intermediateResult = (*ListRunner)->OutputXYZ(&outputFragment)))
@@ -677 +682 @@
     for (int k = 0; k < NDIM; k++) {
       j += k + 1;
-      BoxDimension[k] = 2.5 * (World::getInstance().getConfig()->GetIsAngstroem() ? 1. : 1. / AtomicLengthToAngstroem);
+      BoxDimension[k] = 2.5 * (configuration->GetIsAngstroem() ? 1. : 1. / AtomicLengthToAngstroem);
       cell_size[j] = BoxDimension[k] * 2.;
     }
@@ -684 +689 @@
     // also calculate necessary orbitals
     (*ListRunner)->CountElements(); // this is a bugfix, atoms should shoulds actually be added correctly to this fragment
-    //(*ListRunner)->CalculateOrbitals(*World::getInstance().getConfig);
+    (*ListRunner)->CalculateOrbitals(*configuration);
 
     // change path in config
-    FragmentName = PathBackup;
-    FragmentName += "/";
-    FragmentName += FRAGMENTPREFIX;
-    FragmentName += FragmentNumber;
-    FragmentName += "/";
-    World::getInstance().getConfig()->SetDefaultPath(FragmentName.c_str());
+    //strcpy(PathBackup, configuration->configpath);
+    sprintf(FragmentName, "%s/%s%s/", PathBackup, FRAGMENTPREFIX, FragmentNumber);
+    configuration->SetDefaultPath(FragmentName);
 
     // and save as config
-    FragmentName = prefix + FragmentNumber + ".conf";
+    sprintf(FragmentName, "%s/%s%s.conf", configuration->configpath, FRAGMENTPREFIX, FragmentNumber);
     DoLog(2) && (Log() << Verbose(2) << "Saving bond fragment No. " << FragmentNumber << "/" << FragmentCounter - 1 << " as config ...");
-    if ((intermediateResult = World::getInstance().getConfig()->Save(FragmentName.c_str(), (*ListRunner)->elemente, (*ListRunner))))
+    if ((intermediateResult = configuration->Save(FragmentName, (*ListRunner)->elemente, (*ListRunner))))
       DoLog(0) && (Log() << Verbose(0) << " done." << endl);
     else
@@ -704 +706 @@
 
     // restore old config
-    World::getInstance().getConfig()->SetDefaultPath(PathBackup);
+    configuration->SetDefaultPath(PathBackup);
 
     // and save as mpqc input file
-    FragmentName = prefix + FragmentNumber + ".conf";
+    sprintf(FragmentName, "%s/%s%s.conf", configuration->configpath, FRAGMENTPREFIX, FragmentNumber);
     DoLog(2) && (Log() << Verbose(2) << "Saving bond fragment No. " << FragmentNumber << "/" << FragmentCounter - 1 << " as mpqc input ...");
-    if ((intermediateResult = World::getInstance().getConfig()->SaveMPQC(FragmentName.c_str(), (*ListRunner))))
+    if ((intermediateResult = configuration->SaveMPQC(FragmentName, (*ListRunner))))
       DoLog(2) && (Log() << Verbose(2) << " done." << endl);
     else
@@ -765 +767 @@
 
   // 1. dissect the molecule into connected subgraphs
-  if (configuration->BG != NULL) {
-    if (!configuration->BG->ConstructBondGraph(mol)) {
-      World::getInstance().destroyMolecule(mol);
-      DoeLog(1) && (eLog()<< Verbose(1) << "There are no bonds." << endl);
-      return;
-    }
-  } else {
-    DoeLog(1) && (eLog()<< Verbose(1) << "There is no BondGraph class present to create bonds." << endl);
+  if (!configuration->BG->ConstructBondGraph(mol)) {
+    World::getInstance().destroyMolecule(mol);
+    DoeLog(1) && (eLog()<< Verbose(1) << "There are no bonds." << endl);
     return;
   }