1
/* ========================================================================== */
2
/* === umfpack_qsymbolic ==================================================== */
3
/* ========================================================================== */
5
/* -------------------------------------------------------------------------- */
6
/* UMFPACK Version 4.1 (Apr. 30, 2003), Copyright (c) 2003 by Timothy A. */
7
/* Davis. All Rights Reserved. See ../README for License. */
8
/* email: davis@cise.ufl.edu CISE Department, Univ. of Florida. */
9
/* web: http://www.cise.ufl.edu/research/sparse/umfpack */
10
/* -------------------------------------------------------------------------- */
12
int umfpack_di_qsymbolic
21
const double Control [UMFPACK_CONTROL],
22
double Info [UMFPACK_INFO]
25
long umfpack_dl_qsymbolic
34
const double Control [UMFPACK_CONTROL],
35
double Info [UMFPACK_INFO]
38
int umfpack_zi_qsymbolic
44
const double Ax [ ], const double Az [ ],
47
const double Control [UMFPACK_CONTROL],
48
double Info [UMFPACK_INFO]
51
long umfpack_zl_qsymbolic
57
const double Ax [ ], const double Az [ ],
60
const double Control [UMFPACK_CONTROL],
61
double Info [UMFPACK_INFO]
69
int n_row, n_col, *Ap, *Ai, *Qinit, status ;
70
double Control [UMFPACK_CONTROL], Info [UMFPACK_INFO], *Ax ;
71
status = umfpack_di_qsymbolic (n_row, n_col, Ap, Ai, Ax, Qinit,
72
&Symbolic, Control, Info) ;
78
long n_row, n_col, *Ap, *Ai, *Qinit, status ;
79
double Control [UMFPACK_CONTROL], Info [UMFPACK_INFO], *Ax ;
80
status = umfpack_dl_qsymbolic (n_row, n_col, Ap, Ai, Ax, Qinit,
81
&Symbolic, Control, Info) ;
87
int n_row, n_col, *Ap, *Ai, *Qinit, status ;
88
double Control [UMFPACK_CONTROL], Info [UMFPACK_INFO], *Ax, *Az ;
89
status = umfpack_zi_qsymbolic (n_row, n_col, Ap, Ai, Ax, Az, Qinit,
90
&Symbolic, Control, Info) ;
96
long n_row, n_col, *Ap, *Ai, *Qinit, status ;
97
double Control [UMFPACK_CONTROL], Info [UMFPACK_INFO], *Ax, *Az ;
98
status = umfpack_zl_qsymbolic (n_row, n_col, Ap, Ai, Ax, Az, Qinit,
99
&Symbolic, Control, Info) ;
103
Given the nonzero pattern of a sparse matrix A in column-oriented form, and
104
a sparsity preserving column pre-ordering Qinit, umfpack_*_qsymbolic
105
performs the symbolic factorization of A*Qinit (or A (:,Qinit) in MATLAB
106
notation). This is identical to umfpack_*_symbolic, except that neither
107
COLAMD nor AMD are called and the user input column order Qinit is used
108
instead. Note that in general, the Qinit passed to umfpack_*_qsymbolic
109
can differ from the final Q found in umfpack_*_numeric. The unsymmetric
110
strategy will perform a column etree postordering done in
111
umfpack_*_qsymbolic and sparsity-preserving modifications are made within
112
each frontal matrix during umfpack_*_numeric. The symmetric and 2-by-2
113
strategies will preserve Qinit, unless the matrix is structurally singular.
115
See umfpack_*_symbolic for more information.
117
*** WARNING *** A poor choice of Qinit can easily cause umfpack_*_numeric
118
to use a huge amount of memory and do a lot of work. The "default" symbolic
119
analysis method is umfpack_*_symbolic, not this routine. If you use this
120
routine, the performance of UMFPACK is your responsibility; UMFPACK will
121
not try to second-guess a poor choice of Qinit.
125
The value of Info [UMFPACK_STATUS]; see umfpack_*_symbolic.
126
Also returns UMFPACK_ERROR_invalid_permuation if Qinit is not a valid
131
All arguments are the same as umfpack_*_symbolic, except for the following:
133
Int Qinit [n_col] ; Input argument, not modified.
135
The user's fill-reducing initial column pre-ordering. This must be a
136
permutation of 0..n_col-1. If Qinit [k] = j, then column j is the kth
137
column of the matrix A (:,Qinit) to be factorized. If Qinit is an
138
(Int *) NULL pointer, then COLAMD or AMD are called instead.