-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathschur_product_purge_duplicates_demo_3.wxm
79 lines (60 loc) · 2.35 KB
/
schur_product_purge_duplicates_demo_3.wxm
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 12.09.0 ] */
/* [wxMaxima: input start ] */
file_search_maxima : cons(sconcat("/home/sciamp/Sources/tesi/Young_Tableaux/###.{lisp,mac,mc}"), file_search_maxima)$
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
load (young)$
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
E:new_ydiagram_safe ([3,1]);
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
U:new_ydiagram_safe ([2]);
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
ds:yschur_polynomial_product (E,U)$
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
ds:unique(ds);
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
e:reverse(E@rows)$
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
u:reverse (U@rows)$
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
skews:map (lambda ([x], fill_first_line (remove_tableau_column (generic_skew_tableau (x[1], e, [], 1),[]), u)), ds);
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
l1:map (lambda ([x], if (not emptyp (skews[x])) then
[ds[x][1],length (littlewood_richardson_num (d, e, u,[[skews[x][1],append([skews[x][2]],makelist(0,length (u)-1))]],1,2))]),
makelist (i,i,1,length (skews)))$
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
l2 : delete (false, map (lambda ([x], if (listp (x)) and (x[2] > 0) then x ), l1));
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
poly:map (lambda ([x], x[2]*fl_better_yschur_polynomial_rows (x[1],5)), l2)$
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
prod_lr:sum (poly[i],i,1,length(poly))$
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
polyE : fl_better_yschur_polynomial (E,5);
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
polyU : fl_better_yschur_polynomial (U,5);
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
prod : ratsimp (polyE*polyU);
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
prod_lr_s : ratsimp (prod_lr);
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
is (expand (prod - prod_lr = 0));
/* [wxMaxima: input end ] */
/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$