Changes in src/molecule_fragmentation.cpp [c743f8:83f176]
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/molecule_fragmentation.cpp
rc743f8 r83f176 1 /* 2 * Project: MoleCuilder 3 * Description: creates and alters molecular systems 4 * Copyright (C) 2010 University of Bonn. All rights reserved. 5 * Please see the LICENSE file or "Copyright notice" in builder.cpp for details. 6 */ 7 1 8 /* 2 9 * molecule_fragmentation.cpp … … 456 463 for (molecule::const_iterator iter = begin(); iter != end(); ++iter) { 457 464 #ifdef ADDHYDROGEN 458 if ((*iter)->getType()-> Z!= 1) // skip hydrogen465 if ((*iter)->getType()->getAtomicNumber() != 1) // skip hydrogen 459 466 #endif 460 467 { … … 476 483 for(molecule::const_iterator iter = begin(); iter != end(); ++iter) { 477 484 #ifdef ADDHYDROGEN 478 if ((*iter)->getType()-> Z!= 1) // skip hydrogen485 if ((*iter)->getType()->getAtomicNumber() != 1) // skip hydrogen 479 486 #endif 480 487 { … … 803 810 DoLog(1) && (Log() << Verbose(1) << "Writing OrderAtSite " << ORDERATSITEFILE << " ... " << endl); 804 811 if (file.good()) { 805 for_each(atoms.begin(),atoms.end(),bind2nd(mem_fun(&atom::OutputOrder), &file));812 ActOnAllAtoms( &atom::OutputOrder, &file ); 806 813 file.close(); 807 814 DoLog(1) && (Log() << Verbose(1) << "done." << endl); … … 885 892 for (KeySet::iterator runner = Leaf->begin(); runner != Leaf->end(); runner++) { 886 893 Runner = FindAtom((*runner)); 887 if (Runner->getType()-> Z!= 1) { // skip all those added hydrogens when re-filling snake stack894 if (Runner->getType()->getAtomicNumber() != 1) { // skip all those added hydrogens when re-filling snake stack 888 895 if (ShortestPathList[(*runner)] > SP) { // remove the oldest one with longest shortest path 889 896 SP = ShortestPathList[(*runner)]; … … 974 981 ++iter; 975 982 #ifdef ADDHYDROGEN 976 while ((iter != Leaf->end()) && ((*iter)->getType()-> Z== 1)){ // skip added hydrogen983 while ((iter != Leaf->end()) && ((*iter)->getType()->getAtomicNumber() == 1)){ // skip added hydrogen 977 984 iter++; 978 985 } … … 1352 1359 if ((RestrictedKeySet.find(OtherWalker->nr) != RestrictedKeySet.end()) 1353 1360 #ifdef ADDHYDROGEN 1354 && (OtherWalker->getType()-> Z!= 1)1361 && (OtherWalker->getType()->getAtomicNumber() != 1) 1355 1362 #endif 1356 1363 ) { // skip hydrogens and restrict to fragment
Note:
See TracChangeset
for help on using the changeset viewer.