diff --git a/fsldecorr4d.cc b/fsldecorr4d.cc
index 70530189210615105766b97acd5f4a2caac5c356..feab5913813fd419264e640960b3d023e079f559 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++) {