|
|
last edited 14 years ago by Bill Page |
1 2 3 4 5 6 7 8 | ||
Editor: Bill Page
Time: 2011/06/08 20:40:40 GMT-7 |
||
Note: snails, turtles and deer |
changed: ---eigenvectors(u1::Matrix FRAC POLY INT) -\end{axiom} \end{axiom} \begin{axiom} )set output algebra on )set output tex off \end{axiom} \begin{axiom} eigenvectors(u1::Matrix FRAC POLY INT) \end{axiom} \begin{axiom} )set output algebra off )set output tex on \end{axiom} changed: - ( Y Y ) ( Y Y ) added: Bi-algebra conditions \begin{axiom} ΦΦ:= _ ( λ λ ) / _ ( I X I ) / _ ( Y Y ) test((e,e)/ΦΦ=φφ) \end{axiom} Y-forms Three traces of two graftings of an algebra gives six (2,0)-forms. Left snail and right snail:: LS RS Y /\ /\ Y Y ) ( Y \/ \/ i j j i \/ \/ \ /\ /\ / e f \ / f e \/ \ / \/ \ / \ / f / \ f \/ \/ \begin{axiom} LS:= ( Y Λ )/ _ ( Y I )/ _ V RS:= ( Λ Y )/ _ ( I Y )/ _ V test(LS=RS) \end{axiom} Left and right deer:: RD LD \ /\/ \/\ / Y /\ /\ Y Y ) ( Y \/ \/ i j i j \ /\ / \ /\ / \ f \ / \ / f / \/ \/ \/ \/ \ /\ /\ / e / \ / \ e \/ \ / \/ \ / \ / f / \ f \/ \/ Left and right deer forms are identical but different from snails. \begin{axiom} RD:= ( I Λ I ) / _ ( Y X ) / _ ( Y I ) / _ V LD:= ( I Λ I ) / _ ( X Y ) / _ ( I Y ) / _ V test(LD=RD) test(RD=RS) test(RD=LS) \end{axiom} Left and right turtles:: RT LT /\ / / \ \ /\ ( Y / \ Y ) \ Y Y / \/ \/ i j i j /\ / / \ \ /\ / f / / \ \ f \ / \/ / \ \/ \ \ \ / \ / / \ e / \ e / \ \/ \/ / \ / \ / \ f f / \/ \/ \begin{axiom} RT:= ( Λ I I ) / _ ( I Y I ) / _ ( I Y ) / _ V LT:= ( I I Λ ) / _ ( I Y I ) / _ ( Y I ) / _ V test(LT=RT) \end{axiom} The turles are symmetric \begin{axiom} test(RT = X/RT) test(LT = X/LT) \end{axiom} Five of the six forms are independent. \begin{axiom} test(RT=RS) test(RT=LS) test(RT=RD) test(LT=RS) test(LT=LS) test(LT=RD) \end{axiom}
Grassmann Algebra Is Frobenius In Many Ways
A -dimensional vector space represents Grassmann algebra with
generators
Linear operators over a 4-dimensional vector space representing Grassmann algebra with two generators.
Ref:
-permuted Frobenius Algebras
Zbigniew Oziewicz (UNAM), Gregory Peter Wene (UTSA)
Frobenius algebras and 2D topological quantum field theories
Joachim Kock
We need the Axiom LinearOperator library.
(1) -> )library CARTEN ARITY CMONAL CPROP CLOP CALEY
>> System error: The value 15684 is not of type LIST
Use the following macros for convenient notation
-- summation macro Σ(x,i, n)==reduce(+, [x for i in n])
-- list macro Ξ(f,i, n)==[f for i in n]
-- subscript and superscripts macro sb == subscript
macro sp == superscript
𝐋 is the domain of 4-dimensional linear operators over the rational functions ℚ (Expression Integer), i.e. ratio of polynomials with integer coefficients.
dim:=4
![]() | (1) |
macro ℒ == List
macro ℂ == CaleyDickson
macro ℚ == Expression Integer
𝐋 := ClosedLinearOperator(OVAR ['1,'i, 'j, 'k], ℚ)
There are no library operations named ClosedLinearOperator Use HyperDoc Browse or issue )what op ClosedLinearOperator to learn if there is any operation containing " ClosedLinearOperator " in its name.
Cannot find a definition or applicable library operation named ClosedLinearOperator with argument type(s) Type Type
Perhaps you should use "@" to indicate the required return type,or "$" to specify which version of the function you need.
Generate structure constants for Grassmann Algebra
The structure constants can be obtained by dividing each matrix entry by the list of basis vectors.
Grassmann algebra will be specified by setting the Caley-Dickson parameters (i2, j2) to zero.
i2:=sp('i,[2])
![]() | (2) |
j2:=sp('j,[2])
![]() | (3) |
QQ:=CliffordAlgebra(2,ℚ, matrix [[i2, 0], [0, j2]])
![]() | (4) |
B:ℒ QQ := [monomial(1,[]), monomial(1, [1]), monomial(1, [2]), monomial(1, [1, 2])]
![]() | (5) |
M:Matrix QQ := matrix Ξ(Ξ(B.i*B.j,i, 1..dim), j, 1..dim)
![]() | (6) |
S(y) == map(x +-> coefficient(recip(y)*x,[]), M)
ѕ :=map(S,B)::ℒ ℒ ℒ ℚ
Compiling function S with type CliffordAlgebra(2,Expression(Integer) , [[i[;2], 0], [0, j[;2]]]) -> Matrix(Expression(Integer))
![]() | (7) |
-- structure constants form a tensor operator --Y := Σ(Σ(Σ(ѕ(i)(k)(j)*𝐞.i*𝐝.j*𝐝.k,i, 1..dim), j, 1..dim), k, 1..dim) Y := eval(Σ(Σ(Σ(ѕ(i)(k)(j)*𝐞.i*𝐝.j*𝐝.k, i, 1..dim), j, 1..dim), k, 1..dim), [i2=0, j2=0])
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.
Units
e:=𝐞.1; i:=𝐞.2; j:=𝐞.3; k:=𝐞.4;
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.
Multiplication of arbitrary Grassmann numbers and
a:=Σ(sb('a,[i])*𝐞.i, i, 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.
Multiplication is Associative
test( ( I Y ) / _ ( Y ) = _ ( Y I ) / _ ( Y ) )
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(Y)
Perhaps you should use "@" to indicate the required return type,or "$" to specify which version of the function you need.
A scalar product is denoted by the (2,0)-tensor
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.
We say that the scalar product is associative if the tensor equation holds:
Y = Y U U
In other words, if the (3,0)-tensor:
![]() |
![]() | (8) |
Using the LinearOperator domain in Axiom and some carefully chosen symbols we can easily enter expressions that are both readable and interpreted by Axiom as "graphical calculus" diagrams describing complex products and compositions of linear operators.
ω:𝐋 := ( Y I ) / U - ( I Y ) / U
𝐋 is not a valid type.
An algebra with a non-degenerate associative scalar product is called a [Frobenius Algebra]?.
The Cartan-Killing Trace
Ú:= ( Y Λ ) / _ ( Y I ) / _ V
There are no exposed library operations named Y but there are 2 unexposed operations with that name. Use HyperDoc Browse or issue )display op Y to learn more about the available operations.
Cannot find a definition or applicable library operation named Y 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.
forms is degenerate
Ũ := Ù
![]() | (9) |
test ( Y I ) / Ũ = ( I Y ) / Ũ
There are no exposed library operations named Y but there are 2 unexposed operations with that name. Use HyperDoc Browse or issue )display op Y to learn more about the available operations.
Cannot find a definition or applicable library operation named Y with argument type(s) Variable(I)
Perhaps you should use "@" to indicate the required return type,or "$" to specify which version of the function you need.
General Solution
We may consider the problem where multiplication Y is given,
and look for all associative scalar products
This problem can be solved using linear algebra.
)expose MCALCFN
MultiVariableCalculusFunctions is now explicitly exposed in frame initial J := jacobian(ravel ω,concat map(variables, ravel U)::ℒ Symbol);
There are 1 exposed and 0 unexposed library operations named ravel having 1 argument(s) but none was determined to be applicable. Use HyperDoc Browse,or issue )display op ravel 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 ravel 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 matrix J
transforms the coefficients of the tensor
into coefficients of the tensor
. We are looking for
the general linear family of tensors
such that
J
transforms into
for any such
.
If the null space of the J
matrix is not empty we can use
the basis to find all non-trivial solutions for U:
Ñ:=nullSpace(J)
There are 3 exposed and 3 unexposed library operations named nullSpace having 1 argument(s) but none was determined to be applicable. Use HyperDoc Browse,or issue )display op nullSpace 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 nullSpace with argument type(s) Variable(J)
Perhaps you should use "@" to indicate the required return type,or "$" to specify which version of the function you need.
This defines a family of pre-Frobenius algebras:
zero? eval(ω,ℰ)
There are 10 exposed and 6 unexposed library operations named eval having 2 argument(s) but none was determined to be applicable. Use HyperDoc Browse,or issue )display op eval 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 eval with argument type(s) Variable(ω) Variable(ℰ)
Perhaps you should use "@" to indicate the required return type,or "$" to specify which version of the function you need.
Frobenius Form (co-unit)
d:=ε1*𝐝.1+εi*𝐝.2+εj*𝐝.3+εk*𝐝.4
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.
Express scalar product in terms of Frobenius form
𝔓:=solve(𝔇,Ξ(sb('p, [i]), i, 1..#Ñ)).1
There are 2 exposed and 2 unexposed library operations named # having 1 argument(s) but none was determined to be applicable. Use HyperDoc Browse,or issue )display op # 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 # 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.
In general the pairing is not symmetric!
u1:=matrix Ξ(Ξ(retract((𝐞.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.
)set output algebra on
)set output tex off
eigenvectors(u1::Matrix FRAC POLY INT)
Cannot convert the value from type Variable(u1) to Matrix(Fraction( Polynomial(Integer))) .
)set output algebra off
)set output tex on
The scalar product must be non-degenerate:
Ů:=determinant u1
There are 3 exposed and 1 unexposed library operations named determinant having 1 argument(s) but none was determined to be applicable. Use HyperDoc Browse,or issue )display op determinant 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 determinant with argument type(s) Variable(u1)
Perhaps you should use "@" to indicate the required return type,or "$" to specify which version of the function you need.
Frobenius scalar product of "vectors" and
a:=sb('a,[1])*i+sb('a, [2])*j
![]() | (10) |
b:=sb('b,[1])*i+sb('b, [2])*j
![]() | (11) |
(a,a)/Ų
There are 11 exposed and 15 unexposed library operations named / having 2 argument(s) but none was determined to be applicable. Use HyperDoc Browse,or issue )display op / 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 / with argument type(s) Tuple(Polynomial(Integer)) Variable(Ų)
Perhaps you should use "@" to indicate the required return type,or "$" to specify which version of the function you need.
Co-scalar product
Solve the Snake Relation as a system of linear equations.
Ω:𝐋:=Σ(Σ(script('u,[[i, j]])*𝐞.i*𝐞.j, i, 1..dim), j, 1..dim)
𝐋 is not a valid type.
matrix Ξ(Ξ(retract(Ω/(𝐝.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.
Check "dimension" and the snake relations.
O:𝐋:= Ω / Ų
𝐋 is not a valid type.
Co-algebra
Compute the "three-point" function and use it to define co-multiplication.
W:= (Y I) / Ų
There are no exposed library operations named Y but there are 2 unexposed operations with that name. Use HyperDoc Browse or issue )display op Y to learn more about the available operations.
Cannot find a definition or applicable library operation named Y with argument type(s) Variable(I)
Perhaps you should use "@" to indicate the required return type,or "$" to specify which version of the function you need.
test ( I ΩX ) / ( Y 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(ΩX)
Perhaps you should use "@" to indicate the required return type,or "$" to specify which version of the function you need.
Co-associativity
test( ( λ ) / _ ( 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.
test e / λ = ΩX
![]() | (12) |
Frobenius Condition (fork)
H := Y / λ
![]() | (13) |
test ( λ I ) / ( I Y ) = H
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) Variable(I)
Perhaps you should use "@" to indicate the required return type,or "$" to specify which version of the function you need.
Handle
Φ := λ / Y
![]() | (14) |
Figure 12
φφ:= _ ( Ω Ω ) / _ ( X I I ) / _ ( I X I ) / _ ( I I X ) / _ ( Y Y )
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) Variable(Ω)
Perhaps you should use "@" to indicate the required return type,or "$" to specify which version of the function you need.
Bi-algebra conditions
ΦΦ:= _ ( λ λ ) / _ ( I I X ) / _ ( I X I ) / _ ( I I X ) / _ ( Y Y )
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) Variable(λ)
Perhaps you should use "@" to indicate the required return type,or "$" to specify which version of the function you need.
Bi-algebra conditions
ΦΦ:= _ ( λ λ ) / _ ( I X I ) / _ ( Y Y )
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) Variable(λ)
Perhaps you should use "@" to indicate the required return type,or "$" to specify which version of the function you need.
Three traces of two graftings of an algebra gives six (2,0)-forms.
Left snail and right snail:
LS RS Y /\ /\ Y Y ) ( Y \/ \/ i j j i \/ \/ \ /\ /\ / e f \ / f e \/ \ / \/ \ / \ / f / \ f \/ \/
LS:= ( Y Λ )/ _ ( Y I )/ _ V
There are no exposed library operations named Y but there are 2 unexposed operations with that name. Use HyperDoc Browse or issue )display op Y to learn more about the available operations.
Cannot find a definition or applicable library operation named Y 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.
Left and right deer:
RD LD \ /\/ \/\ / Y /\ /\ Y Y ) ( Y \/ \/ i j i j \ /\ / \ /\ / \ f \ / \ / f / \/ \/ \/ \/ \ /\ /\ / e / \ / \ e \/ \ / \/ \ / \ / f / \ f \/ \/
Left and right deer forms are identical but different from snails.
RD:= ( I Λ I ) / _ ( Y X ) / _ ( Y I ) / _ V
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) Variable(I)
Perhaps you should use "@" to indicate the required return type,or "$" to specify which version of the function you need.
Left and right turtles:
RT LT /\ / / \ \ /\ ( Y / \ Y ) \ Y Y / \/ \/ i j i j /\ / / \ \ /\ / f / / \ \ f \ / \/ / \ \/ \ \ \ / \ / / \ e / \ e / \ \/ \/ / \ / \ / \ f f / \/ \/
RT:= ( Λ I I ) / _ ( I Y I ) / _ ( I Y ) / _ V
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(I)
Perhaps you should use "@" to indicate the required return type,or "$" to specify which version of the function you need.
The turles are symmetric
test(RT = X/RT)
![]() | (15) |
test(LT = X/LT)
![]() | (16) |
Five of the six forms are independent.
test(RT=RS)
![]() | (17) |
test(RT=LS)
![]() | (18) |
test(RT=RD)
![]() | (19) |
test(LT=RS)
![]() | (20) |
test(LT=LS)
![]() | (21) |
test(LT=RD)
![]() | (22) |