- Timestamp:
- Nov 3, 2011, 7:44:01 PM (14 years ago)
- Branches:
- Action_Thermostats, Add_AtomRandomPerturbation, Add_FitFragmentPartialChargesAction, Add_RotateAroundBondAction, Add_SelectAtomByNameAction, Added_ParseSaveFragmentResults, AddingActions_SaveParseParticleParameters, Adding_Graph_to_ChangeBondActions, Adding_MD_integration_tests, Adding_ParticleName_to_Atom, Adding_StructOpt_integration_tests, AtomFragments, Automaking_mpqc_open, AutomationFragmentation_failures, Candidate_v1.5.4, Candidate_v1.6.0, Candidate_v1.6.1, ChangeBugEmailaddress, ChangingTestPorts, ChemicalSpaceEvaluator, CombiningParticlePotentialParsing, Combining_Subpackages, Debian_Package_split, Debian_package_split_molecuildergui_only, Disabling_MemDebug, Docu_Python_wait, EmpiricalPotential_contain_HomologyGraph, EmpiricalPotential_contain_HomologyGraph_documentation, Enable_parallel_make_install, Enhance_userguide, Enhanced_StructuralOptimization, Enhanced_StructuralOptimization_continued, Example_ManyWaysToTranslateAtom, Exclude_Hydrogens_annealWithBondGraph, FitPartialCharges_GlobalError, Fix_BoundInBox_CenterInBox_MoleculeActions, Fix_ChargeSampling_PBC, Fix_ChronosMutex, Fix_FitPartialCharges, Fix_FitPotential_needs_atomicnumbers, Fix_ForceAnnealing, Fix_IndependentFragmentGrids, Fix_ParseParticles, Fix_ParseParticles_split_forward_backward_Actions, Fix_PopActions, Fix_QtFragmentList_sorted_selection, Fix_Restrictedkeyset_FragmentMolecule, Fix_StatusMsg, Fix_StepWorldTime_single_argument, Fix_Verbose_Codepatterns, Fix_fitting_potentials, Fixes, ForceAnnealing_goodresults, ForceAnnealing_oldresults, ForceAnnealing_tocheck, ForceAnnealing_with_BondGraph, ForceAnnealing_with_BondGraph_continued, ForceAnnealing_with_BondGraph_continued_betteresults, ForceAnnealing_with_BondGraph_contraction-expansion, FragmentAction_writes_AtomFragments, FragmentMolecule_checks_bonddegrees, GeometryObjects, Gui_Fixes, Gui_displays_atomic_force_velocity, ImplicitCharges, IndependentFragmentGrids, IndependentFragmentGrids_IndividualZeroInstances, IndependentFragmentGrids_IntegrationTest, IndependentFragmentGrids_Sole_NN_Calculation, JobMarket_RobustOnKillsSegFaults, JobMarket_StableWorkerPool, JobMarket_unresolvable_hostname_fix, MoreRobust_FragmentAutomation, ODR_violation_mpqc_open, PartialCharges_OrthogonalSummation, PdbParser_setsAtomName, PythonUI_with_named_parameters, QtGui_reactivate_TimeChanged_changes, Recreated_GuiChecks, Rewrite_FitPartialCharges, RotateToPrincipalAxisSystem_UndoRedo, SaturateAtoms_findBestMatching, SaturateAtoms_singleDegree, StoppableMakroAction, Subpackage_CodePatterns, Subpackage_JobMarket, Subpackage_LinearAlgebra, Subpackage_levmar, Subpackage_mpqc_open, Subpackage_vmg, Switchable_LogView, ThirdParty_MPQC_rebuilt_buildsystem, TrajectoryDependenant_MaxOrder, TremoloParser_IncreasedPrecision, TremoloParser_MultipleTimesteps, TremoloParser_setsAtomName, Ubuntu_1604_changes, stable
- Children:
- 41a467
- Parents:
- 50e4e5
- git-author:
- Frederik Heber <heber@…> (10/27/11 11:53:58)
- git-committer:
- Frederik Heber <heber@…> (11/03/11 19:44:01)
- Location:
- src/Graph
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
src/Graph/BondGraph.cpp
r50e4e5 r47d041 81 81 LOG(1, "Parsing bond length matrix successful."); 82 82 } else { 83 DoeLog(1) && (eLog()<< Verbose(1) << "Parsing bond length matrix failed." << endl);83 ELOG(1, "Parsing bond length matrix failed."); 84 84 status = false; 85 85 } -
src/Graph/BondGraph.hpp
r50e4e5 r47d041 395 395 AtomSetMixin<container_type,iterator_type,const_iterator_type> &Set) const 396 396 { 397 // DoLog(1) && (Log() << Verbose(1) << "Correcting Bond degree of each bond ... " << endl);397 //LOG(1, "Correcting Bond degree of each bond ... "); 398 398 int No = 0, OldNo = -1; 399 399 do { … … 404 404 } 405 405 } while (OldNo != No); 406 // DoLog(0) && (Log() << Verbose(0) << " done." << endl);406 //LOG(0, " done."); 407 407 return No; 408 408 } -
src/Graph/BreadthFirstSearchAdd.cpp
r50e4e5 r47d041 21 21 22 22 #include "BreadthFirstSearchAdd.hpp" 23 24 #include <sstream> 23 25 24 26 #include "atom.hpp" … … 70 72 PredecessorList[OtherAtom->getNr()] = Walker; // Walker is the predecessor 71 73 ShortestPathList[OtherAtom->getNr()] = ShortestPathList[Walker->getNr()] + 1; 72 DoLog(2) && (Log() << Verbose(2) << "Coloring OtherAtom " << OtherAtom->getName() << " " << GraphEdge::getColorName(ColorList[OtherAtom->getNr()]) << ", its predecessor is " << Walker->getName() << " and its Shortest Path is " << ShortestPathList[OtherAtom->getNr()] << " egde(s) long." << endl);74 LOG(2, "Coloring OtherAtom " << OtherAtom->getName() << " " << GraphEdge::getColorName(ColorList[OtherAtom->getNr()]) << ", its predecessor is " << Walker->getName() << " and its Shortest Path is " << ShortestPathList[OtherAtom->getNr()] << " egde(s) long."); 73 75 if ((((ShortestPathList[OtherAtom->getNr()] < BondOrder) && (Binder != Bond)))) { // Check for maximum distance 74 DoLog(3) && (Log() << Verbose(3));76 std::stringstream output; 75 77 if (AddedAtomList[OtherAtom->getNr()] == NULL) { // add if it's not been so far 76 78 AddedAtomList[OtherAtom->getNr()] = Mol->AddCopyAtom(OtherAtom); 77 DoLog(0) && (Log() << Verbose(0) << "Added OtherAtom " << OtherAtom->getName());79 output << "Added OtherAtom " << OtherAtom->getName(); 78 80 AddedBondList[Binder] = Mol->CopyBond(AddedAtomList[Walker->getNr()], AddedAtomList[OtherAtom->getNr()], Binder); 79 DoLog(0) && (Log() << Verbose(0) << " and bond " << *(AddedBondList[Binder]) << ", ");81 output << " and bond " << *(AddedBondList[Binder]) << ", "; 80 82 } else { // this code should actually never come into play (all GraphEdge::white atoms are not yet present in BondMolecule, that's why they are GraphEdge::white in the first place) 81 DoLog(0) && (Log() << Verbose(0) << "Not adding OtherAtom " << OtherAtom->getName());83 output << "Not adding OtherAtom " << OtherAtom->getName(); 82 84 if (AddedBondList[Binder] == NULL) { 83 85 AddedBondList[Binder] = Mol->CopyBond(AddedAtomList[Walker->getNr()], AddedAtomList[OtherAtom->getNr()], Binder); 84 DoLog(0) && (Log() << Verbose(0) << ", added Bond " << *(AddedBondList[Binder]));86 output << ", added Bond " << *(AddedBondList[Binder]); 85 87 } else 86 DoLog(0) && (Log() << Verbose(0) << ", not added Bond ");88 output << ", not added Bond "; 87 89 } 88 DoLog(0) && (Log() << Verbose(0) << ", putting OtherAtom into queue." << endl); 90 output << ", putting OtherAtom into queue."; 91 LOG(0, output.str()); 89 92 BFSStack.push_front(OtherAtom); 90 93 } else { // out of bond order, then replace 91 94 if ((AddedAtomList[OtherAtom->getNr()] == NULL) && (Binder->Cyclic)) 92 95 ColorList[OtherAtom->getNr()] = GraphEdge::white; // unmark if it has not been queued/added, to make it available via its other bonds (cyclic) 93 if (Binder == Bond) 94 DoLog(3) && (Log() << Verbose(3) << "Not Queueing, is the Root bond"); 95 else if (ShortestPathList[OtherAtom->getNr()] >= BondOrder) 96 DoLog(3) && (Log() << Verbose(3) << "Not Queueing, is out of Bond Count of " << BondOrder); 97 if (!Binder->Cyclic) 98 DoLog(0) && (Log() << Verbose(0) << ", is not part of a cyclic bond, saturating bond with Hydrogen." << endl); 96 { 97 std::stringstream output; 98 if (Binder == Bond) 99 output << "Not Queueing, is the Root bond"; 100 else if (ShortestPathList[OtherAtom->getNr()] >= BondOrder) 101 output << "Not Queueing, is out of Bond Count of " << BondOrder; 102 if (!Binder->Cyclic) 103 output << ", is not part of a cyclic bond, saturating bond with Hydrogen."; 104 LOG(3, output.str()); 105 } 99 106 if (AddedBondList[Binder] == NULL) { 100 107 if ((AddedAtomList[OtherAtom->getNr()] != NULL)) { // .. whether we add or saturate … … 112 119 void BreadthFirstSearchAdd::VisitedNode(molecule *Mol, atom *&Walker, atom *&OtherAtom, bond *&Binder, bond *&Bond) 113 120 { 114 DoLog(3) && (Log() << Verbose(3) << "Not Adding, has already been visited." << endl);121 LOG(3, "Not Adding, has already been visited."); 115 122 // This has to be a cyclic bond, check whether it's present ... 116 123 if (AddedBondList[Binder] == NULL) { … … 147 154 BFSStack.pop_front(); 148 155 const BondList& ListOfBonds = Walker->getListOfBonds(); 149 DoLog(1) && (Log() << Verbose(1) << "Current Walker is: " << Walker->getName() << ", and has " << ListOfBonds.size() << " bonds." << endl);156 LOG(1, "Current Walker is: " << Walker->getName() << ", and has " << ListOfBonds.size() << " bonds."); 150 157 for (BondList::const_iterator Runner = ListOfBonds.begin(); 151 158 Runner != ListOfBonds.end(); … … 154 161 Binder = (*Runner); 155 162 OtherAtom = (*Runner)->GetOtherAtom(Walker); 156 DoLog(2) && (Log() << Verbose(2) << "Current OtherAtom is: " << OtherAtom->getName() << " for bond " << *(*Runner) << "." << endl);163 LOG(2, "Current OtherAtom is: " << OtherAtom->getName() << " for bond " << *(*Runner) << "."); 157 164 if (ColorList[OtherAtom->getNr()] == GraphEdge::white) { 158 165 UnvisitedNode(Mol, Walker, OtherAtom, Binder, Bond); … … 163 170 } 164 171 ColorList[Walker->getNr()] = GraphEdge::black; 165 DoLog(1) && (Log() << Verbose(1) << "Coloring Walker " << Walker->getName() << " GraphEdge::black." << endl);172 LOG(1, "Coloring Walker " << Walker->getName() << " GraphEdge::black."); 166 173 } 167 174 } -
src/Graph/BuildInducedSubgraph.cpp
r50e4e5 r47d041 45 45 { 46 46 // fill parent list with sons 47 DoLog(3) && (Log() << Verbose(3) << "Filling Parent List." << endl);47 LOG(3, "Filling Parent List."); 48 48 for (molecule::const_iterator iter = Son->begin(); iter != Son->end(); ++iter) { 49 49 ParentList[(*iter)->father] = (*iter); -
src/Graph/CyclicStructureAnalysis.cpp
r50e4e5 r47d041 321 321 int MinRingSize = -1; 322 322 323 //Log() << Verbose(1) << "Back edge list - "; 324 //BackEdgeStack->Output(out); 323 //std::stringstream output; 324 //output << "Back edge list - "; 325 //BackEdgeStack->Output(output); 326 //LOG(0, output.str()); 325 327 326 328 LOG(1, "STATUS: Analysing cycles ... " << endl); -
src/Graph/DepthFirstSearchAnalysis.cpp
r50e4e5 r47d041 163 163 const bond *Binder = *BondRunner; 164 164 if (DoLog(2)) { 165 ostream &out = (Log() << Verbose(2)); 166 out << ((Binder->Type == GraphEdge::TreeEdge) ? "TreeEdge " : "BackEdge ") << *Binder << ": <"; 167 out << ((Binder->leftatom->SeparationVertex) ? "SP," : "") << "L" << Binder->leftatom->LowpointNr << " G" << Binder->leftatom->GraphNr << " Comp."; 168 Binder->leftatom->OutputComponentNumber(&out); 169 out << " === "; 170 out << ((Binder->rightatom->SeparationVertex) ? "SP," : "") << "L" << Binder->rightatom->LowpointNr << " G" << Binder->rightatom->GraphNr << " Comp."; 171 Binder->rightatom->OutputComponentNumber(&out); 172 out << ">." << endl; 165 std::stringstream output; 166 output << ((Binder->Type == GraphEdge::TreeEdge) ? "TreeEdge " : "BackEdge ") << *Binder << ": <"; 167 output << ((Binder->leftatom->SeparationVertex) ? "SP," : "") << "L" << Binder->leftatom->LowpointNr << " G" << Binder->leftatom->GraphNr << " Comp."; 168 Binder->leftatom->OutputComponentNumber(&output); 169 output << " === "; 170 output << ((Binder->rightatom->SeparationVertex) ? "SP," : "") << "L" << Binder->rightatom->LowpointNr << " G" << Binder->rightatom->GraphNr << " Comp."; 171 Binder->rightatom->OutputComponentNumber(&output); 172 output << ">."; 173 LOG(2, output.str()); 173 174 } 174 175 if (Binder->Cyclic) // cyclic ??
Note:
See TracChangeset
for help on using the changeset viewer.