1 | /*
|
---|
2 | * atom_trajectoryparticle.hpp
|
---|
3 | *
|
---|
4 | * Created on: Oct 19, 2009
|
---|
5 | * Author: heber
|
---|
6 | */
|
---|
7 |
|
---|
8 | #ifndef ATOM_TRAJECTORYPARTICLE_HPP_
|
---|
9 | #define ATOM_TRAJECTORYPARTICLE_HPP_
|
---|
10 |
|
---|
11 | using namespace std;
|
---|
12 |
|
---|
13 | /*********************************************** includes ***********************************/
|
---|
14 |
|
---|
15 | // include config.h
|
---|
16 | #ifdef HAVE_CONFIG_H
|
---|
17 | #include <config.h>
|
---|
18 | #endif
|
---|
19 |
|
---|
20 | #include <fstream>
|
---|
21 |
|
---|
22 | #include <gsl/gsl_inline.h>
|
---|
23 | #include <gsl/gsl_randist.h>
|
---|
24 |
|
---|
25 | #include "atom_atominfo.hpp"
|
---|
26 | #include "atom_particleinfo.hpp"
|
---|
27 | #include "atom_trajectoryparticleinfo.hpp"
|
---|
28 |
|
---|
29 | /****************************************** forward declarations *****************************/
|
---|
30 |
|
---|
31 | class atom;
|
---|
32 | class config;
|
---|
33 | class ForceMatrix;
|
---|
34 | class Vector;
|
---|
35 |
|
---|
36 | /********************************************** declarations *******************************/
|
---|
37 |
|
---|
38 |
|
---|
39 | class TrajectoryParticle : public TrajectoryParticleInfo, public virtual AtomInfo, public virtual ParticleInfo {
|
---|
40 | public:
|
---|
41 |
|
---|
42 | TrajectoryParticle();
|
---|
43 | virtual ~TrajectoryParticle();
|
---|
44 |
|
---|
45 | // constraint potential and dynamics stuff
|
---|
46 | void AddKineticToTemperature(double *temperature, int step) const;
|
---|
47 | void EvaluateConstrainedForce(int startstep, int endstep, atom **PermutationMap, ForceMatrix *Force) const;
|
---|
48 | void CorrectVelocity(double *ActualTemp, int Step, Vector *CoGVelocity);
|
---|
49 |
|
---|
50 | // trajectory stuff
|
---|
51 | void ResizeTrajectory(int MaxSteps);
|
---|
52 | void CopyStepOnStep(int dest, int src);
|
---|
53 | void VelocityVerletUpdate(int MDSteps, config *configuration, ForceMatrix *Force, const size_t offset);
|
---|
54 | void SumUpKineticEnergy( int Step, double *TotalMass, Vector *TotalVelocity ) const;
|
---|
55 |
|
---|
56 | // thermostats
|
---|
57 | void Thermostat_Woodcock(double ScaleTempFactor, int Step, double *ekin);
|
---|
58 | void Thermostat_Gaussian_init(int Step, double *G, double *E);
|
---|
59 | void Thermostat_Gaussian_least_constraint(int Step, double G_over_E, double *ekin, config *configuration);
|
---|
60 | void Thermostat_Langevin(int Step, gsl_rng * r, double *ekin, config *configuration);
|
---|
61 | void Thermostat_Berendsen(int Step, double ScaleTempFactor, double *ekin, config *configuration);
|
---|
62 | void Thermostat_NoseHoover_init(int Step, double *delta_alpha);
|
---|
63 | void Thermostat_NoseHoover_scale(int Step, double *ekin, config *configuration);
|
---|
64 |
|
---|
65 | std::ostream & operator << (std::ostream &ost) const;
|
---|
66 |
|
---|
67 | private:
|
---|
68 | };
|
---|
69 |
|
---|
70 | ostream & operator << (ostream &ost, const TrajectoryParticle &a);
|
---|
71 |
|
---|
72 | #endif /* ATOM_TRAJECTORYPARTICLE_HPP_ */
|
---|