VasCalc Documentation - An example: Difference between revisions

From Drorbn
Jump to navigationJump to search
No edit summary
 
 
(10 intermediate revisions by 2 users not shown)
Line 1: Line 1:
This is an example of how to use [[VasCalc]]. We check the third Reidemeister move against an almost-invariant (not a technical term).
This is an example of how to use [[VasCalc]]. We check the third Reidemeister move against an almost-invariant (not a technical term).


== Example ==
== The Reidemeister 3 Move ==


We want to use [[VasCalc]] to verify the third Reidemeister move. This is meant
[[Image:reid3.jpg]]
as a small example of how to use the VasCalc package.

The first couple of steps are to load up VasCalc.

{{In|n=1|in=<nowiki><<CDinterface.m </nowiki>}}

{{In|n=2|in=<nowiki>SetVasCalcPath["/home/zavosh/vc"]; </nowiki>}}

Now we need to load the definitions of <math>\Phi</math> and <math>R</math> as
defined in Dror Bar-Natan's paper on Non-Associative Tangles[http://www.math.toronto.edu/~drorbn/papers/nat/nat.pdf]:

{{InOut|n=3|in=<nowiki>Phi = ASeries[1 + (1/24)*CD[Line[1],
Line[2], Line[1, 2]] - (1/24)*CD[Line[2], Line[1], Line[1, 2]], 3, 0, 3]</nowiki>| out=<nowiki>ASeries[3, 0, {«JavaObject[vectorSpace.Coefficient]»,
«JavaObject[vectorSpace.Coefficient]», «JavaObject[ChordVector]»,
«JavaObject[vectorSpace.Coefficient]»}] </nowiki>}}


{{InOut|n=4|in=<nowiki>R = ASeries[1 + (1/2)CD[Line[1], Line[1]]
+ (1/8)CD[Line[1, 2], Line[1, 2]] + (1/48)CD[Line[1, 2, 3], Line[1, 2, 3]]
, 2, 0]</nowiki>|
out=<nowiki>ASeries[2, 0, {«JavaObject[vectorSpace.Coefficient]», «JavaObject[ChordVector]»,
«JavaObject[ChordVector]», «JavaObject[ChordVector]»}] </nowiki>}}

Now we define the definition of <math>Z(B(k),n)</math>. Where <math>Z</math> stands for our (pseudo) Vassiliev invariant and <math>Z(B(k),n)</math> is the invariant of a braid on <math>n</math> strands with a crossing on the strands <math>k</math> and <math>k+1</math>.


{{In|n=5|in=<nowiki>Clear[Z] </nowiki>}}

