From a5cd1457ccbbddc57209d655b4107fcd5374a652 Mon Sep 17 00:00:00 2001 From: Ryan Curtin Date: Mon, 15 Nov 2010 19:05:48 +0000 Subject: [PATCH] Fix "translation" errors in DeflationGaussUpdate_() and fix spelling in comment. --- fastlib/branches/fastlib-stl/mlpack/fastica/fastica.h | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/fastlib/branches/fastlib-stl/mlpack/fastica/fastica.h b/fastlib/branches/fastlib-stl/mlpack/fastica/fastica.h index 7e30acfa94..c1ac74694d 100644 --- a/fastlib/branches/fastlib-stl/mlpack/fastica/fastica.h +++ b/fastlib/branches/fastlib-stl/mlpack/fastica/fastica.h @@ -324,15 +324,16 @@ class FastICA { * Deflation Newton-Raphson using Gaussian contrast function */ void DeflationGaussUpdate_(index_t n, const arma::mat& X, arma::vec& w) { - arma::vec u, u_sq, u_sq_a, ex, temp1; + arma::vec u, u_sq, u_sq_a, ex; u = trans(X) * w; u_sq = pow(u, 2); u_sq_a = -a2_ * u_sq; // u is gauss ex = exp(u_sq_a / 2.0); + u = (ex % u); ex += (ex % u_sq_a); - // ex id dGauss + // ex is dGauss w *= -accu(ex); w += (X * u);