From 49f1a5ef142865e7981525b87b572bd6d0928db1 Mon Sep 17 00:00:00 2001
From: Mark Jenkinson <mark@fmrib.ox.ac.uk>
Date: Tue, 27 Nov 2007 15:03:54 +0000
Subject: [PATCH] Added verboseness, but not working yet
---
fsldecorr4d.cc | 6 ++++++
1 file changed, 6 insertions(+)
diff --git a/fsldecorr4d.cc b/fsldecorr4d.cc
index 7053018..feab591 100644
--- a/fsldecorr4d.cc
+++ b/fsldecorr4d.cc
@@ -107,8 +107,10 @@ int do_work(int argc, char *argv[])
ColumnVector XY(nc);
XY=0.0; XX=0.0;
+ if (verbose.value()) { cout << "Calculating matrix values:" << endl; }
// calculate matrix elements (4D correlations)
for (int n=1; n<=nc; n++) {
+ if (verbose.value()) { cout << "Component #" << n << endl; }
for (int t=0; t<nt; t++) {
for (int z=smaps[n-1].minz(); z<=smaps[n-1].maxz(); z++) {
for (int y=smaps[n-1].miny(); y<=smaps[n-1].maxy(); y++) {
@@ -131,10 +133,14 @@ int do_work(int argc, char *argv[])
}
}
+ if (verbose.value()) { cout << "Finding amplitudes" << endl; }
// find amplitudes for each component
ColumnVector Beta(nc);
Beta = pinv(XX)*XY;
+ if (verbose.value()) { cout << "Amplitudes = " << Beta.t() << endl << endl; }
+
+ if (verbose.value()) { cout << "Removing components" << endl; }
// remove components from input data
volume4D<float> vout(vin);
for (int n=1; n<=nc; n++) {
--
GitLab