{{In|n=6|in=<nowiki>Z[B[k_],n_] /; k > 1 :=
Module[{ser, tbl},
ser = Nest[AddStrand[#, #[[1]]] &, Nest[DoubleStrand[#, 0] &, Phi, k - 2], n - k - 1].
Nest[AddStrand[#, #[[1]]] &, Nest[AddStrand[#, 0] &, R, k - 1], n - k - 1] .
Nest[AddStrand[#, #[[1]]] &, Nest[DoubleStrand[#, 0] &,
PermuteStrand[(Phi)^(-1), {{2, 3}}], k - 2], n - k - 1];
tbl = Table[i, {i, ser[[1]]}];
tbl = ReplacePart[tbl, k + 1, k];
tbl = ReplacePart[tbl, k, k + 1];
ASeries[tbl, ser[[3]]]]</nowiki>}}

Here we need to make an exception for the case <math>k=1</math> because the above will not work:

{{In|n=7|in=<nowiki>Z[B[1], n_] := Module[{tbl, ser},
ser = Nest[AddStrand[#, #[[1]]] &, R, n - 2];
tbl = Table[i, {i, ser[[1]]}];
tbl = ReplacePart[tbl, 2, 1];
tbl = ReplacePart[tbl, 1, 2];
ASeries[tbl, ser[[3]]]]</nowiki>}}

Now we are ready to verify the move. If <math>Z</math> is an invariant of the Reidemeister three move
we must have:

<math>
Z(B(k),n).Z(B(k+1),n).Z(B(k),n) = Z(B(k+1),n).Z(B(k),n).Z(B(k+1),n)
</math>

We calculate the right and left hand sides separately for <math>k=3</math> and <math>n=6</math>.

{{InOut|n=8|in=<nowiki>ll = Z[B[3], 6].Z[B[4], 6].Z[B[3], 6]</nowiki>|
out=<nowiki>ASeries[{1, 2, 5, 4, 3, 6}, {«JavaObject[vectorSpace.Coefficient]»,
«JavaObject[ChordVector]», «JavaObject[ChordVector]», «JavaObject[ChordVector]»}]</nowiki>}}

{{InOut|n=9|in=<nowiki>rr = Z[B[4], 6].Z[B[3], 6].Z[B[4], 6]</nowiki>|
out=<nowiki>ASeries[{1, 2, 5, 4, 3, 6}, {«JavaObject[vectorSpace.Coefficient]»,
«JavaObject[ChordVector]», «JavaObject[ChordVector]», «JavaObject[ChordVector]»}]</nowiki>}}


Now to compare:

{{InOut|n=10|in=<nowiki>CD[Reduce[ll - rr]]</nowiki>|out=0}}

Voila!

==References==
*<nowiki>[1]</nowiki> D. Bar-Natan, [http://www.math.toronto.edu/~drorbn/papers/nat/nat.pdf <i>Non-Associative Tangles</i>]. Geometric Topology (proceedings of the Georgia International Topology Conference, W. H. Kazez ed.), 139-183, Amer. Math. Soc. and International Press, Providence, 1997.

Latest revision as of 06:32, 16 August 2006

This is an example of how to use VasCalc. We check the third Reidemeister move against an almost-invariant (not a technical term).

The Reidemeister 3 Move

We want to use VasCalc to verify the third Reidemeister move. This is meant as a small example of how to use the VasCalc package.

The first couple of steps are to load up VasCalc.

In[1]:= <<CDinterface.m
In[2]:= SetVasCalcPath["/home/zavosh/vc"];

Now we need to load the definitions of and as defined in Dror Bar-Natan's paper on Non-Associative Tangles[1]:

In[3]:= Phi = ASeries[1 + (1/24)*CD[Line[1], Line[2], Line[1, 2]] - (1/24)*CD[Line[2], Line[1], Line[1, 2]], 3, 0, 3]
Out[3]= ASeries[3, 0, {«JavaObject[vectorSpace.Coefficient]», «JavaObject[vectorSpace.Coefficient]», «JavaObject[ChordVector]», «JavaObject[vectorSpace.Coefficient]»}]


In[4]:= R = ASeries[1 + (1/2)CD[Line[1], Line[1]] + (1/8)CD[Line[1, 2], Line[1, 2]] + (1/48)CD[Line[1, 2, 3], Line[1, 2, 3]] , 2, 0]
Out[4]= ASeries[2, 0, {«JavaObject[vectorSpace.Coefficient]», «JavaObject[ChordVector]», «JavaObject[ChordVector]», «JavaObject[ChordVector]»}]

Now we define the definition of . Where stands for our (pseudo) Vassiliev invariant and is the invariant of a braid on strands with a crossing on the strands and .


In[5]:= Clear[Z]
In[6]:= Z[B[k_],n_] /; k > 1 := Module[{ser, tbl}, ser = Nest[AddStrand[#, #[[1]]] &, Nest[DoubleStrand[#, 0] &, Phi, k - 2], n - k - 1]. Nest[AddStrand[#, #[[1]]] &, Nest[AddStrand[#, 0] &, R, k - 1], n - k - 1] . Nest[AddStrand[#, #[[1]]] &, Nest[DoubleStrand[#, 0] &, PermuteStrand[(Phi)^(-1), {{2, 3}}], k - 2], n - k - 1]; tbl = Table[i, {i, ser[[1]]}]; tbl = ReplacePart[tbl, k + 1, k]; tbl = ReplacePart[tbl, k, k + 1]; ASeries[tbl, ser[[3]]]]

Here we need to make an exception for the case because the above will not work:

In[7]:= Z[B[1], n_] := Module[{tbl, ser}, ser = Nest[AddStrand[#, #[[1]]] &, R, n - 2]; tbl = Table[i, {i, ser[[1]]}]; tbl = ReplacePart[tbl, 2, 1]; tbl = ReplacePart[tbl, 1, 2]; ASeries[tbl, ser[[3]]]]

Now we are ready to verify the move. If is an invariant of the Reidemeister three move we must have:

We calculate the right and left hand sides separately for and .

In[8]:= ll = Z[B[3], 6].Z[B[4], 6].Z[B[3], 6]
Out[8]= ASeries[{1, 2, 5, 4, 3, 6}, {«JavaObject[vectorSpace.Coefficient]», «JavaObject[ChordVector]», «JavaObject[ChordVector]», «JavaObject[ChordVector]»}]
In[9]:= rr = Z[B[4], 6].Z[B[3], 6].Z[B[4], 6]
Out[9]= ASeries[{1, 2, 5, 4, 3, 6}, {«JavaObject[vectorSpace.Coefficient]», «JavaObject[ChordVector]», «JavaObject[ChordVector]», «JavaObject[ChordVector]»}]


Now to compare:

In[10]:= CD[Reduce[ll - rr]]
Out[10]= 0

Voila!

References

  • [1] D. Bar-Natan, Non-Associative Tangles. Geometric Topology (proceedings of the Georgia International Topology Conference, W. H. Kazez ed.), 139-183, Amer. Math. Soc. and International Press, Providence, 1997.