aboutsummaryrefslogtreecommitdiff
path: root/lamp.cpp
diff options
context:
space:
mode:
authorSamuel Fadel <samuelfadel@gmail.com>2015-06-01 18:11:33 -0300
committerSamuel Fadel <samuelfadel@gmail.com>2015-06-01 18:11:33 -0300
commitb7e1060b7cb71b30b91cc65c011b719e255c0a06 (patch)
tree5411613d7b0529510dbd94ba26bb4aa2aa7e58be /lamp.cpp
parent4e0b46a727f6ea727b9e7920150609c58ce65fce (diff)
Minor code changes in main/lamp/forceScheme.
Diffstat (limited to 'lamp.cpp')
-rw-r--r--lamp.cpp3
1 files changed, 1 insertions, 2 deletions
diff --git a/lamp.cpp b/lamp.cpp
index b7a1027..7d36877 100644
--- a/lamp.cpp
+++ b/lamp.cpp
@@ -16,7 +16,7 @@ void mp::lamp(const arma::mat &X, const arma::uvec &sampleIndices, const arma::m
const arma::mat &Xs = X.rows(sampleIndices);
arma::uword sampleSize = sampleIndices.n_elem;
- #pragma omp parallel for shared(X, Xs, sampleIndices, Ys, Y)
+ #pragma omp parallel for shared(X, Xs, Ys, Y)
for (arma::uword i = 0; i < X.n_rows; i++) {
const arma::rowvec &point = X.row(i);
@@ -52,7 +52,6 @@ void mp::lamp(const arma::mat &X, const arma::uvec &sampleIndices, const arma::m
arma::svd(U, s, V, At * B);
arma::mat M = U.cols(0, 1) * V.t();
- // the projection of point i
Y(i, arma::span(0, 1)) = (point - Xtil) * M + Ytil;
}
}