diff toolboxes/distance_learning/mlr/distance/setDistanceFullMKL.m @ 0:e9a9cd732c1e tip

first hg version after svn

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/toolboxes/distance_learning/mlr/distance/setDistanceFullMKL.m	Tue Feb 10 15:05:51 2015 +0000
@@ -0,0 +1,32 @@
+function D = setDistanceFullMKL(X, W, Ifrom, Ito)
+%
+% D = setDistanceFullMKL(X, W, Ifrom, Ito)
+%
+%   X       = d-by-n-by-m data matrix
+%   W       = d-by-d-by-m PSD matrix
+%   Ifrom   = k-by-1 vector of source points
+%   Ito     = j-by-1 vector of destination points
+%
+%   D = n-by-n matrix of squared euclidean distances from Ifrom to Ito
+%       D is sparse, and only the rows corresponding to Ifrom and
+%       columns corresponding to Ito are populated.
+
+    [d,n,m]       = size(X);
+
+    D = 0;
+
+    if nargin < 4
+        Ito     = 1:n;
+    end
+
+    parfor i = 1:m
+        [vecs,vals] = eig(0.5 * (W(:,:,i) + W(:,:,i)'));
+        L           = real(abs(vals)).^0.5 * vecs';
+
+        Vfrom   = L * X(:,Ifrom,i);
+
+        Vto     = L * X(:,Ito,i);
+
+        D = D + distToFrom(n, Vto, Vfrom, Ito, Ifrom);
+    end
+end