1 files changed, 35 insertions, 0 deletions
diff --git a/sys/plio/plcompare.x b/sys/plio/plcompare.x
new file mode 100644
index 00000000..59ba9d1f
--- /dev/null
+++ b/sys/plio/plcompare.x
@@ -0,0 +1,35 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+include	<plset.h>
+include	<plio.h>
+
+# PL_COMPARE -- Compare two masks for equality, optionally noting any
+# differences on the output file.  PL_EQUAL is returned if the two masks
+# are equivalent.
+
+int procedure pl_compare (pl_1, pl_2, outfd)
+
+pointer	pl_1, pl_2		#I masks to be compared
+int	outfd			#I file for diagnostic output, or NULL
+
+int	i
+bool	pll_equal()
+
+begin
+	if (PL_NAXES(pl_1) != PL_NAXES(pl_2) || PL_NLP(pl_1) != PL_NLP(pl_2)) {
+	    if (outfd != NULL)
+		call fprintf (outfd, "the masks are not the same size\n")
+	    return (PL_NOTEQUAL)
+	}
+
+	do i = 1, PL_NLP(pl_1)
+	    if (!pll_equal (LL(pl_1,PL_LP(pl_1,i)), LL(pl_2,PL_LP(pl_2,i)))) {
+		if (outfd != NULL) {
+		    call fprintf (outfd, "masks differ at line %d\n")
+			call pargi (i)
+		}
+		return (PL_NOTEQUAL)
+	    }
+
+	return (PL_EQUAL)
+end