34 |
|
* |
35 |
|
* [1] Meineke, et al., J. Comp. Chem. 26, 252-271 (2005). |
36 |
|
* [2] Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006). |
37 |
< |
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008). |
38 |
< |
* [4] Vardeman & Gezelter, in progress (2009). |
37 |
> |
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008). |
38 |
> |
* [4] Kuang & Gezelter, J. Chem. Phys. 133, 164101 (2010). |
39 |
> |
* [4] , Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011). * |
40 |
|
* |
40 |
– |
* |
41 |
|
* AlphaHull.hpp |
42 |
|
* |
43 |
|
* Purpose: To calculate alphahull, hull volume using the QuickHull algorithm provided by QHull. |
60 |
|
#include <vector> |
61 |
|
#include <string> |
62 |
|
|
63 |
< |
|
63 |
> |
using namespace std; |
64 |
|
namespace OpenMD { |
65 |
|
class AlphaHull : public Hull { |
66 |
|
public: |
67 |
< |
|
68 |
< |
AlphaHull(double alpha); |
67 |
> |
|
68 |
> |
AlphaHull(RealType alpha); |
69 |
|
virtual ~AlphaHull(){}; |
70 |
< |
|
71 |
< |
void computeHull( std::vector<StuntDouble*> bodydoubles ); |
72 |
< |
|
70 |
> |
|
71 |
> |
void computeHull( vector<StuntDouble*> bodydoubles ); |
72 |
> |
|
73 |
|
/* Total area of Hull*/ |
74 |
< |
RealType getArea(){return area_;} |
75 |
< |
|
74 |
> |
RealType getArea(){ return area_; } |
75 |
> |
|
76 |
|
/* Total Volume enclosed by Hull */ |
77 |
|
RealType getVolume(){ return volume_; } |
78 |
< |
|
79 |
< |
std::vector<Triangle> getMesh(){return Triangles_;} |
78 |
> |
|
79 |
> |
vector<Triangle> getMesh(){ return Triangles_; } |
80 |
|
void printHull(const std::string& geomFileName); |
81 |
< |
|
81 |
> |
|
82 |
|
protected: |
83 |
|
RealType volume_; |
84 |
|
RealType area_; |
85 |
+ |
Mat3x3d boundingBox_; |
86 |
|
int dim_; |
87 |
< |
double alpha_; |
87 |
> |
RealType alpha_; |
88 |
|
const std::string options_; |
89 |
|
|
90 |
|
private: |
91 |
< |
std::vector<Triangle> Triangles_; |
91 |
> |
vector<Triangle> Triangles_; |
92 |
|
}; |
93 |
|
} |
94 |
< |
#endif /*MATH_CONVEXHULL_HPP_*/ |
94 |
> |
#endif |