Changes in src/atom_bondedparticle.cpp [e138de:791138]
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/atom_bondedparticle.cpp
re138de r791138 44 44 void BondedParticle::OutputBondOfAtom() const 45 45 { 46 Log() << Verbose(4) << "Atom " << Name << "/" << nr << " with " << ListOfBonds.size() << " bonds: " ;46 Log() << Verbose(4) << "Atom " << Name << "/" << nr << " with " << ListOfBonds.size() << " bonds: " << endl; 47 47 int TotalDegree = 0; 48 48 for (BondList::const_iterator Runner = ListOfBonds.begin(); Runner != ListOfBonds.end(); ++Runner) { 49 Log() << Verbose( 0) << **Runner << "\t";49 Log() << Verbose(4) << **Runner << endl; 50 50 TotalDegree += (*Runner)->BondDegree; 51 51 } 52 Log() << Verbose( 0) << " -- TotalDegree: " << TotalDegree << endl;52 Log() << Verbose(4) << " -- TotalDegree: " << TotalDegree << endl; 53 53 }; 54 54 … … 75 75 status = true; 76 76 } else { 77 Log() << Verbose(1) << "ERROR: "<< *Binder << " does not contain " << *this << "." << endl;77 eLog() << Verbose(1) << *Binder << " does not contain " << *this << "." << endl; 78 78 } 79 79 } else { 80 Log() << Verbose(1) << "ERROR:Binder is " << Binder << "." << endl;80 eLog() << Verbose(1) << "Binder is " << Binder << "." << endl; 81 81 } 82 82 return status; … … 94 94 status = true; 95 95 } else { 96 Log() << Verbose(1) << "ERROR: "<< *Binder << " does not contain " << *this << "." << endl;96 eLog() << Verbose(1) << *Binder << " does not contain " << *this << "." << endl; 97 97 } 98 98 } else { 99 Log() << Verbose(1) << "ERROR:Binder is " << Binder << "." << endl;99 eLog() << Verbose(1) << "Binder is " << Binder << "." << endl; 100 100 } 101 101 return status; … … 121 121 bond *CandidateBond = NULL; 122 122 123 NoBonds = CountBonds(); 123 124 //Log() << Verbose(3) << "Walker " << *this << ": " << (int)this->type->NoValenceOrbitals << " > " << NoBonds << "?" << endl; 124 NoBonds = CountBonds();125 125 if ((int)(type->NoValenceOrbitals) > NoBonds) { // we have a mismatch, check all bonding partners for mismatch 126 126 for (BondList::const_iterator Runner = ListOfBonds.begin(); Runner != ListOfBonds.end(); (++Runner)) { 127 127 OtherWalker = (*Runner)->GetOtherAtom(this); 128 128 OtherNoBonds = OtherWalker->CountBonds(); 129 //Log() << Verbose(3) << "OtherWalker " << *OtherWalker << ": " << (int)OtherWalker->type->NoValenceOrbitals << " > " << NoBonds << "?" << endl;130 if ((int)(OtherWalker->type->NoValenceOrbitals) > NoBonds) { // check if possible candidate129 //Log() << Verbose(3) << "OtherWalker " << *OtherWalker << ": " << (int)OtherWalker->type->NoValenceOrbitals << " > " << OtherNoBonds << "?" << endl; 130 if ((int)(OtherWalker->type->NoValenceOrbitals) > OtherNoBonds) { // check if possible candidate 131 131 if ((CandidateBond == NULL) || (ListOfBonds.size() > OtherWalker->ListOfBonds.size())) { // pick the one with fewer number of bonds first 132 132 CandidateBond = (*Runner); … … 137 137 if ((CandidateBond != NULL)) { 138 138 CandidateBond->BondDegree++; 139 Log() << Verbose(2) << "Increased bond degree for bond " << *CandidateBond << "." << endl;139 //Log() << Verbose(2) << "Increased bond degree for bond " << *CandidateBond << "." << endl; 140 140 } else { 141 //Log() << Verbose(2) << "Could not find correct degree for atom " << *this << "." << endl;141 eLog() << Verbose(2) << "Could not find correct degree for atom " << *this << "." << endl; 142 142 FalseBondDegree++; 143 143 }
Note:
See TracChangeset
for help on using the changeset viewer.