add RMSD calculations on nodes positions

TITAN2D · Jul 20, 2015 · acab7cb · acab7cb
1 parent 411d79a
commit acab7cb
Showing 1 changed file with 2 additions and 2 deletions.
diff --git 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"