ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/OpenMD/branches/development/src/selection/HullFinder.cpp
(Generate patch)

Comparing branches/development/src/selection/HullFinder.cpp (file contents):
Revision 1803 by gezelter, Wed Oct 3 14:20:07 2012 UTC vs.
Revision 1866 by gezelter, Thu Apr 25 14:32:56 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 97 | Line 97 | namespace OpenMD {
97      sprintf( painCave.errMsg,
98               "HullFinder : Hull calculation is not possible without libqhull.\n"
99               "\tPlease rebuild OpenMD with qhull enabled.");
100 +    painCave.severity = OPENMD_ERROR;
101 +    painCave.isFatal = 1;
102 +    simError();
103 + #endif
104 +    
105 +    std::vector<Triangle> sMesh = surfaceMesh_->getMesh();
106 +    int nTriangles = sMesh.size();
107 +    // Loop over the mesh faces
108 +    std::vector<Triangle>::iterator face;
109 +    std::vector<StuntDouble*>::iterator vertex;
110 +
111 +    // This will work in parallel because the triangles returned by the mesh
112 +    // have a NULL stuntDouble if this processor doesn't own the
113 +    
114 +    for (face = sMesh.begin(); face != sMesh.end(); ++face) {
115 +      Triangle thisTriangle = *face;
116 +      std::vector<StuntDouble*> vertexSDs = thisTriangle.getVertices();
117 +      for (vertex = vertexSDs.begin(); vertex != vertexSDs.end(); ++vertex) {
118 +        if ((*vertex) != NULL) {
119 +          bsResult.setBitOn((*vertex)->getGlobalIndex());          
120 +        }
121 +      }
122 +    }
123 +    return bsResult;
124 +  }
125 +
126 +  OpenMDBitSet HullFinder::findHull(int frame) {
127 +    Snapshot* currSnapshot = info_->getSnapshotManager()->getSnapshot(frame);
128 +    OpenMDBitSet bsResult(nStuntDoubles_);
129 + #ifdef HAVE_QHULL
130 +    surfaceMesh_->computeHull(localSites_);
131 + #else
132 +    sprintf( painCave.errMsg,
133 +             "HullFinder : Hull calculation is not possible without libqhull.\n"
134 +             "\tPlease rebuild OpenMD with qhull enabled.");
135        painCave.severity = OPENMD_ERROR;
136        painCave.isFatal = 1;
137        simError();

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines