3 int cs_fkeep (
cs *
A,
int (*fkeep) (
int,
int,
double,
void *),
void *other)
5 int j, p, nz = 0, n, *Ap, *Ai ;
7 if (!
CS_CSC (A) || !fkeep)
return (-1) ;
8 n = A->
n ; Ap = A->
p ; Ai = A->
i ; Ax = A->
x ;
9 for (j = 0 ; j < n ; j++)
13 for ( ; p < Ap [j+1] ; p++)
15 if (fkeep (Ai [p], j, Ax ? Ax [p] : 1, other))
17 if (Ax) Ax [nz] = Ax [p] ;
int cs_sprealloc(cs *A, int nzmax)
int cs_fkeep(cs *A, int(*fkeep)(int, int, double, void *), void *other)