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

Edit detail for Snake Relation revision 2 of 11

1 2 3 4 5 6 7 8 9 10 11
Editor: Bill Page
Time: 2011/04/22 07:58:53 GMT-7
Note: dimension

added:

$$
\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} 
}
$$

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

changed:
-d1:=(I*Ω)/(U*I);
-d2:=(Ω*I)/(I*U);
Iright:=(I*Ω)/(U*I);
Ileft:=(Ω*I)/(I*U);

changed:
-eq1:=equate(d1,I);
-eq2:=equate(d2,I);
eq1:=equate(Iright,I);
eq2:=equate(Ileft,I);

added:
\end{axiom}

This is equivalent to a matrix inverse (transposed!)
\begin{axiom}
Um:=matrix Ξ(Ξ((𝐞.i*𝐞.j)/U, i,1..dim), j,1..dim)
mU:=transpose inverse map(retract,Um)
Ωm:=Σ(Σ(mU(i,j)*(𝐞.i*𝐞.j), i,1..dim), j,1..dim)

changed:
-inverse map(retract,Um)
test(Ω=Ωm)

changed:
-Check "dimension": It depends on parameters!
Check that the snake relation holds

removed:
-
-d:𝐋:=
-    (   Ω    )  /
-    (   U    )

added:
The quantity "dimension" depends on $U$!

\begin{axiom}
d:𝐋:=

       Ω      /
       U    

\end{axiom}

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 need the Axiom LinearOperator? library.

fricas
(1) -> )library MONAL PROP LIN
)library cannot find the file LIN. 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

Use the following macros for convenient notation

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

𝐋 is the domain of 4-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(dim, OVAR [], ℚ)
There are no library operations named LinearOperator Use HyperDoc Browse or issue )what op LinearOperator to learn if there is any operation containing " LinearOperator " in its name.
Cannot find a definition or applicable library operation named LinearOperator with argument type(s) PositiveInteger Type Type
Perhaps you should use "@" to indicate the required return type, or "$" to specify which version of the function you need.

A scalar product (pairing) is denoted by

fricas
U:=Σ(Σ(script('u,[[],[i,j]])*𝐝.i*𝐝.j, i,1..dim), j,1..dim)
There are no library operations named 𝐝 Use HyperDoc Browse or issue )what op 𝐝 to learn if there is any operation containing " 𝐝 " in its name. Cannot find a definition or applicable library operation named 𝐝 with argument type(s) PositiveInteger
Perhaps you should use "@" to indicate the required return type, or "$" to specify which version of the function you need. FriCAS will attempt to step through and interpret the code. There are no library operations named 𝐝 Use HyperDoc Browse or issue )what op 𝐝 to learn if there is any operation containing " 𝐝 " in its name.
Cannot find a definition or applicable library operation named 𝐝 with argument type(s) PositiveInteger
Perhaps you should use "@" to indicate the required return type, or "$" to specify which version of the function you need.

Co-pairing

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

fricas
Ω:𝐋:=Σ(Σ(script('u,[[i,j]])*𝐞.i*𝐞.j, i,1..dim), j,1..dim)
𝐋 is not a valid type.

This is equivalent to a matrix inverse (transposed!)

fricas
Um:=matrix Ξ(Ξ((𝐞.i*𝐞.j)/U, i,1..dim), j,1..dim)
There are no library operations named 𝐞 Use HyperDoc Browse or issue )what op 𝐞 to learn if there is any operation containing " 𝐞 " in its name. Cannot find a definition or applicable library operation named 𝐞 with argument type(s) PositiveInteger
Perhaps you should use "@" to indicate the required return type, or "$" to specify which version of the function you need. FriCAS will attempt to step through and interpret the code. There are no library operations named 𝐞 Use HyperDoc Browse or issue )what op 𝐞 to learn if there is any operation containing " 𝐞 " in its name.
Cannot find a definition or applicable library operation named 𝐞 with argument type(s) PositiveInteger
Perhaps you should use "@" to indicate the required return type, or "$" to specify which version of the function you need.

Check that the snake relation holds

fricas
test
    (  I Ω   )  /
    (   U I  )  =  I
There are no exposed library operations named I but there is one unexposed operation with that name. Use HyperDoc Browse or issue )display op I to learn more about the available operation.
Cannot find a definition or applicable library operation named I with argument type(s) Variable(Ω)
Perhaps you should use "@" to indicate the required return type, or "$" to specify which version of the function you need.

The quantity "dimension" depends on U!

fricas
d:𝐋:=
Ω / U
𝐋 is not a valid type.