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

Comparing trunk/src/applications/staticProps/TetrahedralityParamZ.cpp (file contents):
Revision 1846 by gezelter, Thu Jan 31 15:55:47 2013 UTC vs.
Revision 2071 by gezelter, Sat Mar 7 21:41:51 2015 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   * [6]  Kuang & Gezelter, Mol. Phys., 110, 691-701 (2012).
# Line 57 | Line 57 | namespace OpenMD {
57                                               const std::string& sele1,
58                                               const std::string& sele2,
59                                               double rCut, int nzbins)
60 <    : StaticAnalyser(info, filename), selectionScript1_(sele1),
61 <      evaluator1_(info), seleMan1_(info), selectionScript2_(sele2),
62 <      evaluator2_(info), seleMan2_(info), nZBins_(nzbins) {
60 >    : StaticAnalyser(info, filename),
61 >      selectionScript1_(sele1), selectionScript2_(sele2),
62 >      seleMan1_(info), seleMan2_(info), evaluator1_(info), evaluator2_(info),
63 >      nZBins_(nzbins) {
64      
65      evaluator1_.loadScriptString(sele1);
66      if (!evaluator1_.isDynamic()) {
# Line 235 | Line 236 | namespace OpenMD {
236        qZstream << "#z\tQk\n";
237        for (unsigned int i = 0; i < sliceQ_.size(); ++i) {
238          RealType z = zAve * (i+0.5) / sliceQ_.size();
239 <        qZstream << z << "\t" << sliceQ_[i] / sliceCount_[i] << "\n";
239 >        if (sliceCount_[i] != 0) {
240 >          qZstream << z << "\t" << sliceQ_[i] / sliceCount_[i] << "\n";
241 >        }
242        }
243        
244      } else {      

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines