Initial commit

1 files changed, 257 insertions, 0 deletions

diff --git a/pkg/images/imutil/src/generic/imamul.x b/pkg/images/imutil/src/generic/imamul.x
new file mode 100644
index 00000000..05fdf8a4
--- /dev/null
+++ b/pkg/images/imutil/src/generic/imamul.x
@@ -0,0 +1,257 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+include	<imhdr.h>
+
+# IMA_MUL -- Image arithmetic multiplication.
+
+
+procedure ima_muls (im_a, im_b, im_c, a, b)
+
+pointer	im_a, im_b, im_c
+short	a, b
+
+int	len
+pointer	im[3], buf[3]
+long	v[IM_MAXDIM, 3]
+
+int	ima_nls()
+
+begin
+	# Loop through all of the image lines.
+	im[1] = im_c
+	len = IM_LEN (im[1], 1)
+	call amovkl (long(1), v, 3 * IM_MAXDIM)
+
+	# If imagea is constant then read imageb.  If the constant
+	# is 1 do a vector move to imagec otherwise do a vector
+	# multiply to imagec.
+	if (im_a == NULL) {
+	    im[2] = im_b
+	    while (ima_nls (im, buf, v, 2) != EOF) {
+		if (a == 1)
+		    call amovs (Mems[buf[2]], Mems[buf[1]], len)
+		else
+		    call amulks (Mems[buf[2]], a, Mems[buf[1]], len)
+	    }
+
+	# If imageb is constant then read imagea.  If the constant
+	# is 1 do a vector move to imagec otherwise do a vector
+	# multiply to imagec.
+	} else if (im_b == NULL) {
+	    im[2] = im_a
+	    while (ima_nls (im, buf, v, 2) != EOF) {
+		if (b == 1)
+		    call amovs (Mems[buf[2]], Mems[buf[1]], len)
+		else
+		    call amulks (Mems[buf[2]], b, Mems[buf[1]], len)
+	    }
+
+	# Read imagea and imageb and do the vector multiply to imagec.
+	} else {
+	    im[2] = im_a
+	    im[3] = im_b
+	    while (ima_nls (im, buf, v, 3) != EOF)
+		call amuls (Mems[buf[2]], Mems[buf[3]], Mems[buf[1]], len)
+	}
+end
+
+procedure ima_muli (im_a, im_b, im_c, a, b)
+
+pointer	im_a, im_b, im_c
+int	a, b
+
+int	len
+pointer	im[3], buf[3]
+long	v[IM_MAXDIM, 3]
+
+int	ima_nli()
+
+begin
+	# Loop through all of the image lines.
+	im[1] = im_c
+	len = IM_LEN (im[1], 1)
+	call amovkl (long(1), v, 3 * IM_MAXDIM)
+
+	# If imagea is constant then read imageb.  If the constant
+	# is 1 do a vector move to imagec otherwise do a vector
+	# multiply to imagec.
+	if (im_a == NULL) {
+	    im[2] = im_b
+	    while (ima_nli (im, buf, v, 2) != EOF) {
+		if (a == 1)
+		    call amovi (Memi[buf[2]], Memi[buf[1]], len)
+		else
+		    call amulki (Memi[buf[2]], a, Memi[buf[1]], len)
+	    }
+
+	# If imageb is constant then read imagea.  If the constant
+	# is 1 do a vector move to imagec otherwise do a vector
+	# multiply to imagec.
+	} else if (im_b == NULL) {
+	    im[2] = im_a
+	    while (ima_nli (im, buf, v, 2) != EOF) {
+		if (b == 1)
+		    call amovi (Memi[buf[2]], Memi[buf[1]], len)
+		else
+		    call amulki (Memi[buf[2]], b, Memi[buf[1]], len)
+	    }
+
+	# Read imagea and imageb and do the vector multiply to imagec.
+	} else {
+	    im[2] = im_a
+	    im[3] = im_b
+	    while (ima_nli (im, buf, v, 3) != EOF)
+		call amuli (Memi[buf[2]], Memi[buf[3]], Memi[buf[1]], len)
+	}
+end
+
+procedure ima_mull (im_a, im_b, im_c, a, b)
+
+pointer	im_a, im_b, im_c
+long	a, b
+
+int	len
+pointer	im[3], buf[3]
+long	v[IM_MAXDIM, 3]
+
+int	ima_nll()
+
+begin
+	# Loop through all of the image lines.
+	im[1] = im_c
+	len = IM_LEN (im[1], 1)
+	call amovkl (long(1), v, 3 * IM_MAXDIM)
+
+	# If imagea is constant then read imageb.  If the constant
+	# is 1 do a vector move to imagec otherwise do a vector
+	# multiply to imagec.
+	if (im_a == NULL) {
+	    im[2] = im_b
+	    while (ima_nll (im, buf, v, 2) != EOF) {
+		if (a == 1)
+		    call amovl (Meml[buf[2]], Meml[buf[1]], len)
+		else
+		    call amulkl (Meml[buf[2]], a, Meml[buf[1]], len)
+	    }
+
+	# If imageb is constant then read imagea.  If the constant
+	# is 1 do a vector move to imagec otherwise do a vector
+	# multiply to imagec.
+	} else if (im_b == NULL) {
+	    im[2] = im_a
+	    while (ima_nll (im, buf, v, 2) != EOF) {
+		if (b == 1)
+		    call amovl (Meml[buf[2]], Meml[buf[1]], len)
+		else
+		    call amulkl (Meml[buf[2]], b, Meml[buf[1]], len)
+	    }
+
+	# Read imagea and imageb and do the vector multiply to imagec.
+	} else {
+	    im[2] = im_a
+	    im[3] = im_b
+	    while (ima_nll (im, buf, v, 3) != EOF)
+		call amull (Meml[buf[2]], Meml[buf[3]], Meml[buf[1]], len)
+	}
+end
+
+procedure ima_mulr (im_a, im_b, im_c, a, b)
+
+pointer	im_a, im_b, im_c
+real	a, b
+
+int	len
+pointer	im[3], buf[3]
+long	v[IM_MAXDIM, 3]
+
+int	ima_nlr()
+
+begin
+	# Loop through all of the image lines.
+	im[1] = im_c
+	len = IM_LEN (im[1], 1)
+	call amovkl (long(1), v, 3 * IM_MAXDIM)
+
+	# If imagea is constant then read imageb.  If the constant
+	# is 1 do a vector move to imagec otherwise do a vector
+	# multiply to imagec.
+	if (im_a == NULL) {
+	    im[2] = im_b
+	    while (ima_nlr (im, buf, v, 2) != EOF) {
+		if (a == 1.0)
+		    call amovr (Memr[buf[2]], Memr[buf[1]], len)
+		else
+		    call amulkr (Memr[buf[2]], a, Memr[buf[1]], len)
+	    }
+
+	# If imageb is constant then read imagea.  If the constant
+	# is 1 do a vector move to imagec otherwise do a vector
+	# multiply to imagec.
+	} else if (im_b == NULL) {
+	    im[2] = im_a
+	    while (ima_nlr (im, buf, v, 2) != EOF) {
+		if (b == 1.0)
+		    call amovr (Memr[buf[2]], Memr[buf[1]], len)
+		else
+		    call amulkr (Memr[buf[2]], b, Memr[buf[1]], len)
+	    }
+
+	# Read imagea and imageb and do the vector multiply to imagec.
+	} else {
+	    im[2] = im_a
+	    im[3] = im_b
+	    while (ima_nlr (im, buf, v, 3) != EOF)
+		call amulr (Memr[buf[2]], Memr[buf[3]], Memr[buf[1]], len)
+	}
+end
+
+procedure ima_muld (im_a, im_b, im_c, a, b)
+
+pointer	im_a, im_b, im_c
+double	a, b
+
+int	len
+pointer	im[3], buf[3]
+long	v[IM_MAXDIM, 3]
+
+int	ima_nld()
+
+begin
+	# Loop through all of the image lines.
+	im[1] = im_c
+	len = IM_LEN (im[1], 1)
+	call amovkl (long(1), v, 3 * IM_MAXDIM)
+
+	# If imagea is constant then read imageb.  If the constant
+	# is 1 do a vector move to imagec otherwise do a vector
+	# multiply to imagec.
+	if (im_a == NULL) {
+	    im[2] = im_b
+	    while (ima_nld (im, buf, v, 2) != EOF) {
+		if (a == 1.0D0)
+		    call amovd (Memd[buf[2]], Memd[buf[1]], len)
+		else
+		    call amulkd (Memd[buf[2]], a, Memd[buf[1]], len)
+	    }
+
+	# If imageb is constant then read imagea.  If the constant
+	# is 1 do a vector move to imagec otherwise do a vector
+	# multiply to imagec.
+	} else if (im_b == NULL) {
+	    im[2] = im_a
+	    while (ima_nld (im, buf, v, 2) != EOF) {
+		if (b == 1.0D0)
+		    call amovd (Memd[buf[2]], Memd[buf[1]], len)
+		else
+		    call amulkd (Memd[buf[2]], b, Memd[buf[1]], len)
+	    }
+
+	# Read imagea and imageb and do the vector multiply to imagec.
+	} else {
+	    im[2] = im_a
+	    im[3] = im_b
+	    while (ima_nld (im, buf, v, 3) != EOF)
+		call amuld (Memd[buf[2]], Memd[buf[3]], Memd[buf[1]], len)
+	}
+end
+