ViewVC Help

View File | Revision Log | Show Annotations | View Changeset | Root Listing

root/OpenMD/branches/development/src/math/RMSD.cpp

Comparing:
trunk/src/math/RMSD.cpp (file contents), Revision 1349 by cli2, Wed May 20 19:35:05 2009 UTC vs.
branches/development/src/math/RMSD.cpp (file contents), Revision 1825 by gezelter, Wed Jan 9 19:27:52 2013 UTC

#	Line 1 | Line 1
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,
#	Line 72 | Line 72 | RealType RMSD::calculate_rmsd(std::vector<Vector3d> mo
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		}
#	Line 109 | Line 109 | RotMat3x3d RMSD::optimal_superposition(std::vector<Vec
109
110		/* initialize */
111		Mat3x3d R(0.0);
112	–	RealType E0 = 0.0;
112
113		for (n=0; n < n_vec; n++) {
114

Comparing:
trunk/src/math/RMSD.cpp (property svn:keywords), Revision 1349 by cli2, Wed May 20 19:35:05 2009 UTC vs.
branches/development/src/math/RMSD.cpp (property svn:keywords), Revision 1825 by gezelter, Wed Jan 9 19:27:52 2013 UTC

#	Line 0 | Line 1
1	+	Author Id Revision Date

Diff Legend

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