0001 function [GC, spectrum, CTS] = slltsa(X, G, dl, dg)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031
0032
0033
0034
0035
0036
0037
0038
0039
0040
0041
0042
0043
0044
0045
0046
0047
0048 if ~isnumeric(X) || ndims(X) ~= 2
0049 error('X should be a 2D numeric matrix');
0050 end
0051 n = size(X, 2);
0052
0053 gi = slgraphinfo(G, {[n, n]});
0054
0055 if nargin < 4
0056 dg = dl;
0057 end
0058
0059 if ~strcmp(gi.form, 'adjmat')
0060 G = sladjmat(G, 'sparse', true);
0061 end
0062
0063 out_cts = (nargout >= 3);
0064
0065
0066
0067
0068
0069 [LM, LP] = sllocaltanspace(X, G, dl);
0070
0071
0072 [Y, GM] = sllocaltancoords(LM, LP, X, G, 'T');
0073 if ~out_cts
0074 clear LM LP;
0075 end
0076
0077
0078 if out_cts
0079 [GC, spectrum, LT] = sllocalcoordalign(GM, Y, dg);
0080 else
0081 [GC, spectrum] = sllocalcoordalign(GM, Y, dg);
0082 end
0083
0084
0085 if out_cts
0086 CTS.means = LM;
0087 CTS.bases = LP;
0088 CTS.ftrans = LT;
0089 end
0090