1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
|
/*
* CXSPARSE: a Concise Sparse Matrix package - Extended.
* Copyright (c) 2006-2009, Timothy A. Davis.
* http://www.cise.ufl.edu/research/sparse/CXSparse
*
* CXSparse is free software; you can redistribute it and/or
* modify it under the terms of the GNU Lesser General Public
* License as published by the Free Software Foundation; either
* version 2.1 of the License, or (at your option) any later version.
*
* CXSparse is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public
* License along with this Module; if not, write to the Free Software
* Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
*/
#include "cs.h"
/* print a sparse matrix */
/* CS_INT cs_print (const cs *A, CS_INT brief) */
/* { */
/* CS_INT p, j, m, n, nzmax, nz, *Ap, *Ai ; */
/* CS_ENTRY *Ax ; */
/* if (!A) { printf ("(null)\n") ; return (0) ; } */
/* m = A->m ; n = A->n ; Ap = A->p ; Ai = A->i ; Ax = A->x ; */
/* nzmax = A->nzmax ; nz = A->nz ; */
/* printf ("CXSparse Version %d.%d.%d, %s. %s\n", CS_VER, CS_SUBVER, */
/* CS_SUBSUB, CS_DATE, CS_COPYRIGHT) ; */
/* if (nz < 0) */
/* { */
/* printf (""CS_ID"-by-"CS_ID", nzmax: "CS_ID" nnz: "CS_ID", 1-norm: %g\n", m, n, nzmax, */
/* Ap [n], cs_norm (A)) ; */
/* for (j = 0 ; j < n ; j++) */
/* { */
/* printf (" col "CS_ID" : locations "CS_ID" to "CS_ID"\n", j, Ap [j], Ap [j+1]-1); */
/* for (p = Ap [j] ; p < Ap [j+1] ; p++) */
/* { */
/* #ifdef CS_COMPLEX */
/* printf (" "CS_ID" : (%g, %g)\n", Ai [p], */
/* Ax ? CS_REAL (Ax [p]) : 1, Ax ? CS_IMAG (Ax [p]) : 0) ; */
/* #else */
/* printf (" "CS_ID" : %g\n", Ai [p], Ax ? Ax [p] : 1) ; */
/* #endif */
/* if (brief && p > 20) { printf (" ...\n") ; return (1) ; } */
/* } */
/* } */
/* } */
/* else */
/* { */
/* printf ("triplet: "CS_ID"-by-"CS_ID", nzmax: "CS_ID" nnz: "CS_ID"\n", m, n, nzmax, nz) ; */
/* for (p = 0 ; p < nz ; p++) */
/* { */
/* #ifdef CS_COMPLEX */
/* printf (" "CS_ID" "CS_ID" : (%g, %g)\n", Ai [p], Ap [p], */
/* Ax ? CS_REAL (Ax [p]) : 1, Ax ? CS_IMAG (Ax [p]) : 0) ; */
/* #else */
/* printf (" "CS_ID" "CS_ID" : %g\n", Ai [p], Ap [p], Ax ? Ax [p] : 1) ; */
/* #endif */
/* if (brief && p > 20) { printf (" ...\n") ; return (1) ; } */
/* } */
/* } */
/* return (1) ; */
/* } */
|
