← Back to branch summary

~ubuntu-branches/ubuntu/natty/suitesparse/natty

« back to all changes in this revision

Viewing changes to CSparse/MATLAB/Test/chol_updown.m

  • Committer: Bazaar Package Importer
  • Author(s): Christophe Prud'homme
  • Date: 2006-12-22 10:16:15 UTC
  • Revision ID: james.westby@ubuntu.com-20061222101615-2ohaj8902oix2rnk
Tags: upstream-2.3.1
Import upstream version 2.3.1

Show diffs side-by-side

added added

removed removed

Lines of Context:
CSparse/MATLAB/Test/chol_updown.m
 
1
function [L, w] = chol_updown (L, sigma, w)
 
2
%CHOL_UPDOWN update or downdate a Cholesky factorization.
 
3
% Example:
 
4
%   [L, w] = chol_updown (L, sigma, w)
 
5
% See also: cs_demo
 
6
 
 
7
%   Copyright 2006, Timothy A. Davis.
 
8
%   http://www.cise.ufl.edu/research/sparse
 
9
 
 
10
beta = 1 ;
 
11
n = size (L,1) ;
 
12
if (n == 1)
 
13
    L = sqrt (L*L'+sigma*w*w') ;
 
14
    return ;
 
15
end
 
16
for k = 1:n
 
17
    alpha = w(k) / L(k,k) ;
 
18
    beta2 = sqrt (beta^2 + sigma*alpha^2) ;
 
19
    gamma = sigma * alpha / (beta2 * beta) ;
 
20
    if (sigma > 0)
 
21
        % update
 
22
        delta = beta / beta2 ;
 
23
        L (k,k) = delta * L (k,k) + gamma * w (k) ;
 
24
        w1 = w (k+1:n) ;
 
25
        w (k+1:n) = w (k+1:n) - alpha * L (k+1:n,k) ;
 
26
        L (k+1:n,k) = delta * L (k+1:n,k) + gamma * w1 ;
 
27
    else
 
28
        % downdate
 
29
        delta = beta2 / beta ;
 
30
        L (k,k) = delta * L (k,k) ;
 
31
        w (k+1:n) = w (k+1:n) - alpha * L (k+1:n,k) ;
 
32
        L (k+1:n,k) = delta * L (k+1:n,k) + gamma * w (k+1:n) ;
 
33
    end
 
34
    w (k) = alpha ;
 
35
    beta = beta2 ;
 
36
end