diff --git a/meldata.cc b/meldata.cc
index 8a5ead4e8f07089aebb1a9ad2bbe7838544753ba..27956f658a79a89f786486afce48d79360ac4b5e 100644
--- a/meldata.cc
+++ b/meldata.cc
@@ -310,10 +310,12 @@ namespace Melodic{
 						newWM=(newWM*whiteMatrix);
 					}
 					else{
+					  if(opts.debug.value())
+					    message(" --mod_pca ");
 						Matrix tmp1, tmp2;
-						tmp1 = whiteMatrix * alldat;  
-						tmp1 = remmean(tmp1,2) * tmpData.t();  
-						tmp2 = pinv(tmp1);  
+						tmp1 = whiteMatrix * alldat;
+						tmp1 = remmean(tmp1,2) * tmpData.t();
+						tmp2 = pinv(tmp1.t()).t();  
 						std_pca(tmp1 * tmpData, RXweight, Corr, pcaE, pcaD);
 						calc_white(pcaE, pcaD, order, newWM, newDWM);		
 						newDWM=(tmp2*newDWM);