Skip to content

Let users set L/U storage estimates via superlumt_options_t (#27) #37

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
xiaoyeli merged 1 commit into from
Oct 10, 2026
Merged

Let users set L/U storage estimates via superlumt_options_t (#27) #37

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
3 changes: 3 additions & 0 deletions EXAMPLE/pclinsolx.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -146,6 +146,9 @@ int main(int argc, char *argv[])
superlumt_options.diag_pivot_thresh = u;
superlumt_options.SymmetricMode = NO;
superlumt_options.PrintStat = NO;
superlumt_options.fill_lusup = 0; /* use sp_ienv(6) */
superlumt_options.fill_ucol = 0; /* use sp_ienv(7) */
superlumt_options.fill_lsub = 0; /* use sp_ienv(8) */
superlumt_options.perm_c = perm_c;
superlumt_options.perm_r = perm_r;
superlumt_options.work = work;
Expand Down
3 changes: 3 additions & 0 deletions EXAMPLE/pclinsolx1.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -164,6 +164,9 @@ int main(int argc, char *argv[])
superlumt_options.diag_pivot_thresh = u;
superlumt_options.SymmetricMode = NO;
superlumt_options.PrintStat = NO;
superlumt_options.fill_lusup = 0; /* use sp_ienv(6) */
superlumt_options.fill_ucol = 0; /* use sp_ienv(7) */
superlumt_options.fill_lsub = 0; /* use sp_ienv(8) */
superlumt_options.perm_c = perm_c;
superlumt_options.perm_r = perm_r;
superlumt_options.work = work;
Expand Down
3 changes: 3 additions & 0 deletions EXAMPLE/pclinsolx2.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -157,6 +157,9 @@ int main(int argc, char *argv[])
superlumt_options.usepr = usepr;
superlumt_options.drop_tol = drop_tol;
superlumt_options.PrintStat = NO;
superlumt_options.fill_lusup = 0; /* use sp_ienv(6) */
superlumt_options.fill_ucol = 0; /* use sp_ienv(7) */
superlumt_options.fill_lsub = 0; /* use sp_ienv(8) */
superlumt_options.perm_c = perm_c;
superlumt_options.perm_r = perm_r;
superlumt_options.work = work;
Expand Down
3 changes: 3 additions & 0 deletions EXAMPLE/pdlinsolx.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -146,6 +146,9 @@ int main(int argc, char *argv[])
superlumt_options.diag_pivot_thresh = u;
superlumt_options.SymmetricMode = NO;
superlumt_options.PrintStat = NO;
superlumt_options.fill_lusup = 0; /* use sp_ienv(6) */
superlumt_options.fill_ucol = 0; /* use sp_ienv(7) */
superlumt_options.fill_lsub = 0; /* use sp_ienv(8) */
superlumt_options.perm_c = perm_c;
superlumt_options.perm_r = perm_r;
superlumt_options.work = work;
Expand Down
3 changes: 3 additions & 0 deletions EXAMPLE/pdlinsolx1.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -164,6 +164,9 @@ int main(int argc, char *argv[])
superlumt_options.diag_pivot_thresh = u;
superlumt_options.SymmetricMode = NO;
superlumt_options.PrintStat = NO;
superlumt_options.fill_lusup = 0; /* use sp_ienv(6) */
superlumt_options.fill_ucol = 0; /* use sp_ienv(7) */
superlumt_options.fill_lsub = 0; /* use sp_ienv(8) */
superlumt_options.perm_c = perm_c;
superlumt_options.perm_r = perm_r;
superlumt_options.work = work;
Expand Down
3 changes: 3 additions & 0 deletions EXAMPLE/pdlinsolx2.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -157,6 +157,9 @@ int main(int argc, char *argv[])
superlumt_options.usepr = usepr;
superlumt_options.drop_tol = drop_tol;
superlumt_options.PrintStat = NO;
superlumt_options.fill_lusup = 0; /* use sp_ienv(6) */
superlumt_options.fill_ucol = 0; /* use sp_ienv(7) */
superlumt_options.fill_lsub = 0; /* use sp_ienv(8) */
superlumt_options.perm_c = perm_c;
superlumt_options.perm_r = perm_r;
superlumt_options.work = work;
Expand Down
3 changes: 3 additions & 0 deletions EXAMPLE/pslinsolx.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -146,6 +146,9 @@ int main(int argc, char *argv[])
superlumt_options.diag_pivot_thresh = u;
superlumt_options.SymmetricMode = NO;
superlumt_options.PrintStat = NO;
superlumt_options.fill_lusup = 0; /* use sp_ienv(6) */
superlumt_options.fill_ucol = 0; /* use sp_ienv(7) */
superlumt_options.fill_lsub = 0; /* use sp_ienv(8) */
superlumt_options.perm_c = perm_c;
superlumt_options.perm_r = perm_r;
superlumt_options.work = work;
Expand Down
3 changes: 3 additions & 0 deletions EXAMPLE/pslinsolx1.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -164,6 +164,9 @@ int main(int argc, char *argv[])
superlumt_options.diag_pivot_thresh = u;
superlumt_options.SymmetricMode = NO;
superlumt_options.PrintStat = NO;
superlumt_options.fill_lusup = 0; /* use sp_ienv(6) */
superlumt_options.fill_ucol = 0; /* use sp_ienv(7) */
superlumt_options.fill_lsub = 0; /* use sp_ienv(8) */
superlumt_options.perm_c = perm_c;
superlumt_options.perm_r = perm_r;
superlumt_options.work = work;
Expand Down
3 changes: 3 additions & 0 deletions EXAMPLE/pslinsolx2.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -157,6 +157,9 @@ int main(int argc, char *argv[])
superlumt_options.usepr = usepr;
superlumt_options.drop_tol = drop_tol;
superlumt_options.PrintStat = NO;
superlumt_options.fill_lusup = 0; /* use sp_ienv(6) */
superlumt_options.fill_ucol = 0; /* use sp_ienv(7) */
superlumt_options.fill_lsub = 0; /* use sp_ienv(8) */
superlumt_options.perm_c = perm_c;
superlumt_options.perm_r = perm_r;
superlumt_options.work = work;
Expand Down
3 changes: 3 additions & 0 deletions EXAMPLE/pzlinsolx.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -146,6 +146,9 @@ int main(int argc, char *argv[])
superlumt_options.diag_pivot_thresh = u;
superlumt_options.SymmetricMode = NO;
superlumt_options.PrintStat = NO;
superlumt_options.fill_lusup = 0; /* use sp_ienv(6) */
superlumt_options.fill_ucol = 0; /* use sp_ienv(7) */
superlumt_options.fill_lsub = 0; /* use sp_ienv(8) */
superlumt_options.perm_c = perm_c;
superlumt_options.perm_r = perm_r;
superlumt_options.work = work;
Expand Down
3 changes: 3 additions & 0 deletions EXAMPLE/pzlinsolx1.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -164,6 +164,9 @@ int main(int argc, char *argv[])
superlumt_options.diag_pivot_thresh = u;
superlumt_options.SymmetricMode = NO;
superlumt_options.PrintStat = NO;
superlumt_options.fill_lusup = 0; /* use sp_ienv(6) */
superlumt_options.fill_ucol = 0; /* use sp_ienv(7) */
superlumt_options.fill_lsub = 0; /* use sp_ienv(8) */
superlumt_options.perm_c = perm_c;
superlumt_options.perm_r = perm_r;
superlumt_options.work = work;
Expand Down
3 changes: 3 additions & 0 deletions EXAMPLE/pzlinsolx2.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -157,6 +157,9 @@ int main(int argc, char *argv[])
superlumt_options.usepr = usepr;
superlumt_options.drop_tol = drop_tol;
superlumt_options.PrintStat = NO;
superlumt_options.fill_lusup = 0; /* use sp_ienv(6) */
superlumt_options.fill_ucol = 0; /* use sp_ienv(7) */
superlumt_options.fill_lsub = 0; /* use sp_ienv(8) */
superlumt_options.perm_c = perm_c;
superlumt_options.perm_r = perm_r;
superlumt_options.work = work;
Expand Down
9 changes: 9 additions & 0 deletions SRC/pcgssvx.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -202,6 +202,15 @@ pcgssvx(int_t nprocs, superlumt_options_t *superlumt_options, SuperMatrix *A,
* performing the factorization, and returns it in
* superlu_memusage->total_needed; no other side effects.
*
* o fill_lusup, fill_ucol, fill_lsub (int_t)
* Storage estimates for the values of the L supernodes, the
* columns of U, and the row subscripts of L, respectively.
* A negative value is a fill-in growth factor: its magnitude
* times nnz(A) entries are allocated. A positive value is the
* number of entries to allocate. 0 means use the defaults from
* sp_ienv(6), sp_ienv(7) and sp_ienv(8). Increase these if the
* factorization aborts with "Storage for ... exceeded".
*
* -------------------------------------------------------------------------
* ** The following 3 arrays must be allocated before calling this routine.
* -------------------------------------------------------------------------
Expand Down
3 changes: 3 additions & 0 deletions SRC/pcgstrf_init.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -138,6 +138,9 @@ pcgstrf_init(int_t nprocs, fact_t fact, trans_t trans, yes_no_t refact,
superlumt_options->drop_tol = drop_tol;
superlumt_options->SymmetricMode = NO;
superlumt_options->PrintStat = NO;
superlumt_options->fill_lusup = 0; /* use sp_ienv(6) */
superlumt_options->fill_ucol = 0; /* use sp_ienv(7) */
superlumt_options->fill_lsub = 0; /* use sp_ienv(8) */

/*
* The following should be retained for repeated factorizations.
Expand Down
19 changes: 13 additions & 6 deletions SRC/pcmemory.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -269,10 +269,16 @@ pcgstrf_MemInit(int_t n, int_t annz, superlumt_options_t *superlumt_options,
complex *ucol;
int_t *usub, *xusub, *xusub_end;
int_t nzlmax, nzumax, nzlumax;
int_t FILL_LUSUP = sp_ienv(6); /* Guess the fill-in growth for LUSUP */
int_t FILL_UCOL = sp_ienv(7); /* Guess the fill-in growth for UCOL */
int_t FILL_LSUB = sp_ienv(8); /* Guess the fill-in growth for LSUB */
int_t FILL_LUSUP, FILL_UCOL, FILL_LSUB; /* Guess the fill-in growth */

/* Storage estimates for L\U: from the options if set, else sp_ienv(). */
FILL_LUSUP = superlumt_options->fill_lusup;
FILL_UCOL = superlumt_options->fill_ucol;
FILL_LSUB = superlumt_options->fill_lsub;
if ( FILL_LUSUP == 0 ) FILL_LUSUP = sp_ienv(6);
if ( FILL_UCOL == 0 ) FILL_UCOL = sp_ienv(7);
if ( FILL_LSUB == 0 ) FILL_LSUB = sp_ienv(8);

no_expand = 0;
ndim = n;
iword = sizeof(int_t);
Expand Down Expand Up @@ -334,9 +340,10 @@ pcgstrf_MemInit(int_t n, int_t annz, superlumt_options_t *superlumt_options,

while ( !ucol || !lsub || !usub ) {
/*SUPERLU_ABORT("Not enough core in LUMemInit()");*/
#if (PRNTlevel==1)
printf(".. pcgstrf_MemInit(): #retries " IFMT "\n", ++retries);
#endif
++retries;
fprintf(stderr, "** pcgstrf_MemInit(): cannot allocate storage for L\\U"
" (nzlmax " IFMT ", nzumax " IFMT "); halving the request,"
" retry " IFMT "\n", nzlmax, nzumax, retries);
if ( whichspace == SYSTEM ) {
SUPERLU_FREE(ucol);
SUPERLU_FREE(lsub);
Expand Down
9 changes: 9 additions & 0 deletions SRC/pdgssvx.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -202,6 +202,15 @@ pdgssvx(int_t nprocs, superlumt_options_t *superlumt_options, SuperMatrix *A,
* performing the factorization, and returns it in
* superlu_memusage->total_needed; no other side effects.
*
* o fill_lusup, fill_ucol, fill_lsub (int_t)
* Storage estimates for the values of the L supernodes, the
* columns of U, and the row subscripts of L, respectively.
* A negative value is a fill-in growth factor: its magnitude
* times nnz(A) entries are allocated. A positive value is the
* number of entries to allocate. 0 means use the defaults from
* sp_ienv(6), sp_ienv(7) and sp_ienv(8). Increase these if the
* factorization aborts with "Storage for ... exceeded".
*
* -------------------------------------------------------------------------
* ** The following 3 arrays must be allocated before calling this routine.
* -------------------------------------------------------------------------
Expand Down
3 changes: 3 additions & 0 deletions SRC/pdgstrf_init.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -138,6 +138,9 @@ pdgstrf_init(int_t nprocs, fact_t fact, trans_t trans, yes_no_t refact,
superlumt_options->drop_tol = drop_tol;
superlumt_options->SymmetricMode = NO;
superlumt_options->PrintStat = NO;
superlumt_options->fill_lusup = 0; /* use sp_ienv(6) */
superlumt_options->fill_ucol = 0; /* use sp_ienv(7) */
superlumt_options->fill_lsub = 0; /* use sp_ienv(8) */

/*
* The following should be retained for repeated factorizations.
Expand Down
19 changes: 13 additions & 6 deletions SRC/pdmemory.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -269,10 +269,16 @@ pdgstrf_MemInit(int_t n, int_t annz, superlumt_options_t *superlumt_options,
double *ucol;
int_t *usub, *xusub, *xusub_end;
int_t nzlmax, nzumax, nzlumax;
int_t FILL_LUSUP = sp_ienv(6); /* Guess the fill-in growth for LUSUP */
int_t FILL_UCOL = sp_ienv(7); /* Guess the fill-in growth for UCOL */
int_t FILL_LSUB = sp_ienv(8); /* Guess the fill-in growth for LSUB */
int_t FILL_LUSUP, FILL_UCOL, FILL_LSUB; /* Guess the fill-in growth */

/* Storage estimates for L\U: from the options if set, else sp_ienv(). */
FILL_LUSUP = superlumt_options->fill_lusup;
FILL_UCOL = superlumt_options->fill_ucol;
FILL_LSUB = superlumt_options->fill_lsub;
if ( FILL_LUSUP == 0 ) FILL_LUSUP = sp_ienv(6);
if ( FILL_UCOL == 0 ) FILL_UCOL = sp_ienv(7);
if ( FILL_LSUB == 0 ) FILL_LSUB = sp_ienv(8);

no_expand = 0;
ndim = n;
iword = sizeof(int_t);
Expand Down Expand Up @@ -334,9 +340,10 @@ pdgstrf_MemInit(int_t n, int_t annz, superlumt_options_t *superlumt_options,

while ( !ucol || !lsub || !usub ) {
/*SUPERLU_ABORT("Not enough core in LUMemInit()");*/
#if (PRNTlevel==1)
printf(".. pdgstrf_MemInit(): #retries " IFMT "\n", ++retries);
#endif
++retries;
fprintf(stderr, "** pdgstrf_MemInit(): cannot allocate storage for L\\U"
" (nzlmax " IFMT ", nzumax " IFMT "); halving the request,"
" retry " IFMT "\n", nzlmax, nzumax, retries);
if ( whichspace == SYSTEM ) {
SUPERLU_FREE(ucol);
SUPERLU_FREE(lsub);
Expand Down
9 changes: 9 additions & 0 deletions SRC/psgssvx.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -202,6 +202,15 @@ psgssvx(int_t nprocs, superlumt_options_t *superlumt_options, SuperMatrix *A,
* performing the factorization, and returns it in
* superlu_memusage->total_needed; no other side effects.
*
* o fill_lusup, fill_ucol, fill_lsub (int_t)
* Storage estimates for the values of the L supernodes, the
* columns of U, and the row subscripts of L, respectively.
* A negative value is a fill-in growth factor: its magnitude
* times nnz(A) entries are allocated. A positive value is the
* number of entries to allocate. 0 means use the defaults from
* sp_ienv(6), sp_ienv(7) and sp_ienv(8). Increase these if the
* factorization aborts with "Storage for ... exceeded".
*
* -------------------------------------------------------------------------
* ** The following 3 arrays must be allocated before calling this routine.
* -------------------------------------------------------------------------
Expand Down
3 changes: 3 additions & 0 deletions SRC/psgstrf_init.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -138,6 +138,9 @@ psgstrf_init(int_t nprocs, fact_t fact, trans_t trans, yes_no_t refact,
superlumt_options->drop_tol = drop_tol;
superlumt_options->SymmetricMode = NO;
superlumt_options->PrintStat = NO;
superlumt_options->fill_lusup = 0; /* use sp_ienv(6) */
superlumt_options->fill_ucol = 0; /* use sp_ienv(7) */
superlumt_options->fill_lsub = 0; /* use sp_ienv(8) */

/*
* The following should be retained for repeated factorizations.
Expand Down
19 changes: 13 additions & 6 deletions SRC/psmemory.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -269,10 +269,16 @@ psgstrf_MemInit(int_t n, int_t annz, superlumt_options_t *superlumt_options,
float *ucol;
int_t *usub, *xusub, *xusub_end;
int_t nzlmax, nzumax, nzlumax;
int_t FILL_LUSUP = sp_ienv(6); /* Guess the fill-in growth for LUSUP */
int_t FILL_UCOL = sp_ienv(7); /* Guess the fill-in growth for UCOL */
int_t FILL_LSUB = sp_ienv(8); /* Guess the fill-in growth for LSUB */
int_t FILL_LUSUP, FILL_UCOL, FILL_LSUB; /* Guess the fill-in growth */

/* Storage estimates for L\U: from the options if set, else sp_ienv(). */
FILL_LUSUP = superlumt_options->fill_lusup;
FILL_UCOL = superlumt_options->fill_ucol;
FILL_LSUB = superlumt_options->fill_lsub;
if ( FILL_LUSUP == 0 ) FILL_LUSUP = sp_ienv(6);
if ( FILL_UCOL == 0 ) FILL_UCOL = sp_ienv(7);
if ( FILL_LSUB == 0 ) FILL_LSUB = sp_ienv(8);

no_expand = 0;
ndim = n;
iword = sizeof(int_t);
Expand Down Expand Up @@ -334,9 +340,10 @@ psgstrf_MemInit(int_t n, int_t annz, superlumt_options_t *superlumt_options,

while ( !ucol || !lsub || !usub ) {
/*SUPERLU_ABORT("Not enough core in LUMemInit()");*/
#if (PRNTlevel==1)
printf(".. psgstrf_MemInit(): #retries " IFMT "\n", ++retries);
#endif
++retries;
fprintf(stderr, "** psgstrf_MemInit(): cannot allocate storage for L\\U"
" (nzlmax " IFMT ", nzumax " IFMT "); halving the request,"
" retry " IFMT "\n", nzlmax, nzumax, retries);
if ( whichspace == SYSTEM ) {
SUPERLU_FREE(ucol);
SUPERLU_FREE(lsub);
Expand Down
9 changes: 9 additions & 0 deletions SRC/pzgssvx.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -202,6 +202,15 @@ pzgssvx(int_t nprocs, superlumt_options_t *superlumt_options, SuperMatrix *A,
* performing the factorization, and returns it in
* superlu_memusage->total_needed; no other side effects.
*
* o fill_lusup, fill_ucol, fill_lsub (int_t)
* Storage estimates for the values of the L supernodes, the
* columns of U, and the row subscripts of L, respectively.
* A negative value is a fill-in growth factor: its magnitude
* times nnz(A) entries are allocated. A positive value is the
* number of entries to allocate. 0 means use the defaults from
* sp_ienv(6), sp_ienv(7) and sp_ienv(8). Increase these if the
* factorization aborts with "Storage for ... exceeded".
*
* -------------------------------------------------------------------------
* ** The following 3 arrays must be allocated before calling this routine.
* -------------------------------------------------------------------------
Expand Down
3 changes: 3 additions & 0 deletions SRC/pzgstrf_init.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -138,6 +138,9 @@ pzgstrf_init(int_t nprocs, fact_t fact, trans_t trans, yes_no_t refact,
superlumt_options->drop_tol = drop_tol;
superlumt_options->SymmetricMode = NO;
superlumt_options->PrintStat = NO;
superlumt_options->fill_lusup = 0; /* use sp_ienv(6) */
superlumt_options->fill_ucol = 0; /* use sp_ienv(7) */
superlumt_options->fill_lsub = 0; /* use sp_ienv(8) */

/*
* The following should be retained for repeated factorizations.
Expand Down
19 changes: 13 additions & 6 deletions SRC/pzmemory.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -269,10 +269,16 @@ pzgstrf_MemInit(int_t n, int_t annz, superlumt_options_t *superlumt_options,
doublecomplex *ucol;
int_t *usub, *xusub, *xusub_end;
int_t nzlmax, nzumax, nzlumax;
int_t FILL_LUSUP = sp_ienv(6); /* Guess the fill-in growth for LUSUP */
int_t FILL_UCOL = sp_ienv(7); /* Guess the fill-in growth for UCOL */
int_t FILL_LSUB = sp_ienv(8); /* Guess the fill-in growth for LSUB */
int_t FILL_LUSUP, FILL_UCOL, FILL_LSUB; /* Guess the fill-in growth */

/* Storage estimates for L\U: from the options if set, else sp_ienv(). */
FILL_LUSUP = superlumt_options->fill_lusup;
FILL_UCOL = superlumt_options->fill_ucol;
FILL_LSUB = superlumt_options->fill_lsub;
if ( FILL_LUSUP == 0 ) FILL_LUSUP = sp_ienv(6);
if ( FILL_UCOL == 0 ) FILL_UCOL = sp_ienv(7);
if ( FILL_LSUB == 0 ) FILL_LSUB = sp_ienv(8);

no_expand = 0;
ndim = n;
iword = sizeof(int_t);
Expand Down Expand Up @@ -334,9 +340,10 @@ pzgstrf_MemInit(int_t n, int_t annz, superlumt_options_t *superlumt_options,

while ( !ucol || !lsub || !usub ) {
/*SUPERLU_ABORT("Not enough core in LUMemInit()");*/
#if (PRNTlevel==1)
printf(".. pzgstrf_MemInit(): #retries " IFMT "\n", ++retries);
#endif
++retries;
fprintf(stderr, "** pzgstrf_MemInit(): cannot allocate storage for L\\U"
" (nzlmax " IFMT ", nzumax " IFMT "); halving the request,"
" retry " IFMT "\n", nzlmax, nzumax, retries);
if ( whichspace == SYSTEM ) {
SUPERLU_FREE(ucol);
SUPERLU_FREE(lsub);
Expand Down
9 changes: 9 additions & 0 deletions SRC/slu_mt_util.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -292,6 +292,15 @@ typedef struct {
yes_no_t SymmetricMode;
yes_no_t PrintStat;

/* Storage estimates for the L and U factors, used by p?gstrf_MemInit().
Each behaves like the corresponding sp_ienv() parameter: a negative
value is a fill-in growth factor whose magnitude is multiplied by
nnz(A); a positive value is an absolute number of entries; 0 (the
default) means use sp_ienv(6), sp_ienv(7) or sp_ienv(8). */
int_t fill_lusup; /* values of the L supernodes; sp_ienv(6) */
int_t fill_ucol; /* columns of U; sp_ienv(7) */
int_t fill_lsub; /* row subscripts of L; sp_ienv(8) */

/* The following arrays are persistent during repeated factorizations. */
int_t *perm_c;
int_t *perm_r;
Expand Down
3 changes: 3 additions & 0 deletions TESTING/pcdrive.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -204,6 +204,9 @@ int main(int argc, char *argv[])
superlumt_options.relax = relax;
superlumt_options.diag_pivot_thresh = diag_pivot_thresh;
superlumt_options.drop_tol = drop_tol;
superlumt_options.fill_lusup = 0; /* use sp_ienv(6) */
superlumt_options.fill_ucol = 0; /* use sp_ienv(7) */
superlumt_options.fill_lsub = 0; /* use sp_ienv(8) */
superlumt_options.perm_c = perm_c;
superlumt_options.perm_r = perm_r;
superlumt_options.work = work;
Expand Down
3 changes: 3 additions & 0 deletions TESTING/pddrive.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -204,6 +204,9 @@ int main(int argc, char *argv[])
superlumt_options.relax = relax;
superlumt_options.diag_pivot_thresh = diag_pivot_thresh;
superlumt_options.drop_tol = drop_tol;
superlumt_options.fill_lusup = 0; /* use sp_ienv(6) */
superlumt_options.fill_ucol = 0; /* use sp_ienv(7) */
superlumt_options.fill_lsub = 0; /* use sp_ienv(8) */
superlumt_options.perm_c = perm_c;
superlumt_options.perm_r = perm_r;
superlumt_options.work = work;
Expand Down
Loading
Loading