Changeset f5d7e1

Timestamp:

Oct 17, 2008, 1:24:10 PM (17 years ago)

Author:

Frederik Heber <heber@…>

Children:

59f86c

Parents:

fa2d7d

Message:

molecule::PickLocalBackEdges(): BUGFIX: if Walker as not in subgraph, next candidate wasn't popped from stack due to wrong brackets of if clause.

This caused to loop to run indefinitely (hence, the new output messages for the functions in molecule::FragmentMolecule(). It is fixed.

File:

• molecuilder/src/molecules.cpp (modified) (3 diffs)

molecuilder/src/molecules.cpp

@@ -2942 +2942 @@
   while (MolecularWalker->next != NULL) {
     MolecularWalker = MolecularWalker->next;
-    *out << Verbose(0) << "Analysing the cycles of subgraph " << MolecularWalker->Leaf << " with nr. " << FragmentCounter << "." << endl;
     LocalBackEdgeStack = new StackClass<bond *> (MolecularWalker->Leaf->BondCount);
 //    // check the list of local atoms for debugging
@@ -2951 +2950 @@
 //      else
 //        *out << "\t" << ListOfLocalAtoms[FragmentCounter][i]->Name;
+    *out << Verbose(0) << "Gathering local back edges for subgraph " << MolecularWalker->Leaf << " with nr. " << FragmentCounter << "." << endl;
     MolecularWalker->Leaf->PickLocalBackEdges(out, ListOfLocalAtoms[FragmentCounter++], BackEdgeStack, LocalBackEdgeStack);
+    *out << Verbose(0) << "Analysing the cycles of subgraph " << MolecularWalker->Leaf << " with nr. " << FragmentCounter << "." << endl;
     MolecularWalker->Leaf->CyclicStructureAnalysis(out, LocalBackEdgeStack, MinimumRingSize);
+    *out << Verbose(0) << "Done with Analysing the cycles of subgraph " << MolecularWalker->Leaf << " with nr. " << FragmentCounter << "." << endl;
     delete(LocalBackEdgeStack);
   }
@@ -3094 +3096 @@
   do {  // go through all bonds and push local ones
     Walker = ListOfLocalAtoms[Binder->leftatom->nr];  // get one atom in the reference molecule
-    if (Walker == NULL) // if this Walker exists in the subgraph ... 
-      continue;
-    for(int i=0;i<NumberOfBondsPerAtom[Walker->nr];i++) {    // go through the local list of bonds
-      OtherAtom = ListOfBondsPerAtom[Walker->nr][i]->GetOtherAtom(Walker);
-      if (OtherAtom == ListOfLocalAtoms[Binder->rightatom->nr]) { // found the bond
-        LocalStack->Push(ListOfBondsPerAtom[Walker->nr][i]);
-        break;
-      }
-    }
+    if (Walker != NULL) // if this Walker exists in the subgraph ... 
+        for(int i=0;i<NumberOfBondsPerAtom[Walker->nr];i++) {    // go through the local list of bonds
+        OtherAtom = ListOfBondsPerAtom[Walker->nr][i]->GetOtherAtom(Walker);
+              if (OtherAtom == ListOfLocalAtoms[Binder->rightatom->nr]) { // found the bond
+              LocalStack->Push(ListOfBondsPerAtom[Walker->nr][i]);
+                                        *out << Verbose(3) << "Found local edge " << *(ListOfBondsPerAtom[Walker->nr][i]) << "." << endl;
+                break;
+            }
+        }
     Binder = ReferenceStack->PopFirst();  // loop the stack for next item
+                *out << Verbose(3) << "Current candidate edge " << Binder << "." << endl;
     ReferenceStack->Push(Binder);
   } while (FirstBond != Binder);