0byt3m1n1
Path:
C:
/
Program Files
/
R
/
R-4.4.1
/
library
/
Matrix
/
scripts
/
[
Home
]
File: wall.patch
diff -ruN src/SuiteSparse/CCOLAMD/Source/ccolamd.c src/SuiteSparse-patched/CCOLAMD/Source/ccolamd.c --- src/SuiteSparse/CCOLAMD/Source/ccolamd.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CCOLAMD/Source/ccolamd.c 2024-03-16 01:25:07 @@ -1545,7 +1545,9 @@ Int *dead_cols ; Int set1 ; Int set2 ; +#ifndef NDEBUG Int cs ; +#endif int ok ; @@ -1894,7 +1896,9 @@ p [k] = col ; ASSERT (A [col] == EMPTY) ; +#ifndef NDEBUG cs = CMEMBER (col) ; +#endif ASSERT (k >= cset_start [cs] && k < cset_start [cs+1]) ; A [col] = k ; @@ -1911,8 +1915,8 @@ if (A [col] == EMPTY) { k = Col [col].shared2.order ; - cs = CMEMBER (col) ; #ifndef NDEBUG + cs = CMEMBER (col) ; dead_cols [cs]-- ; #endif ASSERT (k >= cset_start [cs] && k < cset_start [cs+1]) ; diff -ruN src/SuiteSparse/CHOLMOD/Check/cholmod_check.c src/SuiteSparse-patched/CHOLMOD/Check/cholmod_check.c --- src/SuiteSparse/CHOLMOD/Check/cholmod_check.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Check/cholmod_check.c 2024-03-16 01:25:07 @@ -216,7 +216,6 @@ //-------------------------------------------------------------------------- RETURN_IF_NULL_COMMON (FALSE) ; - int init_print = print ; P2 ("%s", "\n") ; @@ -965,7 +964,7 @@ ) { void *Xx, *Xz ; - Int i, j, d, nrow, ncol, nzmax, nz, init_print, count ; + Int i, j, d, nrow, ncol, nzmax, init_print, count ; const char *type = "dense" ; //-------------------------------------------------------------------------- @@ -1034,7 +1033,6 @@ { init_print = print ; ETC_START (count, 9) ; - nz = nrow * ncol ; for (j = 0 ; j < ncol ; j++) { ETC (j == ncol-1, count, 5) ; @@ -1487,10 +1485,10 @@ void *Lx, *Lz ; Int *Lp, *Li, *Lnz, *Lnext, *Lprev, *Perm, *ColCount, *Lpi, *Lpx, *Super, *Ls ; - Int n, nzmax, j, p, pend, i, nz, ordering, space, is_monotonic, minor, - count, precise, init_print, ilast, lnz, head, tail, jprev, plast, - jnext, examine_super, nsuper, s, k1, k2, psi, psend, psx, nsrow, nscol, - ps2, psxend, ssize, xsize, maxcsize, maxesize, nsrow2, jj, ii ; + Int n, nzmax, j, p, pend, i, nz, space, is_monotonic, minor, + count, init_print, ilast, lnz, head, tail, jprev, plast, + jnext, examine_super, nsuper, s, k1, k2, psi, psend, psx = 0, nsrow, nscol, + psxend, ssize, xsize, maxcsize, maxesize, nsrow2, jj, ii ; Int check_Lpx ; const char *type = "factor" ; @@ -1512,7 +1510,6 @@ n = L->n ; minor = L->minor ; - ordering = L->ordering ; int xtype = L->xtype ; int dtype = L->dtype ; @@ -1520,8 +1517,6 @@ ColCount = L->ColCount ; lnz = 0 ; - precise = Common->precise ; - P3 (" "ID"", n) ; P3 ("-by-"ID"", n) ; @@ -1896,7 +1891,6 @@ nsrow = psend - psi ; nscol = k2 - k1 ; nsrow2 = nsrow - nscol ; - ps2 = psi + nscol ; if (check_Lpx) { diff -ruN src/SuiteSparse/CHOLMOD/Check/cholmod_write.c src/SuiteSparse-patched/CHOLMOD/Check/cholmod_write.c --- src/SuiteSparse/CHOLMOD/Check/cholmod_write.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Check/cholmod_write.c 2024-03-16 01:25:07 @@ -354,7 +354,7 @@ double x = 0, z = 0 ; void *Ax, *Az ; Int *Ap, *Ai, *Anz, *Zp, *Zi, *Znz ; - Int nrow, ncol, is_complex, symmetry, i, j, q, iz, p, nz, is_binary, stype, + Int nrow, ncol, is_complex, i, j, q, iz, p, nz, is_binary, stype, is_integer, asym, is_sym, apacked, zpacked, pend, qend, zsym ; int ok, xtype, dtype ; @@ -527,7 +527,6 @@ // A is rectangular or unsymmetric ok = ok && (fprintf (f, " general\n") > 0) ; is_sym = FALSE ; - symmetry = CHOLMOD_MM_UNSYMMETRIC ; break ; case CHOLMOD_MM_SYMMETRIC: @@ -535,7 +534,6 @@ // A is symmetric ok = ok && (fprintf (f, " symmetric\n") > 0) ; is_sym = TRUE ; - symmetry = CHOLMOD_MM_SYMMETRIC ; break ; case CHOLMOD_MM_HERMITIAN: @@ -543,14 +541,12 @@ // A is Hermitian ok = ok && (fprintf (f, " Hermitian\n") > 0) ; is_sym = TRUE ; - symmetry = CHOLMOD_MM_HERMITIAN ; break ; case CHOLMOD_MM_SKEW_SYMMETRIC: // A is skew symmetric ok = ok && (fprintf (f, " skew-symmetric\n") > 0) ; is_sym = TRUE ; - symmetry = CHOLMOD_MM_SKEW_SYMMETRIC ; break ; } diff -ruN src/SuiteSparse/CHOLMOD/Cholesky/cholmod_analyze.c src/SuiteSparse-patched/CHOLMOD/Cholesky/cholmod_analyze.c --- src/SuiteSparse/CHOLMOD/Cholesky/cholmod_analyze.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Cholesky/cholmod_analyze.c 2024-03-16 01:25:07 @@ -399,7 +399,7 @@ Int *First, *Level, *Work4n, *Cmember, *CParent, *ColCount, *Lperm, *Parent, *Post, *Perm, *Lparent, *Lcolcount ; cholmod_factor *L ; - Int k, n, method, nmethods, status, default_strategy, ncol, + Int k, n, method, nmethods, status, default_strategy, skip_analysis, skip_best ; Int amd_backup ; int ok = TRUE ; @@ -421,7 +421,6 @@ //-------------------------------------------------------------------------- n = A->nrow ; - ncol = A->ncol ; //-------------------------------------------------------------------------- // set the default strategy diff -ruN src/SuiteSparse/CHOLMOD/Cholesky/cholmod_factorize.c src/SuiteSparse-patched/CHOLMOD/Cholesky/cholmod_factorize.c --- src/SuiteSparse/CHOLMOD/Cholesky/cholmod_factorize.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Cholesky/cholmod_factorize.c 2024-03-16 01:25:07 @@ -104,7 +104,7 @@ { cholmod_sparse *S, *F, *A1, *A2 ; - Int nrow, ncol, stype, convert, n, grow2, status ; + Int nrow, ncol, stype, convert, grow2, status ; int ok = TRUE ; //-------------------------------------------------------------------------- @@ -118,7 +118,6 @@ RETURN_IF_XTYPE_INVALID (L, CHOLMOD_PATTERN, CHOLMOD_ZOMPLEX, FALSE) ; nrow = A->nrow ; ncol = A->ncol ; - n = L->n ; stype = A->stype ; if (L->n != A->nrow) { diff -ruN src/SuiteSparse/CHOLMOD/Cholesky/cholmod_resymbol.c src/SuiteSparse-patched/CHOLMOD/Cholesky/cholmod_resymbol.c --- src/SuiteSparse/CHOLMOD/Cholesky/cholmod_resymbol.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Cholesky/cholmod_resymbol.c 2024-03-16 01:25:07 @@ -102,8 +102,6 @@ //-------------------------------------------------------------------------- Int stype = A->stype ; - Int nrow = A->nrow ; - Int ncol = A->ncol ; // s = 2*nrow + (stype ? 0 : ncol) int ok = TRUE ; diff -ruN src/SuiteSparse/CHOLMOD/Cholesky/cholmod_solve.c src/SuiteSparse-patched/CHOLMOD/Cholesky/cholmod_solve.c --- src/SuiteSparse/CHOLMOD/Cholesky/cholmod_solve.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Cholesky/cholmod_solve.c 2024-03-16 01:25:07 @@ -203,8 +203,8 @@ cholmod_dense *Y = NULL, *X = NULL ; cholmod_sparse *C, *Yset, C_header, Yset_header, *Xset ; Int *Perm = NULL, *IPerm = NULL ; - Int n, nrhs, ncols, k1, nr, ytype, k, blen, p, i, d, nrow ; - Int Cp [2], Ysetp [2], *Ci, *Yseti, ysetlen ; + Int n, nrhs, ncols, k1, nr, ytype, k, blen, p, d, nrow ; + Int Cp [2], Ysetp [2], *Ci, *Yseti ; RETURN_IF_NULL_COMMON (FALSE) ; RETURN_IF_NULL (L, FALSE) ; diff -ruN src/SuiteSparse/CHOLMOD/Cholesky/t_cholmod_psolve_worker.c src/SuiteSparse-patched/CHOLMOD/Cholesky/t_cholmod_psolve_worker.c --- src/SuiteSparse/CHOLMOD/Cholesky/t_cholmod_psolve_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Cholesky/t_cholmod_psolve_worker.c 2024-03-16 01:25:07 @@ -265,12 +265,17 @@ //-------------------------------------------------------------------------- Real *Yx, *Yz, *Xx, *Xz ; - Int k2, nk, p, k, j, nrow, ncol, d, dj, j2 ; + Int k2, p, k, j, nrow, ncol, d, dj, j2 ; +#ifndef NDEBUG + Int nk ; +#endif ncol = X->ncol ; nrow = X->nrow ; k2 = MIN (k1+ncols, ncol) ; +#ifndef NDEBUG nk = MAX (k2 - k1, 0) ; +#endif d = X->d ; Xx = X->x ; Xz = X->z ; diff -ruN src/SuiteSparse/CHOLMOD/Cholesky/t_cholmod_resymbol_worker.c src/SuiteSparse-patched/CHOLMOD/Cholesky/t_cholmod_resymbol_worker.c --- src/SuiteSparse/CHOLMOD/Cholesky/t_cholmod_resymbol_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Cholesky/t_cholmod_resymbol_worker.c 2024-03-16 01:25:07 @@ -29,12 +29,13 @@ Int *Ap = A->p ; Int *Anz = A->nz ; bool apacked = A->packed ; - bool sorted = A->sorted ; int stype = A->stype ; Int *Li = L->i ; Real *Lx = L->x ; +#ifdef ZOMPLEX Real *Lz = L->z ; +#endif Int *Lp = L->p ; Int *Lnz = L->nz ; diff -ruN src/SuiteSparse/CHOLMOD/Cholesky/t_cholmod_rowfac_worker.c src/SuiteSparse-patched/CHOLMOD/Cholesky/t_cholmod_rowfac_worker.c --- src/SuiteSparse/CHOLMOD/Cholesky/t_cholmod_rowfac_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Cholesky/t_cholmod_rowfac_worker.c 2024-03-16 01:25:07 @@ -49,11 +49,14 @@ #ifdef ZOMPLEX Real yz [1], lz [1], fz [1] ; #endif - Real *Ax, *Az, *Lx, *Lz, *Wx, *Wz, *Fx, *Fz ; + Real *Ax, *Lx, *Wx, *Fx ; +#ifdef ZOMPLEX + Real *Az, *Lz, *Wz, *Fz ; +#endif Int *Ap, *Anz, *Ai, *Lp, *Lnz, *Li, *Lnext, *Flag, *Stack, *Fp, *Fi, *Fnz, *Iwork ; Int i, p, k, t, pf, pfend, top, s, mark, pend, n, lnz, is_ll, multadds, - use_bound, packed, stype, Fpacked, sorted, nzmax, len, parent ; + use_bound, packed, stype, Fpacked, sorted, len, parent ; #ifndef REAL Int dk_imaginary ; #endif @@ -71,7 +74,9 @@ Fp = NULL ; Fi = NULL ; Fx = NULL ; +#ifdef ZOMPLEX Fz = NULL ; +#endif Fnz = NULL ; Fpacked = TRUE ; } @@ -81,7 +86,9 @@ Fp = F->p ; Fi = F->i ; Fx = F->x ; +#ifdef ZOMPLEX Fz = F->z ; +#endif Fnz = F->nz ; Fpacked = F->packed ; } @@ -89,7 +96,9 @@ Ap = A->p ; // size A->ncol+1, column pointers of A Ai = A->i ; // size nz = Ap [A->ncol], row indices of A Ax = A->x ; // size nz, numeric values of A +#ifdef ZOMPLEX Az = A->z ; +#endif Anz = A->nz ; packed = A->packed ; sorted = A->sorted ; @@ -153,8 +162,9 @@ Lnext = L->next ; // size n+2 Li = L->i ; // size L->nzmax, can change in size Lx = L->x ; // size L->nzmax or 2*L->nzmax, can change in size +#ifdef ZOMPLEX Lz = L->z ; // size L->nzmax for zomplex case, can change in size - nzmax = L->nzmax ; +#endif ASSERT (Lnz != NULL && Li != NULL && Lx != NULL) ; //-------------------------------------------------------------------------- @@ -166,9 +176,13 @@ Flag = Common->Flag ; // size n, Flag [i] < mark must hold Wx = Common->Xwork ; // size n if real, 2*n if complex or // zomplex. Xwork [i] == 0 must hold. +#ifdef ZOMPLEX Wz = Wx + n ; // size n for zomplex case only +#endif mark = Common->mark ; +#ifndef NDEBUG size_t wsize = (L->xtype == CHOLMOD_REAL ? 1:2) * ((size_t) n) ; +#endif ASSERT (Common->xworkbytes >= wsize * sizeof (Real)) ; //-------------------------------------------------------------------------- @@ -388,7 +402,9 @@ } Li = L->i ; // L->i, L->x, L->z may have moved Lx = L->x ; +#ifdef ZOMPLEX Lz = L->z ; +#endif p = Lp [i] + lnz ; // contents of L->p changed ASSERT (p < Lp [Lnext [i]]) ; } diff -ruN src/SuiteSparse/CHOLMOD/Cholesky/t_cholmod_solve_worker.c src/SuiteSparse-patched/CHOLMOD/Cholesky/t_cholmod_solve_worker.c --- src/SuiteSparse/CHOLMOD/Cholesky/t_cholmod_solve_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Cholesky/t_cholmod_solve_worker.c 2024-03-16 01:25:07 @@ -70,7 +70,9 @@ Int *Lp = L->p ; Real *Lx = L->x ; Real *Yx = Y->x ; +#ifdef ZOMPLEX Real *Yz = Y->z ; +#endif if (Yset) { @@ -209,7 +211,9 @@ //-------------------------------------------------------------------------- Real *Yx = Y->x ; +#ifdef ZOMPLEX Real *Yz = Y->z ; +#endif Int *Ysetp = Yset->p ; Int *Yseti = Yset->i ; Int ysetlen = Ysetp [1] ; @@ -227,7 +231,9 @@ // Y (C) = B (Bset) Real *Bx = B->x ; +#ifdef ZOMPLEX Real *Bz = B->z ; +#endif Int *Bsetp = Bset->p ; Int *Bseti = Bset->i ; Int *Bsetnz = Bset->nz ; @@ -262,12 +268,16 @@ { Real *Xx = X->x ; +#ifdef ZOMPLEX Real *Xz = X->z ; +#endif Int *Xseti = Xset->i ; Int *Xsetp = Xset->p ; Real *Yx = Y->x ; +#ifdef ZOMPLEX Real *Yz = Y->z ; +#endif Int *Ysetp = Yset->p ; Int *Yseti = Yset->i ; diff -ruN src/SuiteSparse/CHOLMOD/Cholesky/t_cholmod_spsolve_worker.c src/SuiteSparse-patched/CHOLMOD/Cholesky/t_cholmod_spsolve_worker.c --- src/SuiteSparse/CHOLMOD/Cholesky/t_cholmod_spsolve_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Cholesky/t_cholmod_spsolve_worker.c 2024-03-16 01:25:07 @@ -30,12 +30,16 @@ Int *Bp = B->p ; Int *Bi = B->i ; Real *Bx = B->x ; +#ifdef ZOMPLEX Real *Bz = B->z ; +#endif Int *Bnz = B->nz ; bool packed = B->packed ; Real *B4x = B4->x ; +#ifdef ZOMPLEX Real *B4z = B4->z ; +#endif Int n = B4->nrow ; @@ -78,13 +82,17 @@ Int *Xp = X->p ; Int *Xi = X->i ; Real *Xx = X->x ; +#ifdef ZOMPLEX Real *Xz = X->z ; +#endif size_t px = (*xnz) ; size_t nzmax = X->nzmax ; Real *X4x = X4->x ; +#ifdef ZOMPLEX Real *X4z = X4->z ; +#endif Int n = X4->nrow ; //-------------------------------------------------------------------------- @@ -139,7 +147,9 @@ } Xi = X->i ; Xx = X->x ; +#ifdef ZOMPLEX Xz = X->z ; +#endif } Xi [px] = i ; ASSIGN (Xx, Xz, px, X4x, X4z, p) ; @@ -180,7 +190,9 @@ bool packed = B->packed ; Real *B4x = B4->x ; +#ifdef ZOMPLEX Real *B4z = B4->z ; +#endif Int n = B4->nrow ; diff -ruN src/SuiteSparse/CHOLMOD/MatrixOps/cholmod_ssmult.c src/SuiteSparse-patched/CHOLMOD/MatrixOps/cholmod_ssmult.c --- src/SuiteSparse/CHOLMOD/MatrixOps/cholmod_ssmult.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/MatrixOps/cholmod_ssmult.c 2024-03-16 01:25:07 @@ -84,7 +84,7 @@ // check inputs //-------------------------------------------------------------------------- - cholmod_sparse *C = NULL, *A2 = NULL, *B2 = NULL, *C2 = NULL ; + cholmod_sparse *C = NULL, *A2 = NULL, *B2 = NULL ; RETURN_IF_NULL_COMMON (NULL) ; RETURN_IF_NULL (A, NULL) ; @@ -196,7 +196,6 @@ Int ncol = B->ncol ; // get workspace - void *W = Common->Xwork ; // size nrow, unused if values is false Int *Flag = Common->Flag ; // size nrow, Flag [0..nrow-1] < mark on input //-------------------------------------------------------------------------- diff -ruN src/SuiteSparse/CHOLMOD/MatrixOps/t_cholmod_drop_worker.c src/SuiteSparse-patched/CHOLMOD/MatrixOps/t_cholmod_drop_worker.c --- src/SuiteSparse/CHOLMOD/MatrixOps/t_cholmod_drop_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/MatrixOps/t_cholmod_drop_worker.c 2024-03-16 01:25:07 @@ -59,7 +59,9 @@ Int *Ap = A->p ; Int *Ai = A->i ; Real *Ax = A->x ; +#ifdef ZOMPLEX Real *Az = A->z ; +#endif Int *Anz = A->nz ; bool packed = A->packed ; Int ncol = A->ncol ; diff -ruN src/SuiteSparse/CHOLMOD/MatrixOps/t_cholmod_horzcat_worker.c src/SuiteSparse-patched/CHOLMOD/MatrixOps/t_cholmod_horzcat_worker.c --- src/SuiteSparse/CHOLMOD/MatrixOps/t_cholmod_horzcat_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/MatrixOps/t_cholmod_horzcat_worker.c 2024-03-16 01:25:07 @@ -25,23 +25,35 @@ Int *Ap = A->p ; Int *Anz = A->nz ; Int *Ai = A->i ; +#ifndef PATTERN Real *Ax = A->x ; +#ifdef ZOMPLEX Real *Az = A->z ; +#endif +#endif bool apacked = A->packed ; Int ancol = A->ncol ; Int *Bp = B->p ; Int *Bnz = B->nz ; Int *Bi = B->i ; +#ifndef PATTERN Real *Bx = B->x ; +#ifdef ZOMPLEX Real *Bz = B->z ; +#endif +#endif bool bpacked = B->packed ; Int bncol = B->ncol ; Int *Cp = C->p ; Int *Ci = C->i ; +#ifndef PATTERN Real *Cx = C->x ; +#ifdef ZOMPLEX Real *Cz = C->z ; +#endif +#endif Int ncol = C->ncol ; //-------------------------------------------------------------------------- diff -ruN src/SuiteSparse/CHOLMOD/MatrixOps/t_cholmod_norm_worker.c src/SuiteSparse-patched/CHOLMOD/MatrixOps/t_cholmod_norm_worker.c --- src/SuiteSparse/CHOLMOD/MatrixOps/t_cholmod_norm_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/MatrixOps/t_cholmod_norm_worker.c 2024-03-16 01:25:07 @@ -36,7 +36,9 @@ Int nrow = X->nrow ; Int d = X->d ; Real *Xx = X->x ; +#ifdef ZOMPLEX Real *Xz = X->z ; +#endif double xnorm = 0 ; @@ -158,8 +160,12 @@ Int *Ap = A->p ; Int *Ai = A->i ; Int *Anz = A->nz ; +#ifndef PATTERN Real *Ax = A->x ; +#ifdef ZOMPLEX Real *Az = A->z ; +#endif +#endif Int ncol = A->ncol ; Int nrow = A->nrow ; bool packed = A->packed ; diff -ruN src/SuiteSparse/CHOLMOD/MatrixOps/t_cholmod_scale_worker.c src/SuiteSparse-patched/CHOLMOD/MatrixOps/t_cholmod_scale_worker.c --- src/SuiteSparse/CHOLMOD/MatrixOps/t_cholmod_scale_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/MatrixOps/t_cholmod_scale_worker.c 2024-03-16 01:25:07 @@ -28,12 +28,16 @@ Int *Anz = A->nz ; Int *Ai = A->i ; Real *Ax = A->x ; +#ifdef ZOMPLEX Real *Az = A->z ; +#endif bool packed = A->packed ; Int ncol = A->ncol ; Real *Sx = S->x ; +#ifdef ZOMPLEX Real *Sz = S->z ; +#endif //-------------------------------------------------------------------------- // scale the matrix @@ -55,7 +59,9 @@ Int i = Ai [p] ; // t = S (i) * A (i,j) Real tx [2] ; +#ifdef ZOMPLEX Real tz [1] ; +#endif MULT (tx, tz, 0, Sx, Sz, i, Ax, Az, p) ; // A (i,j) = t ASSIGN (Ax, Az, p, tx, tz, 0) ; @@ -74,7 +80,9 @@ { // s = S (j) Real sx [2] ; +#ifdef ZOMPLEX Real sz [1] ; +#endif ASSIGN (sx, sz, 0, Sx, Sz, j) ; Int p = Ap [j] ; @@ -83,7 +91,9 @@ { // t = A (i,j) * s Real tx [2] ; +#ifdef ZOMPLEX Real tz [1] ; +#endif MULT (tx, tz, 0, Ax, Az, p, sx, sz, 0) ; // A (i,j) = t ASSIGN (Ax, Az, p, tx, tz, 0) ; @@ -102,7 +112,9 @@ { // s = S (j) Real sx [2] ; +#ifdef ZOMPLEX Real sz [1] ; +#endif ASSIGN (sx, sz, 0, Sx, Sz, j) ; Int p = Ap [j] ; @@ -112,7 +124,9 @@ Int i = Ai [p] ; // t = A (i,j) * S (i) Real tx [2] ; +#ifdef ZOMPLEX Real tz [1] ; +#endif MULT (tx, tz, 0, Ax, Az, p, Sx, Sz, i) ; // A (i,j) = s * t MULT (Ax, Az, p, sx, sz, 0, tx, tz, 0) ; @@ -129,7 +143,9 @@ // s = S (0) Real sx [2] ; +#ifdef ZOMPLEX Real sz [1] ; +#endif ASSIGN (sx, sz, 0, Sx, Sz, 0) ; for (Int j = 0 ; j < ncol ; j++) @@ -140,7 +156,9 @@ { // t = s * A (i,j) Real tx [2] ; +#ifdef ZOMPLEX Real tz [1] ; +#endif MULT (tx, tz, 0, sx, sz, 0, Ax, Az, p) ; // A (i,j) = t ASSIGN (Ax, Az, p, tx, tz, 0) ; diff -ruN src/SuiteSparse/CHOLMOD/MatrixOps/t_cholmod_sdmult_worker.c src/SuiteSparse-patched/CHOLMOD/MatrixOps/t_cholmod_sdmult_worker.c --- src/SuiteSparse/CHOLMOD/MatrixOps/t_cholmod_sdmult_worker.c 2024-03-16 01:30:04 +++ src/SuiteSparse-patched/CHOLMOD/MatrixOps/t_cholmod_sdmult_worker.c 2024-03-16 01:25:07 @@ -55,28 +55,37 @@ #endif size_t ny = transpose ? A->ncol : A->nrow ; // required length of Y +#ifdef ZOMPLEX size_t nx = transpose ? A->nrow : A->ncol ; // required length of X +#endif - Int nrow = A->nrow ; Int ncol = A->ncol ; Int *Ap = A->p ; Int *Anz = A->nz ; Int *Ai = A->i ; Real *Ax = A->x ; +#ifdef ZOMPLEX Real *Az = A->z ; +#endif bool packed = A->packed ; Real *Xx = X->x ; +#ifdef ZOMPLEX Real *Xz = X->z ; +#endif Real *Yx = Y->x ; +#ifdef ZOMPLEX Real *Yz = Y->z ; +#endif Int kcol = X->ncol ; size_t dy = Y->d ; size_t dx = X->d ; Real *w = W ; +#ifdef ZOMPLEX Real *Wz = (W) ? W + 4*nx : W ; +#endif //-------------------------------------------------------------------------- // Y = beta * Y @@ -120,7 +129,9 @@ //-------------------------------------------------------------------------- Yx = Y->x ; +#ifdef ZOMPLEX Yz = Y->z ; +#endif Int k = 0 ; if (A->stype == 0) diff -ruN src/SuiteSparse/CHOLMOD/MatrixOps/t_cholmod_ssmult_worker.c src/SuiteSparse-patched/CHOLMOD/MatrixOps/t_cholmod_ssmult_worker.c --- src/SuiteSparse/CHOLMOD/MatrixOps/t_cholmod_ssmult_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/MatrixOps/t_cholmod_ssmult_worker.c 2024-03-16 01:25:07 @@ -23,35 +23,51 @@ // get inputs //-------------------------------------------------------------------------- - Real *W = Common->Xwork ; - Int *Ap = A->p ; Int *Anz = A->nz ; Int *Ai = A->i ; +#ifndef PATTERN Real *Ax = A->x ; +#ifdef ZOMPLEX Real *Az = A->z ; +#endif +#endif bool apacked = A->packed ; Int *Bp = B->p ; Int *Bnz = B->nz ; Int *Bi = B->i ; +#ifndef PATTERN Real *Bx = B->x ; +#ifdef ZOMPLEX Real *Bz = B->z ; +#endif +#endif bool bpacked = B->packed ; // get the size of C +#ifdef ZOMPLEX Int nrow = A->nrow ; +#endif Int ncol = B->ncol ; // get workspace +#ifndef PATTERN Real *Wx = Common->Xwork ; // size nrow, unused if C is pattern +#ifdef ZOMPLEX Real *Wz = Wx + nrow ; // only used for the zomplex case +#endif +#endif Int *Flag = Common->Flag ; // size nrow, Flag [0..nrow-1] < mark on input Int *Cp = C->p ; Int *Ci = C->i ; +#ifndef PATTERN Real *Cx = C->x ; +#ifdef ZOMPLEX Real *Cz = C->z ; +#endif +#endif //-------------------------------------------------------------------------- // C = A*B @@ -77,8 +93,12 @@ Int k = Bi [pb] ; // b = Bx [pb] ; +#ifndef PATTERN Real bx [2] ; +#ifdef ZOMPLEX Real bz [1] ; +#endif +#endif ASSIGN (bx, bz, 0, Bx, Bz, pb) ; // add the nonzero pattern of A(:,k) to the pattern of C(:,j) diff -ruN src/SuiteSparse/CHOLMOD/MatrixOps/t_cholmod_submatrix_worker.c src/SuiteSparse-patched/CHOLMOD/MatrixOps/t_cholmod_submatrix_worker.c --- src/SuiteSparse/CHOLMOD/MatrixOps/t_cholmod_submatrix_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/MatrixOps/t_cholmod_submatrix_worker.c 2024-03-16 01:25:07 @@ -29,14 +29,22 @@ Int *Ap = A->p ; Int *Ai = A->i ; Int *Anz = A->nz ; +#ifndef PATTERN Real *Ax = A->x ; +#ifdef ZOMPLEX Real *Az = A->z ; +#endif +#endif bool packed = A->packed ; Int *Cp = C->p ; Int *Ci = C->i ; +#ifndef PATTERN Real *Cx = C->x ; +#ifdef ZOMPLEX Real *Cz = C->z ; +#endif +#endif Int cncol = C->ncol ; //-------------------------------------------------------------------------- diff -ruN src/SuiteSparse/CHOLMOD/MatrixOps/t_cholmod_vertcat_worker.c src/SuiteSparse-patched/CHOLMOD/MatrixOps/t_cholmod_vertcat_worker.c --- src/SuiteSparse/CHOLMOD/MatrixOps/t_cholmod_vertcat_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/MatrixOps/t_cholmod_vertcat_worker.c 2024-03-16 01:25:07 @@ -25,22 +25,34 @@ Int *Ap = A->p ; Int *Anz = A->nz ; Int *Ai = A->i ; +#ifndef PATTERN Real *Ax = A->x ; +#ifdef ZOMPLEX Real *Az = A->z ; +#endif +#endif bool apacked = A->packed ; Int anrow = A->nrow ; Int *Bp = B->p ; Int *Bnz = B->nz ; Int *Bi = B->i ; +#ifndef PATTERN Real *Bx = B->x ; +#ifdef ZOMPLEX Real *Bz = B->z ; +#endif +#endif bool bpacked = B->packed ; Int *Cp = C->p ; Int *Ci = C->i ; +#ifndef PATTERN Real *Cx = C->x ; +#ifdef ZOMPLEX Real *Cz = C->z ; +#endif +#endif Int ncol = C->ncol ; //-------------------------------------------------------------------------- diff -ruN src/SuiteSparse/CHOLMOD/Modify/cholmod_rowadd.c src/SuiteSparse-patched/CHOLMOD/Modify/cholmod_rowadd.c --- src/SuiteSparse/CHOLMOD/Modify/cholmod_rowadd.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Modify/cholmod_rowadd.c 2024-03-16 01:25:07 @@ -136,7 +136,6 @@ RETURN_IF_NULL (R, FALSE) ; RETURN_IF_XTYPE_INVALID (L, CHOLMOD_PATTERN, CHOLMOD_REAL, FALSE) ; RETURN_IF_XTYPE_INVALID (R, CHOLMOD_REAL, CHOLMOD_REAL, FALSE) ; - Int n = L->n ; Int k = kadd ; if (kadd >= L->n || k < 0) { diff -ruN src/SuiteSparse/CHOLMOD/Modify/cholmod_rowdel.c src/SuiteSparse-patched/CHOLMOD/Modify/cholmod_rowdel.c --- src/SuiteSparse/CHOLMOD/Modify/cholmod_rowdel.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Modify/cholmod_rowdel.c 2024-03-16 01:25:07 @@ -122,7 +122,6 @@ RETURN_IF_NULL_COMMON (FALSE) ; RETURN_IF_NULL (L, FALSE) ; RETURN_IF_XTYPE_INVALID (L, CHOLMOD_PATTERN, CHOLMOD_REAL, FALSE) ; - Int n = L->n ; Int k = kdel ; if (kdel >= L->n || k < 0) { diff -ruN src/SuiteSparse/CHOLMOD/Modify/t_cholmod_updown_worker.c src/SuiteSparse-patched/CHOLMOD/Modify/t_cholmod_updown_worker.c --- src/SuiteSparse/CHOLMOD/Modify/t_cholmod_updown_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Modify/t_cholmod_updown_worker.c 2024-03-16 01:25:07 @@ -148,7 +148,9 @@ double fl = 0 ; bool use_colmark = (colmark != NULL) ; +#ifndef NDEBUG Int *ps1 = NULL ; +#endif Int *ps2 = NULL ; //-------------------------------------------------------------------------- @@ -699,7 +701,9 @@ if (nsets > 0) { +#ifndef NDEBUG ps1 = Set_ps1 [0] ; +#endif ps2 = Set_ps2 [0] ; } diff -ruN src/SuiteSparse/CHOLMOD/Partition/cholmod_metis.c src/SuiteSparse-patched/CHOLMOD/Partition/cholmod_metis.c --- src/SuiteSparse/CHOLMOD/Partition/cholmod_metis.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Partition/cholmod_metis.c 2024-03-16 01:25:07 @@ -578,7 +578,7 @@ idx_t *Mp, *Mi, *Mperm, *Miperm ; cholmod_sparse *B ; Int i, j, n, nz, p, identity ; - idx_t nn, zero = 0 ; + idx_t nn ; RETURN_IF_NULL_COMMON (FALSE) ; RETURN_IF_NULL (A, FALSE) ; diff -ruN src/SuiteSparse/CHOLMOD/Partition/cholmod_nesdis.c src/SuiteSparse-patched/CHOLMOD/Partition/cholmod_nesdis.c --- src/SuiteSparse/CHOLMOD/Partition/cholmod_nesdis.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Partition/cholmod_nesdis.c 2024-03-16 01:25:07 @@ -681,8 +681,11 @@ ) { - Int n, mark, cj, j, sj, sn, p, i, snode, pstart, pdest, pend, nd_components, + Int mark, cj, j, sj, sn, p, i, snode, pstart, pdest, pend, nd_components, part, first, save_mark ; +#ifndef NDEBUG + Int n ; +#endif Int *Bp, *Bi, *Flag ; //-------------------------------------------------------------------------- @@ -702,7 +705,9 @@ Bp = B->p ; Bi = B->i ; +#ifndef NDEBUG n = B->nrow ; +#endif ASSERT (cnode >= EMPTY && cnode < n) ; ASSERT (IMPLIES (cnode >= 0, Flag [cnode] < EMPTY)) ; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/GKlib/gk_mkpqueue.h src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/GKlib/gk_mkpqueue.h --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/GKlib/gk_mkpqueue.h 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/GKlib/gk_mkpqueue.h 2024-03-16 01:25:07 @@ -11,6 +11,38 @@ #ifndef _GK_MKPQUEUE_H #define _GK_MKPQUEUE_H +#ifndef NDEBUG +#define DO_CHECKHEAP\ +do {\ + gk_idx_t i, j;\ + size_t nnodes;\ + gk_idx_t *locator;\ + KVT *heap;\ +\ + heap = queue->heap;\ + locator = queue->locator;\ + nnodes = queue->nnodes;\ +\ + if (nnodes == 0)\ + return 1;\ +\ + ASSERT(locator[heap[0].val] == 0);\ + for (i=1; i<nnodes; i++) {\ + ASSERT(locator[heap[i].val] == i);\ + ASSERT(!KEY_LT(heap[i].key, heap[(i-1)/2].key));\ + }\ + for (i=1; i<nnodes; i++)\ + ASSERT(!KEY_LT(heap[i].key, heap[0].key));\ +\ + for (j=i=0; i<queue->maxnodes; i++) {\ + if (locator[i] != -1)\ + j++;\ + }\ + ASSERTP(j == nnodes, ("%jd %jd\n", (intmax_t)j, (intmax_t)nnodes));\ +} while (0) +#else +#define DO_CHECKHEAP +#endif #define GK_MKPQUEUE(FPRFX, PQT, KVT, KT, VT, KVMALLOC, KMAX, KEY_LT)\ /*************************************************************************/\ @@ -380,32 +412,7 @@ /**************************************************************************/\ int FPRFX ## CheckHeap(PQT *queue)\ {\ - gk_idx_t i, j;\ - size_t nnodes;\ - gk_idx_t *locator;\ - KVT *heap;\ -\ - heap = queue->heap;\ - locator = queue->locator;\ - nnodes = queue->nnodes;\ -\ - if (nnodes == 0)\ - return 1;\ -\ - ASSERT(locator[heap[0].val] == 0);\ - for (i=1; i<nnodes; i++) {\ - ASSERT(locator[heap[i].val] == i);\ - ASSERT(!KEY_LT(heap[i].key, heap[(i-1)/2].key));\ - }\ - for (i=1; i<nnodes; i++)\ - ASSERT(!KEY_LT(heap[i].key, heap[0].key));\ -\ - for (j=i=0; i<queue->maxnodes; i++) {\ - if (locator[i] != -1)\ - j++;\ - }\ - ASSERTP(j == nnodes, ("%jd %jd\n", (intmax_t)j, (intmax_t)nnodes));\ -\ + DO_CHECKHEAP;\ return 1;\ }\ diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/GKlib/mcore.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/GKlib/mcore.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/GKlib/mcore.c 2024-03-16 01:30:04 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/GKlib/mcore.c 2024-03-16 01:28:20 @@ -71,6 +71,7 @@ if (mcore == NULL) return; +#ifndef NDEBUG if (showstats) printf("\n gk_mcore statistics\n" " coresize: %12zu nmops: %12zu cmop: %6zu\n" @@ -89,6 +90,7 @@ " cur_callocs: %6zu cur_hallocs: %6zu cmop: %6zu\n", mcore->cur_callocs, mcore->cur_hallocs, mcore->cmop); } +#endif gk_free((void **)&mcore->core, &mcore->mops, &mcore, LTERM); @@ -107,6 +109,7 @@ if (mcore == NULL) return; +#ifndef NDEBUG if (showstats) printf("\n gk_mcore statistics\n" " nmops: %12zu cmop: %6zu\n" @@ -125,6 +128,7 @@ " cur_hallocs: %6zu cmop: %6zu\n", mcore->cur_hallocs, mcore->cmop); } +#endif free(mcore->mops); free(mcore); diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/GKlib/memory.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/GKlib/memory.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/GKlib/memory.c 2024-03-16 01:30:04 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/GKlib/memory.c 2024-03-16 01:28:20 @@ -152,8 +152,10 @@ ptr = (void *)malloc(nbytes); if (ptr == NULL) { +#ifndef NDEBUG fprintf(stderr, " Current memory used: %10zu bytes\n", gk_GetCurMemoryUsed()); fprintf(stderr, " Maximum memory used: %10zu bytes\n", gk_GetMaxMemoryUsed()); +#endif gk_errexit(SIGMEM, "***Memory allocation failed for %s. Requested size: %zu bytes", msg, nbytes); return NULL; @@ -187,8 +189,10 @@ ptr = (void *)realloc(oldptr, nbytes); if (ptr == NULL) { +#ifndef NDEBUG fprintf(stderr, " Maximum memory used: %10zu bytes\n", gk_GetMaxMemoryUsed()); fprintf(stderr, " Current memory used: %10zu bytes\n", gk_GetCurMemoryUsed()); +#endif gk_errexit(SIGMEM, "***Memory realloc failed for %s. " "Requested size: %zu bytes", msg, nbytes); return NULL; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/balance.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/balance.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/balance.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/balance.c 2024-03-16 01:25:07 @@ -40,7 +40,7 @@ **************************************************************************/ void Bnd2WayBalance(ctrl_t *ctrl, graph_t *graph, real_t *ntpwgts) { - idx_t i, ii, j, k, kwgt, nvtxs, nbnd, nswaps, from, to, pass, me, tmp; + idx_t i, ii, j, k, kwgt, nvtxs, nbnd, nswaps, from, to, tmp; idx_t *xadj, *vwgt, *adjncy, *adjwgt, *where, *id, *ed, *bndptr, *bndind, *pwgts; idx_t *moved, *perm; rpq_t *queue; @@ -168,7 +168,7 @@ **************************************************************************/ void General2WayBalance(ctrl_t *ctrl, graph_t *graph, real_t *ntpwgts) { - idx_t i, ii, j, k, kwgt, nvtxs, nbnd, nswaps, from, to, pass, me, tmp; + idx_t i, ii, j, k, kwgt, nvtxs, nbnd, nswaps, from, to, tmp; idx_t *xadj, *vwgt, *adjncy, *adjwgt, *where, *id, *ed, *bndptr, *bndind, *pwgts; idx_t *moved, *perm; rpq_t *queue; @@ -280,8 +280,8 @@ **************************************************************************/ void McGeneral2WayBalance(ctrl_t *ctrl, graph_t *graph, real_t *ntpwgts) { - idx_t i, ii, j, k, l, kwgt, nvtxs, ncon, nbnd, nswaps, from, to, pass, - me, limit, tmp, cnum; + idx_t i, ii, j, k, l, kwgt, nvtxs, ncon, nbnd, nswaps, from, to, + limit, tmp, cnum; idx_t *xadj, *adjncy, *vwgt, *adjwgt, *where, *pwgts, *id, *ed, *bndptr, *bndind; idx_t *moved, *swaps, *perm, *qnum, *qsizes; idx_t higain, mincut, newcut, mincutorder; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/coarsen.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/coarsen.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/coarsen.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/coarsen.c 2024-03-16 01:25:07 @@ -21,7 +21,7 @@ /*************************************************************************/ graph_t *CoarsenGraph(ctrl_t *ctrl, graph_t *graph) { - idx_t i, eqewgts, level=0; + idx_t i, eqewgts; IFSET(ctrl->dbglvl, METIS_DBG_TIME, gk_startcputimer(ctrl->CoarsenTmr)); @@ -62,7 +62,6 @@ graph = graph->coarser; eqewgts = 0; - level++; ASSERT(CheckGraph(graph, 0, 1)); @@ -148,8 +147,8 @@ /**************************************************************************/ idx_t Match_RM(ctrl_t *ctrl, graph_t *graph) { - idx_t i, pi, ii, j, jj, jjinc, k, nvtxs, ncon, cnvtxs, maxidx, last_unmatched; - idx_t *xadj, *vwgt, *adjncy, *adjwgt, *maxvwgt; + idx_t i, pi, j, k, nvtxs, ncon, cnvtxs, maxidx, last_unmatched; + idx_t *xadj, *vwgt, *adjncy, *maxvwgt; idx_t *match, *cmap, *perm; size_t nunmatched=0; @@ -162,7 +161,6 @@ xadj = graph->xadj; vwgt = graph->vwgt; adjncy = graph->adjncy; - adjwgt = graph->adjwgt; cmap = graph->cmap; maxvwgt = ctrl->maxvwgt; @@ -275,7 +273,7 @@ /**************************************************************************/ idx_t Match_SHEM(ctrl_t *ctrl, graph_t *graph) { - idx_t i, pi, ii, j, jj, jjinc, k, nvtxs, ncon, cnvtxs, maxidx, maxwgt, + idx_t i, pi, j, k, nvtxs, ncon, cnvtxs, maxidx, maxwgt, last_unmatched, avgdegree; idx_t *xadj, *vwgt, *adjncy, *adjwgt, *maxvwgt; idx_t *match, *cmap, *degrees, *perm, *tperm; @@ -437,7 +435,7 @@ idx_t Match_2HopAny(ctrl_t *ctrl, graph_t *graph, idx_t *perm, idx_t *match, idx_t cnvtxs, size_t *r_nunmatched, size_t maxdegree) { - idx_t i, pi, ii, j, jj, k, nvtxs; + idx_t i, pi, j, jj, nvtxs; idx_t *xadj, *adjncy, *colptr, *rowind; idx_t *cmap; size_t nunmatched; @@ -516,7 +514,7 @@ idx_t Match_2HopAll(ctrl_t *ctrl, graph_t *graph, idx_t *perm, idx_t *match, idx_t cnvtxs, size_t *r_nunmatched, size_t maxdegree) { - idx_t i, pi, pk, ii, j, jj, k, nvtxs, mask, idegree; + idx_t i, pi, pk, j, jj, k, nvtxs, mask, idegree; idx_t *xadj, *adjncy; idx_t *cmap, *mark; ikv_t *keys; @@ -621,7 +619,7 @@ void CreateCoarseGraph(ctrl_t *ctrl, graph_t *graph, idx_t cnvtxs, idx_t *match) { - idx_t j, jj, k, kk, l, m, istart, iend, nvtxs, nedges, ncon, cnedges, + idx_t j, jj, k, kk, m, istart, iend, nvtxs, nedges, ncon, cnedges, v, u, mask, dovsize; idx_t *xadj, *vwgt, *vsize, *adjncy, *adjwgt; idx_t *cmap, *htable; @@ -932,7 +930,7 @@ void CreateCoarseGraphPerm(ctrl_t *ctrl, graph_t *graph, idx_t cnvtxs, idx_t *match, idx_t *perm) { - idx_t i, j, jj, k, kk, l, m, istart, iend, nvtxs, nedges, ncon, cnedges, + idx_t i, j, jj, k, kk, m, istart, iend, nvtxs, nedges, ncon, cnedges, v, u, mask, dovsize; idx_t *xadj, *vwgt, *vsize, *adjncy, *adjwgt; idx_t *cmap, *htable; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/compress.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/compress.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/compress.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/compress.c 2024-03-16 01:25:07 @@ -151,7 +151,7 @@ idx_t *vwgt, idx_t *iperm, real_t factor) { idx_t i, j, k, l, nlarge, pnvtxs, pnedges; - idx_t *pxadj, *padjncy, *padjwgt, *pvwgt; + idx_t *pxadj, *padjncy, *pvwgt; idx_t *perm; graph_t *graph=NULL; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/contig.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/contig.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/contig.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/contig.c 2024-03-16 01:25:07 @@ -32,7 +32,7 @@ idx_t FindPartitionInducedComponents(graph_t *graph, idx_t *where, idx_t *cptr, idx_t *cind) { - idx_t i, ii, j, jj, k, me=0, nvtxs, first, last, nleft, ncmps; + idx_t i, j, k, me=0, nvtxs, first, last, nleft, ncmps; idx_t *xadj, *adjncy; idx_t *touched, *perm, *todo; idx_t mustfree_ccsr=0, mustfree_where=0; @@ -267,7 +267,7 @@ idx_t FindSepInducedComponents(ctrl_t *ctrl, graph_t *graph, idx_t *cptr, idx_t *cind) { - idx_t i, j, k, nvtxs, first, last, nleft, ncmps, wgt; + idx_t i, j, k, nvtxs, first, last, nleft, ncmps; idx_t *xadj, *adjncy, *where, *touched, *queue; nvtxs = graph->nvtxs; @@ -335,14 +335,13 @@ /*************************************************************************/ void EliminateComponents(ctrl_t *ctrl, graph_t *graph) { - idx_t i, ii, j, jj, k, me, nparts, nvtxs, ncon, ncmps, other, + idx_t i, ii, j, jj, me, nparts, nvtxs, ncon, ncmps, ncand, target; idx_t *xadj, *adjncy, *vwgt, *adjwgt, *where, *pwgts; idx_t *cptr, *cind, *cpvec, *pcptr, *pcind, *cwhere; idx_t cid, bestcid, *cwgt, *bestcwgt; idx_t ntodo, oldntodo, *todo; rkv_t *cand; - real_t *tpwgts; idx_t *vmarker=NULL, *pmarker=NULL, *modind=NULL; /* volume specific work arrays */ WCOREPUSH; @@ -358,7 +357,6 @@ pwgts = graph->pwgts; nparts = ctrl->nparts; - tpwgts = ctrl->tpwgts; cptr = iwspacemalloc(ctrl, nvtxs+1); cind = iwspacemalloc(ctrl, nvtxs); @@ -531,12 +529,11 @@ void MoveGroupContigForCut(ctrl_t *ctrl, graph_t *graph, idx_t to, idx_t gid, idx_t *ptr, idx_t *ind) { - idx_t i, ii, iii, j, jj, k, l, nvtxs, nbnd, from, me; + idx_t i, ii, iii, j, k, nbnd, from, me; idx_t *xadj, *adjncy, *adjwgt, *where, *bndptr, *bndind; ckrinfo_t *myrinfo; cnbr_t *mynbrs; - nvtxs = graph->nvtxs; xadj = graph->xadj; adjncy = graph->adjncy; adjwgt = graph->adjwgt; @@ -602,12 +599,11 @@ idx_t *ptr, idx_t *ind, idx_t *vmarker, idx_t *pmarker, idx_t *modind) { - idx_t i, ii, iii, j, jj, k, l, nvtxs, from, me, other, xgain; + idx_t i, ii, iii, j, k, l, from, other, xgain; idx_t *xadj, *vsize, *adjncy, *where; vkrinfo_t *myrinfo, *orinfo; vnbr_t *mynbrs, *onbrs; - nvtxs = graph->nvtxs; xadj = graph->xadj; vsize = graph->vsize; adjncy = graph->adjncy; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/debug.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/debug.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/debug.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/debug.c 2024-03-16 01:28:20 @@ -47,7 +47,7 @@ /*************************************************************************/ idx_t ComputeVolume(graph_t *graph, idx_t *where) { - idx_t i, j, k, me, nvtxs, nparts, totalv; + idx_t i, j, k, nvtxs, nparts, totalv; idx_t *xadj, *adjncy, *vsize, *marker; @@ -106,7 +106,9 @@ maxcut = cuts[iargmax(nparts, cuts)]; +#ifndef NDEBUG printf("%zu => %"PRIDX"\n", iargmax(nparts, cuts), maxcut); +#endif gk_free((void **)&cuts, LTERM); @@ -120,6 +122,7 @@ /*************************************************************************/ idx_t CheckBnd(graph_t *graph) { +#ifndef NDEBUG idx_t i, j, nvtxs, nbnd; idx_t *xadj, *adjncy, *where, *bndptr, *bndind; @@ -145,6 +148,7 @@ } ASSERTP(nbnd == graph->nbnd, ("%"PRIDX" %"PRIDX"\n", nbnd, graph->nbnd)); +#endif return 1; } @@ -157,6 +161,7 @@ /*************************************************************************/ idx_t CheckBnd2(graph_t *graph) { +#ifndef NDEBUG idx_t i, j, nvtxs, nbnd, id, ed; idx_t *xadj, *adjncy, *where, *bndptr, *bndind; @@ -183,6 +188,7 @@ } ASSERTP(nbnd == graph->nbnd, ("%"PRIDX" %"PRIDX"\n", nbnd, graph->nbnd)); +#endif return 1; } @@ -194,6 +200,7 @@ /*************************************************************************/ idx_t CheckNodeBnd(graph_t *graph, idx_t onbnd) { +#ifndef NDEBUG idx_t i, j, nvtxs, nbnd; idx_t *xadj, *adjncy, *where, *bndptr, *bndind; @@ -219,6 +226,7 @@ ASSERTP(bndptr[i] != -1, ("%"PRIDX" %"PRIDX"\n", i, bndptr[i])); } } +#endif return 1; } @@ -231,6 +239,7 @@ /*************************************************************************/ idx_t CheckRInfo(ctrl_t *ctrl, ckrinfo_t *rinfo) { +#ifndef NDEBUG idx_t i, j; cnbr_t *nbrs; @@ -242,6 +251,7 @@ ("%"PRIDX" %"PRIDX" %"PRIDX" %"PRIDX"\n", i, j, nbrs[i].pid, nbrs[j].pid)); } +#endif return 1; } @@ -254,15 +264,14 @@ /*************************************************************************/ idx_t CheckNodePartitionParams(graph_t *graph) { - idx_t i, j, k, l, nvtxs, me, other; - idx_t *xadj, *adjncy, *adjwgt, *vwgt, *where; + idx_t i, j, nvtxs, me, other; + idx_t *xadj, *adjncy, *vwgt, *where; idx_t edegrees[2], pwgts[3]; nvtxs = graph->nvtxs; xadj = graph->xadj; vwgt = graph->vwgt; adjncy = graph->adjncy; - adjwgt = graph->adjwgt; where = graph->where; /*------------------------------------------------------------ @@ -308,6 +317,7 @@ /*************************************************************************/ idx_t IsSeparable(graph_t *graph) { +#ifndef NDEBUG idx_t i, j, nvtxs, other; idx_t *xadj, *adjncy, *where; @@ -327,6 +337,7 @@ xadj[adjncy[j]+1]-xadj[adjncy[j]])); } } +#endif return 1; } @@ -338,8 +349,8 @@ /*************************************************************************/ void CheckKWayVolPartitionParams(ctrl_t *ctrl, graph_t *graph) { - idx_t i, ii, j, k, kk, l, nvtxs, nbnd, mincut, minvol, me, other, pid; - idx_t *xadj, *vsize, *adjncy, *pwgts, *where, *bndind, *bndptr; + idx_t i, ii, j, k, kk, nvtxs, me, other, pid; + idx_t *xadj, *vsize, *adjncy, *where; vkrinfo_t *rinfo, *myrinfo, *orinfo, tmprinfo; vnbr_t *mynbrs, *onbrs, *tmpnbrs; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/fm.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/fm.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/fm.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/fm.c 2024-03-16 01:25:07 @@ -28,7 +28,7 @@ /*************************************************************************/ void FM_2WayCutRefine(ctrl_t *ctrl, graph_t *graph, real_t *ntpwgts, idx_t niter) { - idx_t i, ii, j, k, kwgt, nvtxs, nbnd, nswaps, from, to, pass, me, limit, tmp; + idx_t i, ii, j, k, kwgt, nvtxs, nbnd, nswaps, from, to, pass, limit, tmp; idx_t *xadj, *vwgt, *adjncy, *adjwgt, *where, *id, *ed, *bndptr, *bndind, *pwgts; idx_t *moved, *swaps, *perm; rpq_t *queues[2]; @@ -207,7 +207,7 @@ void FM_Mc2WayCutRefine(ctrl_t *ctrl, graph_t *graph, real_t *ntpwgts, idx_t niter) { idx_t i, ii, j, k, l, kwgt, nvtxs, ncon, nbnd, nswaps, from, to, pass, - me, limit, tmp, cnum; + limit, tmp, cnum; idx_t *xadj, *adjncy, *vwgt, *adjwgt, *pwgts, *where, *id, *ed, *bndptr, *bndind; idx_t *moved, *swaps, *perm, *qnum; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/graph.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/graph.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/graph.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/graph.c 2024-03-16 01:25:07 @@ -17,8 +17,7 @@ graph_t *SetupGraph(ctrl_t *ctrl, idx_t nvtxs, idx_t ncon, idx_t *xadj, idx_t *adjncy, idx_t *vwgt, idx_t *vsize, idx_t *adjwgt) { - idx_t i, j, k, sum; - real_t *nvwgt; + idx_t i, j; graph_t *graph; /* allocate the graph and fill in the fields */ diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/initpart.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/initpart.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/initpart.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/initpart.c 2024-03-16 01:25:07 @@ -114,18 +114,15 @@ void RandomBisection(ctrl_t *ctrl, graph_t *graph, real_t *ntpwgts, idx_t niparts) { - idx_t i, ii, j, k, nvtxs, pwgts[2], zeromaxpwgt, from, me, - bestcut=0, icut, mincut, inbfs; - idx_t *xadj, *vwgt, *adjncy, *adjwgt, *where; + idx_t i, ii, nvtxs, pwgts[2], zeromaxpwgt, + bestcut=0, inbfs; + idx_t *vwgt, *where; idx_t *perm, *bestwhere; WCOREPUSH; nvtxs = graph->nvtxs; - xadj = graph->xadj; vwgt = graph->vwgt; - adjncy = graph->adjncy; - adjwgt = graph->adjwgt; Allocate2WayPartitionMemory(ctrl, graph); where = graph->where; @@ -191,9 +188,9 @@ { idx_t i, j, k, nvtxs, drain, nleft, first, last, pwgts[2], oneminpwgt, onemaxpwgt, - from, me, bestcut=0, icut, mincut, inbfs; - idx_t *xadj, *vwgt, *adjncy, *adjwgt, *where; - idx_t *queue, *touched, *gain, *bestwhere; + bestcut=0, inbfs; + idx_t *xadj, *vwgt, *adjncy, *where; + idx_t *queue, *touched, *bestwhere; WCOREPUSH; @@ -201,7 +198,6 @@ xadj = graph->xadj; vwgt = graph->vwgt; adjncy = graph->adjncy; - adjwgt = graph->adjwgt; Allocate2WayPartitionMemory(ctrl, graph); where = graph->where; @@ -325,7 +321,7 @@ void McRandomBisection(ctrl_t *ctrl, graph_t *graph, real_t *ntpwgts, idx_t niparts) { - idx_t i, ii, j, k, nvtxs, ncon, from, bestcut=0, mincut, inbfs, qnum; + idx_t i, ii, nvtxs, ncon, bestcut=0, inbfs, qnum; idx_t *bestwhere, *where, *perm, *counts; idx_t *vwgt; @@ -385,7 +381,7 @@ void McGrowBisection(ctrl_t *ctrl, graph_t *graph, real_t *ntpwgts, idx_t niparts) { - idx_t i, j, k, nvtxs, ncon, from, bestcut=0, mincut, inbfs; + idx_t nvtxs, bestcut=0, inbfs; idx_t *bestwhere, *where; WCOREPUSH; @@ -434,9 +430,9 @@ idx_t niparts) { idx_t i, j, k, nvtxs, drain, nleft, first, last, pwgts[2], oneminpwgt, - onemaxpwgt, from, me, bestcut=0, icut, mincut, inbfs; - idx_t *xadj, *vwgt, *adjncy, *adjwgt, *where, *bndind; - idx_t *queue, *touched, *gain, *bestwhere; + onemaxpwgt, bestcut=0, inbfs; + idx_t *xadj, *vwgt, *adjncy, *where, *bndind; + idx_t *queue, *touched, *bestwhere; WCOREPUSH; @@ -444,7 +440,6 @@ xadj = graph->xadj; vwgt = graph->vwgt; adjncy = graph->adjncy; - adjwgt = graph->adjwgt; bestwhere = iwspacemalloc(ctrl, nvtxs); queue = iwspacemalloc(ctrl, nvtxs); @@ -570,7 +565,7 @@ void GrowBisectionNode2(ctrl_t *ctrl, graph_t *graph, real_t *ntpwgts, idx_t niparts) { - idx_t i, j, k, nvtxs, bestcut=0, mincut, inbfs; + idx_t i, j, nvtxs, bestcut=0, inbfs; idx_t *xadj, *where, *bndind, *bestwhere; WCOREPUSH; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/kmetis.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/kmetis.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/kmetis.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/kmetis.c 2024-03-16 01:25:07 @@ -102,10 +102,9 @@ /*************************************************************************/ idx_t MlevelKWayPartitioning(ctrl_t *ctrl, graph_t *graph, idx_t *part) { - idx_t i, j, objval=0, curobj=0, bestobj=0; + idx_t i, objval=0, curobj=0, bestobj=0; real_t curbal=0.0, bestbal=0.0; graph_t *cgraph; - int status; for (i=0; i<ctrl->ncuts; i++) { @@ -171,7 +170,10 @@ /*************************************************************************/ void InitKWayPartitioning(ctrl_t *ctrl, graph_t *graph) { - idx_t i, ntrials, options[METIS_NOPTIONS], curobj=0, bestobj=0; + idx_t i, options[METIS_NOPTIONS], curobj=0; +#ifdef XXX + idx_t ntrials, bestobj=0; +#endif idx_t *bestwhere=NULL; real_t *ubvec=NULL; int status; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/kwayfm.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/kwayfm.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/kwayfm.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/kwayfm.c 2024-03-16 01:25:07 @@ -61,12 +61,12 @@ real_t ffactor, idx_t omode) { /* Common variables to all types of kway-refinement/balancing routines */ - idx_t i, ii, iii, j, k, l, pass, nvtxs, nparts, gain; + idx_t i, ii, iii, j, k, pass, nvtxs, nparts, gain; idx_t from, me, to, oldcut, vwgt; idx_t *xadj, *adjncy, *adjwgt; idx_t *where, *pwgts, *perm, *bndptr, *bndind, *minwgt, *maxwgt, *itpwgts; idx_t nmoved, nupd, *vstatus, *updptr, *updind; - idx_t maxndoms, *safetos=NULL, *nads=NULL, *doms=NULL, **adids=NULL, **adwgts=NULL; + idx_t maxndoms, *safetos=NULL, *nads=NULL, *doms=NULL, **adids=NULL; idx_t *bfslvl=NULL, *bfsind=NULL, *bfsmrk=NULL; idx_t bndtype = (omode == OMODE_REFINE ? BNDTYPE_REFINE : BNDTYPE_BALANCE); @@ -118,7 +118,6 @@ nads = ctrl->nads; adids = ctrl->adids; - adwgts = ctrl->adwgts; doms = iset(nparts, 0, ctrl->pvec1); } @@ -371,12 +370,12 @@ real_t ffactor, idx_t omode) { /* Common variables to all types of kway-refinement/balancing routines */ - idx_t i, ii, iii, j, k, l, pass, nvtxs, nparts, gain; + idx_t i, ii, iii, j, k, pass, nvtxs, nparts, gain; idx_t from, me, to, oldcut, vwgt; idx_t *xadj, *adjncy; - idx_t *where, *pwgts, *perm, *bndptr, *bndind, *minwgt, *maxwgt, *itpwgts; + idx_t *where, *pwgts, *perm, *bndind, *minwgt, *maxwgt, *itpwgts; idx_t nmoved, nupd, *vstatus, *updptr, *updind; - idx_t maxndoms, *safetos=NULL, *nads=NULL, *doms=NULL, **adids=NULL, **adwgts=NULL; + idx_t maxndoms, *safetos=NULL, *nads=NULL, *doms=NULL, **adids=NULL; idx_t *bfslvl=NULL, *bfsind=NULL, *bfsmrk=NULL; idx_t bndtype = (omode == OMODE_REFINE ? BNDTYPE_REFINE : BNDTYPE_BALANCE); @@ -393,7 +392,6 @@ nvtxs = graph->nvtxs; xadj = graph->xadj; adjncy = graph->adjncy; - bndptr = graph->bndptr; bndind = graph->bndind; where = graph->where; pwgts = graph->pwgts; @@ -425,7 +423,6 @@ nads = ctrl->nads; adids = ctrl->adids; - adwgts = ctrl->adwgts; doms = iset(nparts, 0, ctrl->pvec1); } @@ -685,12 +682,12 @@ real_t ffactor, idx_t omode) { /* Common variables to all types of kway-refinement/balancing routines */ - idx_t i, ii, iii, j, k, l, pass, nvtxs, ncon, nparts, gain; + idx_t i, ii, iii, j, k, pass, nvtxs, ncon, nparts, gain; idx_t from, me, to, cto, oldcut; idx_t *xadj, *vwgt, *adjncy, *adjwgt; idx_t *where, *pwgts, *perm, *bndptr, *bndind, *minwgt, *maxwgt; idx_t nmoved, nupd, *vstatus, *updptr, *updind; - idx_t maxndoms, *safetos=NULL, *nads=NULL, *doms=NULL, **adids=NULL, **adwgts=NULL; + idx_t maxndoms, *safetos=NULL, *nads=NULL, *doms=NULL, **adids=NULL; idx_t *bfslvl=NULL, *bfsind=NULL, *bfsmrk=NULL; idx_t bndtype = (omode == OMODE_REFINE ? BNDTYPE_REFINE : BNDTYPE_BALANCE); real_t *ubfactors, *pijbm; @@ -765,7 +762,6 @@ nads = ctrl->nads; adids = ctrl->adids; - adwgts = ctrl->adwgts; doms = iset(nparts, 0, ctrl->pvec1); } @@ -1027,12 +1023,12 @@ real_t ffactor, idx_t omode) { /* Common variables to all types of kway-refinement/balancing routines */ - idx_t i, ii, iii, j, k, l, pass, nvtxs, ncon, nparts, gain; + idx_t i, ii, iii, j, k, pass, nvtxs, ncon, nparts, gain; idx_t from, me, to, cto, oldcut; idx_t *xadj, *vwgt, *adjncy; - idx_t *where, *pwgts, *perm, *bndptr, *bndind, *minwgt, *maxwgt; + idx_t *where, *pwgts, *perm, *bndind, *minwgt, *maxwgt; idx_t nmoved, nupd, *vstatus, *updptr, *updind; - idx_t maxndoms, *safetos=NULL, *nads=NULL, *doms=NULL, **adids=NULL, **adwgts=NULL; + idx_t maxndoms, *safetos=NULL, *nads=NULL, *doms=NULL, **adids=NULL; idx_t *bfslvl=NULL, *bfsind=NULL, *bfsmrk=NULL; idx_t bndtype = (omode == OMODE_REFINE ? BNDTYPE_REFINE : BNDTYPE_BALANCE); real_t *ubfactors, *pijbm; @@ -1053,7 +1049,6 @@ xadj = graph->xadj; vwgt = graph->vwgt; adjncy = graph->adjncy; - bndptr = graph->bndptr; bndind = graph->bndind; where = graph->where; pwgts = graph->pwgts; @@ -1104,7 +1099,6 @@ nads = ctrl->nads; adids = ctrl->adids; - adwgts = ctrl->adwgts; doms = iset(nparts, 0, ctrl->pvec1); } @@ -1462,7 +1456,7 @@ idx_t *updind, idx_t bndtype, idx_t *vmarker, idx_t *pmarker, idx_t *modind) { - idx_t i, ii, iii, j, jj, k, kk, l, u, nmod, other, me, myidx; + idx_t i, ii, iii, j, jj, k, kk, u, nmod, other, me, myidx; idx_t *xadj, *vsize, *adjncy, *where; vkrinfo_t *myrinfo, *orinfo; vnbr_t *mynbrs, *onbrs; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/kwayrefine.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/kwayrefine.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/kwayrefine.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/kwayrefine.c 2024-03-16 01:25:07 @@ -148,7 +148,7 @@ /**************************************************************************/ void ComputeKWayPartitionParams(ctrl_t *ctrl, graph_t *graph) { - idx_t i, j, k, l, nvtxs, ncon, nparts, nbnd, mincut, me, other; + idx_t i, j, k, nvtxs, ncon, nparts, nbnd, mincut, me, other; idx_t *xadj, *vwgt, *adjncy, *adjwgt, *pwgts, *where, *bndind, *bndptr; nparts = ctrl->nparts; @@ -561,8 +561,8 @@ /*************************************************************************/ void ComputeKWayVolGains(ctrl_t *ctrl, graph_t *graph) { - idx_t i, ii, j, k, l, nvtxs, nparts, me, other, pid; - idx_t *xadj, *vsize, *adjncy, *adjwgt, *where, + idx_t i, ii, j, k, nvtxs, nparts, me, other; + idx_t *xadj, *vsize, *adjncy, *where, *bndind, *bndptr, *ophtable; vkrinfo_t *myrinfo, *orinfo; vnbr_t *mynbrs, *onbrs; @@ -575,7 +575,6 @@ xadj = graph->xadj; vsize = graph->vsize; adjncy = graph->adjncy; - adjwgt = graph->adjwgt; where = graph->where; bndind = graph->bndind; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/minconn.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/minconn.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/minconn.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/minconn.c 2024-03-16 01:28:20 @@ -18,16 +18,13 @@ void ComputeSubDomainGraph(ctrl_t *ctrl, graph_t *graph) { idx_t i, ii, j, pid, other, nparts, nvtxs, nnbrs; - idx_t *xadj, *adjncy, *adjwgt, *where; + idx_t *where; idx_t *pptr, *pind; idx_t nads=0, *vadids, *vadwgts; WCOREPUSH; nvtxs = graph->nvtxs; - xadj = graph->xadj; - adjncy = graph->adjncy; - adjwgt = graph->adjwgt; where = graph->where; nparts = ctrl->nparts; @@ -264,10 +261,12 @@ avg = total/nparts; max = nads[iargmax(nparts, nads)]; +#ifndef NDEBUG IFSET(ctrl->dbglvl, METIS_DBG_CONNINFO, printf("Adjacent Subdomain Stats: Total: %3"PRIDX", " "Max: %3"PRIDX"[%zu], Avg: %3"PRIDX"\n", total, max, iargmax(nparts, nads), avg)); +#endif if (max < badfactor*avg) break; @@ -477,12 +476,11 @@ void MoveGroupMinConnForCut(ctrl_t *ctrl, graph_t *graph, idx_t to, idx_t nind, idx_t *ind) { - idx_t i, ii, j, jj, k, l, nvtxs, nbnd, from, me; + idx_t i, ii, j, k, nbnd, from, me; idx_t *xadj, *adjncy, *adjwgt, *where, *bndptr, *bndind; ckrinfo_t *myrinfo; cnbr_t *mynbrs; - nvtxs = graph->nvtxs; xadj = graph->xadj; adjncy = graph->adjncy; adjwgt = graph->adjwgt; @@ -561,12 +559,11 @@ void MoveGroupMinConnForVol(ctrl_t *ctrl, graph_t *graph, idx_t to, idx_t nind, idx_t *ind, idx_t *vmarker, idx_t *pmarker, idx_t *modind) { - idx_t i, ii, j, jj, k, l, nvtxs, from, me, other, xgain, ewgt; + idx_t i, ii, j, k, l, from, me, other, xgain, ewgt; idx_t *xadj, *vsize, *adjncy, *where; vkrinfo_t *myrinfo, *orinfo; vnbr_t *mynbrs, *onbrs; - nvtxs = graph->nvtxs; xadj = graph->xadj; vsize = graph->vsize; adjncy = graph->adjncy; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/ometis.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/ometis.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/ometis.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/ometis.c 2024-03-16 01:25:07 @@ -48,7 +48,6 @@ graph_t *graph=NULL; ctrl_t *ctrl; idx_t *cptr, *cind, *piperm; - int numflag = 0; /* set up malloc cleaning code and signal catchers */ if (!gk_malloc_init()) @@ -183,12 +182,10 @@ void MlevelNestedDissection(ctrl_t *ctrl, graph_t *graph, idx_t *order, idx_t lastvtx) { - idx_t i, j, nvtxs, nbnd; + idx_t i, nbnd; idx_t *label, *bndind; graph_t *lgraph, *rgraph; - nvtxs = graph->nvtxs; - MlevelNodeBisectionMultiple(ctrl, graph); IFSET(ctrl->dbglvl, METIS_DBG_SEPINFO, @@ -236,7 +233,7 @@ void MlevelNestedDissectionCC(ctrl_t *ctrl, graph_t *graph, idx_t *order, idx_t lastvtx) { - idx_t i, j, nvtxs, nbnd, ncmps, rnvtxs, snvtxs; + idx_t i, nvtxs, nbnd, ncmps, rnvtxs, snvtxs; idx_t *label, *bndind; idx_t *cptr, *cind; graph_t **sgraphs; @@ -423,7 +420,7 @@ graph_t **r_rgraph) { idx_t i, ii, j, k, l, istart, iend, mypart, nvtxs, snvtxs[3], snedges[3]; - idx_t *xadj, *vwgt, *adjncy, *adjwgt, *label, *where, *bndptr, *bndind; + idx_t *xadj, *vwgt, *adjncy, *label, *where, *bndptr, *bndind; idx_t *sxadj[2], *svwgt[2], *sadjncy[2], *sadjwgt[2], *slabel[2]; idx_t *rename; idx_t *auxadjncy; @@ -437,7 +434,6 @@ xadj = graph->xadj; vwgt = graph->vwgt; adjncy = graph->adjncy; - adjwgt = graph->adjwgt; label = graph->label; where = graph->where; bndptr = graph->bndptr; @@ -552,8 +548,8 @@ graph_t **SplitGraphOrderCC(ctrl_t *ctrl, graph_t *graph, idx_t ncmps, idx_t *cptr, idx_t *cind) { - idx_t i, ii, iii, j, k, l, istart, iend, mypart, nvtxs, snvtxs, snedges; - idx_t *xadj, *vwgt, *adjncy, *adjwgt, *label, *where, *bndptr, *bndind; + idx_t i, ii, iii, j, k, l, istart, iend, nvtxs, snvtxs, snedges; + idx_t *xadj, *vwgt, *adjncy, *label, *where, *bndptr, *bndind; idx_t *sxadj, *svwgt, *sadjncy, *sadjwgt, *slabel; idx_t *rename; idx_t *auxadjncy; @@ -567,7 +563,6 @@ xadj = graph->xadj; vwgt = graph->vwgt; adjncy = graph->adjncy; - adjwgt = graph->adjwgt; label = graph->label; where = graph->where; bndptr = graph->bndptr; @@ -654,7 +649,7 @@ /*************************************************************************/ void MMDOrder(ctrl_t *ctrl, graph_t *graph, idx_t *order, idx_t lastvtx) { - idx_t i, j, k, nvtxs, nofsub, firstvtx; + idx_t i, k, nvtxs, nofsub, firstvtx; idx_t *xadj, *adjncy, *label; idx_t *perm, *iperm, *head, *qsize, *list, *marker; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/parmetis.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/parmetis.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/parmetis.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/parmetis.c 2024-03-16 01:25:07 @@ -105,7 +105,7 @@ void MlevelNestedDissectionP(ctrl_t *ctrl, graph_t *graph, idx_t *order, idx_t lastvtx, idx_t npes, idx_t cpos, idx_t *sizes) { - idx_t i, j, nvtxs, nbnd; + idx_t i, nvtxs, nbnd; idx_t *label, *bndind; graph_t *lgraph, *rgraph; @@ -161,7 +161,6 @@ int METIS_ComputeVertexSeparator(idx_t *nvtxs, idx_t *xadj, idx_t *adjncy, idx_t *vwgt, idx_t *options, idx_t *r_sepsize, idx_t *part) { - idx_t i, j; graph_t *graph; ctrl_t *ctrl; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/pmetis.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/pmetis.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/pmetis.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/pmetis.c 2024-03-16 01:25:07 @@ -157,7 +157,7 @@ idx_t MlevelRecursiveBisection(ctrl_t *ctrl, graph_t *graph, idx_t nparts, idx_t *part, real_t *tpwgts, idx_t fpart) { - idx_t i, j, nvtxs, ncon, objval; + idx_t i, nvtxs, ncon, objval; idx_t *label, *where; graph_t *lgraph, *rgraph; real_t wsum, *tpwgts2; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/separator.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/separator.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/separator.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/separator.c 2024-03-16 01:25:07 @@ -20,7 +20,7 @@ **************************************************************************/ void ConstructSeparator(ctrl_t *ctrl, graph_t *graph) { - idx_t i, j, k, nvtxs, nbnd; + idx_t i, j, nvtxs, nbnd; idx_t *xadj, *where, *bndind; WCOREPUSH; @@ -70,7 +70,10 @@ { idx_t i, ii, j, jj, k, l, nvtxs, nbnd, bnvtxs[3], bnedges[2], csize; idx_t *xadj, *adjncy, *bxadj, *badjncy; - idx_t *where, *bndind, *bndptr, *vmap, *ivmap, *cover; + idx_t *where, *bndind, *vmap, *ivmap, *cover; +#ifndef NDEBUG + idx_t *bndptr; +#endif WCOREPUSH; @@ -80,7 +83,9 @@ nbnd = graph->nbnd; bndind = graph->bndind; +#ifndef NDEBUG bndptr = graph->bndptr; +#endif where = graph->where; vmap = iwspacemalloc(ctrl, nvtxs); diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/sfm.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/sfm.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/sfm.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/sfm.c 2024-03-16 01:25:07 @@ -476,7 +476,7 @@ void FM_2WayNodeBalance(ctrl_t *ctrl, graph_t *graph) { idx_t i, ii, j, k, jj, kk, nvtxs, nbnd, nswaps, gain; - idx_t badmaxpwgt, higain, oldgain, pass, to, other; + idx_t badmaxpwgt, higain, oldgain, to, other; idx_t *xadj, *vwgt, *adjncy, *where, *pwgts, *edegrees, *bndind, *bndptr; idx_t *perm, *moved; rpq_t *queue; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/srefine.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/srefine.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/srefine.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/srefine.c 2024-03-16 01:25:07 @@ -136,7 +136,7 @@ /*************************************************************************/ void Project2WayNodePartition(ctrl_t *ctrl, graph_t *graph) { - idx_t i, j, nvtxs; + idx_t i, nvtxs; idx_t *cmap, *where, *cwhere; graph_t *cgraph; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/stat.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/stat.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/stat.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/stat.c 2024-03-16 01:25:07 @@ -20,8 +20,8 @@ **************************************************************************/ void ComputePartitionInfoBipartite(graph_t *graph, idx_t nparts, idx_t *where) { - idx_t i, j, k, nvtxs, ncon, mustfree=0; - idx_t *xadj, *adjncy, *vwgt, *vsize, *adjwgt, *kpwgts, *tmpptr; + idx_t i, j, nvtxs, ncon, mustfree=0; + idx_t *xadj, *adjncy, *vwgt, *vsize, *adjwgt, *kpwgts; idx_t *padjncy, *padjwgt, *padjcut; nvtxs = graph->nvtxs; @@ -126,7 +126,6 @@ { idx_t i, j, nvtxs, ncon; idx_t *kpwgts, *vwgt; - real_t balance; nvtxs = graph->nvtxs; ncon = graph->ncon; diff -ruN src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/wspace.c src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/wspace.c --- src/SuiteSparse/CHOLMOD/SuiteSparse_metis/libmetis/wspace.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/SuiteSparse_metis/libmetis/wspace.c 2024-03-16 01:28:20 @@ -81,12 +81,14 @@ { gk_mcoreDestroy(&ctrl->mcore, ctrl->dbglvl&METIS_DBG_INFO); +#ifndef NDEBUG IFSET(ctrl->dbglvl, METIS_DBG_INFO, printf(" nbrpool statistics\n" " nbrpoolsize: %12zu nbrpoolcpos: %12zu\n" " nbrpoolreallocs: %12zu\n\n", ctrl->nbrpoolsize, ctrl->nbrpoolcpos, ctrl->nbrpoolreallocs)); +#endif gk_free((void **)&ctrl->cnbrpool, &ctrl->vnbrpool, LTERM); ctrl->nbrpoolsize = 0; diff -ruN src/SuiteSparse/CHOLMOD/Supernodal/cholmod_super_numeric.c src/SuiteSparse-patched/CHOLMOD/Supernodal/cholmod_super_numeric.c --- src/SuiteSparse/CHOLMOD/Supernodal/cholmod_super_numeric.c 2024-03-16 01:30:04 +++ src/SuiteSparse-patched/CHOLMOD/Supernodal/cholmod_super_numeric.c 2024-03-16 01:25:07 @@ -281,10 +281,6 @@ // supernodal numerical factorization, using template routine //-------------------------------------------------------------------------- - float s_beta [2] ; - s_beta [0] = (float) beta [0] ; - s_beta [1] = (float) beta [1] ; - switch ((A->xtype + A->dtype) % 8) { case CHOLMOD_REAL + CHOLMOD_SINGLE: diff -ruN src/SuiteSparse/CHOLMOD/Supernodal/cholmod_super_symbolic.c src/SuiteSparse-patched/CHOLMOD/Supernodal/cholmod_super_symbolic.c --- src/SuiteSparse/CHOLMOD/Supernodal/cholmod_super_symbolic.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Supernodal/cholmod_super_symbolic.c 2024-03-16 01:25:07 @@ -173,11 +173,13 @@ csize, maxcsize, ss, nscol0, nscol1, ns, nfsuper, newzeros, totzeros, merge, snext, esize, maxesize, nrelax0, nrelax1, nrelax2, Asorted ; int find_xsize ; +#if (defined(CHOLMOD_HAS_CUDA) && defined(CHOLMOD_INT64)) const char* env_use_gpu; const char* env_max_bytes; size_t max_bytes; const char* env_max_fraction; double max_fraction ; +#endif RETURN_IF_NULL_COMMON (FALSE) ; RETURN_IF_NULL (A, FALSE) ; @@ -243,10 +245,11 @@ // real/complex/zomples only), and only when the GPU is requested and // available. + #ifdef CHOLMOD_INT64 + max_bytes = 0; max_fraction = 0; - #ifdef CHOLMOD_INT64 if ( Common->useGPU == EMPTY ) { // useGPU not explicitly requested by the user, but not explicitly @@ -552,9 +555,6 @@ double xns = (double) ns ; double xtotsize = (xns * (xns+1) / 2) + xns * (lnz1 - nscol1) ; double z = xtotzeros / xtotsize ; - - Int totsize ; - totsize = (ns * (ns+1) / 2) + ns * (Snz [s+1] - nscol1) ; PRINT2 (("oldzeros "ID" newzeros "ID" xtotsize %g z %g\n", Zeros [s+1], newzeros, xtotsize, z)) ; diff -ruN src/SuiteSparse/CHOLMOD/Supernodal/t_cholmod_super_numeric_worker.c src/SuiteSparse-patched/CHOLMOD/Supernodal/t_cholmod_super_numeric_worker.c --- src/SuiteSparse/CHOLMOD/Supernodal/t_cholmod_super_numeric_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Supernodal/t_cholmod_super_numeric_worker.c 2024-03-16 01:25:07 @@ -152,15 +152,21 @@ #ifdef BLAS_TIMER double tstart, blas_time ; #endif - Real *Lx, *Ax, *Fx, *Az, *Fz, *C ; + Real *Lx, *Ax, *Fx, *C ; +#ifdef ZOMPLEX + Real *Az, *Fz ; +#endif Int *Super, *Head, *Ls, *Lpi, *Lpx, *Map, *SuperMap, *RelativeMap, *Next, - *Lpos, *Fp, *Fi, *Fnz, *Ap, *Ai, *Anz, *Iwork, *Next_save, *Lpos_save, - *Previous; + *Lpos, *Fp, *Fi, *Fnz, *Ap, *Ai, *Anz, *Iwork, *Next_save, *Lpos_save; Int nsuper, n, s, k1, k2, nscol, psi, psx, psend, nsrow, - pj, d, kd1, kd2, info, ndcol, ndrow, pdi, pdx, pdend, pdi1, pdi2, + d, kd1, kd2, info, ndcol, ndrow, pdi, pdx, pdend, pdi2, ndrow1, ndrow2, dancestor, sparent, dnext, nsrow2, ndrow3, stype, Apacked, Fpacked, repeat_supernode, nscol2, ss, - tail, nscol_new = 0; + nscol_new = 0; +#if (defined (CHOLMOD_HAS_CUDA) && defined (DOUBLE)) + Int *Previous ; + Int tail ; +#endif info = 0 ; ASSERT (L->dtype == A->dtype) ; @@ -215,7 +221,9 @@ Lpos = Iwork + 2*((size_t) n) + nsuper ; // size nsuper Next_save = Iwork + 2*((size_t) n) + 2*((size_t) nsuper) ;// size nsuper Lpos_save = Iwork + 2*((size_t) n) + 3*((size_t) nsuper) ;// size nsuper +#if (defined (CHOLMOD_HAS_CUDA) && defined (DOUBLE)) Previous = Iwork + 2*((size_t) n) + 4*((size_t) nsuper) ;// size nsuper +#endif Map = Common->Flag ; // size n, use Flag as workspace for Map array Head = Common->Head ; // size n+1, only Head [0..nsuper-1] used @@ -273,7 +281,9 @@ Fp = NULL ; Fi = NULL ; Fx = NULL ; +#ifdef ZOMPLEX Fz = NULL ; +#endif Fnz = NULL ; Fpacked = TRUE ; } @@ -284,7 +294,9 @@ Fp = F->p ; Fi = F->i ; Fx = F->x ; +#ifdef ZOMPLEX Fz = F->z ; +#endif Fnz = F->nz ; Fpacked = F->packed ; } @@ -292,7 +304,9 @@ Ap = A->p ; Ai = A->i ; Ax = A->x ; +#ifdef ZOMPLEX Az = A->z ; +#endif Anz = A->nz ; Apacked = A->packed ; @@ -303,8 +317,10 @@ #endif Int ii ; +#ifdef _OPENMP #pragma omp parallel for num_threads(nthreads) \ if ( n > 128 ) schedule (static) +#endif for (ii = 0 ; ii < n ; ii++) { Map [ii] = EMPTY ; @@ -370,8 +386,10 @@ #endif Int p ; +#ifdef _OPENMP #pragma omp parallel for num_threads(nthreads) \ schedule (static) if ( pend - psx > 1024 ) +#endif for (p = psx ; p < pend ; p++) { L_CLEAR (Lx,p) ; @@ -390,8 +408,10 @@ #endif Int k ; +#ifdef _OPENMP #pragma omp parallel for num_threads(nthreads) \ if ( nsrow > 128 ) +#endif for (k = 0 ; k < nsrow ; k++) { PRINT1 ((" "ID" map "ID"\n", Ls [psi+k], k)) ; @@ -416,8 +436,6 @@ // copy matrix into supernode s (lower triangular part only) //---------------------------------------------------------------------- - Int pk = psx ; - #ifdef _OPENMP double work ; if (stype != 0) @@ -435,8 +453,10 @@ nthreads = cholmod_nthreads (work, Common) ; #endif +#ifdef _OPENMP #pragma omp parallel for num_threads(nthreads) \ if ( k2-k1 > 64 ) +#endif for (k = k1 ; k < k2 ; k++) { if (stype != 0) @@ -894,8 +914,10 @@ #endif Int i ; +#ifdef _OPENMP #pragma omp parallel for num_threads(nthreads) \ if ( ndrow2 > 64 ) +#endif for (i = 0 ; i < ndrow2 ; i++) { RelativeMap [i] = Map [Ls [pdi1 + i]] ; @@ -912,8 +934,10 @@ #endif Int j ; +#ifdef _OPENMP #pragma omp parallel for num_threads(nthreads) \ if (ndrow1 > 64 ) +#endif for (j = 0 ; j < ndrow1 ; j++) // cols k1:k2-1 { ASSERT (RelativeMap [j] == Map [Ls [pdi1 + j]]) ; diff -ruN src/SuiteSparse/CHOLMOD/Supernodal/t_cholmod_super_solve_worker.c src/SuiteSparse-patched/CHOLMOD/Supernodal/t_cholmod_super_solve_worker.c --- src/SuiteSparse/CHOLMOD/Supernodal/t_cholmod_super_solve_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Supernodal/t_cholmod_super_solve_worker.c 2024-03-16 01:25:07 @@ -37,12 +37,11 @@ Real minus_one [2], one [2] ; Int *Lpi, *Lpx, *Ls, *Super ; Int nsuper, k1, k2, psi, psend, psx, nsrow, nscol, ii, s, - nsrow2, n, ps2, j, i, d, nrhs ; + nsrow2, ps2, j, i, d, nrhs ; nrhs = X->ncol ; Ex = E->x ; Xx = X->x ; - n = L->n ; d = X->d ; nsuper = L->nsuper ; @@ -326,12 +325,11 @@ Real minus_one [2], one [2] ; Int *Lpi, *Lpx, *Ls, *Super ; Int nsuper, k1, k2, psi, psend, psx, nsrow, nscol, ii, s, - nsrow2, n, ps2, j, i, d, nrhs ; + nsrow2, ps2, j, i, d, nrhs ; nrhs = X->ncol ; Ex = E->x ; Xx = X->x ; - n = L->n ; d = X->d ; nsuper = L->nsuper ; diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_aat.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_aat.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_aat.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_aat.c 2024-03-16 01:25:07 @@ -106,7 +106,6 @@ bool ignore_diag = (mode < 0) ; bool values = (mode > 0) && (A->xtype != CHOLMOD_PATTERN) ; int axtype = (values) ? A->xtype : CHOLMOD_PATTERN ; - bool conj = (mode >= 2) ; Int nrow = A->nrow ; Int ncol = A->ncol ; @@ -121,12 +120,8 @@ // get the sizes of the entries of C and A //-------------------------------------------------------------------------- - size_t ei = sizeof (Int) ; - size_t e = (A->dtype == CHOLMOD_SINGLE) ? sizeof (float) : sizeof (double) ; size_t ew = ((axtype == CHOLMOD_PATTERN) ? 0 : ((axtype == CHOLMOD_COMPLEX) ? 2 : 1)) ; - size_t ex = e * ew ; - size_t ez = e * ((axtype == CHOLMOD_ZOMPLEX) ? 1 : 0) ; //-------------------------------------------------------------------------- // allocate workspace diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_aat_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_aat_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_aat_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_aat_worker.c 2024-03-16 01:25:07 @@ -35,20 +35,32 @@ Int *Ap = (Int *) A->p ; Int *Anz = (Int *) A->nz ; Int *Ai = (Int *) A->i ; +#ifndef PATTERN Real *Ax = (Real *) A->x ; +#ifdef ZOMPLEX Real *Az = (Real *) A->z ; +#endif +#endif bool packed = A->packed ; Int *Fp = (Int *) F->p ; Int *Fi = (Int *) F->i ; +#ifndef PATTERN Real *Fx = (Real *) F->x ; +#ifdef ZOMPLEX Real *Fz = (Real *) F->z ; +#endif +#endif ASSERT (F->packed) ; Int *Cp = (Int *) C->p ; Int *Ci = (Int *) C->i ; +#ifndef PATTERN Real *Cx = (Real *) C->x ; +#ifdef ZOMPLEX Real *Cz = (Real *) C->z ; +#endif +#endif ASSERT (C->packed) ; //-------------------------------------------------------------------------- @@ -93,8 +105,12 @@ //------------------------------------------------------------------ Int t = Fi [pf] ; +#ifndef PATTERN Real fx [2] = {0,0} ; +#ifdef ZOMPLEX Real fz [1] = {0} ; +#endif +#endif ASSIGN (fx, fz, 0, Fx, Fz, pf) ; //------------------------------------------------------------------ diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_add.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_add.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_add.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_add.c 2024-03-16 01:25:07 @@ -118,12 +118,6 @@ // get the sizes of the entries of C, A, and B //-------------------------------------------------------------------------- - size_t ei = sizeof (Int) ; - size_t e = (dtype == CHOLMOD_SINGLE) ? sizeof (float) : sizeof (double) ; - size_t ex = e * ((xtype == CHOLMOD_PATTERN) ? 0 : - ((xtype == CHOLMOD_COMPLEX) ? 2 : 1)) ; - size_t ez = e * ((xtype == CHOLMOD_ZOMPLEX) ? 1 : 0) ; - //-------------------------------------------------------------------------- // convert/sort A and/or B, if needed //-------------------------------------------------------------------------- diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_band_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_band_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_band_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_band_worker.c 2024-03-16 01:25:07 @@ -30,14 +30,22 @@ Int *Ap = (Int *) A->p ; Int *Anz = (Int *) A->nz ; Int *Ai = (Int *) A->i ; +#ifndef PATTERN Real *Ax = (Real *) A->x ; +#ifdef ZOMPLEX Real *Az = (Real *) A->z ; +#endif +#endif bool packed = A->packed ; Int *Cp = (Int *) C->p ; Int *Ci = (Int *) C->i ; +#ifndef PATTERN Real *Cx = (Real *) C->x ; +#ifdef ZOMPLEX Real *Cz = (Real *) C->z ; +#endif +#endif // columns outside of j1:j2 have no entries in diagonals k1:k2 Int j1 = MAX (k1, 0) ; diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_change_factor.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_change_factor.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_change_factor.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_change_factor.c 2024-03-16 01:25:07 @@ -698,8 +698,6 @@ bool make_ldl = (!to_ll && L->is_ll) ; Int n = L->n ; - Int *Lp = (Int *) L->p ; - Int *Li = (Int *) L->i ; Int *Lnz = (Int *) L->nz ; bool grow = false ; @@ -900,11 +898,8 @@ // get inputs //-------------------------------------------------------------------------- - Int n = L->n ; Int nsuper = L->nsuper ; Int *Lpi = (Int *) L->pi ; - Int *Lpx = (Int *) L->px ; - Int *Ls = (Int *) L->s ; Int *Super = (Int *) L->super ; size_t ei = sizeof (Int) ; @@ -1088,7 +1083,6 @@ // get the sizes of the entries //-------------------------------------------------------------------------- - size_t ei = sizeof (Int) ; size_t e = (L->dtype == CHOLMOD_SINGLE) ? sizeof (float) : sizeof (double) ; size_t ex = e * ((to_xtype == CHOLMOD_COMPLEX) ? 2 : 1) ; size_t xs = L->xsize ; diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_change_factor_1_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_change_factor_1_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_change_factor_1_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_change_factor_1_worker.c 2024-03-16 01:25:07 @@ -36,7 +36,9 @@ Int *Lp = (Int *) L->p ; Int *Li = (Int *) L->i ; Real *Lx = (Real *) L->x ; +#ifdef ZOMPLEX Real *Lz = (Real *) L->z ; +#endif Int n = L->n ; //-------------------------------------------------------------------------- @@ -44,7 +46,9 @@ //-------------------------------------------------------------------------- Real onex [2] = {1,0} ; +#ifdef ZOMPLEX Real onez [1] = {0} ; +#endif for (Int j = 0 ; j < n ; j++) { diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_change_factor_2_template.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_change_factor_2_template.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_change_factor_2_template.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_change_factor_2_template.c 2024-03-16 01:25:07 @@ -34,7 +34,9 @@ //------------------------------------------------------------------------------ { +#ifndef IN_PLACE Int pnew = 0 ; +#endif if (make_ll) { diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_change_factor_2_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_change_factor_2_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_change_factor_2_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_change_factor_2_worker.c 2024-03-16 01:25:07 @@ -50,7 +50,9 @@ Int *Lp = (Int *) L->p ; Int *Li = (Int *) L->i ; Real *Lx = (Real *) L->x ; +#ifdef ZOMPLEX Real *Lz = (Real *) L->z ; +#endif Int *Lnz = (Int *) L->nz ; //-------------------------------------------------------------------------- diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_change_factor_3_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_change_factor_3_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_change_factor_3_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_change_factor_3_worker.c 2024-03-16 01:25:07 @@ -100,7 +100,6 @@ Int *Lp = (Int *) L->p ; // simplicial col pointers Int *Li = (Int *) L->i ; // simplicial row indices Int *Lnz = (Int *) L->nz ; // simplicial column counts - Int lnz = L->nzmax ; // size of Li, Lp, and Lx //---------------------------------------------------------------------- // convert supernodal LL' to simplicial LL' or LDL' (packed/unpacked) diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_change_xdtype.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_change_xdtype.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_change_xdtype.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_change_xdtype.c 2024-03-16 01:25:07 @@ -214,7 +214,7 @@ int output_dtype = to_xdtype & 4 ; // double or single if (output_xtype <= CHOLMOD_PATTERN || - L->is_super && output_xtype == CHOLMOD_ZOMPLEX) + (L->is_super && output_xtype == CHOLMOD_ZOMPLEX)) { // output_xtype not supported ERROR (CHOLMOD_INVALID, "invalid xtype") ; diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_copy_dense2_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_copy_dense2_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_copy_dense2_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_copy_dense2_worker.c 2024-03-16 01:25:07 @@ -36,9 +36,13 @@ //-------------------------------------------------------------------------- Real *Xx = (Real *) X->x ; +#ifdef ZOMPLEX Real *Xz = (Real *) X->z ; +#endif Real *Yx = (Real *) Y->x ; +#ifdef ZOMPLEX Real *Yz = (Real *) Y->z ; +#endif size_t nrow = X->nrow ; size_t ncol = X->ncol ; size_t xd = X->d ; @@ -50,9 +54,13 @@ size_t e = (X->dtype == CHOLMOD_SINGLE) ? sizeof (float) : sizeof (double) ; size_t fx = ((X->xtype == CHOLMOD_COMPLEX) ? 2 : 1) ; +#ifdef ZOMPLEX size_t fz = ((X->xtype == CHOLMOD_ZOMPLEX) ? 1 : 0) ; +#endif size_t e_fx_nrow = e * fx * nrow ; +#ifdef ZOMPLEX size_t e_fz_nrow = e * fz * nrow ; +#endif //-------------------------------------------------------------------------- // copy X = Y diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_copy_factor.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_copy_factor.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_copy_factor.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_copy_factor.c 2024-03-16 01:25:07 @@ -69,7 +69,6 @@ size_t e = (L->dtype == CHOLMOD_SINGLE) ? sizeof (float) : sizeof (double) ; size_t ex = e * ((L->xtype == CHOLMOD_PATTERN) ? 0 : ((L->xtype == CHOLMOD_COMPLEX) ? 2 : 1)) ; - size_t ez = e * ((L->xtype == CHOLMOD_ZOMPLEX) ? 1 : 0) ; //-------------------------------------------------------------------------- // allocate the new factor H, H->Perm, and H->ColCount diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_copy_factor_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_copy_factor_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_copy_factor_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_copy_factor_worker.c 2024-03-16 01:25:07 @@ -29,11 +29,15 @@ Int *Li = (Int *) L->i ; Int *Lnz = (Int *) L->nz ; Real *Lx = (Real *) L->x ; +#ifdef ZOMPLEX Real *Lz = (Real *) L->z ; +#endif Int *Hi = (Int *) H->i ; Real *Hx = (Real *) H->x ; +#ifdef ZOMPLEX Real *Hz = (Real *) H->z ; +#endif //-------------------------------------------------------------------------- // copy each column diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_copy_sparse_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_copy_sparse_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_copy_sparse_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_copy_sparse_worker.c 2024-03-16 01:25:07 @@ -31,13 +31,21 @@ Int *Ap = (Int *) A->p ; Int *Anz = (Int *) A->nz ; Int *Ai = (Int *) A->i ; +#ifndef PATTERN Real *Ax = (Real *) A->x ; +#ifdef ZOMPLEX Real *Az = (Real *) A->z ; +#endif +#endif size_t ncol = A->ncol ; Int *Ci = (Int *) C->i ; +#ifndef PATTERN Real *Cx = (Real *) C->x ; +#ifdef ZOMPLEX Real *Cz = (Real *) C->z ; +#endif +#endif //-------------------------------------------------------------------------- // copy the contents from A to C diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_copy_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_copy_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_copy_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_copy_worker.c 2024-03-16 01:25:07 @@ -45,14 +45,21 @@ Int *Ap = (Int *) A->p ; Int *Anz = (Int *) A->nz ; Int *Ai = (Int *) A->i ; +#ifndef PATTERN Real *Ax = (Real *) A->x ; +#ifdef ZOMPLEX Real *Az = (Real *) A->z ; +#endif +#endif bool packed = A->packed ; - Int *Cp = (Int *) C->p ; Int *Ci = (Int *) C->i ; +#ifndef PATTERN Real *Cx = (Real *) C->x ; +#ifdef ZOMPLEX Real *Cz = (Real *) C->z ; +#endif +#endif bool keep_diag = !ignore_diag ; //-------------------------------------------------------------------------- diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_dense_nnz_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_dense_nnz_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_dense_nnz_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_dense_nnz_worker.c 2024-03-16 01:25:07 @@ -21,7 +21,9 @@ //-------------------------------------------------------------------------- Real *Xx = (Real *) X->x ; +#ifdef ZOMPLEX Real *Xz = (Real *) X->z ; +#endif Int nrow = (Int) X->nrow ; Int ncol = (Int) X->ncol ; diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_dense_to_sparse_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_dense_to_sparse_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_dense_to_sparse_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_dense_to_sparse_worker.c 2024-03-16 01:25:07 @@ -34,7 +34,9 @@ //-------------------------------------------------------------------------- Real *Xx = (Real *) X->x ; +#ifdef ZOMPLEX Real *Xz = (Real *) X->z ; +#endif Int nrow = X->nrow ; Int ncol = X->ncol ; Int d = X->d ; @@ -42,7 +44,9 @@ Int *Cp = (Int *) C->p ; Int *Ci = (Int *) C->i ; Real *Cx = (Real *) C->x ; +#ifdef ZOMPLEX Real *Cz = (Real *) C->z ; +#endif bool pattern = (C->xtype == CHOLMOD_PATTERN) ; //-------------------------------------------------------------------------- diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_eye_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_eye_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_eye_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_eye_worker.c 2024-03-16 01:25:07 @@ -21,14 +21,18 @@ //-------------------------------------------------------------------------- Real *Xx = (Real *) X->x ; +#ifdef ZOMPLEX Real *Xz = (Real *) X->z ; +#endif Int nrow = (Int) X->nrow ; Int ncol = (Int) X->ncol ; Int n = MIN (nrow, ncol) ; Real onex [2] = {1,0} ; +#ifdef ZOMPLEX Real onez [1] = {0} ; +#endif for (Int k = 0 ; k < n ; k++) { diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_ones_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_ones_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_ones_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_ones_worker.c 2024-03-16 01:25:07 @@ -21,11 +21,15 @@ //-------------------------------------------------------------------------- Real *Xx = (Real *) X->x ; +#ifdef ZOMPLEX Real *Xz = (Real *) X->z ; +#endif size_t nzmax = X->nzmax ; Real onex [2] = {1,0} ; +#ifdef ZOMPLEX Real onez [1] = {0} ; +#endif for (Int k = 0 ; k < nzmax ; k++) { diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_pack_factor_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_pack_factor_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_pack_factor_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_pack_factor_worker.c 2024-03-16 01:25:07 @@ -29,7 +29,9 @@ Int *Lp = (Int *) L->p ; Int *Li = (Int *) L->i ; Real *Lx = (Real *) L->x ; +#ifdef ZOMPLEX Real *Lz = (Real *) L->z ; +#endif Int *Lnz = (Int *) L->nz ; Int *Lnext = (Int *) L->next ; diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_realloc.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_realloc.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_realloc.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_realloc.c 2024-03-16 01:25:07 @@ -34,7 +34,9 @@ int ok ; bool newly_allocated = (p == NULL) ; +#ifndef NDEBUG void *pold = p ; +#endif size_t nold = (*n) ; p = SuiteSparse_realloc (nnew, *n, size, p, &ok) ; diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_reallocate_column.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_reallocate_column.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_reallocate_column.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_reallocate_column.c 2024-03-16 01:25:07 @@ -84,7 +84,6 @@ need = MAX (need, 1) ; double slack = MAX (Common->grow1, 1.0) * ((double) need) + Common->grow2 ; slack = MIN (slack, (double) (n-j)) ; - size_t nslack = (size_t) floor (slack) ; need = MAX (need, slack) ; need = MAX (need, 1) ; need = MIN (need, n-j) ; diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_reallocate_column_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_reallocate_column_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_reallocate_column_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_reallocate_column_worker.c 2024-03-16 01:25:07 @@ -26,7 +26,9 @@ Int *Lnz = (Int *) L->nz ; Int *Li = (Int *) L->i ; Real *Lx = (Real *) L->x ; +#ifdef ZOMPLEX Real *Lz = (Real *) L->z ; +#endif Int len = Lnz [j] ; //-------------------------------------------------------------------------- diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_sort_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_sort_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_sort_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_sort_worker.c 2024-03-16 01:25:07 @@ -185,8 +185,12 @@ Int *Ap = (Int *) A->p ; Int *Ai = (Int *) A->i ; Int *Anz = (Int *) A->nz ; +#ifndef PATTERN Real *Ax = (Real *) A->x ; +#ifdef ZOMPLEX Real *Az = (Real *) A->z ; +#endif +#endif Int ncol = A->ncol ; bool packed = A->packed ; uint64_t seed = 42 ; @@ -207,7 +211,6 @@ if (i < ilast) { // sort Ai, Ax, Ax [pa:pend-1] according to row index Ai - Int jnz = pend - pa ; CM_QSRT (Ai, Ax, Az, pa, pend - pa, &seed) ; break ; } diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_sparse_to_dense_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_sparse_to_dense_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_sparse_to_dense_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_sparse_to_dense_worker.c 2024-03-16 01:25:07 @@ -30,13 +30,19 @@ //-------------------------------------------------------------------------- Real *Xx = (Real *) X->x ; +#ifdef ZOMPLEX Real *Xz = (Real *) X->z ; +#endif Int *Ap = (Int *) A->p ; Int *Ai = (Int *) A->i ; Int *Anz = (Int *) A->nz ; +#ifndef PATTERN Real *Ax = (Real *) A->x ; +#ifdef ZOMPLEX Real *Az = (Real *) A->z ; +#endif +#endif Int nrow = (Int) A->nrow ; Int ncol = (Int) A->ncol ; diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_sparse_to_triplet_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_sparse_to_triplet_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_sparse_to_triplet_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_sparse_to_triplet_worker.c 2024-03-16 01:25:07 @@ -21,19 +21,26 @@ // get inputs //-------------------------------------------------------------------------- - Int nrow = A->nrow ; Int ncol = A->ncol ; bool packed = (bool) A->packed ; Int *Ap = (Int *) A->p ; Int *Ai = (Int *) A->i ; +#ifndef PATTERN Real *Ax = (Real *) A->x ; +#ifdef ZOMPLEX Real *Az = (Real *) A->z ; +#endif +#endif Int *Anz = (Int *) A->nz ; Int *Ti = (Int *) T->i ; Int *Tj = (Int *) T->j ; +#ifndef PATTERN Real *Tx = (Real *) T->x ; +#ifdef ZOMPLEX Real *Tz = (Real *) T->z ; +#endif +#endif Int k = 0 ; diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_speye_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_speye_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_speye_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_speye_worker.c 2024-03-16 01:25:07 @@ -22,15 +22,23 @@ Int *Ap = (Int *) A->p ; Int *Ai = (Int *) A->i ; +#ifndef PATTERN Real *Ax = (Real *) A->x ; +#ifdef ZOMPLEX Real *Az = (Real *) A->z ; +#endif +#endif Int ncol = (Int) A->ncol ; Int nrow = (Int) A->nrow ; Int n = MIN (nrow, ncol) ; +#ifndef PATTERN Real onex [2] = {1,0} ; +#ifdef ZOMPLEX Real onez [1] = {0} ; +#endif +#endif for (Int k = 0 ; k < n ; k++) { diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_transpose_sym_permuted.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_transpose_sym_permuted.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_transpose_sym_permuted.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_transpose_sym_permuted.c 2024-03-16 01:25:07 @@ -44,7 +44,10 @@ #endif { // C(jnew,inew) = conj (A(iold,jold)) - Int pc = Wi [inew]++ ; +#ifdef NUMERIC + Int pc = +#endif + Wi [inew]++ ; #ifdef NUMERIC ASSIGN_CONJ_OR_NCONJ (Cx, Cz, pc, Ax, Az, pa) ; Ci [pc] = jnew ; @@ -53,7 +56,10 @@ else { // C(inew,jnew) = A(iold,jold) - Int pc = Wi [jnew]++ ; +#ifdef NUMERIC + Int pc = +#endif + Wi [jnew]++ ; #ifdef NUMERIC ASSIGN (Cx, Cz, pc, Ax, Az, pa) ; Ci [pc] = inew ; diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_transpose_sym_unpermuted.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_transpose_sym_unpermuted.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_transpose_sym_unpermuted.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_transpose_sym_unpermuted.c 2024-03-16 01:25:07 @@ -40,7 +40,10 @@ if (i > j) continue ; #endif // C(j,i) = conj (A(i,j)) - Int pc = Wi [i]++ ; +#ifdef NUMERIC + Int pc = +#endif + Wi [i]++ ; #ifdef NUMERIC ASSIGN_CONJ_OR_NCONJ (Cx, Cz, pc, Ax, Az, pa) ; Ci [pc] = j ; diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_transpose_sym_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_transpose_sym_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_transpose_sym_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_transpose_sym_worker.c 2024-03-16 01:25:07 @@ -28,13 +28,20 @@ Int *Ap = (Int *) A->p ; Int *Ai = (Int *) A->i ; Int *Anz = (Int *) A->nz ; +#ifndef PATTERN Real *Ax = (Real *) A->x ; +#ifdef ZOMPLEX Real *Az = (Real *) A->z ; +#endif +#endif - Int *Cp = (Int *) C->p ; Int *Ci = (Int *) C->i ; +#ifndef PATTERN Real *Cx = (Real *) C->x ; +#ifdef ZOMPLEX Real *Cz = (Real *) C->z ; +#endif +#endif //-------------------------------------------------------------------------- // compute pattern and values of C diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_transpose_unsym_template.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_transpose_unsym_template.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_transpose_unsym_template.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_transpose_unsym_template.c 2024-03-16 01:25:07 @@ -59,7 +59,10 @@ for ( ; p < pend ; p++) { // get A(i,j) and count it or get its place in C - Int pc = Wi [Ai [p]]++ ; +#ifdef NUMERIC + Int pc = +#endif + Wi [Ai [p]]++ ; #ifdef NUMERIC // C(j,i) = conj (A(i,j)) ASSIGN_CONJ_OR_NCONJ (Cx, Cz, pc, Ax, Az, p) ; diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_transpose_unsym_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_transpose_unsym_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_transpose_unsym_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_transpose_unsym_worker.c 2024-03-16 01:25:07 @@ -27,14 +27,21 @@ Int *Ap = (Int *) A->p ; Int *Ai = (Int *) A->i ; Int *Anz = (Int *) A->nz ; +#ifndef PATTERN Real *Ax = (Real *) A->x ; +#ifdef ZOMPLEX Real *Az = (Real *) A->z ; +#endif +#endif Int ncol = A->ncol ; - Int *Cp = (Int *) C->p ; Int *Ci = (Int *) C->i ; +#ifndef PATTERN Real *Cx = (Real *) C->x ; +#ifdef ZOMPLEX Real *Cz = (Real *) C->z ; +#endif +#endif //-------------------------------------------------------------------------- // compute the pattern and values of C diff -ruN src/SuiteSparse/CHOLMOD/Utility/t_cholmod_triplet_to_sparse_worker.c src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_triplet_to_sparse_worker.c --- src/SuiteSparse/CHOLMOD/Utility/t_cholmod_triplet_to_sparse_worker.c 2024-01-24 19:45:24 +++ src/SuiteSparse-patched/CHOLMOD/Utility/t_cholmod_triplet_to_sparse_worker.c 2024-03-16 01:25:07 @@ -25,13 +25,21 @@ Int *Rp = (Int *) R->p ; Int *Ri = (Int *) R->i ; Int *Rnz = (Int *) R->nz ; +#ifndef PATTERN Real *Rx = (Real *) R->x ; +#ifdef ZOMPLEX Real *Rz = (Real *) R->z ; +#endif +#endif Int *Ti = (Int *) T->i ; Int *Tj = (Int *) T->j ; +#ifndef PATTERN Real *Tx = (Real *) T->x ; +#ifdef ZOMPLEX Real *Tz = (Real *) T->z ; +#endif +#endif size_t nrow = T->nrow ; size_t ncol = T->ncol ; Int nz = T->nnz ;