← Back to branch summary

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

« back to all changes in this revision

Viewing changes to CXSparse/Demo/cs_demo.out

  • 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/Demo/cs_demo.out
1
 
make[1]: Entering directory `/amd/netapp3/vol/research0a/research18/sparse/SuiteSparse/CXSparse/Demo'
 
1
make[1]: Entering directory `/amd/birch/export/research07/sparse/SuiteSparse/CXSparse/Demo'
2
2
./cs_idemo < ../Matrix/t2
3
3
 
4
 
--- cs_idemo
 
4
--- cs_idemo, size of CS_INT: 4
5
5
 
6
6
T:
7
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
7
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
8
8
triplet: 4-by-4, nzmax: 16 nnz: 10
9
9
    2 2 : (3, 3.14159)
10
10
    1 0 : (3.1, 42)
18
18
    2 1 : (1.7, 1)
19
19
 
20
20
Treal:
21
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
21
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
22
22
triplet: 4-by-4, nzmax: 16 nnz: 10
23
23
    2 2 : 3
24
24
    1 0 : 3.1
32
32
    2 1 : 1.7
33
33
 
34
34
Timag:
35
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
35
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
36
36
triplet: 4-by-4, nzmax: 16 nnz: 10
37
37
    2 2 : 3.14159
38
38
    1 0 : 42
46
46
    2 1 : 1
47
47
 
48
48
A:
49
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
49
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
50
50
4-by-4, nzmax: 10 nnz: 10, 1-norm: 106.075
51
51
    col 0 : locations 0 to 2
52
52
      1 : (3.1, 42)
64
64
      1 : (0.9, 99)
65
65
 
66
66
C1 = real(A):
67
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
67
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
68
68
4-by-4, nzmax: 10 nnz: 10, 1-norm: 11.1
69
69
    col 0 : locations 0 to 2
70
70
      1 : 3.1
82
82
      1 : 0.9
83
83
 
84
84
C2 = imag(A):
85
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
85
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
86
86
4-by-4, nzmax: 10 nnz: 10, 1-norm: 106
87
87
    col 0 : locations 0 to 2
88
88
      1 : 42
100
100
      1 : 99
101
101
 
102
102
A1:
103
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
103
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
104
104
4-by-4, nzmax: 10 nnz: 10, 1-norm: 11.1
105
105
    col 0 : locations 0 to 2
106
106
      1 : (3.1, 0)
118
118
      1 : (0.9, 0)
119
119
 
120
120
A2:
121
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
121
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
122
122
4-by-4, nzmax: 10 nnz: 10, 1-norm: 106
123
123
    col 0 : locations 0 to 2
124
124
      1 : (0, 42)
136
136
      1 : (0, 99)
137
137
 
138
138
B = conj(A):
139
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
139
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
140
140
4-by-4, nzmax: 10 nnz: 10, 1-norm: 106.075
141
141
    col 0 : locations 0 to 2
142
142
      1 : (3.1, -42)
154
154
      1 : (0.9, -99)
155
155
./cs_ldemo < ../Matrix/t2
156
156
 
157
 
--- cs_ldemo
 
157
--- cs_ldemo, size of CS_INT: 8
158
158
 
159
159
T:
160
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
160
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
161
161
triplet: 4-by-4, nzmax: 16 nnz: 10
162
162
    2 2 : (3, 3.14159)
163
163
    1 0 : (3.1, 42)
171
171
    2 1 : (1.7, 1)
172
172
 
173
173
Treal:
174
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
174
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
175
175
triplet: 4-by-4, nzmax: 16 nnz: 10
176
176
    2 2 : 3
177
177
    1 0 : 3.1
185
185
    2 1 : 1.7
186
186
 
187
187
Timag:
188
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
188
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
189
189
triplet: 4-by-4, nzmax: 16 nnz: 10
190
190
    2 2 : 3.14159
191
191
    1 0 : 42
199
199
    2 1 : 1
200
200
 
201
201
A:
202
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
202
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
203
203
4-by-4, nzmax: 10 nnz: 10, 1-norm: 106.075
204
204
    col 0 : locations 0 to 2
205
205
      1 : (3.1, 42)
217
217
      1 : (0.9, 99)
218
218
 
219
219
C1 = real(A):
220
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
220
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
221
221
4-by-4, nzmax: 10 nnz: 10, 1-norm: 11.1
222
222
    col 0 : locations 0 to 2
223
223
      1 : 3.1
235
235
      1 : 0.9
236
236
 
237
237
C2 = imag(A):
238
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
238
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
239
239
4-by-4, nzmax: 10 nnz: 10, 1-norm: 106
240
240
    col 0 : locations 0 to 2
241
241
      1 : 42
253
253
      1 : 99
254
254
 
255
255
A1:
256
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
256
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
257
257
4-by-4, nzmax: 10 nnz: 10, 1-norm: 11.1
258
258
    col 0 : locations 0 to 2
259
259
      1 : (3.1, 0)
271
271
      1 : (0.9, 0)
272
272
 
273
273
A2:
274
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
274
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
275
275
4-by-4, nzmax: 10 nnz: 10, 1-norm: 106
276
276
    col 0 : locations 0 to 2
277
277
      1 : (0, 42)
289
289
      1 : (0, 99)
290
290
 
291
291
B = conj(A):
292
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
292
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
293
293
4-by-4, nzmax: 10 nnz: 10, 1-norm: 106.075
294
294
    col 0 : locations 0 to 2
295
295
      1 : (3.1, -42)
307
307
      1 : (0.9, -99)
308
308
./cs_demo1 < ../Matrix/t1
309
309
T:
310
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
310
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
311
311
triplet: 4-by-4, nzmax: 16 nnz: 10
312
312
    2 2 : 3
313
313
    1 0 : 3.1
320
320
    0 0 : 4.5
321
321
    2 1 : 1.7
322
322
A:
323
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
323
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
324
324
4-by-4, nzmax: 10 nnz: 10, 1-norm: 11.1
325
325
    col 0 : locations 0 to 2
326
326
      1 : 3.1
337
337
      3 : 1
338
338
      1 : 0.9
339
339
AT:
340
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
340
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
341
341
4-by-4, nzmax: 10 nnz: 10, 1-norm: 7.7
342
342
    col 0 : locations 0 to 1
343
343
      0 : 4.5
354
354
      1 : 0.4
355
355
      3 : 1
356
356
D:
357
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
357
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
358
358
4-by-4, nzmax: 16 nnz: 16, 1-norm: 139.58
359
359
    col 0 : locations 0 to 3
360
360
      1 : 13.95
380
380
 
381
381
--- Matrix: 4-by-4, nnz: 10 (sym: 0: nnz 0), norm: 1.11e+01
382
382
blocks: 1 singletons: 0 structural rank: 4
383
 
QR   natural    time:     0.00 resid: 1.04e-17
384
 
