1 files changed, 76 insertions, 0 deletions
diff --git a/arraymap.scm b/arraymap.scm
new file mode 100644
index 0000000..18ee64a
--- /dev/null
+++ b/arraymap.scm
@@ -0,0 +1,76 @@
+;;;; "arraymap.scm", applicative routines for arrays in Scheme.  
+;;; Copyright (c) 1993 Aubrey Jaffer
+;
+;Permission to copy this software, to redistribute it, and to use it
+;for any purpose is granted, subject to the following restrictions and
+;understandings.
+;
+;1.  Any copy made of this software must include this copyright notice
+;in full.
+;
+;2.  I have made no warrantee or representation that the operation of
+;this software will be error-free, and I am under no obligation to
+;provide any services, by way of maintenance, update, or otherwise.
+;
+;3.  In conjunction with products arising from the use of this
+;material, there shall be no use of my name in any advertising,
+;promotional, or sales literature without prior written consent in
+;each case.
+
+(require 'array)
+
+(define (array-map! ra0 proc . ras)
+  (define (ramap rshape inds)
+    (if (null? (cdr rshape))
+	(do ((i (cadar rshape) (+ -1 i))
+	     (is (cons (cadar rshape) inds)
+		 (cons (+ -1 i) inds)))
+	    ((< i (caar rshape)))
+	  (apply array-set! ra0
+		 (apply proc (map (lambda (ra) (apply array-ref ra is))
+				  ras))
+		 is))
+	(let ((crshape (cdr rshape))
+	      (ll (caar rshape)))
+	  (do ((i (cadar rshape) (+ -1 i)))
+	      ((< i ll))
+	    (ramap crshape (cons i inds))))))
+  (ramap (reverse (array-shape ra0)) '()))
+
+(define (array-for-each proc . ras)
+  (define (rafe rshape inds)
+    (if (null? (cdr rshape))
+	(do ((i (caar rshape) (+ 1 i)))
+	    ((> i (cadar rshape)))
+	  (apply proc
+		 (map (lambda (ra)
+			(apply array-ref ra (reverse (cons i inds)))) ras)))
+	(let ((crshape (cdr rshape))
+	      (ll (cadar rshape)))
+	  (do ((i (caar rshape) (+ 1 i)))
+	      ((> i ll))
+	    (rafe crshape (cons i inds))))))
+  (rafe (array-shape (car ras)) '()))
+
+(define (shape->indexes shape)
+  (define ra0 (apply make-array '() shape))
+  (define (ramap rshape inds)
+    (if (null? (cdr rshape))
+	(do ((i (cadar rshape) (+ -1 i))
+	     (is (cons (cadar rshape) inds)
+		 (cons (+ -1 i) inds)))
+	    ((< i (caar rshape)))
+	  (apply array-set! ra0 is is))
+	(let ((crshape (cdr rshape))
+	      (ll (caar rshape)))
+	  (do ((i (cadar rshape) (+ -1 i)))
+	      ((< i ll))
+	    (ramap crshape (cons i inds))))))
+  (ramap (reverse shape) '())
+  ra0)
+
+(define (array-indexes ra)
+  (shape->indexes (array-shape ra)))
+
+(define (array-copy! source dest)
+  (array-map! dest identity source))