Как правильно использовать умножение матриц в clojure?

Question

[void_phoenix]

(ns my-namespace
(:use clojure.core.matrix)
(:require [clojure.core.matrix.operators :as M]))
=> nil
(def a (matrix [[1 2 3] [4 5 6]]))
=> (var my-namespace/a)
(def b (matrix [[10 20] [20 30] [30 40]]))
=> (var my-namespace/b)
(M/* a b)
RuntimeException Incompatible shapes, cannot broadcast [3 2] to [2 3] clojure.core.matrix.impl.persistent-vector/eval5700/fn--5701 (persistent_vector.clj:156)

Не могу понять, в чем проблема. В теории две матрицы должны нормально перемножаться...

Answer 1

[Александр Мелеховец]

А вы правильное умножение используете?

Comments

[void_phoenix]

Ну, насколько я понял конструкцию :require :as, я назначил аллиас M на clojure.core.matrix.operators через
(:require [clojure.core.matrix.operators :as M])
то есть используется * именно из clojure.core.matrix.operators.
Если использовать простое умножение, получается другая ошибка, которую и можно ожидать при применении оператора, на рассчитанного на матрицы.

[Александр Мелеховец]

M/* = clojure.core.matrix.mul - поэлементное умножение
clojure.core.matrix.mmul - матричное умножение

[void_phoenix]

Спасибо!

Answer 2

[void_phoenix]

Вообще все выглядит так, как будто переопределили оператор умножения с математического, на простое перемножение элементов матриц, как в сложении. Только я не очень понимаю, зачем это было сделано.