provide tensor function, just creating arrays

decons.lisp

@@ -3,9 +3,9 @@
 (defpackage :decons
   (:use :common-lisp)
   (:local-nicknames (:util :scopes/util))
-  (:export #:*pi*
-           #:area #:circle
+  (:export #:*pi* #:area #:circle
            #:absv #:double #:remainder
+           #:scalar-p #:tensor #:at
            #:line
            ))
 
@@ -34,6 +34,22 @@
     v
     (remainder (- v d) d)))
 
+;;;; library stuff
+
+(defgeneric scalar-p (x)
+  (:method (x) t)
+  (:method ((x list)) nil)
+  (:method ((x array)) nil))
+
+(defun tensor (s v)
+  (make-array s :initial-contents v))
+
+(defun at (a &rest subs)
+  (apply #'aref a subs))
+
+(defun (setf at) (v a &rest subs)
+  (setf (apply #'aref a subs) v))
+
 ;;;; parameterized functions
 
 (defun line (x)

test-decons.lisp

@@ -11,6 +11,7 @@
 (defun run ()
   (let ((t:*test-suite* (t:test-suite "decons")))
     (test-basic)
+    (test-lib)
     (test-line)
     (t:show-result)))
 
@@ -24,6 +25,17 @@
   (== (decons:remainder 7 4) 3)
 )
 
+(deftest test-lib ()
+  (== (decons:scalar-p 7) t)
+  (== (decons:scalar-p '(a b)) nil)
+  (let ((t1 (decons:tensor '(4) '(0 1 2 3))))
+    (== (decons:at t1 0) 0)
+    (setf (decons:at t1 0) 5)
+    (== (decons:at t1 0) 5)
+    (== (decons:scalar-p t1) nil)
+    )
+)
+
 (deftest test-line ()
   (let ((p1 (decons:line 0.0))
         (p2 (decons:line 1.0))