*** empty log message ***

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);