Threading a link by a polynomial: Difference between revisions
From Knot Atlas
Jump to navigationJump to search
Sam.panitch (talk | contribs) No edit summary |
Sam.panitch (talk | contribs) No edit summary |
||
| Line 2: | Line 2: | ||
As an example, we can verify some formulas from Mausbaum: |
As an example, we can verify some formulas from Mausbaum: |
||
{{Startup Note}} |
|||
<!--$$Import["http://katlas.org/w/index.php?title=CableLink.m&action=raw"];$$--> |
<!--$$Import["http://katlas.org/w/index.php?title=CableLink.m&action=raw"];$$--> |
||
<!--$$Import["http://katlas.org/w/index.php?title=CableLink.m&action=raw"];$$--> |
<!--$$Import["http://katlas.org/w/index.php?title=CableLink.m&action=raw"];$$--> |
||
<!--Robot Land, no human edits to "END"--> |
<!--Robot Land, no human edits to "END"--> |
||
{{In| |
{{In| |
||
n = |
n = 2 | |
||
in = <nowiki>Import["http://katlas.org/w/index.php?title=CableComponent.m&action=raw"];</nowiki>}} |
in = <nowiki>Import["http://katlas.org/w/index.php?title=CableComponent.m&action=raw"];</nowiki>}} |
||
<!--END--> |
<!--END--> |
||
| Line 19: | Line 20: | ||
<!--Robot Land, no human edits to "END"--> |
<!--Robot Land, no human edits to "END"--> |
||
{{In| |
{{In| |
||
n = |
n = 3 | |
||
in = <nowiki>hopfLink=PD[X[3,1,4,2],X[2,4,1,3]]; |
in = <nowiki>hopfLink=PD[X[3,1,4,2],X[2,4,1,3]]; |
||
bracket[n_]:=a^n-a^(-n); |
bracket[n_]:=a^n-a^(-n); |
||
| Line 37: | Line 38: | ||
<!--Robot Land, no human edits to "END"--> |
<!--Robot Land, no human edits to "END"--> |
||
{{InOut| |
{{InOut| |
||
n = |
n = 4 | |
||
in = <nowiki>Expand[CableLink[hopfLink, |
in = <nowiki>Expand[CableLink[hopfLink, |
||
R[Subscript[z, 1], 1]*cheb[2, Subscript[z, 2]], {1, 3}, {Subscript[ |
R[Subscript[z, 1], 1]*cheb[2, Subscript[z, 2]], {1, 3}, {Subscript[ |
||
| Line 51: | Line 52: | ||
$$--> |
$$--> |
||
{{InOut| |
{{InOut| |
||
n = |
n = 5 | |
||
in = <nowiki>Expand[CableLink[hopfLink, |
in = <nowiki>Expand[CableLink[hopfLink, |
||
R[Subscript[z, 1], 2]*cheb[4, Subscript[z, 2]], {1, 3}, {Subscript[ |
R[Subscript[z, 1], 2]*cheb[4, Subscript[z, 2]], {1, 3}, {Subscript[ |
||
Revision as of 18:27, 5 August 2025
CableLink[link,poly,strandList,vars], whose code is available here, computes the Kauffman bracket of link (given as a PD) with components L1,L2,...,Ln, cabled by the polynomial poly in the variables z1,z2,...,zn. strandList is a list of strand labels of length n, where the ith element is the first strand label corresponding to component Li.
As an example, we can verify some formulas from Mausbaum:
(For In[1] see Setup)
In[2]:=
|
Import["http://katlas.org/w/index.php?title=CableComponent.m&action=raw"];
|
In[3]:=
|
hopfLink=PD[X[3,1,4,2],X[2,4,1,3]];
bracket[n_]:=a^n-a^(-n);
bracketFact[n_]:=Product[bracket[i],{i,1,n}];
R[z_, n_] := Product[z + lambda[2*i], {i, 0, n - 1}];
cheb[0, z_] = 1;
cheb[1, z_] = z;
cheb[n_, z_] := cheb[n, z] = z*cheb[n - 1, z] - cheb[n - 2, z];
|
In[4]:=
|
Expand[CableLink[hopfLink,
R[Subscript[z, 1], 1]*cheb[2, Subscript[z, 2]], {1, 3}, {Subscript[
z, 1], Subscript[z, 2]}] /. {A -> a^(1/2)}]
Expand[(-1)^1*bracketFact[3]/bracket[1]]
|
Out[4]=
|
-1/a^5 + 1/a + a - a^5
-1/a^5 + 1/a + a - a^5
|
In[5]:=
|
Expand[CableLink[hopfLink,
R[Subscript[z, 1], 2]*cheb[4, Subscript[z, 2]], {1, 3}, {Subscript[
z, 1], Subscript[z, 2]}] /. {A -> a^(1/2)}]
Expand[(-1)^1*bracketFact[3]/bracket[1]]
|
Out[5]=
|
2 + 1/a^14 - 1/a^10 - 1/a^8 - 1/a^6 + 1/a^2 + a^2 - a^6 - a^8 - a^10 + a^14
2 + 1/a^14 - 1/a^10 - 1/a^8 - 1/a^6 + 1/a^2 + a^2 - a^6 - a^8 - a^10 + a^14
|
