| 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 | ]
|