factor/library/test/math/matrices.factor

IN: temporary
USING: kernel lists math matrices namespaces sequences test ;

[ [ [ 1 4 ] [ 2 5 ] [ 3 6 ] ] ]
[ M[ [ 1 4 ] [ 2 5 ] [ 3 6 ] ]M row-list ] unit-test

[
    M[ [ 0 ] [ 0 ] [ 0 ] ]M
] [
    3 1 <zero-matrix>
] unit-test

[
    M[ [ 1 ] [ 2 ] [ 3 ] ]M
] [
    { 1 2 3 } <col-vector>
] unit-test

[
    M[ [ 1 0 0 ]
       [ 0 1 0 ]
       [ 0 0 1 ] ]M
] [
    3 <identity-matrix>
] unit-test

[
    M[ [ 1 0 4 ]
       [ 0 7 0 ]
       [ 6 0 3 ] ]M
] [
    M[ [ 1 0 0 ]
       [ 0 2 0 ]
       [ 0 0 3 ] ]M
       
    M[ [ 0 0 4 ]
       [ 0 5 0 ]
       [ 6 0 0 ] ]M

    m+
] unit-test

[
    M[ [ 1 0 4 ]
       [ 0 7 0 ]
       [ 6 0 3 ] ]M
] [
    M[ [ 1 0 0 ]
       [ 0 2 0 ]
       [ 0 0 3 ] ]M
       
    M[ [ 0 0 -4 ]
       [ 0 -5 0 ]
       [ -6 0 0 ] ]M

    m-
] unit-test

[
    { 10 20 30 }
] [
    10 { 1 2 3 } n*v
] unit-test

[
    M[ [ 6 ] ]M
] [
    M[ [ 3 ] ]M M[ [ 2 ] ]M m.
] unit-test

[
    M[ [ 11 ] ]M
] [
    M[ [ 1 3 ] ]M M[ [ 5 ] [ 2 ] ]M m.
] unit-test

[
    [ [[ 0 0 ]] [[ 1 0 ]] ]
] [
    [ 2 1 [ 2dup cons , ] 2repeat ] make-list
] unit-test

[
    { 3 4 }
] [
    M[ [ 1 0 ]
       [ 0 1 ] ]M

    { 3 4 }

    m.v
] unit-test

[
    { 4 3 }
] [
    M[ [ 0 1 ]
       [ 1 0 ] ]M

    { 3 4 }

    m.v
] unit-test

[ { 0 0 1 } ] [ { 1 0 0 } { 0 1 0 } cross ] unit-test
[ { 1 0 0 } ] [ { 0 1 0 } { 0 0 1 } cross ] unit-test
[ { 0 1 0 } ] [ { 0 0 1 } { 1 0 0 } cross ] unit-test

[ M[ [ 1 2 ] [ 3 4 ] [ 5 6 ] ]M ]
[ M[ [ 1 2 ] [ 3 4 ] [ 5 6 ] ]M transpose transpose ]
unit-test

[ M[ [ 1 3 5 ] [ 2 4 6 ] ]M ]
[ M[ [ 1 3 5 ] [ 2 4 6 ] ]M transpose transpose ]
unit-test

[ M[ [ 1 3 5 ] [ 2 4 6 ] ]M ]
[ M[ [ 1 2 ] [ 3 4 ] [ 5 6 ] ]M transpose ]
unit-test

[
    M[ [ 28 ] ]M
] [
    M[ [ 2 4 6 ] ]M

    M[ [ 1 ]
       [ 2 ]
       [ 3 ] ]M
    
    m.
] unit-test

[
    [ [ 7 ] [ 4 8 ] [ 1 5 9 ] [ 2 6 ] [ 3 ] ]
] [
    M[ [ 1 2 3 ] [ 4 5 6 ] [ 7 8 9 ] ]M
    5 [ 2 - swap <diagonal> ] project-with [ >list ] map
] unit-test