Skip to content
Permalink
Branch: master
Find file Copy path
Find file Copy path
Fetching contributors…
Cannot retrieve contributors at this time
395 lines (338 sloc) 6 KB
USING: math.matrices math.vectors tools.test math kernel ;
IN: math.matrices.tests
{
{ { 0 } { 0 } { 0 } }
} [
3 1 zero-matrix
] unit-test
{
{ { 1 0 0 }
{ 0 1 0 }
{ 0 0 1 } }
} [
3 identity-matrix
] unit-test
{
{ { 1 0 0 }
{ 0 2 0 }
{ 0 0 3 } }
} [
{ 1 2 3 } diagonal-matrix
] unit-test
{
{ { 1 1 1 }
{ 4 2 1 }
{ 9 3 1 }
{ 25 5 1 } }
} [
{ 1 2 3 5 } 3 vandermonde-matrix
] unit-test
{
{
{ 1 0 0 }
{ 0 1 0 }
{ 0 0 1 }
}
} [
3 3 0 eye
] unit-test
{
{
{ 0 1 0 }
{ 0 0 1 }
{ 0 0 0 }
}
} [
3 3 1 eye
] unit-test
{
{
{ 0 0 0 }
{ 1 0 0 }
{ 0 1 0 }
}
} [
3 3 -1 eye
] unit-test
{
{
{ 1 0 0 0 }
{ 0 1 0 0 }
{ 0 0 1 0 }
}
} [
3 4 0 eye
] unit-test
{
{
{ 0 1 0 }
{ 0 0 1 }
{ 0 0 0 }
{ 0 0 0 }
}
} [
4 3 1 eye
] unit-test
{
{
{ 0 0 0 }
{ 1 0 0 }
{ 0 1 0 }
{ 0 0 1 }
}
} [
4 3 -1 eye
] unit-test
{
{ { 1 1/2 1/3 1/4 }
{ 1/2 1/3 1/4 1/5 }
{ 1/3 1/4 1/5 1/6 }
}
} [ 3 4 hilbert-matrix ] unit-test
{
{ { 1 2 3 4 }
{ 2 1 2 3 }
{ 3 2 1 2 }
{ 4 3 2 1 } }
} [ 4 toeplitz-matrix ] unit-test
{
{ { 1 2 3 4 }
{ 2 3 4 0 }
{ 3 4 0 0 }
{ 4 0 0 0 } }
} [ 4 hankel-matrix ] unit-test
{
{ { 1 0 4 }
{ 0 7 0 }
{ 6 0 3 } }
} [
{ { 1 0 0 }
{ 0 2 0 }
{ 0 0 3 } }
{ { 0 0 4 }
{ 0 5 0 }
{ 6 0 0 } }
m+
] unit-test
{
{ { 1 0 4 }
{ 0 7 0 }
{ 6 0 3 } }
} [
{ { 1 0 0 }
{ 0 2 0 }
{ 0 0 3 } }
{ { 0 0 -4 }
{ 0 -5 0 }
{ -6 0 0 } }
m-
] unit-test
{
{ 10 20 30 }
} [
10 { 1 2 3 } n*v
] unit-test
{
{ 3 4 }
} [
{ { 1 0 }
{ 0 1 } }
{ 3 4 }
m.v
] unit-test
{
{ 4 3 }
} [
{ { 0 1 }
{ 1 0 } }
{ 3 4 }
m.v
] unit-test
{
{ { 6 } }
} [
{ { 3 } } { { 2 } } m.
] unit-test
{
{ { 11 } }
} [
{ { 1 3 } } { { 5 } { 2 } } m.
] unit-test
{
{ { 28 } }
} [
{ { 2 4 6 } }
{ { 1 }
{ 2 }
{ 3 } }
m.
] 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
{ { 0.0 -0.707 0.707 } } [ { 1.0 0.0 0.0 } { 0.0 0.707 0.707 } cross ] unit-test
{ { 0 -2 2 } } [ { -1 -1 -1 } { 1 -1 -1 } cross ] unit-test
{ { 1 0 0 } } [ { 1 1 0 } { 1 0 0 } proj ] unit-test
{ { { 4181 6765 } { 6765 10946 } } }
[ { { 0 1 } { 1 1 } } 20 m^n ] unit-test
[ { { 0 1 } { 1 1 } } -20 m^n ] [ negative-power-matrix? ] must-fail-with
{
{ { 0 5 0 10 } { 6 7 12 14 } { 0 15 0 20 } { 18 21 24 28 } }
}
[ { { 1 2 } { 3 4 } } { { 0 5 } { 6 7 } } kron ] unit-test
{
{
{ 1 1 1 1 }
{ 1 -1 1 -1 }
{ 1 1 -1 -1 }
{ 1 -1 -1 1 }
}
} [ { { 1 1 } { 1 -1 } } dup kron ] unit-test
{
{
{ 1 1 1 1 1 1 1 1 }
{ 1 -1 1 -1 1 -1 1 -1 }
{ 1 1 -1 -1 1 1 -1 -1 }
{ 1 -1 -1 1 1 -1 -1 1 }
{ 1 1 1 1 -1 -1 -1 -1 }
{ 1 -1 1 -1 -1 1 -1 1 }
{ 1 1 -1 -1 -1 -1 1 1 }
{ 1 -1 -1 1 -1 1 1 -1 }
}
} [ { { 1 1 } { 1 -1 } } dup dup kron kron ] unit-test
{
{
{ 1 1 1 1 1 1 1 1 }
{ 1 -1 1 -1 1 -1 1 -1 }
{ 1 1 -1 -1 1 1 -1 -1 }
{ 1 -1 -1 1 1 -1 -1 1 }
{ 1 1 1 1 -1 -1 -1 -1 }
{ 1 -1 1 -1 -1 1 -1 1 }
{ 1 1 -1 -1 -1 -1 1 1 }
{ 1 -1 -1 1 -1 1 1 -1 }
}
} [ { { 1 1 } { 1 -1 } } dup dup kron swap kron ] unit-test
! kron is not generally commutative, make sure we have the right order
{
{
{ 1 2 3 4 5 1 2 3 4 5 }
{ 6 7 8 9 10 6 7 8 9 10 }
{ 1 2 3 4 5 -1 -2 -3 -4 -5 }
{ 6 7 8 9 10 -6 -7 -8 -9 -10 }
}
}
[
{ { 1 1 } { 1 -1 } }
{ { 1 2 3 4 5 } { 6 7 8 9 10 } } kron
] unit-test
{
{
{ 1 1 2 2 3 3 4 4 5 5 }
{ 1 -1 2 -2 3 -3 4 -4 5 -5 }
{ 6 6 7 7 8 8 9 9 10 10 }
{ 6 -6 7 -7 8 -8 9 -9 10 -10 }
}
}
[
{ { 1 1 } { 1 -1 } }
{ { 1 2 3 4 5 } { 6 7 8 9 10 } } swap kron
] unit-test
{
{ { 5 10 15 }
{ 6 12 18 }
{ 7 14 21 } }
} [ { 5 6 7 } { 1 2 3 } outer ] unit-test
CONSTANT: test-points {
{ 80 27 89 } { 80 27 88 } { 75 25 90 }
{ 62 24 87 } { 62 22 87 } { 62 23 87 }
{ 62 24 93 } { 62 24 93 } { 58 23 87 }
{ 58 18 80 } { 58 18 89 } { 58 17 88 }
{ 58 18 82 } { 58 19 93 } { 50 18 89 }
{ 50 18 86 } { 50 19 72 } { 50 19 79 }
{ 50 20 80 } { 56 20 82 } { 70 20 91 }
}
{
{
{ 84+2/35 22+23/35 24+4/7 }
{ 22+23/35 9+104/105 6+87/140 }
{ 24+4/7 6+87/140 28+5/7 }
}
} [
test-points sample-cov-matrix
] unit-test
{
{
{ 80+8/147 21+85/147 23+59/147 }
{ 21+85/147 9+227/441 6+15/49 }
{ 23+59/147 6+15/49 27+17/49 }
}
} [
test-points population-cov-matrix
] unit-test
{
{
{ 5 5 }
{ 5 5 }
}
} [
2 2 5 <matrix>
] unit-test
{
{
{ 5 5 }
{ 5 5 }
}
} [
2 2 [ 5 ] make-matrix
] unit-test
{
{
{ 0 1 2 }
{ 1 2 3 }
}
} [
2 3 [ + ] make-matrix-with-indices
] unit-test
{
{
{ 0 1 }
{ 0 1 }
}
} [
2 square-rows
] unit-test
{
{
{ 0 0 }
{ 1 1 }
}
} [
2 square-cols
] unit-test
{
{
{ 5 6 }
{ 5 6 }
}
} [
{ 5 6 } square-rows
] unit-test
{
{
{ 5 5 }
{ 6 6 }
}
} [
{ 5 6 } square-cols
] unit-test
{ t } [ { } square-matrix? ] unit-test
{ t } [ { { 1 } } square-matrix? ] unit-test
{ t } [ { { 1 2 } { 3 4 } } square-matrix? ] unit-test
{ f } [ { { 1 } { 2 3 } } square-matrix? ] unit-test
{ f } [ { { 1 2 } } square-matrix? ] unit-test
{ 9 }
[ { { 2 -2 1 } { 1 3 -1 } { 2 -4 2 } } m-1norm ] unit-test
{ 8 }
[ { { 2 -2 1 } { 1 3 -1 } { 2 -4 2 } } m-infinity-norm ] unit-test
{ 2.0 }
[ { { 1 1 } { 1 1 } } frobenius-norm ] unit-test
You can’t perform that action at this time.