Changes in src/molecule_fragmentation.cpp [35b698:112b09]

File:

• src/molecule_fragmentation.cpp (modified) (20 diffs)

src/molecule_fragmentation.cpp

@@ -82 +82 @@
  * -# Scans TEFactors file and sets the TEFactor of each key set in the temporary graph accordingly
  * Finally, the temporary graph is inserted into the given \a FragmentList for return.
- * \param &path path to file
+ * \param *out output stream for debugging
+ * \param *path path to file
  * \param *FragmentList empty, filled on return
  * \return true - parsing successfully, false - failure on parsing (FragmentList will be NULL)
  */
-bool ParseKeySetFile(std::string &path, Graph *&FragmentList)
+bool ParseKeySetFile(char *path, Graph *&FragmentList)
 {
   bool status = true;
@@ -93 +94 @@
   GraphTestPair testGraphInsert;
   int NumberOfFragments = 0;
-  string filename;
+  char filename[MAXSTRINGSIZE];
 
   if (FragmentList == NULL) { // check list pointer
@@ -101 +102 @@
   // 1st pass: open file and read
   DoLog(1) && (Log() << Verbose(1) << "Parsing the KeySet file ... " << endl);
-  filename = path + KEYSETFILE;
-  InputFile.open(filename.c_str());
-  if (InputFile.good()) {
+  sprintf(filename, "%s/%s%s", path, FRAGMENTPREFIX, KEYSETFILE);
+  InputFile.open(filename);
+  if (InputFile != NULL) {
     // each line represents a new fragment
     char buffer[MAXSTRINGSIZE];
@@ -180 +181 @@
 
 /** Stores key sets to file.
+ * \param *out output stream for debugging
  * \param KeySetList Graph with Keysets
- * \param &path path to file
+ * \param *path path to file
  * \return true - file written successfully, false - writing failed
  */
-bool StoreKeySetFile(Graph &KeySetList, std::string &path)
-{
+bool StoreKeySetFile(Graph &KeySetList, char *path)
+{
+  ofstream output;
   bool status =  true;
-  string line = path + KEYSETFILE;
-  ofstream output(line.c_str());
+  string line;
 
   // open KeySet file
+  line = path;
+  line.append("/");
+  line += FRAGMENTPREFIX;
+  line += KEYSETFILE;
+  output.open(line.c_str(), ios::out);
   DoLog(1) && (Log() << Verbose(1) << "Saving key sets of the total graph ... ");
-  if(output.good()) {
+  if(output != NULL) {
     for(Graph::iterator runner = KeySetList.begin(); runner != KeySetList.end(); runner++) {
       for (KeySet::iterator sprinter = (*runner).first.begin();sprinter != (*runner).first.end(); sprinter++) {
@@ -295 +302 @@
 
 /** Scans the adaptive order file and insert (index, value) into map.
- * \param &path path to ENERGYPERFRAGMENT file (may be NULL if Order is non-negative)
+ * \param *out output stream for debugging
+ * \param *path path to ENERGYPERFRAGMENT file (may be NULL if Order is non-negative)
  * \param &IndexedKeySetList list to find key set for a given index \a No
  * \return adaptive criteria list from file
  */
-map<int, pair<double,int> > * ScanAdaptiveFileIntoMap(std::string &path, map<int,KeySet> &IndexKeySetList)
+map<int, pair<double,int> > * ScanAdaptiveFileIntoMap(char *path, map<int,KeySet> &IndexKeySetList)
 {
   map<int, pair<double,int> > *AdaptiveCriteriaList = new map<int, pair<double,int> >;
@@ -305 +313 @@
   double Value = 0.;
   char buffer[MAXSTRINGSIZE];
-  string filename = path + ENERGYPERFRAGMENT;
-  ifstream InputFile(filename.c_str());
-
-  if (InputFile.fail()) {
-    DoeLog(1) && (eLog() << Verbose(1) << "Cannot find file " << filename << "." << endl);
-    return AdaptiveCriteriaList;
-  }
+  sprintf(buffer, "%s/%s%s.dat", path, FRAGMENTPREFIX, ENERGYPERFRAGMENT);
+  ifstream InputFile(buffer, ios::in);
 
   if (CountLinesinFile(InputFile) > 0) {
@@ -416 +419 @@
 
 /** Checks whether the OrderAtSite is still below \a Order at some site.
+ * \param *out output stream for debugging
  * \param *AtomMask defines true/false per global Atom::nr to mask in/out each nuclear site, used to activate given number of site to increment order adaptively
  * \param *GlobalKeySetList list of keysets with global ids (valid in "this" molecule) needed for adaptive increase
  * \param Order desired Order if positive, desired exponent in threshold criteria if negative (0 is single-step)
  * \param *MinimumRingSize array of max. possible order to avoid loops
- * \param path path to ENERGYPERFRAGMENT file (may be NULL if Order is non-negative)
+ * \param *path path to ENERGYPERFRAGMENT file (may be NULL if Order is non-negative)
  * \return true - needs further fragmentation, false - does not need fragmentation
  */
-bool molecule::CheckOrderAtSite(bool *AtomMask, Graph *GlobalKeySetList, int Order, int *MinimumRingSize, std::string path)
+bool molecule::CheckOrderAtSite(bool *AtomMask, Graph *GlobalKeySetList, int Order, int *MinimumRingSize, char *path)
 {
   bool status = false;
@@ -581 +585 @@
  * of vertex indices: Global always means the index in "this" molecule, whereas local refers to the molecule or
  * subgraph in the MoleculeListClass.
+ * \param *out output stream for debugging
  * \param Order up to how many neighbouring bonds a fragment contains in BondOrderScheme::BottumUp scheme
- * \param &prefix path and prefix of the bond order configs to be written
+ * \param *configuration configuration for writing config files for each fragment
  * \return 1 - continue, 2 - stop (no fragmentation occured)
  */
-int molecule::FragmentMolecule(int Order, std::string &prefix)
+int molecule::FragmentMolecule(int Order, config *configuration)
 {
   MoleculeListClass *BondFragments = NULL;
@@ -619 +624 @@
 
   // === compare it with adjacency file ===
-  FragmentationToDo = FragmentationToDo && CheckAdjacencyFileAgainstMolecule(prefix, ListOfAtoms);
+  FragmentationToDo = FragmentationToDo && CheckAdjacencyFileAgainstMolecule(configuration->configpath, ListOfAtoms);
   delete[](ListOfAtoms);
 
@@ -653 +658 @@
 
   // ===== 3. if structure still valid, parse key set file and others =====
-  FragmentationToDo = FragmentationToDo && ParseKeySetFile(prefix, ParsedFragmentList);
+  FragmentationToDo = FragmentationToDo && ParseKeySetFile(configuration->configpath, ParsedFragmentList);
 
   // ===== 4. check globally whether there's something to do actually (first adaptivity check)
-  FragmentationToDo = FragmentationToDo && ParseOrderAtSiteFromFile(prefix);
+  FragmentationToDo = FragmentationToDo && ParseOrderAtSiteFromFile(configuration->configpath);
 
   // =================================== Begin of FRAGMENTATION ===============================
@@ -667 +672 @@
   AtomMask[getAtomCount()] = false;
   FragmentationToDo = false;  // if CheckOrderAtSite just ones recommends fragmentation, we will save fragments afterwards
-  while ((CheckOrder = CheckOrderAtSite(AtomMask, ParsedFragmentList, Order, MinimumRingSize, prefix))) {
+  while ((CheckOrder = CheckOrderAtSite(AtomMask, ParsedFragmentList, Order, MinimumRingSize, configuration->configpath))) {
     FragmentationToDo = FragmentationToDo || CheckOrder;
     AtomMask[getAtomCount()] = true;   // last plus one entry is used as marker that we have been through this loop once already in CheckOrderAtSite()
@@ -722 +727 @@
     KeySet test = (*runner).first;
     DoLog(0) && (Log() << Verbose(0) << "Fragment No." << (*runner).second.first << " with TEFactor " << (*runner).second.second << "." << endl);
-    BondFragments->insert(StoreFragmentFromKeySet(test, World::getInstance().getConfig()));
+    BondFragments->insert(StoreFragmentFromKeySet(test, configuration));
     k++;
   }
@@ -734 +739 @@
 
     DoLog(1) && (Log() << Verbose(1) << "Writing " << BondFragments->ListOfMolecules.size() << " possible bond fragmentation configs" << endl);
-    if (BondFragments->OutputConfigForListOfFragments(prefix, SortIndex))
+    if (BondFragments->OutputConfigForListOfFragments(configuration, SortIndex))
       DoLog(1) && (Log() << Verbose(1) << "All configs written." << endl);
     else
@@ -740 +745 @@
 
     // store force index reference file
-    BondFragments->StoreForcesFile(prefix, SortIndex);
+    BondFragments->StoreForcesFile(configuration->configpath, SortIndex);
 
     // store keysets file
-    StoreKeySetFile(TotalGraph, prefix);
+    StoreKeySetFile(TotalGraph, configuration->configpath);
 
     {
       // store Adjacency file
-      std::string filename = prefix + ADJACENCYFILE;
-      StoreAdjacencyToFile(filename);
+      char filename[MAXSTRINGSIZE];
+      strcpy(filename, FRAGMENTPREFIX);
+      strcat(filename, ADJACENCYFILE);
+      StoreAdjacencyToFile(configuration->configpath, filename);
     }
 
     // store Hydrogen saturation correction file
-    BondFragments->AddHydrogenCorrection(prefix);
+    BondFragments->AddHydrogenCorrection(configuration->configpath);
 
     // store adaptive orders into file
-    StoreOrderAtSiteFile(prefix);
+    StoreOrderAtSiteFile(configuration->configpath);
 
     // restore orbital and Stop values
-    //CalculateOrbitals(*configuration);
+    CalculateOrbitals(*configuration);
 
     // free memory for bond part
@@ -775 +782 @@
 /** Stores pairs (Atom::nr, Atom::AdaptiveOrder) into file.
  * Atoms not present in the file get "-1".
- * \param &path path to file ORDERATSITEFILE
+ * \param *out output stream for debugging
+ * \param *path path to file ORDERATSITEFILE
  * \return true - file writable, false - not writable
  */
-bool molecule::StoreOrderAtSiteFile(std::string &path)
-{
-  string line;
+bool molecule::StoreOrderAtSiteFile(char *path)
+{
+  stringstream line;
   ofstream file;
 
-  line = path + ORDERATSITEFILE;
-  file.open(line.c_str());
+  line << path << "/" << FRAGMENTPREFIX << ORDERATSITEFILE;
+  file.open(line.str().c_str());
   DoLog(1) && (Log() << Verbose(1) << "Writing OrderAtSite " << ORDERATSITEFILE << " ... " << endl);
-  if (file.good()) {
+  if (file != NULL) {
     ActOnAllAtoms( &atom::OutputOrder, &file );
     file.close();
@@ -792 +800 @@
     return true;
   } else {
-    DoLog(1) && (Log() << Verbose(1) << "failed to open file " << line << "." << endl);
+    DoLog(1) && (Log() << Verbose(1) << "failed to open file " << line.str() << "." << endl);
     return false;
   }
@@ -799 +807 @@
 /** Parses pairs(Atom::nr, Atom::AdaptiveOrder) from file and stores in molecule's Atom's.
  * Atoms not present in the file get "0".
- * \param &path path to file ORDERATSITEFILEe
+ * \param *out output stream for debugging
+ * \param *path path to file ORDERATSITEFILEe
  * \return true - file found and scanned, false - file not found
  * \sa ParseKeySetFile() and CheckAdjacencyFileAgainstMolecule() as this is meant to be used in conjunction with the two
  */
-bool molecule::ParseOrderAtSiteFromFile(std::string &path)
+bool molecule::ParseOrderAtSiteFromFile(char *path)
 {
   unsigned char *OrderArray = new unsigned char[getAtomCount()];
@@ -809 +818 @@
   bool status;
   int AtomNr, value;
-  string line;
+  stringstream line;
   ifstream file;
 
@@ -818 +827 @@
 
   DoLog(1) && (Log() << Verbose(1) << "Begin of ParseOrderAtSiteFromFile" << endl);
-  line = path + ORDERATSITEFILE;
-  file.open(line.c_str());
-  if (file.good()) {
+  line << path << "/" << FRAGMENTPREFIX << ORDERATSITEFILE;
+  file.open(line.str().c_str());
+  if (file != NULL) {
     while (!file.eof()) { // parse from file
       AtomNr = -1;
@@ -841 +850 @@
     status = true;
   } else {
-    DoLog(1) && (Log() << Verbose(1) << "\t ... failed to open file " << line << "." << endl);
+    DoLog(1) && (Log() << Verbose(1) << "\t ... failed to open file " << line.str() << "." << endl);
     status = false;
   }