comparison toolboxes/FullBNT-1.0.7/netlab3.3/gtmmag.m @ 0:e9a9cd732c1e tip

first hg version after svn

+function mags = gtmmag(net, latent_data)
+%GTMMAG	Magnification factors for a GTM
+%
+%	Description
+%	 MAGS = GTMMAG(NET, LATENTDATA) takes a GTM structure NET, and
+%	computes the magnification factors for each point the latent space
+%	contained in LATENTDATA.
+%
+%	See also
+%	GTM, GTMPOST, GTMLMEAN
+%
+
+%	Copyright (c) Ian T Nabney (1996-2001)
+
+errstring = consist(net, 'gtm');
+if ~isempty(errstring)
+  error(errstring);
+end
+
+Jacs = rbfjacob(net.rbfnet, latent_data);
+nlatent = size(latent_data, 1);
+mags = zeros(nlatent, 1);
+temp = zeros(net.rbfnet.nin, net.rbfnet.nout);
+for m = 1:nlatent
+  temp = squeeze(Jacs(m, :, :));  % Turn into a 2d matrix
+  mags(m) = sqrt(det(temp*temp'));
+end