1 | test_outer = ->
|
2 | run_test [
|
3 | "outer(a,b)",
|
4 | "a*b",
|
5 |
|
6 | "outer(a,[b1,b2])",
|
7 | "[a*b1,a*b2]",
|
8 |
|
9 | "outer([a1,a2],b)",
|
10 | "[a1*b,a2*b]",
|
11 |
|
12 | "H33=hilbert(3)",
|
13 | "",
|
14 |
|
15 | "H44=hilbert(4)",
|
16 | "",
|
17 |
|
18 | "H55=hilbert(5)",
|
19 | "",
|
20 |
|
21 | "H3344=outer(H33,H44)",
|
22 | "",
|
23 |
|
24 | "H4455=outer(H44,H55)",
|
25 | "",
|
26 |
|
27 | "H33444455=outer(H33,H44,H44,H55)",
|
28 | "",
|
29 |
|
30 | "simplify(inner(H3344,H4455)-contract(H33444455,4,5))",
|
31 | "0",
|
32 | ]
|