← Back to branch summary

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

« back to all changes in this revision

Viewing changes to CXSparse/MATLAB/Test/dmperm_test.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/MATLAB/Test/dmperm_test.m
 
1
function dmperm_test
 
2
%DMPERM_TEST test cs_dmperm
 
3
%
 
4
% Example:
 
5
%   dmperm_test
 
6
% See also: testall
 
7
 
 
8
%   Copyright 2006-2007, Timothy A. Davis.
 
9
%   http://www.cise.ufl.edu/research/sparse
 
10
 
 
11
index = UFget ;
 
12
 
 
13
f = find (index.nrows ~= index.ncols) ;
 
14
[ignore i] = sort (index.nrows(f) ./ index.ncols(f)) ;
 
15
f = [209:211 f(i)] ;
 
16
 
 
17
nmat = length(f) ;
 
18
 
 
19
tt1 = zeros (1,nmat) ;
 
20
tt2 = zeros (1,nmat) ;
 
21
tt3 = zeros (1,nmat) ;
 
22
tt4 = zeros (1,nmat) ;
 
23
mm = zeros (1,nmat) ;
 
24
nn = zeros (1,nmat) ;
 
25
ss = zeros (1,nmat) ;
 
26
me = zeros (1,nmat) ;
 
27
ne = zeros (1,nmat) ;
 
28
 
 
29
p = cs_dmperm (sparse (1)) ;
 
30
 
 
31
for k = 1:length(f)
 
32
    
 
33
    i = f(k) ;
 
34
 
 
35
    Prob = UFget (i)                                                        %#ok
 
36
 
 
37
    A = Prob.A ;
 
38
    [m n] = size (A) ;
 
39
    if (m > n)
 
40
        % make sure A is short and fat
 
41
        A = A' ;
 
42
    end
 
43
    % C is tall and thin
 
44
    C = A' ;
 
45
    [m n] = size (A) ;
 
46
 
 
47
    k1 = 0 ;
 
48
    t1 = 0 ;
 
49
    while (t1 < 1)
 
50
        tic
 
51
        p = cs_dmperm (A) ;
 
52
        t = toc ;
 
53
        t1 = t1 + t ;
 
54
        k1 = k1 + 1 ;
 
55
    end
 
56
    t1 = t1 / k1 ;
 
57
    s1 = sum (p > 0) ;
 
58
 
 
59
    k2 = 0 ;
 
60
    t2 = 0 ;
 
61
    while (t2 < 1)
 
62
        tic
 
63
        p = cs_dmperm (C) ;
 
64
        t = toc ;
 
65
        t2 = t2 + t ;
 
66
        k2 = k2 + 1 ;
 
67
    end
 
68
    t2 = t2 / k2 ;
 
69
    s2 = sum (p > 0) ;
 
70
 
 
71
    k3 = 0 ;
 
72
    t3 = 0 ;
 
73
    while (t3 < 1)
 
74
        tic
 
75
        p = cs_dmperm_orig (A) ;
 
76
        t = toc ;
 
77
        t3 = t3 + t ;
 
78
        k3 = k3 + 1 ;
 
79
    end
 
80
    t3 = t3 / k3 ;
 
81
 
 
82
    k4 = 0 ;
 
83
    t4 = 0 ;
 
84
    while (t4 < 1)
 
85
        tic
 
86
        p = cs_dmperm_orig (A') ;
 
87
        t = toc ;
 
88
        t4 = t4 + t ;
 
89
        k4 = k4 + 1 ;
 
90
    end
 
91
    t4 = t4 / k4 ;
 
92
 
 
93
    sprnk = sum (p > 0) ;
 
94
    nempty = full (sum (sum (spones (A)) == 0)) ;
 
95
    mempty = full (sum (sum (spones (C)) == 0)) ;
 
96
 
 
97
    fprintf ('[m %d:%d n %d:%d (%d)]:\n', m, mempty, n, nempty, sprnk) ;
 
98
    fprintf ('    A: t1 %10.6f (%6d)   C: t2 %10.6f (%6d) new\n', ...
 
99
        t1, k1, t2, k2) ;
 
100
    fprintf ('    A: t3 %10.6f (%6d)   C: t4 %10.6f (%6d) orig\n', ...
 
101
        t3, k3, t4, k4) ;
 
102
 
 
103
    if (s1 ~= sprnk | s2 ~= sprnk)                                          %#ok
 
104
        s1                                                                  %#ok
 
105
        s2                                                                  %#ok
 
106
        sprnk                                                               %#ok
 
107
        error ('!') ;
 
108
    end
 
109
 
 
110
    tt1 (k) = t1 ;
 
111
    tt2 (k) = t2 ;
 
112
    tt3 (k) = t3 ;
 
113
    tt4 (k) = t4 ;
 
114
    mm (k) = m ;
 
115
    nn (k) = n ;
 
116
    ss (k) = sprnk ;
 
117
    me (k) = mempty ;
 
118
    ne (k) = nempty ;
 
119
 
 
120
    clear A C
 
121
 
 
122
    semilogy (ss(1:k) ./ nn(1:k), tt1(1:k) ./ tt3(1:k), 'o', ...
 
123
        [0 1], [1 1], 'r-') ;
 
124
    drawnow
 
125
 
 
126
end