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

Comparing trunk/src/applications/staticProps/GofR.cpp (file contents):
Revision 507 by gezelter, Fri Apr 15 22:04:00 2005 UTC vs.
Revision 963 by tim, Wed May 17 21:51:42 2006 UTC

# Line 46 | Line 46 | namespace oopse {
46  
47   namespace oopse {
48  
49 <  GofR::GofR(SimInfo* info, const std::string& filename, const std::string& sele1, const std::string& sele2, double len, int nrbins)
49 >  GofR::GofR(SimInfo* info, const std::string& filename, const std::string& sele1, const std::string& sele2, RealType len, int nrbins)
50      : RadialDistrFunc(info, filename, sele1, sele2), len_(len), nRBins_(nrbins){
51  
52        deltaR_ = len_ /nRBins_;
# Line 70 | Line 70 | namespace oopse {
70    void GofR::processHistogram() {
71  
72      int nPairs = getNPairs();
73 <    double volume = info_->getSnapshotManager()->getCurrentSnapshot()->getVolume();
74 <    double pairDensity = nPairs /volume * 2.0;
75 <    double pairConstant = ( 4.0 * NumericConstant::PI * pairDensity ) / 3.0;
73 >    RealType volume = info_->getSnapshotManager()->getCurrentSnapshot()->getVolume();
74 >    RealType pairDensity = nPairs /volume * 2.0;
75 >    RealType pairConstant = ( 4.0 * NumericConstant::PI * pairDensity ) / 3.0;
76  
77      for(int i = 0 ; i < histogram_.size(); ++i){
78  
79 <      double rLower = i * deltaR_;
80 <      double rUpper = rLower + deltaR_;
81 <      double volSlice = ( rUpper * rUpper * rUpper ) - ( rLower * rLower * rLower );
82 <      double nIdeal = volSlice * pairConstant;
79 >      RealType rLower = i * deltaR_;
80 >      RealType rUpper = rLower + deltaR_;
81 >      RealType volSlice = ( rUpper * rUpper * rUpper ) - ( rLower * rLower * rLower );
82 >      RealType nIdeal = volSlice * pairConstant;
83  
84        avgGofr_[i] += histogram_[i] / nIdeal;    
85      }
# Line 97 | Line 97 | namespace oopse {
97      Vector3d r12 = pos2 - pos1;
98      currentSnapshot_->wrapVector(r12);
99  
100 <    double distance = r12.length();
100 >    RealType distance = r12.length();
101  
102      if (distance < len_) {
103        int whichBin = distance / deltaR_;
# Line 114 | Line 114 | namespace oopse {
114        rdfStream << "selection2: (" << selectionScript2_ << ")\n";
115        rdfStream << "#r\tcorrValue\n";
116        for (int i = 0; i < avgGofr_.size(); ++i) {
117 <        double r = deltaR_ * (i + 0.5);
117 >        RealType r = deltaR_ * (i + 0.5);
118          rdfStream << r << "\t" << avgGofr_[i]/nProcessed_ << "\n";
119        }
120          

Diff Legend

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