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

Comparing trunk/src/brains/DataStorage.hpp (file contents):
Revision 323 by tim, Sun Feb 13 15:56:10 2005 UTC vs.
Revision 1245 by chuckv, Tue May 27 16:39:06 2008 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 54 | Line 54 | namespace oopse {
54   #include <math/SquareMatrix3.hpp>
55  
56   namespace oopse {
57 <    /**
58 <     * @class DataStorage
59 <     * @warning do not try to insert element into (or ease element from) private member data
60 <     * of DataStorage directly.
61 <     * @todo DataStorage may need refactorying. Every std::vector can inherit from the same base class
62 <     * which will make it easy to maintain
63 <     */
64 <    class DataStorage {
65 <        public:
57 >  /**
58 >   * @class DataStorage
59 >   * @warning do not try to insert element into (or ease element from) private member data
60 >   * of DataStorage directly.
61 >   * @todo DataStorage may need refactoring. Every std::vector can inherit from the same base class
62 >   * which will make it easy to maintain
63 >   */
64 >  class DataStorage {
65 >  public:
66  
67 <            enum{
68 <                dslPosition = 1,
69 <                dslVelocity = 2,
70 <                dslAmat = 4,
71 <                dslAngularMomentum = 8,
72 <                dslElectroFrame = 16,
73 <                dslZAngle = 32,
74 <                dslForce = 64,
75 <                dslTorque = 128
76 <            };
67 >    enum{
68 >      dslPosition = 1,
69 >      dslVelocity = 2,
70 >      dslAmat = 4,
71 >      dslAngularMomentum = 8,
72 >      dslElectroFrame = 16,
73 >      dslZAngle = 32,
74 >      dslForce = 64,
75 >      dslTorque = 128,
76 >      dslParticlePot = 256
77 >    };
78  
79  
80 <            DataStorage();
81 <            DataStorage(int size, int storageLayout = 255);
82 <            /** return the size of this DataStorage. */
83 <            int getSize();
84 <            /**
85 <             * Changes the size of this DataStorage.
86 <             * @param size new size of this DataStorage
87 <             */
88 <            void resize(int newSize);
89 <            /**
90 <             * Reallocates memory manually. The main reason for using reserve() is efficiency
91 <             * if you know the capacity to which your std::vector must eventually grow, then it is usually more
92 <             * efficient to allocate that memory all at once.
93 <             */
94 <            void reserve(int size);
95 <            /**
96 <             * Copies data inside DataStorage class.
97 <             * Copy function actually call std::copy for every std::vector in DataStorage class.
98 <             * One Precondition of std::copy is that target is not within the range [soruce, soruce + num]
99 <             * @param souce
100 <             * @param num number of element to be moved
101 <             * @param target
102 <             */
103 <            void copy(int source, int num, int target);
104 <            /** Returns the storage layout  */
105 <            int getStorageLayout();
106 <            /** Sets the storage layout  */
107 <            void setStorageLayout(int layout);
108 <            /** Returns the pointer of internal array */
109 <            double *getArrayPointer(int whichArray);
80 >    DataStorage();
81 >    DataStorage(int size, int storageLayout = 511);
82 >    /** return the size of this DataStorage. */
83 >    int getSize();
84 >    /**
85 >     * Changes the size of this DataStorage.
86 >     * @param size new size of this DataStorage
87 >     */
88 >    void resize(int newSize);
89 >    /**
90 >     * Reallocates memory manually. The main reason for using reserve() is efficiency
91 >     * if you know the capacity to which your std::vector must eventually grow, then it is usually more
92 >     * efficient to allocate that memory all at once.
93 >     */
94 >    void reserve(int size);
95 >    /**
96 >     * Copies data inside DataStorage class.
97 >     * Copy function actually call std::copy for every std::vector in DataStorage class.
98 >     * One Precondition of std::copy is that target is not within the range [soruce, soruce + num]
99 >     * @param souce
100 >     * @param num number of element to be moved
101 >     * @param target
102 >     */
103 >    void copy(int source, int num, int target);
104 >    /** Returns the storage layout  */
105 >    int getStorageLayout();
106 >    /** Sets the storage layout  */
107 >    void setStorageLayout(int layout);
108 >    /** Returns the pointer of internal array */
109 >    RealType *getArrayPointer(int whichArray);
110  
111 <            std::vector<Vector3d> position;               /** position array */
112 <            std::vector<Vector3d> velocity;               /** velocity array */
113 <            std::vector<RotMat3x3d> aMat;            /** rotation matrix array */
114 <            std::vector<Vector3d> angularMomentum;/** angular momentum array (body-fixed) */
115 <            std::vector<Mat3x3d> electroFrame;                /** the lab frame unit std::vector array*/
116 <            std::vector<double> zAngle;              /** z -angle array */        
117 <            std::vector<Vector3d> force;               /** force array */
118 <            std::vector<Vector3d> torque;               /** torque array */
111 >    std::vector<Vector3d> position;               /** position array */
112 >    std::vector<Vector3d> velocity;               /** velocity array */
113 >    std::vector<RotMat3x3d> aMat;            /** rotation matrix array */
114 >    std::vector<Vector3d> angularMomentum;/** angular momentum array (body-fixed) */
115 >    std::vector<Mat3x3d> electroFrame;                /** the lab frame unit std::vector array*/
116 >    std::vector<RealType> zAngle;              /** z -angle array */        
117 >    std::vector<Vector3d> force;               /** force array */
118 >    std::vector<Vector3d> torque;               /** torque array */
119 >    std::vector<RealType> particlePot;         /** pair potential arrray */
120  
121 <            static int getBytesPerStuntDouble(int layout);
121 >    static int getBytesPerStuntDouble(int layout);
122  
123 <        private:
123 >  private:
124  
125 <            double* internalGetArrayPointer(std::vector<Vector3d>& v);
125 >    RealType* internalGetArrayPointer(std::vector<Vector3d>& v);
126              
127 <            double* internalGetArrayPointer(std::vector<RotMat3x3d>& v);
128 <            double* internalGetArrayPointer(std::vector<double>& v);
127 >    RealType* internalGetArrayPointer(std::vector<RotMat3x3d>& v);
128 >    RealType* internalGetArrayPointer(std::vector<RealType>& v);
129 >
130              
131 <            template<typename T>
132 <            void internalResize(std::vector<T>& v, int newSize);
131 >    template<typename T>
132 >    void internalResize(std::vector<T>& v, int newSize);
133  
134 <            template<typename T>
135 <            void internalCopy(std::vector<T>& v, int source,  int num, int target);
134 >    template<typename T>
135 >    void internalCopy(std::vector<T>& v, int source,  int num, int target);
136              
137 <            int size_;
138 <            int storageLayout_;
137 >    int size_;
138 >    int storageLayout_;
139  
140 <    };
140 >  };
141  
142   }
143   #endif //BRAINS_DATASTORAGE_HPP

Diff Legend

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