Repatch (from linux) of OSX IRAF

1 files changed, 82 insertions, 0 deletions

diff --git a/unix/as.vax/amapr.s b/unix/as.vax/amapr.s
new file mode 100644
index 00000000..5ba41092
--- /dev/null
+++ b/unix/as.vax/amapr.s
@@ -0,0 +1,82 @@
+# AMAPR -- Linear transformation, type real.  The range of pixel values
+# A1 to A2 are mapped into the range B1 to B2 using a linear transformation.
+# Values less than A1 or greater than A2 are mapped into the values B1 and
+# B2 upon output.
+
+	.data	0
+
+	.set	A,	4
+	.set	B,	8
+	.set	NPIX,	12
+	.set	A1,	16
+	.set	A2,	20
+	.set	B1,	24
+	.set	B2,	28
+
+	.align	2
+	.globl	_amapr_
+	.text
+
+	# AMAPR (a, b, npix, a1, a2, b1, b2)
+	#
+	# scalar = real (b2 - b1) / real (a2 - a1)
+	# minout = min (b1, b2)
+	# maxout = max (b1, b2)
+	# 
+	# do i = 1, npix
+	#     b[i] = max(minout, min(maxout,
+	#	  PIXEL((a[i] - a1) * scalar) + b1))
+	#
+	# Registers:
+	#	r0	last_a
+	#	r1	a
+	#	r2	b
+	#	r3	scalar
+	#	r4	minout
+	#	r5	maxout
+	#	r6	a1
+	#	r7	b1
+
+_amapr_:
+	.word	01774			# save r2-r9
+	movl	A(ap), r1
+	movl	B(ap), r2
+	mull3	$4, *NPIX(ap), r0
+	addl2	r1, r0
+	movf	*A1(ap), r6
+	movf	*B1(ap), r7
+	movf	*A2(ap), r8
+	movf	*B2(ap), r9
+
+	subf3	r7, r9, r3		# r3 = (b2 - b1) / (a2 - a1)
+	subf3	r6, r8, r4
+	divf2	r4, r3
+
+	cmpf	r7, r9			# b1 <= b2
+	bleq	L1
+	movf	r9, r4			# no, min=b2, max=b1
+	movf	r7, r5
+	jbr	L2
+L1:	movf	r7, r4			# yes, min=b1, max=b2
+	movf	r9, r5
+L2:
+	subf3	r6, (r1)+, r8		# r8 = a[i] - a1
+	mulf2	r3, r8			# (..) * scalar
+	addf2	r7, r8			# (..) + b1
+
+	cmpf	r8, r4			# r8 < minout?
+	bgtr	L3
+	movf	r4, (r2)+
+	jbr	L5
+L3:
+	cmpf	r8, r5			# r8 > maxout?
+	blss	L4
+	movf	r5, (r2)+
+	jbr	L5
+L4:
+	movf	r8, (r2)+		# new value in range
+L5:
+	cmpl	r1, r0
+	blssu	L2			# loop again
+
+	ret