[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