← Back to branch summary

~logan/ubuntu/trusty/suitesparse/4.2.1-3ubuntu1

« back to all changes in this revision

Viewing changes to CXSparse_newfiles/MATLAB/Test/chol_updown2.m

  • Committer: Bazaar Package Importer
  • Author(s): Christophe Prud'homme
  • Date: 2007-05-29 09:36:29 UTC
  • mfrom: (1.1.1 upstream)
  • Revision ID: james.westby@ubuntu.com-20070529093629-zowquo0b7slkk6nc
Tags: 3.0.0-2
http://bugs.debian.org/426349
http://bugs.debian.org/426354
* suitesparse builds properly twice in a row
* Bug fix: "suitesparse - FTBFS: Broken build depens: libgfortran1-dev",
  thanks to Bastian Blank (Closes: #426349).
* Bug fix: "suitesparse_3.0.0-1: FTBFS: build-depends on
  libgfortran1-dev", thanks to Steve Langasek (Closes: #426354).

Show diffs side-by-side

added added

removed removed

Lines of Context:
CXSparse_newfiles/MATLAB/Test/chol_updown2.m
 
1
function [L, w] = chol_updown2 (L, sigma, w)
 
2
%CHOL_UPDOWN2 Cholesky update/downdate (real and complex)
 
3
% (real or complex)
 
4
% Example:
 
5
%   [L, w] = chol_updown2 (L, sigma, w)
 
6
% See also: cs_demo
 
7
 
 
8
%   Copyright 2006-2007, Timothy A. Davis, William W. Hager
 
9
%   http://www.cise.ufl.edu/research/sparse
 
10
 
 
11
beta = 1 ;
 
12
n = size (L,1) ;
 
13
if (n == 1)
 
14
    wnew = L\w ;
 
15
    L = sqrt (L*L'+sigma*w*w') ;
 
16
    w = wnew ;
 
17
    return ;
 
18
end
 
19
for k = 1:n
 
20
    alpha = w(k) / L(k,k) ;
 
21
    beta2 = sqrt (beta*beta + sigma*alpha*conj(alpha)) ;
 
22
    gamma = sigma * conj(alpha) / (beta2 * beta) ;
 
23
    if (sigma > 0)
 
24
        % update
 
25
        delta = beta / beta2 ;
 
26
        L (k,k) = delta * L (k,k) + gamma * w (k) ;
 
27
        phase = abs (L (k, k))/L (k, k) ;
 
28
        L (k, k) = phase*L (k, k) ;
 
29
 
 
30
        w1 = w (k+1:n) ;
 
31
        w (k+1:n) = w (k+1:n) - alpha * L (k+1:n,k) ;
 
32
        L (k+1:n,k) = phase * (delta * L (k+1:n,k) + gamma * w1) ;
 
33
 
 
34
    else
 
35
        % downdate
 
36
        delta = beta2 / beta ;
 
37
        L (k,k) = delta * L (k,k) ;
 
38
        phase = abs (L (k, k))/L (k, k) ;
 
39
        L (k, k) = phase*L (k, k) ;
 
40
 
 
41
        w (k+1:n) = w (k+1:n) - alpha * L (k+1:n,k) ;
 
42
        L (k+1:n,k) = phase * (delta * L (k+1:n,k) + gamma * w (k+1:n)) ;
 
43
 
 
44
    end
 
45
    w (k) = alpha ;
 
46
    beta = beta2 ;
 
47
end