Kronecker product: Difference between revisions
Content added Content deleted
m (→JS ES6) |
m (→{{header|REXX}}: added whitespace, added a couple of comments.) |
||
Line 2,721:
A little extra coding was added to make the matrix glyphs and elements alignment look nicer.
<lang rexx>/*REXX program calculates the Kronecker product of two arbitrary size matrices. */
w=
aMat= 2x2 1 2 3 4 /*define A matrix size and elements.*/
bMat= 2x2 0 5 6 7 /* " B " " " " */
Line 2,727:
call makeMat 'B', bMat /* " B " " " */
call KronMat 'Kronecker product' /*calculate the Kronecker product. */
w= 0;
aMat= 3x3 0 1 0 1 1 1 0 1 0 /*define A matrix size and elements.*/
bMat= 3x4 1 1 1 1 1 0 0 1 1 1 1 1 /* " B " " " " */
Line 2,736:
/*──────────────────────────────────────────────────────────────────────────────────────*/
KronMat: parse arg what; parse var @.a.shape aRows aCols
#= 0;
do rA=1 for aRows
do rB=1 for bRows; #= #+1; ##= 0;
do cA=1 for aCols; x= @.a.rA.cA
do cB=1 for bCols; y= @.b.rB.cB; ##= ##+1;
@.what.#.##=xy; w= max(w, length(xy) )
end /*cB*/
end /*cA*/
Line 2,750:
makeMat: parse arg what, size elements; arg , row 'X' col .; @.what.shape=row col
#=0; do r=1 for row /* [↓] bump item#; get item; max width*/
do c=1 for col;
@.what.r.c=_
end /*c*/ /* [↑] define an element of WHAT matrix*/
Line 2,756:
call showMat what, size; return
/*──────────────────────────────────────────────────────────────────────────────────────*/
showMat: parse arg what, size .; z= '┌'; parse var size row "X" col;
say; say $ copies('═',7) "matrix" what copies('═',7)
do r=1 for row; _= '│' /*start with long vertical bar*/
do c=1 for col; _=_ right(@.what.r.c, w); if r==1 then z=z left('',w)
end /*c*/
if r==1 then do; z=z '┐'; say $ $ z; end /*show the top part of matrix.*/
say $ $ _ '│' /*append a long vertical bar. */
end /*r*/
say $ $ translate(z, '└┘', "┌┐"); return /*show the bot part of matrix.*/</lang>
| |||