login  home  contents  what's new  discussion  bug reports     help  links  subscribe  changes  refresh  edit

Edit detail for Snake Relation revision 11 of 11

1 2 3 4 5 6 7 8 9 10 11
Editor: Bill Page
Time: 2013/04/23 21:47:01 GMT+0
Note: missing CARTEN library

changed:
-)library MONAL PROP LOP
)library CARTEN MONAL PROP LOP

Non-degeneracy of the pairing

Ref:


\scalebox{1} % Change this value to rescale the drawing.
{
\begin{pspicture}(0,-0.82)(5.22,0.82)
\psbezier[linewidth=0.04](0.0,0.0)(0.0,-0.8)(0.6,-0.8)(0.6,0.0)
\psbezier[linewidth=0.04](0.6,0.0)(0.6,0.8)(1.2,0.8)(1.2,0.0)
\psline[linewidth=0.04cm](0.0,0.0)(0.0,0.6)
\psline[linewidth=0.04cm](1.2,0.0)(1.2,-0.6)
\psline[linewidth=0.04cm](2.6,0.6)(2.6,-0.6)
\psbezier[linewidth=0.04](5.2,0.0)(5.2,-0.8)(4.6,-0.8)(4.6,0.0)
\psbezier[linewidth=0.04](4.6,0.0)(4.6,0.8)(4.0,0.8)(4.0,0.0)
\psline[linewidth=0.04cm](5.2,0.0)(5.2,0.6)
\psline[linewidth=0.04cm](4.0,0.0)(4.0,-0.6)
\usefont{T1}{ptm}{m}{n}
\rput(1.8948437,0.105){=}
\usefont{T1}{ptm}{m}{n}
\rput(3.2948437,0.105){=}
\end{pspicture} 
}
 

We use the Axiom LinearOperator library

fricas
(1) -> )library CARTEN MONAL PROP LOP
CartesianTensor is now explicitly exposed in frame initial CartesianTensor will be automatically loaded when needed from /var/aw/var/LatexWiki/CARTEN.NRLIB/CARTEN Monoidal is now explicitly exposed in frame initial Monoidal will be automatically loaded when needed from /var/aw/var/LatexWiki/MONAL.NRLIB/MONAL Prop is now explicitly exposed in frame initial Prop will be automatically loaded when needed from /var/aw/var/LatexWiki/PROP.NRLIB/PROP LinearOperator is now explicitly exposed in frame initial LinearOperator will be automatically loaded when needed from /var/aw/var/LatexWiki/LOP.NRLIB/LOP

and convenient notation

fricas
macro Σ(x,i,n)==reduce(+,[x for i in n])
Type: Void
fricas
macro Ξ(f,i,n)==[f for i in n]
Type: Void
fricas
macro sb == subscript
Type: Void
fricas
macro sp == superscript
Type: Void

Let 𝐋 be the domain of 2-dimensional linear operators

fricas
dim:=2

\label{eq1}2(1)
Type: PositiveInteger?
fricas
macro ℒ == List
Type: Void
fricas
macro ℚ == Expression Integer
Type: Void
fricas
𝐋 := LinearOperator(OVAR ['1, '2], ℚ)

\label{eq2}\hbox{\axiomType{LinearOperator}\ } \left({{\hbox{\axiomType{OrderedVariableList}\ } \left({\left[ 1, \: 2 \right]}\right)}, \:{\hbox{\axiomType{Expression}\ } \left({\hbox{\axiomType{Integer}\ }}\right)}}\right)(2)
Type: Type
fricas
𝐞:ℒ 𝐋      := basisOut()

\label{eq3}\left[{|_{1}}, \:{|_{2}}\right](3)
Type: List(LinearOperator(OrderedVariableList([1,2]),Expression(Integer)))
fricas
𝐝:ℒ 𝐋      := basisIn()

\label{eq4}\left[{|_{\ }^{1}}, \:{|_{\ }^{2}}\right](4)
Type: List(LinearOperator(OrderedVariableList([1,2]),Expression(Integer)))
fricas
I:𝐋:=[1]   -- identity for composition

\label{eq5}{|_{1}^{1}}+{|_{2}^{2}}(5)
Type: LinearOperator(OrderedVariableList([1,2]),Expression(Integer))
fricas
X:𝐋:=[2,1] -- twist

