From 076008fd0e1be7f4b8d4f58a5568ea5e0a5b3817 Mon Sep 17 00:00:00 2001
From: Mark Jenkinson <mark@fmrib.ox.ac.uk>
Date: Thu, 12 Apr 2012 16:24:03 +0000
Subject: [PATCH] Put back in the multiplication by 100 to get percentages

---
 calc_grad_perc_dev.cc | 18 +++++++++---------
 1 file changed, 9 insertions(+), 9 deletions(-)

diff --git a/calc_grad_perc_dev.cc b/calc_grad_perc_dev.cc
index 86e13fa..5cd738b 100644
--- a/calc_grad_perc_dev.cc
+++ b/calc_grad_perc_dev.cc
@@ -64,27 +64,27 @@ int do_work(int argc, char* argv[])
   xpd.addvolume((tmp[0]+tmp[1])*0.5f);
   tmp=lrzgrad(fw[0],mask);
   xpd.addvolume((tmp[0]+tmp[1])*0.5f);
-  xpd[0]/=fw.xdim();
-  xpd[1]/=fw.ydim();
-  xpd[2]/=fw.zdim();
+  xpd[0]*=100.0/fw.xdim();
+  xpd[1]*=100.0/fw.ydim();
+  xpd[2]*=100.0/fw.zdim();
   tmp=lrxgrad(fw[1],mask);
   ypd.addvolume((tmp[0]+tmp[1])*0.5f);
   tmp=lrygrad(fw[1],mask);
   ypd.addvolume((tmp[0]+tmp[1])*0.5f);
   tmp=lrzgrad(fw[1],mask);
   ypd.addvolume((tmp[0]+tmp[1])*0.5f);
-  ypd[0]/=fw.xdim();
-  ypd[1]/=fw.ydim();
-  ypd[2]/=fw.zdim();
+  ypd[0]*=100.0/fw.xdim();
+  ypd[1]*=100.0/fw.ydim();
+  ypd[2]*=100.0/fw.zdim();
   tmp=lrxgrad(fw[2],mask);
   zpd.addvolume((tmp[0]+tmp[1])*0.5f);
   tmp=lrygrad(fw[2],mask);
   zpd.addvolume((tmp[0]+tmp[1])*0.5f);
   tmp=lrzgrad(fw[2],mask);
   zpd.addvolume((tmp[0]+tmp[1])*0.5f);
-  zpd[0]/=fw.xdim();
-  zpd[1]/=fw.ydim();
-  zpd[2]/=fw.zdim();
+  zpd[0]*=100.0/fw.xdim();
+  zpd[1]*=100.0/fw.ydim();
+  zpd[2]*=100.0/fw.zdim();
   string bname=fslbasename(outname.value());
   save_volume4D(xpd,bname+"_x");
   save_volume4D(ypd,bname+"_y");
-- 
GitLab