5 int i, j, p, q, i2, j2, n, *Ap, *Ai, *
Cp, *Ci, *w ;
8 if (!
CS_CSC (A))
return (NULL) ;
9 n = A->
n ; Ap = A->
p ; Ai = A->
i ; Ax = A->
x ;
10 C =
cs_spalloc (n, n, Ap [n], values && (Ax != NULL), 0) ;
12 if (!C || !w)
return (
cs_done (C, w, NULL, 0)) ;
13 Cp = C->
p ; Ci = C->
i ; Cx = C->
x ;
14 for (j = 0 ; j < n ; j++)
16 j2 = pinv ? pinv [j] : j ;
17 for (p = Ap [j] ; p < Ap [j+1] ; p++)
21 i2 = pinv ? pinv [i] : i ;
26 for (j = 0 ; j < n ; j++)
28 j2 = pinv ? pinv [j] : j ;
29 for (p = Ap [j] ; p < Ap [j+1] ; p++)
33 i2 = pinv ? pinv [i] : i ;
35 if (Cx) Cx [q] = Ax [p] ;
38 return (
cs_done (C, w, NULL, 1)) ;
cs * cs_spalloc(int m, int n, int nzmax, int values, int triplet)
cs * cs_done(cs *C, void *w, void *x, int ok)
void * cs_calloc(int n, size_t size)
double cs_cumsum(int *p, int *c, int n)
cs * cs_symperm(const cs *A, const int *pinv, int values)