38 |
|
#include "config.h" |
39 |
|
#include "math/SquareMatrix3.hpp" |
40 |
|
|
41 |
< |
namespace oopse { |
41 |
> |
namespace OpenMD { |
42 |
|
|
43 |
|
class RMSD{ |
44 |
|
public: |
52 |
|
ref_ = ref; |
53 |
|
ref_com = V3Zero; |
54 |
|
|
55 |
< |
for (int n=0; n<ref_.size(); n++) { |
55 |
> |
for (unsigned int n = 0; n < ref_.size(); n++) { |
56 |
|
ref_com += ref_[n]; |
57 |
|
} |
58 |
|
ref_com /= (RealType)ref.size(); |
71 |
|
*/ |
72 |
|
RealType calculate_rmsd(std::vector<Vector3d> mov, |
73 |
|
Vector3d mov_com, |
74 |
< |
Vector3d mov_to_ref, |
75 |
< |
RotMat3x3d U); |
74 |
> |
Vector3d mov_to_ref); |
75 |
|
|
77 |
– |
|
76 |
|
/* |
77 |
|
* optimal_superposition() |
78 |
|
* |
79 |
|
* Returns best-fit rotation matrix |
80 |
|
*/ |
81 |
< |
RotMat3x3d optimal_superposition(std::vector<Vector3d> mov); |
81 |
> |
RotMat3x3d optimal_superposition(std::vector<Vector3d> mov, |
82 |
> |
Vector3d mov_com, |
83 |
> |
Vector3d mov_to_ref); |
84 |
|
|
85 |
< |
|
85 |
> |
|
86 |
|
protected: |
87 |
|
std::vector<Vector3d> ref_; |
88 |
|
Vector3d ref_com; |