Changeset 329d0f for molecuilder/src/molecules.cpp

Timestamp:

Jun 23, 2008, 11:23:50 AM (17 years ago)

Author:

Frederik Heber <heber@…>

Children:

d7b1bf

Parents:

7a40c8

Message:

StoreFragmentFromKeySet(): BUGFIX - commented out the Leaf->AddBond for non-hydrogen atoms by accident

This happened because the adding was hidden in the "*out << ... << endl;" directive. I put it onto an extra line and surrounded it with "*out <<" and so on so that it may still be used.

File:

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

molecuilder/src/molecules.cpp

@@ -3123 +3123 @@
 //          *out << ", whose son is " << *SonList[OtherFather->nr] << "." << endl; 
           if (OtherFather->nr > FatherOfRunner->nr) { // add bond (nr check is for adding only one of both variants: ab, ba)
-//            *out << Verbose(3) << "Adding Bond: " << Leaf->AddBond(Runner, SonList[OtherFather->nr], ListOfBondsPerAtom[FatherOfRunner->nr][i]->BondDegree) << "." << endl;
+//            *out << Verbose(3) << "Adding Bond: ";
+//            *out << 
+            Leaf->AddBond(Runner, SonList[OtherFather->nr], ListOfBondsPerAtom[FatherOfRunner->nr][i]->BondDegree);
+//            *out << "." << endl;
             //NumBonds[Runner->nr]++;
           } else { 
@@ -3743 +3746 @@
   bool flag = true;
 
-//  *out << Verbose(1) << "Begin of ScanForPeriodicCorrection." << endl;
+  *out << Verbose(1) << "Begin of ScanForPeriodicCorrection." << endl;
 
   ColorList = (enum Shading *) Malloc(sizeof(enum Shading)*AtomCount, "molecule::ScanForPeriodicCorrection: *ColorList");
@@ -3757 +3760 @@
       for (int i=NDIM;i--;) {
         tmp = fabs(Binder->leftatom->x.x[i] - Binder->rightatom->x.x[i]);
-        //*out << Verbose(3) << "Checking " << i << "th distance of " << *Binder->leftatom << " to " << *Binder->rightatom << ": " << tmp << "." << endl;
+        *out << Verbose(3) << "Checking " << i << "th distance of " << *Binder->leftatom << " to " << *Binder->rightatom << ": " << tmp << "." << endl;
         if (tmp > BondDistance) {
           OtherBinder = Binder->next; // note down binding partner for later re-insertion
           unlink(Binder);   // unlink bond
-//          *out << Verbose(2) << "Correcting at bond " << *Binder << "." << endl;
+          *out << Verbose(2) << "Correcting at bond " << *Binder << "." << endl;
           flag = true;
           break;
@@ -3775 +3778 @@
       }
       TranslationVector.MatrixMultiplication(matrix);
-      //*out << "Translation vector is ";
-      //TranslationVector.Output(out);
-      //*out << endl;
+      *out << "Translation vector is ";
+      TranslationVector.Output(out);
+      *out << endl;
       // apply to all atoms of first component via BFS
       for (int i=AtomCount;i--;)
@@ -3799 +3802 @@
       link(Binder, OtherBinder);
     } else {
-//      *out << Verbose(2) << "No corrections for this fragment." << endl;
+      *out << Verbose(2) << "No corrections for this fragment." << endl;
     }
     //delete(CompStack);
@@ -3808 +3811 @@
   Free((void **)&ColorList, "molecule::ScanForPeriodicCorrection: *ColorList");
   Free((void **)&matrix, "molecule::ScanForPeriodicCorrection: *matrix");
-//  *out << Verbose(1) << "End of ScanForPeriodicCorrection." << endl;
+  *out << Verbose(1) << "End of ScanForPeriodicCorrection." << endl;
 };