Написать функцию, которая найдет произведение двух матриц, заданных в виде массивов. Предусмотреть проверку корректности

(defun mx-product (w v
                   &aux
                   (a (array-dimension w 0))
                   (b (array-dimension w 1))
                   (c (array-dimension v 1))
                   (z (make-array `(,a ,c) :initial-element 0)))
  (when (= b (array-dimension v 0))
    (loop for i below a do
          (loop for k below c do
                (setf (aref z i k)
                      (loop for j below b
                            sum (* (aref w i j)
                                   (aref v j k)))))
          finally (return z))))
 
> (mx-product #2a((1 2) (3 4)) #2a((1 2 3) (4 5 6)))
#2A((9 12 15) (19 26 33))
> (mx-product #2a((1 2) (3 4) (5 6)) #2a((1 2 3) (4 5 6)))
#2A((9 12 15) (19 26 33) (29 40 51))
> (mx-product #2a((1 2) (3 4) (5 6)) #2a((1 2 3) (4 5 6) (7 8 9)))
NIL

(defun mx-product (w v
                   &aux
                   (a (array-dimension w 0))
                   (b (array-dimension w 1))
                   (c (array-dimension v 1))
                   (z (make-array `(,a ,c) :initial-element 0)))
  (when (= b (array-dimension v 0))
    (loop for i below a do
          (loop for k below c do
                (setf (aref z i k)
                      (loop for j below b
                            sum (* (aref w i j)
                                   (aref v j k)))))
          finally (return z))))
 
> (mx-product #2a((1 2) (3 4)) #2a((1 2 3) (4 5 6)))
#2A((9 12 15) (19 26 33))
> (mx-product #2a((1 2) (3 4) (5 6)) #2a((1 2 3) (4 5 6)))
#2A((9 12 15) (19 26 33) (29 40 51))
> (mx-product #2a((1 2) (3 4) (5 6)) #2a((1 2 3) (4 5 6) (7 8 9)))
NIL

(defun matr-prod (matr1 matr2)
  (let ((matr2* (apply 'mapcar (cons 'list matr2))))
    (mapcar (lambda (r1) (mapcar (lambda (r2) (apply '+ (mapcar '* r1 r2))) matr2*)) matr1)))
 
(matr-prod '((1 2) (3 4)) '((5 6) (7 8)))
 
==> ((19 22) (43 50))

(defun mx-product (w v)
  (when (= (length (car w)) (length v))
    (mapcar #'(lambda (r)
                (apply #'mapcar #'(lambda (&rest c)
                                    (reduce #'+ (mapcar #'* r c)))
                       v))
            w)))
 
> (mx-product '((1 2) (3 4)) '((1 2 3) (4 5 6)))
((9 12 15) (19 26 33))
> (mx-product '((1 2) (3 4) (5 6)) '((1 2 3) (4 5 6)))
((9 12 15) (19 26 33) (29 40 51))
> (mx-product '((1 2) (3 4) (5 6)) '((1 2 3) (4 5 6) (7 8 9)))
NIL

(assert (progn
          (check-type a (array number (* *)))
          (check-type b (array number (* *)))
          (= (array-dimension a 1) (array-dimension b 0)))
        (a b)
        "Cannot multiply ~A by ~A." a b)