ITE / code / shared / embedded / ITL / incompleteCholeskySigma.m

% The following function computes the Incomplete Cholesky decomposition of
% the Gram matrix generated by the random VECTOR samples when the kernel is
% normalized.
% Input:   "samples" should be COLUMN vectors of length (nxm)
%				'kernelSize' is a scalar for the kernel size.
%               "precision" is a scalar for the precision
% Output: 'choleskyMarix' contains nxd Cholesky matrix
% Default:  precision = 10^-6.
% The code uses incompleteCholeskyMulti.
% Author: Sohan Seth (	Date: 01.06.2009

function choleskyMatrix  = incompleteCholeskySigma(samples,kernelSize,precision)

if nargin == 2
    precision = 10^-6;

c = 1/(sqrt(2*pi)*kernelSize)^(size(samples,2));
choleskyMatrix  = incompleteCholeskyMulti(samples,kernelSize,precision);
choleskyMatrix = choleskyMatrix * sqrt(c);
Tip: Filter by directory path e.g. /media app.js to search for public/media/app.js.
Tip: Use camelCasing e.g. ProjME to search for
Tip: Filter by extension type e.g. /repo .js to search for all .js files in the /repo directory.
Tip: Separate your search with spaces e.g. /ssh pom.xml to search for src/ssh/pom.xml.
Tip: Use ↑ and ↓ arrow keys to navigate and return to view the file.
Tip: You can also navigate files with Ctrl+j (next) and Ctrl+k (previous) and view the file with Ctrl+o.
Tip: You can also navigate files with Alt+j (next) and Alt+k (previous) and view the file with Alt+o.