--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/liboctave/UMFPACK/UMFPACK/Source/umf_report_perm.c	Fri Feb 25 19:55:28 2005 +0000
@@ -0,0 +1,85 @@
+/* ========================================================================== */
+/* === UMF_report_perm ====================================================== */
+/* ========================================================================== */
+
+/* -------------------------------------------------------------------------- */
+/* UMFPACK Version 4.4, Copyright (c) 2005 by Timothy A. Davis.  CISE Dept,   */
+/* Univ. of Florida.  All Rights Reserved.  See ../Doc/License for License.   */
+/* web: http://www.cise.ufl.edu/research/sparse/umfpack                       */
+/* -------------------------------------------------------------------------- */
+
+#include "umf_internal.h"
+
+#define PRINTF4U(params) { if (user || prl >= 4) PRINTF (params) ; }
+
+GLOBAL Int UMF_report_perm
+(
+    Int n,
+    const Int P [ ],
+    Int W [ ],		/* workspace of size n */
+    Int prl,
+    Int user
+)
+{
+    Int i, k, valid, prl1 ;
+
+    ASSERT (prl >= 3) ;
+
+    PRINTF4U (("permutation vector, n = "ID". ", n)) ;
+
+    if (n <= 0)
+    {
+	PRINTF (("ERROR: length of permutation is <= 0\n\n")) ;
+	return (UMFPACK_ERROR_n_nonpositive) ;
+    }
+
+    if (!P)
+    {
+	/* if P is (Int *) NULL, this is the identity permutation */
+	PRINTF (("(not present)\n\n")) ;
+	return (UMFPACK_OK) ;
+    }
+
+    if (!W)
+    {
+	PRINTF (("ERROR: out of memory\n\n")) ;
+	return (UMFPACK_ERROR_out_of_memory) ;
+    }
+
+    PRINTF4 (("\n")) ;
+
+    for (i = 0 ; i < n ; i++)
+    {
+	W [i] = TRUE ;
+    }
+
+    prl1 = prl ;
+    for (k = 0 ; k < n ; k++)
+    {
+	i = P [k] ;
+	PRINTF4 (("    "ID" : "ID" ", INDEX (k), INDEX (i))) ;
+	valid = (i >= 0 && i < n) ;
+	if (valid)
+	{
+	    valid = W [i] ;
+	    W [i] = FALSE ;
+	}
+	if (!valid)
+	{
+	    /* out of range or duplicate entry */
+	    PRINTF (("ERROR: invalid\n\n")) ;
+	    return (UMFPACK_ERROR_invalid_permutation) ;
+	}
+	PRINTF4 (("\n")) ;
+	if (prl == 4 && k == 9 && n > 10)
+	{
+	    PRINTF (("    ...\n")) ;
+	    prl-- ;
+	}
+    }
+    prl = prl1 ;
+
+    PRINTF4 (("    permutation vector ")) ;
+    PRINTF4U (("OK\n\n")) ;
+    return (UMFPACK_OK) ;
+}