From 87fbe0df94afc6a57e23eadda041356fe15515d2 Mon Sep 17 00:00:00 2001
From: Christian Beckmann <c.beckmann@donders.ru.nl>
Date: Tue, 3 Nov 2009 09:01:22 +0000
Subject: [PATCH] *** empty log message ***
---
meldata.cc | 8 +++++---
1 file changed, 5 insertions(+), 3 deletions(-)
diff --git a/meldata.cc b/meldata.cc
index 8a5ead4..27956f6 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);
--
GitLab