QR   amd(A'*A)  time:     0.00 resid: 2.03e-17
385
 
LU   natural    time:     0.00 resid: 1.04e-17
386
 
LU   amd(A+A')  time:     0.00 resid: 4.94e-18
387
 
LU   amd(S'*S)  time:     0.00 resid: 4.94e-18
388
 
LU   amd(A'*A)  time:     0.00 resid: 4.94e-18
 
383
QR   natural    time:     0.00 resid: 3.06e-17
 
384
QR   amd(A'*A)  time:     0.00 resid: 2.30e-17
 
385
LU   natural    time:     0.00 resid: 1.53e-17
 
386
LU   amd(A+A')  time:     0.00 resid: 1.53e-17
 
387
LU   amd(S'*S)  time:     0.00 resid: 0.00e+00
 
388
LU   amd(A'*A)  time:     0.00 resid: 1.53e-17
389
389
./cs_demo2 < ../Matrix/fs_183_1
390
390
 
391
391
--- Matrix: 183-by-183, nnz: 988 (sym: 0: nnz 0), norm: 1.70e+09
392
392
zero entries dropped: 71
393
393
tiny entries dropped: 10
394
394
blocks: 38 singletons: 37 structural rank: 183
395
 
QR   natural    time:     0.01 resid: 3.15e-28
396
 
QR   amd(A'*A)  time:     0.00 resid: 2.60e-28
397
 
LU   natural    time:     0.01 resid: 3.08e-28
398
 
LU   amd(A+A')  time:     0.00 resid: 1.42e-27
399
 
LU   amd(S'*S)  time:     0.00 resid: 7.11e-28
400
 
LU   amd(A'*A)  time:     0.00 resid: 7.11e-28
 
395
QR   natural    time:     0.00 resid: 1.42e-27
 
396
QR   amd(A'*A)  time:     0.00 resid: 3.35e-28
 
397
LU   natural    time:     0.01 resid: 6.20e-28
 
398
LU   amd(A+A')  time:     0.00 resid: 1.55e-27
 
399
LU   amd(S'*S)  time:     0.00 resid: 6.98e-28
 
400
LU   amd(A'*A)  time:     0.00 resid: 6.98e-28
401
401
./cs_demo2 < ../Matrix/west0067
402
402
 
403
403
--- Matrix: 67-by-67, nnz: 294 (sym: 0: nnz 0), norm: 6.14e+00
404
404
blocks: 2 singletons: 1 structural rank: 67
405
 
QR   natural    time:     0.00 resid: 2.81e-17
406
 
QR   amd(A'*A)  time:     0.00 resid: 1.30e-17
407
 
LU   natural    time:     0.00 resid: 3.85e-17
408
 
LU   amd(A+A')  time:     0.00 resid: 1.95e-17
409
 
LU   amd(S'*S)  time:     0.00 resid: 2.60e-17
 
405
QR   natural    time:     0.00 resid: 5.19e-17
 
406
QR   amd(A'*A)  time:     0.00 resid: 3.25e-17
 
407
LU   natural    time:     0.00 resid: 3.89e-17
 
408
LU   amd(A+A')  time:     0.00 resid: 2.27e-17
 
409
LU   amd(S'*S)  time:     0.00 resid: 1.95e-17
410
410
LU   amd(A'*A)  time:     0.00 resid: 2.60e-17
411
411
./cs_demo2 < ../Matrix/lp_afiro
412
412
 
413
413
--- Matrix: 27-by-51, nnz: 102 (sym: 0: nnz 0), norm: 3.43e+00
414
414
blocks: 1 singletons: 0 structural rank: 27
415
 
QR   natural    time:     0.00 resid: 4.30e-16
416
 
QR   amd(A'*A)  time:     0.00 resid: 1.92e-16
 
415
QR   natural    time:     0.00 resid: 3.85e-16
 
416
QR   amd(A'*A)  time:     0.00 resid: 1.50e-16
417
417
./cs_demo2 < ../Matrix/ash219
418
418
 
419
419
--- Matrix: 219-by-85, nnz: 438 (sym: 0: nnz 0), norm: 9.00e+00
424
424
 
425
425
--- Matrix: 492-by-490, nnz: 49920 (sym: 0: nnz 0), norm: 9.29e-01
426
426
blocks: 10 singletons: 8 structural rank: 448
427
 
QR   natural    time:     0.26 resid:      nan
428
 
QR   amd(A'*A)  time:     0.31 resid:      nan
 
427
QR   natural    time:     0.13 resid:      nan
 
428
QR   amd(A'*A)  time:     0.19 resid:      nan
429
429
./cs_demo2 < ../Matrix/bcsstk01
430
430
 
431
431
--- Matrix: 48-by-48, nnz: 224 (sym: -1: nnz 400), norm: 3.57e+09
432
432
blocks: 1 singletons: 0 structural rank: 48
433
 
QR   natural    time:     0.00 resid: 3.97e-19
434
 
QR   amd(A'*A)  time:     0.00 resid: 2.36e-19
435
 
LU   natural    time:     0.00 resid: 2.63e-19
436
 
LU   amd(A+A')  time:     0.00 resid: 8.63e-20
437
 
LU   amd(S'*S)  time:     0.00 resid: 2.04e-19
438
 
LU   amd(A'*A)  time:     0.00 resid: 2.04e-19
439
 
Chol natural    time:     0.00 resid: 1.90e-19
440
 
Chol amd(A+A')  time:     0.00 resid: 2.01e-19
 
433
QR   natural    time:     0.00 resid: 3.65e-19
 
434
QR   amd(A'*A)  time:     0.00 resid: 4.02e-19
 
435
LU   natural    time:     0.00 resid: 2.17e-19
 
436
LU   amd(A+A')  time:     0.00 resid: 1.87e-19
 
437
LU   amd(S'*S)  time:     0.00 resid: 2.38e-19
 
438
LU   amd(A'*A)  time:     0.00 resid: 2.38e-19
 
439
Chol natural    time:     0.00 resid: 2.64e-19
 
440
Chol amd(A+A')  time:     0.00 resid: 2.55e-19
441
441
./cs_demo3 < ../Matrix/bcsstk01
442
442
 
443
443
--- Matrix: 48-by-48, nnz: 224 (sym: -1: nnz 400), norm: 3.57e+09
446
446
symbolic chol time     0.00
447
447
numeric  chol time     0.00
448
448
solve    chol time     0.00
449
 
original: resid: 2.01e-19
 
449
original: resid: 2.55e-19
450
450
update:   time:     0.00
451
 
update:   time:     0.00 (incl solve) resid: 8.39e-19
452
 
rechol:   time:     0.00 (incl solve) resid: 1.20e-18
 
451
update:   time:     0.00 (incl solve) resid: 9.66e-19
 
452
rechol:   time:     0.00 (incl solve) resid: 1.55e-18
453
453
downdate: time:     0.00
454
 
downdate: time:     0.00 (incl solve) resid: 3.30e-17
 
454
downdate: time:     0.00 (incl solve) resid: 3.74e-17
455
455
./cs_demo2 < ../Matrix/bcsstk16
456
456
 
457
457
--- Matrix: 4884-by-4884, nnz: 147631 (sym: -1: nnz 290378), norm: 7.01e+09
458
458
blocks: 75 singletons: 74 structural rank: 4884
459
 
QR   amd(A'*A)  time:     2.44 resid: 1.28e-22
460
 
LU   amd(A+A')  time:     1.46 resid: 1.06e-22
461
 
LU   amd(S'*S)  time:     1.41 resid: 1.21e-22
462
 
LU   amd(A'*A)  time:     1.48 resid: 1.70e-22
463
 
Chol amd(A+A')  time:     0.65 resid: 9.28e-23
 
459
QR   amd(A'*A)  time:     1.86 resid: 2.01e-22
 
460
LU   amd(A+A')  time:     1.12 resid: 1.10e-22
 
461
LU   amd(S'*S)  time:     1.12 resid: 1.28e-22
 
462
LU   amd(A'*A)  time:     1.17 resid: 1.78e-22
 
463
Chol amd(A+A')  time:     0.36 resid: 1.19e-22
464
464
./cs_demo3 < ../Matrix/bcsstk16
465
465
 
466
466
--- Matrix: 4884-by-4884, nnz: 147631 (sym: -1: nnz 290378), norm: 7.01e+09
467
467
 
468
468
chol then update/downdate amd(A+A')  
469
469
symbolic chol time     0.03
470
 
numeric  chol time     0.62
 
470
numeric  chol time     0.33
471
471
solve    chol time     0.01
472
 
original: resid: 9.28e-23
 
472
original: resid: 1.19e-22
473
473
update:   time:     0.00
474
 
update:   time:     0.01 (incl solve) resid: 8.62e-24
475
 
rechol:   time:     0.63 (incl solve) resid: 8.72e-24
 
474
update:   time:     0.01 (incl solve) resid: 1.12e-23
 
475
rechol:   time:     0.33 (incl solve) resid: 1.17e-23
476
476
downdate: time:     0.00
477
 
downdate: time:     0.01 (incl solve) resid: 3.60e-22
 
477
downdate: time:     0.00 (incl solve) resid: 4.09e-22
478
478
./cs_di_demo1 < ../Matrix/t1
479
479
T:
480
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
480
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
481
481
triplet: 4-by-4, nzmax: 16 nnz: 10
482
482
    2 2 : 3
483
483
    1 0 : 3.1
490
490
    0 0 : 4.5
491
491
    2 1 : 1.7
492
492
A:
493
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
493
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
494
494
4-by-4, nzmax: 10 nnz: 10, 1-norm: 11.1
495
495
    col 0 : locations 0 to 2
496
496
      1 : 3.1
507
507
      3 : 1
508
508
      1 : 0.9
509
509
AT:
510
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
510
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
511
511
4-by-4, nzmax: 10 nnz: 10, 1-norm: 7.7
512
512
    col 0 : locations 0 to 1
513
513
      0 : 4.5
524
524
      1 : 0.4
525
525
      3 : 1
526
526
D:
527
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
527
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
528
528
4-by-4, nzmax: 16 nnz: 16, 1-norm: 139.58
529
529
    col 0 : locations 0 to 3
530
530
      1 : 13.95
550
550
 
551
551
--- Matrix: 4-by-4, nnz: 10 (sym: 0: nnz 0), norm: 1.11e+01
552
552
blocks: 1 singletons: 0 structural rank: 4
553
 
QR   natural    time:     0.00 resid: 1.04e-17
554
 
QR   amd(A'*A)  time:     0.00 resid: 2.03e-17
555
 
LU   natural    time:     0.00 resid: 1.04e-17
556
 
LU   amd(A+A')  time:     0.00 resid: 4.94e-18
557
 
LU   amd(S'*S)  time:     0.00 resid: 4.94e-18
558
 
LU   amd(A'*A)  time:     0.00 resid: 4.94e-18
 
553
QR   natural    time:     0.00 resid: 3.06e-17
 
554
QR   amd(A'*A)  time:     0.00 resid: 2.30e-17
 
555
LU   natural    time:     0.00 resid: 1.53e-17
 
556
LU   amd(A+A')  time:     0.00 resid: 1.53e-17
 
557
LU   amd(S'*S)  time:     0.00 resid: 0.00e+00
 
558
LU   amd(A'*A)  time:     0.00 resid: 1.53e-17
559
559
./cs_di_demo2 < ../Matrix/fs_183_1
560
560
 
561
561
--- Matrix: 183-by-183, nnz: 988 (sym: 0: nnz 0), norm: 1.70e+09
562
562
zero entries dropped: 71
563
563
tiny entries dropped: 10
564
564
blocks: 38 singletons: 37 structural rank: 183
565
 
QR   natural    time:     0.00 resid: 3.15e-28
566
 
QR   amd(A'*A)  time:     0.01 resid: 2.60e-28
567
 
LU   natural    time:     0.00 resid: 3.08e-28
568
 
LU   amd(A+A')  time:     0.01 resid: 1.42e-27
569
 
LU   amd(S'*S)  time:     0.00 resid: 7.11e-28
570
 
LU   amd(A'*A)  time:     0.00 resid: 7.11e-28
 
565
QR   natural    time:     0.00 resid: 1.42e-27
 
566
QR   amd(A'*A)  time:     0.00 resid: 3.35e-28
 
567
LU   natural    time:     0.00 resid: 6.20e-28
 
568
LU   amd(A+A')  time:     0.00 resid: 1.55e-27
 
569
LU   amd(S'*S)  time:     0.01 resid: 6.98e-28
 
570
LU   amd(A'*A)  time:     0.00 resid: 6.98e-28
571
571
./cs_di_demo2 < ../Matrix/west0067
572
572
 
573
573
--- Matrix: 67-by-67, nnz: 294 (sym: 0: nnz 0), norm: 6.14e+00
574
574
blocks: 2 singletons: 1 structural rank: 67
575
 
QR   natural    time:     0.00 resid: 2.81e-17
576
 
QR   amd(A'*A)  time:     0.00 resid: 1.30e-17
577
 
LU   natural    time:     0.00 resid: 3.85e-17
578
 
LU   amd(A+A')  time:     0.00 resid: 1.95e-17
579
 
LU   amd(S'*S)  time:     0.00 resid: 2.60e-17
 
575
QR   natural    time:     0.00 resid: 5.19e-17
 
576
QR   amd(A'*A)  time:     0.00 resid: 3.25e-17
 
577
LU   natural    time:     0.00 resid: 3.89e-17
 
578
LU   amd(A+A')  time:     0.00 resid: 2.27e-17
 
579
LU   amd(S'*S)  time:     0.00 resid: 1.95e-17
580
580
LU   amd(A'*A)  time:     0.00 resid: 2.60e-17
581
581
./cs_di_demo2 < ../Matrix/lp_afiro
582
582
 
583
583
--- Matrix: 27-by-51, nnz: 102 (sym: 0: nnz 0), norm: 3.43e+00
584
584
blocks: 1 singletons: 0 structural rank: 27
585
 
QR   natural    time:     0.00 resid: 4.30e-16
586
 
QR   amd(A'*A)  time:     0.00 resid: 1.92e-16
 
585
QR   natural    time:     0.00 resid: 3.85e-16
 
586
QR   amd(A'*A)  time:     0.00 resid: 1.50e-16
587
587
./cs_di_demo2 < ../Matrix/ash219
588
588
 
589
589
--- Matrix: 219-by-85, nnz: 438 (sym: 0: nnz 0), norm: 9.00e+00
594
594
 
595
595
--- Matrix: 492-by-490, nnz: 49920 (sym: 0: nnz 0), norm: 9.29e-01
596
596
blocks: 10 singletons: 8 structural rank: 448
597
 
QR   natural    time:     0.24 resid:      nan
598
 
QR   amd(A'*A)  time:     0.31 resid:      nan
 
597
QR   natural    time:     0.13 resid:      nan
 
598
QR   amd(A'*A)  time:     0.18 resid:      nan
599
599
./cs_di_demo2 < ../Matrix/bcsstk01
600
600
 
601
601
--- Matrix: 48-by-48, nnz: 224 (sym: -1: nnz 400), norm: 3.57e+09
602
602
blocks: 1 singletons: 0 structural rank: 48
603
 
QR   natural    time:     0.00 resid: 3.97e-19
604
 
QR   amd(A'*A)  time:     0.00 resid: 2.36e-19
605
 
LU   natural    time:     0.00 resid: 2.63e-19
606
 
LU   amd(A+A')  time:     0.00 resid: 8.63e-20
607
 
LU   amd(S'*S)  time:     0.00 resid: 2.04e-19
608
 
LU   amd(A'*A)  time:     0.00 resid: 2.04e-19
609
 
Chol natural    time:     0.00 resid: 1.90e-19
610
 
Chol amd(A+A')  time:     0.00 resid: 2.01e-19
 
603
QR   natural    time:     0.00 resid: 3.65e-19
 
604
QR   amd(A'*A)  time:     0.00 resid: 4.02e-19
 
605
LU   natural    time:     0.00 resid: 2.17e-19
 
606
LU   amd(A+A')  time:     0.00 resid: 1.87e-19
 
607
LU   amd(S'*S)  time:     0.00 resid: 2.38e-19
 
608
LU   amd(A'*A)  time:     0.00 resid: 2.38e-19
 
609
Chol natural    time:     0.00 resid: 2.64e-19
 
610
Chol amd(A+A')  time:     0.00 resid: 2.55e-19
611
611
./cs_di_demo3 < ../Matrix/bcsstk01
612
612
 
613
613
--- Matrix: 48-by-48, nnz: 224 (sym: -1: nnz 400), norm: 3.57e+09
616
616
symbolic chol time     0.00
617
617
numeric  chol time     0.00
618
618
solve    chol time     0.00
619
 
original: resid: 2.01e-19
 
619
original: resid: 2.55e-19
620
620
update:   time:     0.00
621
 
update:   time:     0.00 (incl solve) resid: 8.39e-19
622
 
rechol:   time:     0.00 (incl solve) resid: 1.20e-18
 
621
update:   time:     0.00 (incl solve) resid: 9.66e-19
 
622
rechol:   time:     0.00 (incl solve) resid: 1.55e-18
623
623
downdate: time:     0.00
624
 
downdate: time:     0.00 (incl solve) resid: 3.30e-17
 
624
downdate: time:     0.00 (incl solve) resid: 3.74e-17
625
625
./cs_di_demo2 < ../Matrix/bcsstk16
626
626
 
627
627
--- Matrix: 4884-by-4884, nnz: 147631 (sym: -1: nnz 290378), norm: 7.01e+09
628
628
blocks: 75 singletons: 74 structural rank: 4884
629
 
QR   amd(A'*A)  time:     2.41 resid: 1.28e-22
630
 
LU   amd(A+A')  time:     1.49 resid: 1.06e-22
631
 
LU   amd(S'*S)  time:     1.40 resid: 1.21e-22
632
 
LU   amd(A'*A)  time:     1.48 resid: 1.70e-22
633
 
Chol amd(A+A')  time:     0.69 resid: 9.28e-23
 
629
QR   amd(A'*A)  time:     1.86 resid: 2.01e-22
 
630
LU   amd(A+A')  time:     1.12 resid: 1.10e-22
 
631
LU   amd(S'*S)  time:     1.12 resid: 1.28e-22
 
632
LU   amd(A'*A)  time:     1.16 resid: 1.78e-22
 
633
Chol amd(A+A')  time:     0.37 resid: 1.19e-22
634
634
./cs_di_demo3 < ../Matrix/bcsstk16
635
635
 
636
636
--- Matrix: 4884-by-4884, nnz: 147631 (sym: -1: nnz 290378), norm: 7.01e+09
637
637
 
638
638
chol then update/downdate amd(A+A')  
639
639
symbolic chol time     0.03
640
 
numeric  chol time     0.61
 
640
numeric  chol time     0.33
641
641
solve    chol time     0.01
642
 
original: resid: 9.28e-23
643
 
update:   time:     0.00
644
 
update:   time:     0.01 (incl solve) resid: 8.62e-24
645
 
rechol:   time:     0.62 (incl solve) resid: 8.72e-24
 
642
original: resid: 1.19e-22
 
643
update:   time:     0.01
 
644
update:   time:     0.01 (incl solve) resid: 1.12e-23
 
645
rechol:   time:     0.34 (incl solve) resid: 1.17e-23
646
646
downdate: time:     0.00
647
 
downdate: time:     0.01 (incl solve) resid: 3.60e-22
 
647
downdate: time:     0.01 (incl solve) resid: 4.09e-22
648
648
./cs_dl_demo1 < ../Matrix/t1
649
649
T:
650
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
650
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
651
651
triplet: 4-by-4, nzmax: 16 nnz: 10
652
652
    2 2 : 3
653
653
    1 0 : 3.1
660
660
    0 0 : 4.5
661
661
    2 1 : 1.7
662
662
A:
663
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
663
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
664
664
4-by-4, nzmax: 10 nnz: 10, 1-norm: 11.1
665
665
    col 0 : locations 0 to 2
666
666
      1 : 3.1
677
677
      3 : 1
678
678
      1 : 0.9
679
679
AT:
680
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
680
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
681
681
4-by-4, nzmax: 10 nnz: 10, 1-norm: 7.7
682
682
    col 0 : locations 0 to 1
683
683
      0 : 4.5
694
694
      1 : 0.4
695
695
      3 : 1
696
696
D:
697
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
697
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
698
698
4-by-4, nzmax: 16 nnz: 16, 1-norm: 139.58
699
699
    col 0 : locations 0 to 3
700
700
      1 : 13.95
720
720
 
721
721
--- Matrix: 4-by-4, nnz: 10 (sym: 0: nnz 0), norm: 1.11e+01
722
722
blocks: 1 singletons: 0 structural rank: 4
723
 
QR   natural    time:     0.00 resid: 1.04e-17
724
 
QR   amd(A'*A)  time:     0.00 resid: 2.03e-17
725
 
LU   natural    time:     0.00 resid: 1.04e-17
726
 
LU   amd(A+A')  time:     0.00 resid: 4.94e-18
727
 
LU   amd(S'*S)  time:     0.00 resid: 4.94e-18
728
 
LU   amd(A'*A)  time:     0.00 resid: 4.94e-18
 
723
QR   natural    time:     0.00 resid: 3.06e-17
 
724
QR   amd(A'*A)  time:     0.00 resid: 2.30e-17
 
725
LU   natural    time:     0.00 resid: 1.53e-17
 
726
LU   amd(A+A')  time:     0.00 resid: 1.53e-17
 
727
LU   amd(S'*S)  time:     0.00 resid: 0.00e+00
 
728
LU   amd(A'*A)  time:     0.00 resid: 1.53e-17
729
729
./cs_dl_demo2 < ../Matrix/fs_183_1
730
730
 
731
731
--- Matrix: 183-by-183, nnz: 988 (sym: 0: nnz 0), norm: 1.70e+09
732
732
zero entries dropped: 71
733
733
tiny entries dropped: 10
734
734
blocks: 38 singletons: 37 structural rank: 183
735
 
QR   natural    time:     0.00 resid: 3.15e-28
736
 
QR   amd(A'*A)  time:     0.01 resid: 2.60e-28
737
 
LU   natural    time:     0.00 resid: 3.08e-28
738
 
LU   amd(A+A')  time:     0.00 resid: 1.42e-27
739
 
LU   amd(S'*S)  time:     0.00 resid: 7.11e-28
740
 
LU   amd(A'*A)  time:     0.00 resid: 7.11e-28
 
735
QR   natural    time:     0.00 resid: 1.42e-27
 
736
QR   amd(A'*A)  time:     0.01 resid: 3.35e-28
 
737
LU   natural    time:     0.00 resid: 6.20e-28
 
738
LU   amd(A+A')  time:     0.00 resid: 1.55e-27
 
739
LU   amd(S'*S)  time:     0.00 resid: 6.98e-28
 
740
LU   amd(A'*A)  time:     0.00 resid: 6.98e-28
741
741
./cs_dl_demo2 < ../Matrix/west0067
742
742
 
743
743
--- Matrix: 67-by-67, nnz: 294 (sym: 0: nnz 0), norm: 6.14e+00
744
744
blocks: 2 singletons: 1 structural rank: 67
745
 
QR   natural    time:     0.00 resid: 2.81e-17
746
 
QR   amd(A'*A)  time:     0.00 resid: 1.30e-17
747
 
LU   natural    time:     0.00 resid: 3.85e-17
748
 
LU   amd(A+A')  time:     0.00 resid: 1.95e-17
749
 
LU   amd(S'*S)  time:     0.00 resid: 2.60e-17
 
745
QR   natural    time:     0.00 resid: 5.19e-17
 
746
QR   amd(A'*A)  time:     0.00 resid: 3.25e-17
 
747
LU   natural    time:     0.00 resid: 3.89e-17
 
748
LU   amd(A+A')  time:     0.00 resid: 2.27e-17
 
749
LU   amd(S'*S)  time:     0.00 resid: 1.95e-17
750
750
LU   amd(A'*A)  time:     0.00 resid: 2.60e-17
751
751
./cs_dl_demo2 < ../Matrix/lp_afiro
752
752
 
753
753
--- Matrix: 27-by-51, nnz: 102 (sym: 0: nnz 0), norm: 3.43e+00
754
754
blocks: 1 singletons: 0 structural rank: 27
755
 
QR   natural    time:     0.00 resid: 4.30e-16
756
 
QR   amd(A'*A)  time:     0.00 resid: 1.92e-16
 
755
QR   natural    time:     0.00 resid: 3.85e-16
 
756
QR   amd(A'*A)  time:     0.00 resid: 1.50e-16
757
757
./cs_dl_demo2 < ../Matrix/ash219
758
758
 
759
759
--- Matrix: 219-by-85, nnz: 438 (sym: 0: nnz 0), norm: 9.00e+00
764
764
 
765
765
--- Matrix: 492-by-490, nnz: 49920 (sym: 0: nnz 0), norm: 9.29e-01
766
766
blocks: 10 singletons: 8 structural rank: 448
767
 
QR   natural    time:     0.26 resid:      nan
768
 
QR   amd(A'*A)  time:     0.33 resid:      nan
 
767
QR   natural    time:     0.13 resid:      nan
 
768
QR   amd(A'*A)  time:     0.20 resid:      nan
769
769
./cs_dl_demo2 < ../Matrix/bcsstk01
770
770
 
771
771
--- Matrix: 48-by-48, nnz: 224 (sym: -1: nnz 400), norm: 3.57e+09
772
772
blocks: 1 singletons: 0 structural rank: 48
773
 
QR   natural    time:     0.00 resid: 3.97e-19
774
 
QR   amd(A'*A)  time:     0.00 resid: 2.36e-19
775
 
LU   natural    time:     0.00 resid: 2.63e-19
776
 
LU   amd(A+A')  time:     0.00 resid: 8.63e-20
777
 
LU   amd(S'*S)  time:     0.00 resid: 2.04e-19
778
 
LU   amd(A'*A)  time:     0.00 resid: 2.04e-19
779
 
Chol natural    time:     0.00 resid: 1.90e-19
780
 
Chol amd(A+A')  time:     0.00 resid: 2.01e-19
 
773
QR   natural    time:     0.00 resid: 3.65e-19
 
774
QR   amd(A'*A)  time:     0.00 resid: 4.02e-19
 
775
LU   natural    time:     0.00 resid: 2.17e-19
 
776
LU   amd(A+A')  time:     0.00 resid: 1.87e-19
 
777
LU   amd(S'*S)  time:     0.00 resid: 2.38e-19
 
778
LU   amd(A'*A)  time:     0.00 resid: 2.38e-19
 
779
Chol natural    time:     0.00 resid: 2.64e-19
 
780
Chol amd(A+A')  time:     0.00 resid: 2.55e-19
781
781
./cs_dl_demo3 < ../Matrix/bcsstk01
782
782
 
783
783
--- Matrix: 48-by-48, nnz: 224 (sym: -1: nnz 400), norm: 3.57e+09
786
786
symbolic chol time     0.00
787
787
numeric  chol time     0.00
788
788
solve    chol time     0.00
789
 
original: resid: 2.01e-19
 
789
original: resid: 2.55e-19
790
790
update:   time:     0.00
791
 
update:   time:     0.00 (incl solve) resid: 8.39e-19
792
 
rechol:   time:     0.00 (incl solve) resid: 1.20e-18
 
791
update:   time:     0.00 (incl solve) resid: 9.66e-19
 
792
rechol:   time:     0.00 (incl solve) resid: 1.55e-18
793
793
downdate: time:     0.00
794
 
downdate: time:     0.00 (incl solve) resid: 3.30e-17
 
794
downdate: time:     0.00 (incl solve) resid: 3.74e-17
795
795
./cs_dl_demo2 < ../Matrix/bcsstk16
796
796
 
797
797
--- Matrix: 4884-by-4884, nnz: 147631 (sym: -1: nnz 290378), norm: 7.01e+09
798
798
blocks: 75 singletons: 74 structural rank: 4884
799
 
QR   amd(A'*A)  time:     2.42 resid: 1.28e-22
800
 
LU   amd(A+A')  time:     1.39 resid: 1.06e-22
801
 
LU   amd(S'*S)  time:     1.37 resid: 1.21e-22
802
 
LU   amd(A'*A)  time:     1.41 resid: 1.70e-22
803
 
Chol amd(A+A')  time:     0.65 resid: 9.28e-23
 
799
QR   amd(A'*A)  time:     2.06 resid: 2.01e-22
 
800
LU   amd(A+A')  time:     1.15 resid: 1.10e-22
 
801
LU   amd(S'*S)  time:     1.10 resid: 1.28e-22
 
802
LU   amd(A'*A)  time:     1.17 resid: 1.78e-22
 
803
Chol amd(A+A')  time:     0.41 resid: 1.19e-22
804
804
./cs_dl_demo3 < ../Matrix/bcsstk16
805
805
 
806
806
--- Matrix: 4884-by-4884, nnz: 147631 (sym: -1: nnz 290378), norm: 7.01e+09
807
807
 
808
808
chol then update/downdate amd(A+A')  
809
809
symbolic chol time     0.03
810
 
numeric  chol time     0.61
811
 
solve    chol time     0.01
812
 
original: resid: 9.28e-23
 
810
numeric  chol time     0.38
 
811
solve    chol time     0.00
 
812
original: resid: 1.19e-22
813
813
update:   time:     0.00
814
 
update:   time:     0.02 (incl solve) resid: 8.62e-24
815
 
rechol:   time:     0.63 (incl solve) resid: 8.72e-24
 
814
update:   time:     0.01 (incl solve) resid: 1.12e-23
 
815
rechol:   time:     0.38 (incl solve) resid: 1.17e-23
816
816
downdate: time:     0.00
817
 
downdate: time:     0.01 (incl solve) resid: 3.60e-22
 
817
downdate: time:     0.01 (incl solve) resid: 4.09e-22
818
818
./cs_ci_demo1 < ../Matrix/t2
819
819
T:
820
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
820
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
821
821
triplet: 4-by-4, nzmax: 16 nnz: 10
822
822
    2 2 : (3, 3.14159)
823
823
    1 0 : (3.1, 42)
830
830
    0 0 : (4.5, 6)
831
831
    2 1 : (1.7, 1)
832
832
A:
833
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
833
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
834
834
4-by-4, nzmax: 10 nnz: 10, 1-norm: 106.075
835
835
    col 0 : locations 0 to 2
836
836
      1 : (3.1, 42)
847
847
      3 : (1, 7)
848
848
      1 : (0.9, 99)
849
849
AT:
850
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
850
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
851
851
4-by-4, nzmax: 10 nnz: 10, 1-norm: 144.296
852
852
    col 0 : locations 0 to 1
853
853
      0 : (4.5, -6)
864
864
      1 : (0.4, -2.71828)
865
865
      3 : (1, -7)
866
866
D:
867
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
867
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
868
868
4-by-4, nzmax: 16 nnz: 16, 1-norm: 25308.3
869
869
    col 0 : locations 0 to 3
870
870
      1 : (265.95, 170.4)
890
890
 
891
891
--- Matrix: 4-by-4, nnz: 10 (sym: 0: nnz 0), norm: 1.06e+02
892
892
blocks: 1 singletons: 0 structural rank: 4
893
 
QR   natural    time:     0.00 resid: 1.07e-17
894
 
QR   amd(A'*A)  time:     0.00 resid: 3.09e-18
895
 
LU   natural    time:     0.00 resid: 2.67e-18
896
 
LU   amd(A+A')  time:     0.00 resid: 2.12e-18
897
 
LU   amd(S'*S)  time:     0.00 resid: 2.93e-18
898
 
LU   amd(A'*A)  time:     0.00 resid: 2.12e-18
 
893
QR   natural    time:     0.00 resid: 2.06e-17
 
894
QR   amd(A'*A)  time:     0.00 resid: 6.36e-18
 
895
LU   natural    time:     0.00 resid: 4.88e-18
 
896
LU   amd(A+A')  time:     0.00 resid: 2.11e-18
 
897
LU   amd(S'*S)  time:     0.00 resid: 6.36e-18
 
898
LU   amd(A'*A)  time:     0.00 resid: 2.11e-18
899
899
./cs_ci_demo2 < ../Matrix/t3
900
900
 
901
901
--- Matrix: 3-by-4, nnz: 12 (sym: 0: nnz 0), norm: 3.06e+00
902
902
blocks: 1 singletons: 0 structural rank: 3
903
 
QR   natural    time:     0.00 resid: 1.21e-16
904
 
QR   amd(A'*A)  time:     0.00 resid: 1.21e-16
 
903
QR   natural    time:     0.00 resid: 1.05e-16
 
904
QR   amd(A'*A)  time:     0.00 resid: 1.05e-16
905
905
./cs_ci_demo2 < ../Matrix/t4
906
906
 
907
907
--- Matrix: 2-by-2, nnz: 3 (sym: 1: nnz 4), norm: 2.83e+00
908
908
blocks: 1 singletons: 0 structural rank: 2
909
 
QR   natural    time:     0.00 resid: 2.74e-17
910
 
QR   amd(A'*A)  time:     0.00 resid: 2.74e-17
 
909
QR   natural    time:     0.00 resid: 5.65e-17
 
910
QR   amd(A'*A)  time:     0.00 resid: 5.65e-17
911
911
LU   natural    time:     0.00 resid: 0.00e+00
912
912
LU   amd(A+A')  time:     0.00 resid: 0.00e+00
913
913
LU   amd(S'*S)  time:     0.00 resid: 0.00e+00
918
918
 
919
919
--- Matrix: 67-by-67, nnz: 294 (sym: 0: nnz 0), norm: 6.17e+00
920
920
blocks: 2 singletons: 1 structural rank: 67
921
 
QR   natural    time:     0.00 resid: 3.47e-17
922
 
QR   amd(A'*A)  time:     0.00 resid: 6.79e-17
923
 
LU   natural    time:     0.00 resid: 5.08e-17
924
 
LU   amd(A+A')  time:     0.00 resid: 7.10e-17
925
 
LU   amd(S'*S)  time:     0.00 resid: 3.28e-17
926
 
LU   amd(A'*A)  time:     0.00 resid: 4.32e-17
 
921
QR   natural    time:     0.00 resid: 8.16e-17
 
922
QR   amd(A'*A)  time:     0.00 resid: 5.34e-17
 
923
LU   natural    time:     0.00 resid: 4.54e-17
 
924
LU   amd(A+A')  time:     0.00 resid: 5.90e-17
 
925
LU   amd(S'*S)  time:     0.00 resid: 5.28e-17
 
926
LU   amd(A'*A)  time:     0.00 resid: 4.76e-17
927
927
./cs_ci_demo2 < ../Matrix/c_mbeacxc
928
928
 
929
929
--- Matrix: 492-by-490, nnz: 49920 (sym: 0: nnz 0), norm: 9.29e-01
930
930
blocks: 10 singletons: 8 structural rank: 448
931
 
QR   natural    time:     0.62 resid:      nan
932
 
QR   amd(A'*A)  time:     0.68 resid:      nan
 
931
QR   natural    time:     0.30 resid:      nan
 
932
QR   amd(A'*A)  time:     0.37 resid:      nan
933
933
./cs_ci_demo2 < ../Matrix/young1c
934
934
 
935
935
--- Matrix: 841-by-841, nnz: 4089 (sym: 0: nnz 0), norm: 7.30e+02
936
936
blocks: 1 singletons: 0 structural rank: 841
937
 
QR   natural    time:     0.02 resid: 9.10e-17
938
 
QR   amd(A'*A)  time:     0.02 resid: 8.11e-17
939
 
LU   natural    time:     0.01 resid: 1.27e-16
940
 
LU   amd(A+A')  time:     0.02 resid: 3.45e-16
941
 
LU   amd(S'*S)  time:     0.02 resid: 2.31e-16
942
 
LU   amd(A'*A)  time:     0.01 resid: 2.31e-16
 
937
QR   natural    time:     0.02 resid: 1.51e-16
 
938
QR   amd(A'*A)  time:     0.01 resid: 1.53e-16
 
939
LU   natural    time:     0.01 resid: 1.39e-16
 
940
LU   amd(A+A')  time:     0.01 resid: 2.95e-16
 
941
LU   amd(S'*S)  time:     0.01 resid: 3.37e-16
 
942
LU   amd(A'*A)  time:     0.01 resid: 3.37e-16
943
943
./cs_ci_demo2 < ../Matrix/qc324
944
944
 
945
945
--- Matrix: 324-by-324, nnz: 26730 (sym: 0: nnz 0), norm: 1.71e+00
946
946
blocks: 1 singletons: 0 structural rank: 324
947
 
QR   natural    time:     0.04 resid: 3.27e-17
948
 
QR   amd(A'*A)  time:     0.06 resid: 3.70e-17
949
 
LU   natural    time:     0.03 resid: 4.19e-17
950
 
LU   amd(A+A')  time:     0.04 resid: 3.73e-17
951
 
LU   amd(S'*S)  time:     0.03 resid: 4.19e-17
952
 
LU   amd(A'*A)  time:     0.04 resid: 4.19e-17
 
947
QR   natural    time:     0.04 resid: 9.08e-17
 
948
QR   amd(A'*A)  time:     0.05 resid: 8.71e-17
 
949
LU   natural    time:     0.02 resid: 6.01e-17
 
950
LU   amd(A+A')  time:     0.03 resid: 4.05e-17
 
951
LU   amd(S'*S)  time:     0.03 resid: 4.71e-17
 
952
LU   amd(A'*A)  time:     0.03 resid: 4.71e-17
953
953
./cs_ci_demo2 < ../Matrix/neumann
954
954
 
955
955
--- Matrix: 1600-by-1600, nnz: 7840 (sym: 0: nnz 0), norm: 1.41e+01
956
956
blocks: 1 singletons: 0 structural rank: 1600
957
 
QR   amd(A'*A)  time:     0.03 resid: 4.89e-16
958
 
LU   amd(A+A')  time:     0.02 resid: 4.32e-16
959
 
LU   amd(S'*S)  time:     0.02 resid: 4.63e-16
960
 
LU   amd(A'*A)  time:     0.02 resid: 4.63e-16
 
957
QR   amd(A'*A)  time:     0.04 resid: 1.04e-15
 
958
LU   amd(A+A')  time:     0.00 resid: 3.55e-16
 
959
LU   amd(S'*S)  time:     0.02 resid: 4.03e-16
 
960
LU   amd(A'*A)  time:     0.02 resid: 4.03e-16
 
961
./cs_ci_demo2 < ../Matrix/c4
 
962
 
 
963
--- Matrix: 4-by-4, nnz: 10 (sym: -1: nnz 16), norm: 7.37e+01
 
964
blocks: 1 singletons: 0 structural rank: 4
 
965
QR   natural    time:     0.00 resid: 5.85e-17
 
966
QR   amd(A'*A)  time:     0.00 resid: 5.85e-17
 
967
LU   natural    time:     0.00 resid: 2.29e-17
 
968
LU   amd(A+A')  time:     0.00 resid: 2.29e-17
 
969
LU   amd(S'*S)  time:     0.00 resid: 2.29e-17
 
970
LU   amd(A'*A)  time:     0.00 resid: 2.29e-17
 
971
Chol natural    time:     0.00 resid: 6.88e-17
 
972
Chol amd(A+A')  time:     0.00 resid: 6.88e-17
 
973
./cs_ci_demo3 < ../Matrix/c4
 
974
 
 
975
--- Matrix: 4-by-4, nnz: 10 (sym: -1: nnz 16), norm: 7.37e+01
 
976
 
 
977
chol then update/downdate amd(A+A')  
 
978
symbolic chol time     0.00
 
979
numeric  chol time     0.00
 
980
solve    chol time     0.00
 
981
original: resid: 6.88e-17
 
982
update:   time:     0.00
 
983
update:   time:     0.00 (incl solve) resid: 6.49e-17
 
984
rechol:   time:     0.00 (incl solve) resid: 6.49e-17
 
985
downdate: time:     0.00
 
986
downdate: time:     0.00 (incl solve) resid: 5.85e-17
961
987
./cs_ci_demo2 < ../Matrix/mhd1280b
962
988
 
963
989
--- Matrix: 1280-by-1280, nnz: 11963 (sym: -1: nnz 22646), norm: 8.00e+01
964
990
tiny entries dropped: 66
965
991
blocks: 20 singletons: 14 structural rank: 1280
966
 
QR   amd(A'*A)  time:     0.01 resid: 1.98e-25
967
 
LU   amd(A+A')  time:     0.00 resid: 2.07e-25
968
 
LU   amd(S'*S)  time:     0.01 resid: 1.53e-25
969
 
LU   amd(A'*A)  time:     0.02 resid: 1.53e-25
970
 
Chol amd(A+A')  time:     0.01 resid: 3.63e-16
 
992
QR   amd(A'*A)  time:     0.01 resid: 6.15e-25
 
993
LU   amd(A+A')  time:     0.00 resid: 2.33e-25
 
994
LU   amd(S'*S)  time:     0.00 resid: 3.96e-25
 
995
LU   amd(A'*A)  time:     0.01 resid: 3.96e-25
 
996
Chol amd(A+A')  time:     0.01 resid: 1.58e-25
971
997
./cs_ci_demo3 < ../Matrix/mhd1280b
972
998
 
973
999
--- Matrix: 1280-by-1280, nnz: 12029 (sym: -1: nnz 22778), norm: 8.00e+01
974
1000
 
975
1001
chol then update/downdate amd(A+A')  
976
 
symbolic chol time     0.01
977
 
numeric  chol time     0.00
 
1002
symbolic chol time     0.00
 
1003
numeric  chol time     0.01
978
1004
solve    chol time     0.00
979
 
original: resid: 3.64e-16
 
1005
original: resid: 1.51e-25
980
1006
update:   time:     0.00
981
 
update:   time:     0.00 (incl solve) resid: 3.64e-16
982
 
rechol:   time:     0.01 (incl solve) resid: 3.64e-16
 
1007
update:   time:     0.00 (incl solve) resid: 1.75e-25
 
1008
rechol:   time:     0.00 (incl solve) resid: 1.71e-25
983
1009
downdate: time:     0.00
984
 
downdate: time:     0.00 (incl solve) resid: 1.89e-15
 
1010
downdate: time:     0.00 (incl solve) resid: 5.85e-25
985
1011
./cs_cl_demo1 < ../Matrix/t2
986
1012
T:
987
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
1013
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
988
1014
triplet: 4-by-4, nzmax: 16 nnz: 10
989
1015
    2 2 : (3, 3.14159)
990
1016
    1 0 : (3.1, 42)
997
1023
    0 0 : (4.5, 6)
998
1024
    2 1 : (1.7, 1)
999
1025
A:
1000
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
1026
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
1001
1027
4-by-4, nzmax: 10 nnz: 10, 1-norm: 106.075
1002
1028
    col 0 : locations 0 to 2
1003
1029
      1 : (3.1, 42)
1014
1040
      3 : (1, 7)
1015
1041
      1 : (0.9, 99)
1016
1042
AT:
1017
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
1043
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
1018
1044
4-by-4, nzmax: 10 nnz: 10, 1-norm: 144.296
1019
1045
    col 0 : locations 0 to 1
1020
1046
      0 : (4.5, -6)
1031
1057
      1 : (0.4, -2.71828)
1032
1058
      3 : (1, -7)
1033
1059
D:
1034
 
CXSparse Version 2.0.6, Dec 7, 2006.  Copyright (c) Timothy A. Davis, 2006
 
1060
CXSparse Version 2.2.0, May 31, 2007.  Copyright (c) Timothy A. Davis, 2006-2007
1035
1061
4-by-4, nzmax: 16 nnz: 16, 1-norm: 25308.3
1036
1062
    col 0 : locations 0 to 3
1037
1063
      1 : (265.95, 170.4)
1057
1083
 
1058
1084
--- Matrix: 4-by-4, nnz: 10 (sym: 0: nnz 0), norm: 1.06e+02
1059
1085
blocks: 1 singletons: 0 structural rank: 4
1060
 
QR   natural    time:     0.00 resid: 1.07e-17
1061
 
QR   amd(A'*A)  time:     0.00 resid: 3.09e-18
1062
 
LU   natural    time:     0.00 resid: 2.67e-18
1063
 
LU   amd(A+A')  time:     0.00 resid: 2.12e-18
1064
 
LU   amd(S'*S)  time:     0.00 resid: 2.93e-18
1065
 
LU   amd(A'*A)  time:     0.00 resid: 2.12e-18
 
1086
QR   natural    time:     0.00 resid: 2.06e-17
 
1087
QR   amd(A'*A)  time:     0.00 resid: 6.36e-18
 
1088
LU   natural    time:     0.00 resid: 4.88e-18
 
1089
LU   amd(A+A')  time:     0.00 resid: 2.11e-18
 
1090
LU   amd(S'*S)  time:     0.00 resid: 6.36e-18
 
1091
LU   amd(A'*A)  time:     0.00 resid: 2.11e-18
1066
1092
./cs_cl_demo2 < ../Matrix/t3
1067
1093
 
1068
1094
--- Matrix: 3-by-4, nnz: 12 (sym: 0: nnz 0), norm: 3.06e+00
1069
1095
blocks: 1 singletons: 0 structural rank: 3
1070
 
QR   natural    time:     0.00 resid: 1.21e-16
1071
 
QR   amd(A'*A)  time:     0.00 resid: 1.21e-16
 
1096
QR   natural    time:     0.00 resid: 1.05e-16
 
1097
QR   amd(A'*A)  time:     0.00 resid: 1.05e-16
1072
1098
./cs_cl_demo2 < ../Matrix/t4
1073
1099
 
1074
1100
--- Matrix: 2-by-2, nnz: 3 (sym: 1: nnz 4), norm: 2.83e+00
1075
1101
blocks: 1 singletons: 0 structural rank: 2
1076
 
QR   natural    time:     0.00 resid: 2.74e-17
1077
 
QR   amd(A'*A)  time:     0.00 resid: 2.74e-17
 
1102
QR   natural    time:     0.00 resid: 5.65e-17
 
1103
QR   amd(A'*A)  time:     0.00 resid: 5.65e-17
1078
1104
LU   natural    time:     0.00 resid: 0.00e+00
1079
1105
LU   amd(A+A')  time:     0.00 resid: 0.00e+00
1080
1106
LU   amd(S'*S)  time:     0.00 resid: 0.00e+00
1085
1111
 
1086
1112
--- Matrix: 67-by-67, nnz: 294 (sym: 0: nnz 0), norm: 6.17e+00
1087
1113
blocks: 2 singletons: 1 structural rank: 67
1088
 
QR   natural    time:     0.00 resid: 3.47e-17
1089
 
QR   amd(A'*A)  time:     0.00 resid: 6.79e-17
1090
 
LU   natural    time:     0.00 resid: 5.08e-17
1091
 
LU   amd(A+A')  time:     0.00 resid: 7.10e-17
1092
 
LU   amd(S'*S)  time:     0.00 resid: 3.28e-17
1093
 
LU   amd(A'*A)  time:     0.00 resid: 4.32e-17
 
1114
QR   natural    time:     0.00 resid: 8.16e-17
 
1115
QR   amd(A'*A)  time:     0.00 resid: 5.34e-17
 
1116
LU   natural    time:     0.00 resid: 4.54e-17
 
1117
LU   amd(A+A')  time:     0.00 resid: 5.90e-17
 
1118
LU   amd(S'*S)  time:     0.00 resid: 5.28e-17
 
1119
LU   amd(A'*A)  time:     0.00 resid: 4.76e-17
1094
1120
./cs_cl_demo2 < ../Matrix/c_mbeacxc
1095
1121
 
1096
1122
--- Matrix: 492-by-490, nnz: 49920 (sym: 0: nnz 0), norm: 9.29e-01
1097
1123
blocks: 10 singletons: 8 structural rank: 448
1098
 
QR   natural    time:     0.63 resid:      nan
1099
 
QR   amd(A'*A)  time:     0.67 resid:      nan
 
1124
QR   natural    time:     0.30 resid:      nan
 
1125
QR   amd(A'*A)  time:     0.35 resid:      nan
1100
1126
./cs_cl_demo2 < ../Matrix/young1c
1101
1127
 
1102
1128
--- Matrix: 841-by-841, nnz: 4089 (sym: 0: nnz 0), norm: 7.30e+02
1103
1129
blocks: 1 singletons: 0 structural rank: 841
1104
 
QR   natural    time:     0.03 resid: 9.10e-17
1105
 
QR   amd(A'*A)  time:     0.01 resid: 8.11e-17
1106
 
LU   natural    time:     0.02 resid: 1.27e-16
1107
 
LU   amd(A+A')  time:     0.02 resid: 3.45e-16
1108
 
LU   amd(S'*S)  time:     0.01 resid: 2.31e-16
1109
 
LU   amd(A'*A)  time:     0.01 resid: 2.31e-16
 
1130
QR   natural    time:     0.01 resid: 1.51e-16
 
1131
QR   amd(A'*A)  time:     0.00 resid: 1.53e-16
 
1132
LU   natural    time:     0.02 resid: 1.39e-16
 
1133
LU   amd(A+A')  time:     0.01 resid: 2.95e-16
 
1134
LU   amd(S'*S)  time:     0.01 resid: 3.37e-16
 
1135
LU   amd(A'*A)  time:     0.01 resid: 3.37e-16
1110
1136
./cs_cl_demo2 < ../Matrix/qc324
1111
1137
 
1112
1138
--- Matrix: 324-by-324, nnz: 26730 (sym: 0: nnz 0), norm: 1.71e+00
1113
1139
blocks: 1 singletons: 0 structural rank: 324
1114
 
QR   natural    time:     0.04 resid: 3.27e-17
1115
 
QR   amd(A'*A)  time:     0.06 resid: 3.70e-17
1116
 
LU   natural    time:     0.03 resid: 4.19e-17
1117
 
LU   amd(A+A')  time:     0.03 resid: 3.73e-17
1118
 
LU   amd(S'*S)  time:     0.04 resid: 4.19e-17
1119
 
LU   amd(A'*A)  time:     0.03 resid: 4.19e-17
 
1140
QR   natural    time:     0.03 resid: 9.08e-17
 
1141
QR   amd(A'*A)  time:     0.04 resid: 8.71e-17
 
1142
LU   natural    time:     0.02 resid: 6.01e-17
 
1143
LU   amd(A+A')  time:     0.02 resid: 4.05e-17
 
1144
LU   amd(S'*S)  time:     0.04 resid: 4.71e-17
 
1145
LU   amd(A'*A)  time:     0.03 resid: 4.71e-17
1120
1146
./cs_cl_demo2 < ../Matrix/neumann
1121
1147
 
1122
1148
--- Matrix: 1600-by-1600, nnz: 7840 (sym: 0: nnz 0), norm: 1.41e+01
1123
1149
blocks: 1 singletons: 0 structural rank: 1600
1124
 
QR   amd(A'*A)  time:     0.03 resid: 4.89e-16
1125
 
LU   amd(A+A')  time:     0.02 resid: 4.32e-16
1126
 
LU   amd(S'*S)  time:     0.02 resid: 4.63e-16
1127
 
LU   amd(A'*A)  time:     0.03 resid: 4.63e-16
 
1150
QR   amd(A'*A)  time:     0.03 resid: 1.04e-15
 
1151
LU   amd(A+A')  time:     0.01 resid: 3.55e-16
 
1152
LU   amd(S'*S)  time:     0.02 resid: 4.03e-16
 
1153
LU   amd(A'*A)  time:     0.02 resid: 4.03e-16
 
1154
./cs_cl_demo2 < ../Matrix/c4
 
1155
 
 
1156
--- Matrix: 4-by-4, nnz: 10 (sym: -1: nnz 16), norm: 7.37e+01
 
1157
blocks: 1 singletons: 0 structural rank: 4
 
1158
QR   natural    time:     0.00 resid: 5.85e-17
 
1159
QR   amd(A'*A)  time:     0.00 resid: 5.85e-17
 
1160
LU   natural    time:     0.00 resid: 2.29e-17
 
1161
LU   amd(A+A')  time:     0.00 resid: 2.29e-17
 
1162
LU   amd(S'*S)  time:     0.00 resid: 2.29e-17
 
1163
LU   amd(A'*A)  time:     0.00 resid: 2.29e-17
 
1164
Chol natural    time:     0.00 resid: 6.88e-17
 
1165
Chol amd(A+A')  time:     0.00 resid: 6.88e-17
 
1166
./cs_cl_demo3 < ../Matrix/c4
 
1167
 
 
1168
--- Matrix: 4-by-4, nnz: 10 (sym: -1: nnz 16), norm: 7.37e+01
 
1169
 
 
1170
chol then update/downdate amd(A+A')  
 
1171
symbolic chol time     0.00
 
1172
numeric  chol time     0.00
 
1173
solve    chol time     0.00
 
1174
original: resid: 6.88e-17
 
1175
update:   time:     0.00
 
1176
update:   time:     0.00 (incl solve) resid: 6.49e-17
 
1177
rechol:   time:     0.00 (incl solve) resid: 6.49e-17
 
1178
downdate: time:     0.00
 
1179
downdate: time:     0.00 (incl solve) resid: 5.85e-17
1128
1180
./cs_cl_demo2 < ../Matrix/mhd1280b
1129
1181
 
1130
1182
--- Matrix: 1280-by-1280, nnz: 11963 (sym: -1: nnz 22646), norm: 8.00e+01
1131
1183
tiny entries dropped: 66
1132
1184
blocks: 20 singletons: 14 structural rank: 1280
1133
 
QR   amd(A'*A)  time:     0.01 resid: 1.98e-25
1134
 
LU   amd(A+A')  time:     0.01 resid: 2.07e-25
1135
 
LU   amd(S'*S)  time:     0.01 resid: 1.53e-25
1136
 
LU   amd(A'*A)  time:     0.00 resid: 1.53e-25
1137
 
Chol amd(A+A')  time:     0.01 resid: 3.63e-16
 
1185
QR   amd(A'*A)  time:     0.01 resid: 6.15e-25
 
1186
LU   amd(A+A')  time:     0.01 resid: 2.33e-25
 
1187
LU   amd(S'*S)  time:     0.01 resid: 3.96e-25
 
1188
LU   amd(A'*A)  time:     0.00 resid: 3.96e-25
 
1189
Chol amd(A+A')  time:     0.01 resid: 1.58e-25
1138
1190
./cs_cl_demo3 < ../Matrix/mhd1280b
1139
1191
 
1140
1192
--- Matrix: 1280-by-1280, nnz: 12029 (sym: -1: nnz 22778), norm: 8.00e+01
1141
1193
 
1142
1194
chol then update/downdate amd(A+A')  
1143
 
symbolic chol time     0.01
 
1195
symbolic chol time     0.00
1144
1196
numeric  chol time     0.00
1145
1197
solve    chol time     0.00
1146
 
original: resid: 3.64e-16
 
1198
original: resid: 1.51e-25
1147
1199
update:   time:     0.00
1148
 
update:   time:     0.00 (incl solve) resid: 3.64e-16
1149
 
rechol:   time:     0.00 (incl solve) resid: 3.64e-16
 
1200
update:   time:     0.00 (incl solve) resid: 1.75e-25
 
1201
rechol:   time:     0.00 (incl solve) resid: 1.71e-25
1150
1202
downdate: time:     0.00
1151
 
downdate: time:     0.00 (incl solve) resid: 1.89e-15
1152
 
make[1]: Leaving directory `/amd/netapp3/vol/research0a/research18/sparse/SuiteSparse/CXSparse/Demo'
 
1203
downdate: time:     0.00 (incl solve) resid: 5.85e-25
 
1204
make[1]: Leaving directory `/amd/birch/export/research07/sparse/SuiteSparse/CXSparse/Demo'