\label{eq6}{|_{1 \  1}^{1 \  1}}+{|_{2 \  1}^{1 \  2}}+{|_{1 \  2}^{2 \  1}}+{|_{2 \  2}^{2 \  2}}(6)
Type: LinearOperator(OrderedVariableList([1,2]),Expression(Integer))

Pairing

A scalar product (pairing) is represented by

fricas
U:=Σ(Σ(sp('u,[i,j])*𝐝.i*𝐝.j, i,1..dim), j,1..dim)

\label{eq7}{{u^{1, \: 1}}\ {|_{\ }^{1 \  1}}}+{{u^{1, \: 2}}\ {|_{\ }^{1 \  2}}}+{{u^{2, \: 1}}\ {|_{\ }^{2 \  1}}}+{{u^{2, \: 2}}\ {|_{\ }^{2 \  2}}}(7)
Type: LinearOperator(OrderedVariableList([1,2]),Expression(Integer))

In general we do not require that it be symmetric.

Co-pairing

Solve the "snake relation" as a system of linear equations.

fricas
Ω:𝐋:=Σ(Σ(sb('u,[i,j])*𝐞.i*𝐞.j, i,1..dim), j,1..dim)

\label{eq8}{{u_{1, \: 1}}\ {|_{1 \  1}}}+{{u_{1, \: 2}}\ {|_{1 \  2}}}+{{u_{2, \: 1}}\ {|_{2 \  1}}}+{{u_{2, \: 2}}\ {|_{2 \  2}}}(8)
Type: LinearOperator(OrderedVariableList([1,2]),Expression(Integer))
fricas
Í:=(I*Ω)/(U*I);
Type: LinearOperator(OrderedVariableList([1,2]),Expression(Integer))
fricas
Ì:=(Ω*I)/(I*U);
Type: LinearOperator(OrderedVariableList([1,2]),Expression(Integer))
fricas
equate(f,g)==map((x,y)+->(x=y),ravel f, ravel g);
Type: Void
fricas
eq1:=equate(Í,I)
fricas
Compiling function equate with type (LinearOperator(
      OrderedVariableList([1,2]),Expression(Integer)), LinearOperator(
      OrderedVariableList([1,2]),Expression(Integer))) -> List(Equation
      (Expression(Integer)))

\label{eq9}\begin{array}{@{}l}
\displaystyle
\left[{{{{u^{1, \: 2}}\ {u_{2, \: 1}}}+{{u^{1, \: 1}}\ {u_{1, \: 1}}}}= 1}, \:{{{{u^{1, \: 2}}\ {u_{2, \: 2}}}+{{u^{1, \: 1}}\ {u_{1, \: 2}}}}= 0}, \: \right.
\
\
\displaystyle
\left.{{{{u^{2, \: 2}}\ {u_{2, \: 1}}}+{{u^{2, \: 1}}\ {u_{1, \: 1}}}}= 0}, \:{{{{u^{2, \: 2}}\ {u_{2, \: 2}}}+{{u^{2, \: 1}}\ {u_{1, \: 2}}}}= 1}\right] 
(9)
Type: List(Equation(Expression(Integer)))
fricas
eq2:=equate(Ì,I)

\label{eq10}\begin{array}{@{}l}
\displaystyle
\left[{{{{u^{2, \: 1}}\ {u_{1, \: 2}}}+{{u^{1, \: 1}}\ {u_{1, \: 1}}}}= 1}, \:{{{{u^{2, \: 1}}\ {u_{2, \: 2}}}+{{u^{1, \: 1}}\ {u_{2, \: 1}}}}= 0}, \: \right.
\
\
\displaystyle
\left.{{{{u^{2, \: 2}}\ {u_{1, \: 2}}}+{{u^{1, \: 2}}\ {u_{1, \: 1}}}}= 0}, \:{{{{u^{2, \: 2}}\ {u_{2, \: 2}}}+{{u^{1, \: 2}}\ {u_{2, \: 1}}}}= 1}\right] 
(10)
Type: List(Equation(Expression(Integer)))
fricas
snake:=solve(concat(eq1,eq2),concat Ξ(Ξ(sb('u,[i,j]), i,1..dim), j,1..dim));
Type: List(List(Equation(Expression(Integer))))
fricas
if #snake ~= 1 then error "no solution"
Type: Void
fricas
Ω:=eval(Ω,snake(1))

\label{eq11}\begin{array}{@{}l}
\displaystyle
{{\frac{u^{2, \: 2}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}\ {|_{1 \  1}}}- 
\
\
\displaystyle
{{\frac{u^{1, \: 2}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}\ {|_{1 \  2}}}- 
\
\
\displaystyle
{{\frac{u^{2, \: 1}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}\ {|_{2 \  1}}}+ 
\
\
\displaystyle
{{\frac{u^{1, \: 1}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}\ {|_{2 \  2}}}
(11)
Type: LinearOperator(OrderedVariableList([1,2]),Expression(Integer))
fricas
matrix Ξ(Ξ(Ω/(𝐝.i*𝐝.j), i,1..dim), j,1..dim)

\label{eq12}\left[ 
\begin{array}{cc}
{\frac{u^{2, \: 2}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}& -{\frac{u^{2, \: 1}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}
\
-{\frac{u^{1, \: 2}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}&{\frac{u^{1, \: 1}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}
(12)
Type: Matrix(LinearOperator(OrderedVariableList([1,2]),Expression(Integer)))

This is equivalent to a matrix inverse (transposed!)

fricas
Um:=matrix Ξ(Ξ((𝐞.i*𝐞.j)/U, i,1..dim), j,1..dim)

\label{eq13}\left[ 
\begin{array}{cc}
{u^{1, \: 1}}&{u^{2, \: 1}}
\
{u^{1, \: 2}}&{u^{2, \: 2}}
(13)
Type: Matrix(LinearOperator(OrderedVariableList([1,2]),Expression(Integer)))
fricas
mU:=transpose inverse map(retract,Um)

\label{eq14}\left[ 
\begin{array}{cc}
{\frac{u^{2, \: 2}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}& -{\frac{u^{1, \: 2}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}
\
-{\frac{u^{2, \: 1}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}&{\frac{u^{1, \: 1}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}
(14)
Type: Matrix(Expression(Integer))
fricas
Ωm:=Σ(Σ(mU(i,j)*(𝐞.i*𝐞.j), i,1..dim), j,1..dim)

\label{eq15}\begin{array}{@{}l}
\displaystyle
{{\frac{u^{2, \: 2}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}\ {|_{1 \  1}}}- 
\
\
\displaystyle
{{\frac{u^{1, \: 2}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}\ {|_{1 \  2}}}- 
\
\
\displaystyle
{{\frac{u^{2, \: 1}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}\ {|_{2 \  1}}}+ 
\
\
\displaystyle
{{\frac{u^{1, \: 1}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}\ {|_{2 \  2}}}
(15)
Type: LinearOperator(OrderedVariableList([1,2]),Expression(Integer))
fricas
-- compare
test(Ω=Ωm)

\label{eq16} \mbox{\rm true} (16)
Type: Boolean

Check that the snake relation holds

fricas
test
    (  I Ω   )  /
    (   U I  )  =  I

\label{eq17} \mbox{\rm true} (17)
Type: Boolean
fricas
test
    (   Ω I  )  /
    (  I U   )  =  I

\label{eq18} \mbox{\rm true} (18)
Type: Boolean

Dimension

This quantity depends on U!


\scalebox{1} % Change this value to rescale the drawing.
{
\begin{pspicture}(0,-0.92)(0.82,0.92)
\psbezier[linewidth=0.04](0.0,-0.1)(0.0,-0.9)(0.8,-0.9)(0.8,-0.1)
\psbezier[linewidth=0.04](0.0,0.1)(0.0,0.9)(0.8,0.9)(0.8,0.1)
\psline[linewidth=0.04cm](0.0,0.1)(0.0,-0.1)
\psline[linewidth=0.04cm](0.8,0.1)(0.8,-0.1)
\end{pspicture} 
}
 

fricas
d:=
    Ω /
    U

\label{eq19}\frac{{2 \ {u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{2, \: 1}}^{2}}-{{u^{1, \: 2}}^{2}}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}(19)
Type: LinearOperator(OrderedVariableList([1,2]),Expression(Integer))

This "twisted" quantity does not.


\scalebox{1} % Change this value to rescale the drawing.
{
\begin{pspicture}(0,-1.22)(0.82,1.22)
\psbezier[linewidth=0.04](0.0,-0.4)(0.0,-1.2)(0.8,-1.2)(0.8,-0.4)
\psbezier[linewidth=0.04](0.0,0.4)(0.0,1.2)(0.8,1.2)(0.8,0.4)
\psbezier[linewidth=0.04,linestyle=dashed,dash=0.16cm 0.16cm](0.0,0.4)(0.0,-0.2)(0.8,0.2)(0.8,-0.4)
\psbezier[linewidth=0.04](0.8,0.4)(0.8,-0.2)(0.0,0.2)(0.0,-0.4)
\end{pspicture} 
}
 
fricas
d':=
     Ω /
     X /
     U

\label{eq20}2(20)
Type: LinearOperator(OrderedVariableList([1,2]),Expression(Integer))

Symmetric Pairing

Repeat the calculation, assuming that U is symmetric.

fricas
sym:=groebner ravel(U-X/U)
There are 2 exposed and 6 unexposed library operations named groebner having 1 argument(s) but none was determined to be applicable. Use HyperDoc Browse, or issue )display op groebner to learn more about the available operations. Perhaps package-calling the operation or using coercions on the arguments will allow you to apply the operation.
Cannot find a definition or applicable library operation named groebner with argument type(s) List(Expression(Integer))
Perhaps you should use "@" to indicate the required return type, or "$" to specify which version of the function you need.

fricas
Um:=matrix Ξ(Ξ((𝐞.i*𝐞.j)/U, i,1..dim), j,1..dim)

\label{eq21}\left[ 
\begin{array}{cc}
{u^{1, \: 1}}&{u^{2, \: 1}}
\
{u^{1, \: 2}}&{u^{2, \: 2}}
(21)
Type: Matrix(LinearOperator(OrderedVariableList([1,2]),Expression(Integer)))
fricas
mU:=transpose inverse map(retract,Um)

\label{eq22}\left[ 
\begin{array}{cc}
{\frac{u^{2, \: 2}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}& -{\frac{u^{1, \: 2}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}
\
-{\frac{u^{2, \: 1}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}&{\frac{u^{1, \: 1}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}
(22)
Type: Matrix(Expression(Integer))
fricas
Ω:=Σ(Σ(mU(i,j)*(𝐞.i*𝐞.j), i,1..dim), j,1..dim)

\label{eq23}\begin{array}{@{}l}
\displaystyle
{{\frac{u^{2, \: 2}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}\ {|_{1 \  1}}}- 
\
\
\displaystyle
{{\frac{u^{1, \: 2}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}\ {|_{1 \  2}}}- 
\
\
\displaystyle
{{\frac{u^{2, \: 1}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}\ {|_{2 \  1}}}+ 
\
\
\displaystyle
{{\frac{u^{1, \: 1}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}}\ {|_{2 \  2}}}
(23)
Type: LinearOperator(OrderedVariableList([1,2]),Expression(Integer))

Check that the snake relation holds

fricas
test
    (  I Ω   )  /
    (   U I  )  =  I

\label{eq24} \mbox{\rm true} (24)
Type: Boolean
fricas
test
    (   Ω I  )  /
    (  I U   )  =  I

\label{eq25} \mbox{\rm true} (25)
Type: Boolean

These quantities no longer depends on U!

fricas
d:=
    Ω /
    U

\label{eq26}\frac{{2 \ {u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{2, \: 1}}^{2}}-{{u^{1, \: 2}}^{2}}}{{{u^{1, \: 1}}\ {u^{2, \: 2}}}-{{u^{1, \: 2}}\ {u^{2, \: 1}}}}(26)
Type: LinearOperator(OrderedVariableList([1,2]),Expression(Integer))

fricas
d':=
     Ω /
     X /
     U

\label{eq27}2(27)
Type: LinearOperator(OrderedVariableList([1,2]),Expression(Integer))

Twist dimension or twist snake? --Bill Page, Sun, 08 May 2011 14:16:39 -0700 reply
TwistedSnakeRelation