In the ILU setting, allow users to allocate and fill the supno/xsup arrays, defining the supernode parition.

Author: xiaoyeli

SRC/prec-independent/ilu_level_symbfact.c

@@ -76,14 +76,16 @@ int_t ilu_level_symbfact
     int_t nnzL = 0, nnzU = 0;
 
     /* Set up supernode partition */
-    Glu_persist->supno = (int_t *) SUPERLU_MALLOC(n * sizeof(int_t));
-    Glu_persist->xsup = (int_t *) SUPERLU_MALLOC((n+1) * sizeof(int_t));
     if ( options->UserDefineSupernode == NO ) {
+	Glu_persist->supno = (int_t *) SUPERLU_MALLOC(n * sizeof(int_t));
+	Glu_persist->xsup = (int_t *) SUPERLU_MALLOC((n+1) * sizeof(int_t));
 	for (i = 0; i < n; i++) { /* set up trivial supernode for now. */
 	    Glu_persist->supno[i] = i;
 	    Glu_persist->xsup[i] = i;
 	}
 	Glu_persist->xsup[n] = n;
+    } else {
+	/* User alreay allocated and defined the above supernode partition */
     }
 
     /* Count nonzeros per column for L & U */

SRC/prec-independent/trfAux.c

@@ -2579,9 +2579,9 @@ void permCol_SymbolicFact3d(superlu_dist_options_t *options, int n, SuperMatrix
 
     t = SuperLU_timer_();
 
-    if ( options->ILU_level != SLU_EMPTY ) {
+    if ( options->ILU_level != SLU_EMPTY ) { /* for any level-based ILU */
 	iinfo = ilu_level_symbfact(options, &GAC, perm_c, etree, Glu_persist, Glu_freeable);
-    } else {
+    } else { /* for complete LU */
 	iinfo = symbfact(options, iam, &GAC, perm_c, etree, Glu_persist, Glu_freeable);
     }