From acab7cb1d11c9c669afa07e2c26d1aa6ab4ab81b Mon Sep 17 00:00:00 2001
From: Nikolay Simakov <ns@compmodel.org>
Date: Sun, 19 Jul 2015 21:25:53 -0400
Subject: [PATCH] add RMSD calculations on nodes positions

---
 src/testutil/compare_hdf5.py | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/testutil/compare_hdf5.py b/src/testutil/compare_hdf5.py
index e382bfe..006a151 100644
--- a/src/testutil/compare_hdf5.py
+++ b/src/testutil/compare_hdf5.py
@@ -26,9 +26,9 @@
 if ds1.shape!=ds2.shape:
     print "\tDatasets have different size!"
 else:
-    ds=ds1-ds2
+    ds=(ds1-ds2)*(ds1-ds2)
     rmsd=math.sqrt(ds.mean())
-    print "\tRMSD=%.3e"%(rmsd,)
+    print "\tRMSD=%.3e Mean1=%.3e RMSD/Mean1=%.3e"%(rmsd,ds1.mean(),rmsd/ds1.mean())
 
 #compare Mash
 print "Comparing Mesh/Connections"