From 40e5a5811c6ffce9b0974e93cdd927cbcf60c157 Mon Sep 17 00:00:00 2001
From: Joe Hunkeler <jhunkeler@gmail.com>
Date: Tue, 11 Aug 2015 16:51:37 -0400
Subject: Repatch (from linux) of OSX IRAF

---
 unix/as.vax/aluir.s | 54 +++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 54 insertions(+)
 create mode 100644 unix/as.vax/aluir.s

(limited to 'unix/as.vax/aluir.s')

diff --git a/unix/as.vax/aluir.s b/unix/as.vax/aluir.s
new file mode 100644
index 00000000..30a37d0d
--- /dev/null
+++ b/unix/as.vax/aluir.s
@@ -0,0 +1,54 @@
+# ALUIR -- Lookup an interpolate a vector of type real onto a real grid.
+# [OBSOLETE - NO LONGER USED.  5/27/87]
+
+	.data	0
+	.align	2
+	.text
+	.globl	_aluir_
+
+	.set	A,	4
+	.set	B,	8
+	.set	X,	12
+	.set	NPIX,	16
+
+	# ALUIR (a, b, x, npix)
+	#
+	# left = int (x[i])
+	# tau = x[i] - left
+	# b[i] = (a[left] * (1-tau)) + (a[left+1] * tau)
+	#
+	# registers:
+	#	r0	max_b
+	#	r1	a
+	#	r2	b
+	#	r3	x
+	#	r4	x[i], tau
+	#	r5	left
+	#	r6
+
+_aluir_:
+	.word	0374			# save r2-r7
+	subl3	$4, A(ap), r1
+	movl	B(ap), r2
+	movl	X(ap), r3
+	mull3	$4, *NPIX(ap), r0
+	addl2	r2, r0
+L1:
+	movf	(r3)+, r4		# get X into r4
+	cvtfl	r4, r5			# r5 = left
+	cvtlf	r5, r6
+	subf2	r6, r4			# r4 = tau = (x[i] - left)
+
+	movf	(r1)[r5], r6
+	mulf3	r4, r6, r7
+	subf2	r7, r6			# r6 = (a[left] * (1-tau))
+
+	incl	r5
+	mulf3	r4, (r1)[r5], r7	# r7 = (a[left+1] * tau)
+
+	addf3	r6, r7, (r2)+		# output result to B
+
+	cmpl	r2, r0
+	blssu	L1
+
+	ret
-- 
cgit