72 |
|
dslAmat = 4, |
73 |
|
dslAngularMomentum = 8, |
74 |
|
dslElectroFrame = 16, |
75 |
< |
dslZAngle = 32, |
76 |
< |
dslForce = 64, |
77 |
< |
dslTorque = 128, |
78 |
< |
dslParticlePot = 256, |
79 |
< |
dslDensity = 512, |
80 |
< |
dslFunctional = 1024, |
81 |
< |
dslFunctionalDerivative = 2048, |
82 |
< |
dslElectricField = 4096, |
83 |
< |
dslSkippedCharge = 8192, |
84 |
< |
dslFlucQPosition = 16384, |
85 |
< |
dslFlucQVelocity = 32768, |
86 |
< |
dslFlucQForce = 65536 |
75 |
> |
dslForce = 32, |
76 |
> |
dslTorque = 64, |
77 |
> |
dslParticlePot = 128, |
78 |
> |
dslDensity = 256, |
79 |
> |
dslFunctional = 512, |
80 |
> |
dslFunctionalDerivative = 1024, |
81 |
> |
dslElectricField = 2048, |
82 |
> |
dslSkippedCharge = 4096, |
83 |
> |
dslFlucQPosition = 8192, |
84 |
> |
dslFlucQVelocity = 16384, |
85 |
> |
dslFlucQForce = 32768 |
86 |
|
}; |
87 |
|
|
88 |
|
DataStorage(); |
127 |
|
vector<RotMat3x3d> aMat; /** rotation matrix array */ |
128 |
|
vector<Vector3d> angularMomentum; /** angular momentum array (body-fixed) */ |
129 |
|
vector<Mat3x3d> electroFrame; /** the lab frame unit vector array*/ |
131 |
– |
vector<RealType> zAngle; /** z-angle array */ |
130 |
|
vector<Vector3d> force; /** force array */ |
131 |
|
vector<Vector3d> torque; /** torque array */ |
132 |
|
vector<RealType> particlePot; /** particle potential arrray */ |