ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/OpenMD/trunk/src/integrators/NPT.hpp
(Generate patch)

Comparing trunk/src/integrators/NPT.hpp (file contents):
Revision 506 by gezelter, Wed Jan 12 22:41:40 2005 UTC vs.
Revision 507 by gezelter, Fri Apr 15 22:04:00 2005 UTC

# Line 1 | Line 1
1 < /*
1 > /*
2   * Copyright (c) 2005 The University of Notre Dame. All Rights Reserved.
3   *
4   * The University of Notre Dame grants you ("Licensee") a
# Line 53 | Line 53 | namespace oopse {
53   #include "integrators/VelocityVerletIntegrator.hpp"
54  
55   namespace oopse {
56 <    class NPT : public VelocityVerletIntegrator {
57 <        public:
56 >  class NPT : public VelocityVerletIntegrator {
57 >  public:
58  
59 <            NPT(SimInfo * info);
60 <            virtual ~NPT();
59 >    NPT(SimInfo * info);
60 >    virtual ~NPT();
61  
62  
63 <            int getMaxIterationNumber() {
64 <                return maxIterNum_;
65 <            }
63 >    int getMaxIterationNumber() {
64 >      return maxIterNum_;
65 >    }
66  
67 <            void setMaxIterationNumber(int maxIter) {
68 <                maxIterNum_ = maxIter;
69 <            }
70 <            double getTauThermostat() {
71 <                return tauThermostat;
72 <            }
67 >    void setMaxIterationNumber(int maxIter) {
68 >      maxIterNum_ = maxIter;
69 >    }
70 >    double getTauThermostat() {
71 >      return tauThermostat;
72 >    }
73  
74 <            void setTauThermostat(double tt) {
75 <                tauThermostat = tt;
76 <            }
74 >    void setTauThermostat(double tt) {
75 >      tauThermostat = tt;
76 >    }
77  
78 <            double getTauBarostat() {
79 <                return tauBarostat;
80 <            }
81 <            void setTauBarostat(double tb) {
82 <                tauBarostat = tb;
83 <            }
78 >    double getTauBarostat() {
79 >      return tauBarostat;
80 >    }
81 >    void setTauBarostat(double tb) {
82 >      tauBarostat = tb;
83 >    }
84  
85 <            double getTargetTemp() {
86 <                return targetTemp;
87 <            }
85 >    double getTargetTemp() {
86 >      return targetTemp;
87 >    }
88              
89 <            void setTargetTemp(double tt) {
90 <                targetTemp = tt;
91 <            }
89 >    void setTargetTemp(double tt) {
90 >      targetTemp = tt;
91 >    }
92  
93 <            double getTargetPressure() {
94 <                return targetTemp;
95 <            }
93 >    double getTargetPressure() {
94 >      return targetTemp;
95 >    }
96              
97 <            void setTargetPressure(double tp) {
98 <                targetPressure = tp;
99 <            }
97 >    void setTargetPressure(double tp) {
98 >      targetPressure = tp;
99 >    }
100  
101 <            double getChiTolerance() {
102 <                return chiTolerance;
103 <            }
101 >    double getChiTolerance() {
102 >      return chiTolerance;
103 >    }
104              
105 <            void setChiTolerance(double tol) {
106 <                chiTolerance = tol;
107 <            }
105 >    void setChiTolerance(double tol) {
106 >      chiTolerance = tol;
107 >    }
108  
109 <            double getEtaTolerance() {
110 <                return etaTolerance;
111 <            }
109 >    double getEtaTolerance() {
110 >      return etaTolerance;
111 >    }
112              
113 <            void setEtaTolerance(double tol) {
114 <                etaTolerance = tol;
115 <            }
113 >    void setEtaTolerance(double tol) {
114 >      etaTolerance = tol;
115 >    }
116  
117 <        protected:
117 >  protected:
118  
119 <            virtual void integrateStep() {
120 <                needStress= true;
121 <                VelocityVerletIntegrator::integrateStep();
122 <            }
119 >    virtual void integrateStep() {
120 >      needStress= true;
121 >      VelocityVerletIntegrator::integrateStep();
122 >    }
123  
124 <            virtual void doUpdate();
124 >    virtual void doUpdate();
125              
126 <            double NkBT;
127 <            double fkBT;
126 >    double NkBT;
127 >    double fkBT;
128  
129 <            double tt2;            
130 <            double tb2;
129 >    double tt2;            
130 >    double tb2;
131              
132 <            double instaTemp;
133 <            double instaPress;
134 <            double instaVol;
132 >    double instaTemp;
133 >    double instaPress;
134 >    double instaVol;
135  
136  
137 <            // targetTemp, targetPressure, and tauBarostat must be set.
138 <            // One of qmass or tauThermostat must be set;
137 >    // targetTemp, targetPressure, and tauBarostat must be set.
138 >    // One of qmass or tauThermostat must be set;
139  
140 <            double targetTemp;
141 <            double targetPressure;
142 <            double tauThermostat;
143 <            double tauBarostat;
140 >    double targetTemp;
141 >    double targetPressure;
142 >    double tauThermostat;
143 >    double tauBarostat;
144  
145 <            std::vector<Vector3d> oldPos;
146 <            std::vector<Vector3d> oldVel;
147 <            std::vector<Vector3d> oldJi;
145 >    std::vector<Vector3d> oldPos;
146 >    std::vector<Vector3d> oldVel;
147 >    std::vector<Vector3d> oldJi;
148  
149 <            double etaTolerance;
149 >    double etaTolerance;
150        
151 <            double chi;
152 <            double integralOfChidt;
153 <            Mat3x3d press;
151 >    double chi;
152 >    double integralOfChidt;
153 >    Mat3x3d press;
154              
155 <        private:
155 >  private:
156  
157 <            virtual void moveA();
158 <            virtual void moveB();
157 >    virtual void moveA();
158 >    virtual void moveB();
159              
160 <            virtual void getVelScaleA(Vector3d& sc, const Vector3d& vel) = 0;
160 >    virtual void getVelScaleA(Vector3d& sc, const Vector3d& vel) = 0;
161              
162 <            virtual void getVelScaleB(Vector3d& sc, int index) = 0;
162 >    virtual void getVelScaleB(Vector3d& sc, int index) = 0;
163              
164 <            virtual void getPosScale(const Vector3d& pos, const Vector3d& COM,
165 <                                                            int index, Vector3d& sc) = 0;
164 >    virtual void getPosScale(const Vector3d& pos, const Vector3d& COM,
165 >                             int index, Vector3d& sc) = 0;
166  
167 <            virtual void calcVelScale() = 0;
167 >    virtual void calcVelScale() = 0;
168  
169 <            virtual bool etaConverged() = 0;
169 >    virtual bool etaConverged() = 0;
170  
171 <            virtual void evolveEtaA() = 0;
171 >    virtual void evolveEtaA() = 0;
172  
173 <            virtual void evolveEtaB() = 0;
173 >    virtual void evolveEtaB() = 0;
174  
175 <            virtual void scaleSimBox() = 0;
175 >    virtual void scaleSimBox() = 0;
176              
177 <            virtual double calcConservedQuantity() = 0;      
177 >    virtual double calcConservedQuantity() = 0;      
178  
179 <            virtual void loadEta() = 0;
180 <            virtual void saveEta() = 0;
179 >    virtual void loadEta() = 0;
180 >    virtual void saveEta() = 0;
181              
182 <            int maxIterNum_;
182 >    int maxIterNum_;
183  
184 <            double chiTolerance;    
185 <    };
184 >    double chiTolerance;    
185 >  };
186  
187   }      //end namespace oopse
188  

Diff Legend

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