[Scipy-svn] r2665 - trunk/Lib/linsolve/SuperLU/SRC
scipy-svn@scip...
scipy-svn@scip...
Wed Jan 31 22:10:05 CST 2007
Author: wnbell
Date: 2007-01-31 22:10:04 -0600 (Wed, 31 Jan 2007)
New Revision: 2665
Modified:
trunk/Lib/linsolve/SuperLU/SRC/get_perm_c.c
Log:
SuperLU workaround by mauger, fixes ticket #311
Modified: trunk/Lib/linsolve/SuperLU/SRC/get_perm_c.c
===================================================================
--- trunk/Lib/linsolve/SuperLU/SRC/get_perm_c.c 2007-02-01 03:48:41 UTC (rev 2664)
+++ trunk/Lib/linsolve/SuperLU/SRC/get_perm_c.c 2007-02-01 04:10:04 UTC (rev 2665)
@@ -402,7 +402,7 @@
t = SuperLU_timer_();
/* Initialize and allocate storage for GENMMD. */
- delta = 1; /* DELTA is a parameter to allow the choice of nodes
+ delta = 0; /* DELTA is a parameter to allow the choice of nodes
whose degree <= min-degree + DELTA. */
maxint = 2147483647; /* 2**31 - 1 */
invp = (int *) SUPERLU_MALLOC((n+delta)*sizeof(int));
@@ -420,7 +420,7 @@
for (i = 0; i <= n; ++i) ++b_colptr[i];
for (i = 0; i < bnz; ++i) ++b_rowind[i];
- genmmd_(&n, b_colptr, b_rowind, perm_c, invp, &delta, dhead,
+ genmmd_(&n, b_colptr, b_rowind, invp, perm_c, &delta, dhead,
qsize, llist, marker, &maxint, &nofsub);
/* Transform perm_c into 0-based indexing. */
More information about the Scipy-svn
mailing list