1 files changed, 25 insertions, 0 deletions
diff --git a/noao/astutil/pdm/pdmcompare.x b/noao/astutil/pdm/pdmcompare.x
new file mode 100644
index 00000000..b2d60978
--- /dev/null
+++ b/noao/astutil/pdm/pdmcompare.x
@@ -0,0 +1,25 @@
+include <mach.h>
+include "pdm.h"
+
+# Compare procedure for qsort.
+
+int procedure pdm_compare (item1, item2)
+
+int	item1		# index of first phase
+int	item2		# index of second phase
+
+pointer	comarray
+common /sortcom/ comarray
+double	p1, p2
+
+begin
+	p1 = Memd[comarray+item1-1]
+	p2 = Memd[comarray+item2-1]
+
+	if (p1 > p2)
+	    return (1)
+	if (p1 == p2)
+	    return (0)
+	if (p1 < p2)
+	    return (-1)
+end