Kronecker product: Difference between revisions

The Frink library [https://frinklang.org/fsp/colorize.fsp?f=Matrix.frink Matrix.frink] contains an implementation of Kronecker product. TheHowever, the following example demonstrates calculating the Kronecker product and typesetting the equations using Frink'smultidimensional arrays and no layoutexternal routineslibraries.

<lang frink>usea Matrix.frink= [[1,2],[3,4]]

ab = new Matrix[ [[10,25],[36,4] 7]]

c = new Matrix[ [[0,1,0],[1,1,1],[0,1,0] ]]

d = new Matrix[ [[1,1,1,1],[1,0,0,1],[1,1,1,1]] ]

formatProd[a,b] := formatTable[[[a.formatMatrix[a], "\u2297", b.formatMatrix[b], "=", a.formatMatrix[KroneckerProduct[a,b].formatMatrix[]]]]</lang>

┌ ┐

┌ ┐ ┌ ┐ │ 0 5 0 10│

│1┌ 2│ │0┐ 5│ ┌ ┐ │ 6 7 12 14│ │

│3│1 4│ ⊗2│ │6 7│ =│0 5│ │ 06 7 15 12 0 20│14│

└│ ┘ │ ⊗ └│ ┘ │ = │18│ 21 24 28│ │

│3 4│ │6 7│ │ 0 15 └0 ┘20│

└ ┘ └ ┘ │ ┌ ┐│

│18 21 │024 0 0 0 1 1 1 1 0 0 0 0│28│

└ │0 0 0 0 1 0 0 1 0 0 0 0│┘

┌ ┐ ┌ ┐ │0 0 0 0 1 1 1 1 0 0 0 0│ ┐

│0 1 0│ │1 1 1 1│ │1 1 1 1 1 1 1 │0 0 0 0 1 1 1 1 1│ 0 0 0 0│

│0 1 0│ │1 1 1 1│ │1 1 1 1 1 1 1 1 1 │0 0 0 0 1 0 0 1 1│ 0 0 0 0│

└ ┘ └ ┘ │0 0 0 0 1 1 1 1 0 0│ 0 0│ │

│0 0 0 0 1 0 01 1 1 0 0 0 0│

┌ ┐ ┌ ┐ │0 0│ 0 0 1 1 1 1 0 0 0 0│ │

│0 1 0│ │1 1 1 1│ │1 1 1 └1 1 1 1 1 1 1 1 ┘1│