ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/OpenMD/trunk/src/brains/SimInfo.cpp
(Generate patch)

Comparing trunk/src/brains/SimInfo.cpp (file contents):
Revision 1796 by gezelter, Mon Sep 10 18:38:44 2012 UTC vs.
Revision 1908 by gezelter, Fri Jul 19 21:25:45 2013 UTC

# Line 35 | Line 35
35   *                                                                      
36   * [1]  Meineke, et al., J. Comp. Chem. 26, 252-271 (2005).            
37   * [2]  Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006).          
38 < * [3]  Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008).          
38 > * [3]  Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008).          
39   * [4]  Kuang & Gezelter,  J. Chem. Phys. 133, 164101 (2010).
40   * [5]  Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011).
41   */
# Line 91 | Line 91 | namespace OpenMD {
91      for (vector<Component*>::iterator i = components.begin();
92           i !=components.end(); ++i) {
93        molStamp = (*i)->getMoleculeStamp();
94 +      if ( (*i)->haveRegion() ) {        
95 +        molStamp->setRegion( (*i)->getRegion() );
96 +      } else {
97 +        // set the region to a disallowed value:
98 +        molStamp->setRegion( -1 );
99 +      }
100 +
101        nMolWithSameStamp = (*i)->getNMol();
102        
103        addMoleculeStamp(molStamp, nMolWithSameStamp);
# Line 780 | Line 787 | namespace OpenMD {
787   #endif
788  
789      return atomTypes;        
790 +  }
791 +
792 +
793 +  int getGlobalCountOfType(AtomType* atype) {
794 +    /*
795 +    set<AtomType*> atypes = getSimulatedAtomTypes();
796 +    map<AtomType*, int> counts_;
797 +
798 +    for(mol = beginMolecule(mi); mol != NULL; mol = nextMolecule(mi)) {
799 +      for(atom = mol->beginAtom(ai); atom != NULL;
800 +          atom = mol->nextAtom(ai)) {
801 +        atom->getAtomType();
802 +      }      
803 +    }    
804 +    */
805 +    return 0;
806    }
807  
808    void SimInfo::setupSimVariables() {
# Line 925 | Line 948 | namespace OpenMD {
948        }
949      }    
950      
928    //scan topology
929
930    int* excludeList = excludedInteractions_.getPairList();
931    int* oneTwoList = oneTwoInteractions_.getPairList();
932    int* oneThreeList = oneThreeInteractions_.getPairList();
933    int* oneFourList = oneFourInteractions_.getPairList();
934
951      topologyDone_ = true;
952    }
953  
# Line 1003 | Line 1019 | namespace OpenMD {
1019    
1020    
1021    StuntDouble* SimInfo::getIOIndexToIntegrableObject(int index) {
1022 <    if (index >= IOIndexToIntegrableObject.size()) {
1022 >    if (index >= int(IOIndexToIntegrableObject.size())) {
1023        sprintf(painCave.errMsg,
1024                "SimInfo::getIOIndexToIntegrableObject Error: Integrable Object\n"
1025                "\tindex exceeds number of known objects!\n");

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines