1 |
|
#include "math/RMSD.hpp" |
2 |
|
#include "math/SVD.hpp" |
3 |
|
|
4 |
< |
using namespace oopse; |
4 |
> |
using namespace OpenMD; |
5 |
|
using namespace JAMA; |
6 |
|
|
7 |
|
RealType RMSD::calculate_rmsd(std::vector<Vector3d> mov, |
72 |
|
stmp.getSubVector(0,s); |
73 |
|
wtmp.getSubMatrix(0,0,w); |
74 |
|
|
75 |
< |
int is_reflection = (v.determinant() * w.determinant()) < 0.0; |
75 |
> |
int is_reflection = (v.determinant() * w.determinant()) < RealType(0.0); |
76 |
|
if (is_reflection) |
77 |
|
s(2) = -s(2); |
78 |
|
|
79 |
|
RealType rmsd_sq = (E0 - 2.0 * s.sum() )/ (RealType)n_vec; |
80 |
< |
rmsd_sq = max(rmsd_sq,0.0); |
80 |
> |
rmsd_sq = max(rmsd_sq, RealType(0.0)); |
81 |
|
RealType rmsd = sqrt(rmsd_sq); |
82 |
|
return rmsd; |
83 |
|
} |