1 /* PR target/56348 */
2 /* { dg-do compile { target ia32 } } */
3 /* { dg-require-effective-target fpic } */
4 /* { dg-options "-O2 -fPIC -mfpmath=sse -msse2" } */
5
6 typedef unsigned int size_t;
7
8 extern double fabs (double __x) __attribute__ ((__nothrow__, __leaf__))
9 __attribute__ ((__const__));
10
11 typedef struct cholmod_sparse_struct
12 {
13 size_t ncol;
14 void *p;
15 } cholmod_sparse;
16
17 int cholmod_l_reallocate_sparse (size_t, cholmod_sparse *, void *);
18
19 int
cholmod_l_drop(double tol,cholmod_sparse * A)20 cholmod_l_drop (double tol, cholmod_sparse * A)
21 {
22 double aij;
23 double *Ax;
24 long long *Ap, *Ai, *Anz;
25 long long packed, i, j, nrow, ncol, p, pend, nz, values;
26 Ap = A->p;
27 ncol = A->ncol;
28 nz = 0;
29 for (j = 0; j < ncol; j++)
30 for (; p < pend; p++)
31 {
32 i = Ai[p];
33 aij = Ax[p];
34 if (i <= j && (fabs (aij) > tol || ((aij) != (aij))))
35 nz++;
36 }
37 Ap[ncol] = nz;
38 cholmod_l_reallocate_sparse (nz, A, 0);
39 }
40