help links subscribe changes refresh edit

	SandBox
SandBoxNijenhuis ←	↑	→ SandBoxNonZeroInteger

SandBoxNonAssociativeAlgebra

last edited 15 years ago by Bill Page

Edit detail for SandBoxNonAssociativeAlgebra revision 24 of 25

	1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
Editor: Bill Page Time: 2008/05/08 07:25:59 GMT-7
Note: cleanup

added:
Multiplication

added:
\end{axiom}

Commutator and Associator
\begin{axiom}
matrix([[commutator(a.x,a.y) for x in 1..4] for y in 1..4])$OutputForm
\end{axiom}

\begin{axiom}
[matrix([[associator(a.x,a.y,a.z) for x in 1..4] for y in 1..4])$OutputForm for z in 1..4]
for x in 1..4 repeat
  for y in 1..4 repeat
    for z in 1..4 repeat
      if associator(a.x,a.y,a.z) ~= 0$V then
        output([[a.x,a.y,a.z],"=",associator(a.x,a.y,a.z)])
\end{axiom}

Volume form?
\begin{axiom}
a.2 * (a.3 * a.4) = (a.2 * a.3) * a.4
\end{axiom}

Check standard properties
\begin{axiom}

added:
Commuting elements

added:
Symbolic computations

changed:
-leftRankPolynomial()$G
q:=leftRankPolynomial()$G
map(factor,coefficients q)

removed:
-From BillPage Fri May 2 08:38:36 -0700 2008
-From: Bill Page
-Date: Fri, 02 May 2008 08:38:36 -0700
-Subject: commutator and associator
-Message-ID: <20080502083836-0700@axiom-wiki.newsynthesis.org>
-
-\begin{axiom}
-matrix([[commutator(a.x,a.y) for x in 1..4] for y in 1..4])$OutputForm
-\end{axiom}
-
-\begin{axiom}
-[matrix([[associator(a.x,a.y,a.z) for x in 1..4] for y in 1..4])$OutputForm for z in 1..4]
-for x in 1..4 repeat
-  for y in 1..4 repeat
-    for z in 1..4 repeat
-      if associator(a.x,a.y,a.z) ~= 0$V then
-        output([[a.x,a.y,a.z],"=",associator(a.x,a.y,a.z)])
-\end{axiom}
-
-Volume form?
-\begin{axiom}
-a.2 * (a.3 * a.4) = (a.2 * a.3) * a.4
-\end{axiom}

Ref: http://arxiv.org/abs/0711.3220

Fourvector algebra

Author: Diego Saa (Submitted on 20 Nov 2007)

Abstract: The algebra of fourvectors is described. The fourvectors are more appropriate than the Hamilton quaternions for its use in Physics and the sciences in general. The fourvectors embrace the 3D vectors in a natural form. It is shown the excellent ability to perform rotations with the use of fourvectors, as well as their use in relativity for producing Lorentz boosts, which are understood as simple rotations.

axiom
_*_*(x,y)==concat(x(1) * y(1) + dot(x(2..), y(2..)), x(1) * y(2..) - x(2..) * y(1) + cross(x(2..), y(2..)))

Type: Void

axiom
e:Vector INT:=[1,0,0,0]

(1)

Type: Vector Integer

axiom
i:Vector INT:=[0,1,0,0]

(2)

Type: Vector Integer

axiom
j:Vector INT:=[0,0,1,0]

(3)

Type: Vector Integer

axiom
k:Vector INT:=[0,0,0,1]

(4)

Type: Vector Integer

axiomtest(e**e=e)  and _
test(i**i=e)  and _
test(j**j=e)  and _
test(k**k=e)  and _
test(e**i=i)  and _
test(e**j=j)  and _
test(e**k=k)  and _
test(i**e=-i) and _
test(j**e=-j) and _
test(k**e=-k) and _
test(i**j=k)  and _
test(j**i=-k) and _
test(k**i=j)  and _
test(i**k=-j) and _
test(j**k=i)  and _
test(k**j=-i)

axiom

Compiling function ** with type (Vector Integer,Vector Integer) -> 
      Vector Integer

(5)

Type: Boolean

Axiom has a domain for NonAssociative? Algebra

This is documented in the article: Computations in Algebras of Fixed Rank by Johannes Grabmeir and Robert Wisbauer, from the book "Computational Algebra" By Klaus G. Fischer, Philippe Loustaunau, Jay Shapiro.

The algebra above can be given by structural constants.

axiom)clear all
   All user variables and function definitions have been cleared.
sc:Vector Matrix Fraction Integer := [ _
[[ 1, 0, 0, 0], _
 [ 0, 1, 0, 0], _
 [ 0, 0, 1, 0], _
 [ 0, 0, 0, 1]], _
[[ 0, 1, 0, 0], _
 [-1, 0, 0, 0], _
 [ 0, 0, 0, 1], _
 [ 0, 0,-1, 0]], _
[[ 0, 0, 1, 0], _
 [ 0, 0, 0,-1], _
 [-1, 0, 0, 0], _
 [ 0, 1, 0, 0]], _
[[ 0, 0, 0, 1], _
 [ 0, 0, 1, 0], _
 [ 0,-1, 0, 0], _
 [-1, 0, 0, 0]]];

Type: Vector Matrix Fraction Integer

axiom
V:=AlgebraGivenByStructuralConstants(Fraction Integer, 4, [e,i,j,k],sc)

(6)

Type: Domain

Multiplication

axiom
a:=basis()$V

(7)

Type: Vector AlgebraGivenByStructuralConstants?(Fraction Integer,4,[e,i,j,k]?,[MATRIX,MATRIX,MATRIX,MATRIX]?)

axiom
matrix([[(a.i * a.j) for j in 1..4] for i in 1..4])$OutputForm

(8)

Type: OutputForm?

Commutator and Associator

axiom
matrix([[commutator(a.x,a.y) for x in 1..4] for y in 1..4])$OutputForm

(9)

Type: OutputForm?

axiom
[matrix([[associator(a.x,a.y,a.z) for x in 1..4] for y in 1..4])$OutputForm for z in 1..4]

(10)

Type: List OutputForm?

axiomfor x in 1..4 repeat
  for y in 1..4 repeat
    for z in 1..4 repeat
      if associator(a.x,a.y,a.z) ~= 0$V then
        output([[a.x,a.y,a.z],"=",associator(a.x,a.y,a.z)])
   [[i,e,e],"=",2i]
   [[i,e,i],"=",- 2e]
   [[i,e,j],"=",- 2k]
   [[i,e,k],"=",2j]
   [[i,i,e],"=",2e]
   [[i,i,i],"=",2i]
   [[i,i,j],"=",2j]
   [[i,i,k],"=",2k]
   [[j,e,e],"=",2j]
   [[j,e,i],"=",2k]
   [[j,e,j],"=",- 2e]
   [[j,e,k],"=",- 2i]
   [[j,j,e],"=",2e]
   [[j,j,i],"=",2i]
   [[j,j,j],"=",2j]
   [[j,j,k],"=",2k]
   [[k,e,e],"=",2k]
   [[k,e,i],"=",- 2j]
   [[k,e,j],"=",2i]
   [[k,e,k],"=",- 2e]
   [[k,k,e],"=",2e]
   [[k,k,i],"=",2i]
   [[k,k,j],"=",2j]
   [[k,k,k],"=",2k]

Type: Void

Volume form?

axiom
a.2 * (a.3 * a.4) = (a.2 * a.3) * a.4

(11)

Type: Equation AlgebraGivenByStructuralConstants?(Fraction Integer,4,[e,i,j,k]?,[MATRIX,MATRIX,MATRIX,MATRIX]?)

Check standard properties

axiom
leftUnit()$V

(12)

Type: Union(AlgebraGivenByStructuralConstants?(Fraction Integer,4,[e,i,j,k]?,[MATRIX,MATRIX,MATRIX,MATRIX]?),...)

axiomrightUnit()$V
   this algebra has no right unit

(13)

Type: Union("failed",...)

axiomalternative?()$V
   algebra is not left alternative

(14)

Type: Boolean

axiomleftAlternative?()$V
   algebra is not left alternative

(15)

Type: Boolean

axiomrightAlternative?()$V
   algebra is not right alternative

(16)

Type: Boolean

axiomassociative?()$V
   algebra is not associative

(17)

Type: Boolean

axiomantiAssociative?()$V
   algebra is not anti-associative

(18)

Type: Boolean

axiom--powerAssociative?()$V
commutative?()$V
   algebra is not commutative

(19)

Type: Boolean

axiomjordanAlgebra?()$V
   algebra is not commutative
   this is not a Jordan algebra

(20)

Type: Boolean

axiomjordanAdmissible?()$V
   algebra is not Jordan admissible

(21)

Type: Boolean

axiomnoncommutativeJordanAlgebra?()$V
   algebra is not flexible
   this is not a noncommutative Jordan algebra, as it is not flexible

(22)

Type: Boolean

axiomlieAlgebra?()$V
   algebra is not anti-commutative
   this is not a Lie algebra

(23)

Type: Boolean

axiomlieAdmissible?()$V
   algebra is not Lie admissible

(24)

Type: Boolean

axiomjacobiIdentity?()$V
   Jacobi identity does not hold

(25)

Type: Boolean

Commuting elements

axiom
V has FramedNonAssociativeAlgebra(Fraction Integer)

(26)

Type: Boolean

axiom
basisOfCommutingElements()$AlgebraPackage(Fraction Integer,V)

(27)

Type: List AlgebraGivenByStructuralConstants?(Fraction Integer,4,[e,i,j,k]?,[MATRIX,MATRIX,MATRIX,MATRIX]?)

axiom
basisOfCenter()$AlgebraPackage(Fraction Integer,V)

(28)

Type: List AlgebraGivenByStructuralConstants?(Fraction Integer,4,[e,i,j,k]?,[MATRIX,MATRIX,MATRIX,MATRIX]?)

axiom
basisOfCentroid()$AlgebraPackage(Fraction Integer,V)

(29)

Type: List Matrix Fraction Integer

axiom
basisOfNucleus()$AlgebraPackage(Fraction Integer,V)

(30)

Type: List AlgebraGivenByStructuralConstants?(Fraction Integer,4,[e,i,j,k]?,[MATRIX,MATRIX,MATRIX,MATRIX]?)

axiom
basisOfLeftNucloid()$AlgebraPackage(Fraction Integer,V)

(31)

Type: List Matrix Fraction Integer

Symbolic computations

axiom
G:=GenericNonAssociativeAlgebra(Fraction Integer, 4, [e,i,j,k],sc)

(32)

Type: Domain

Look for Idempotents

axiom
conditionsForIdempotents()$G

(33)

Type: List Polynomial Fraction Integer

axiom
gb:=groebnerFactorize %

(34)

Type: List List Polynomial Fraction Integer

axiom
associatorDependence()$G

(35)

Type: List Vector Fraction Polynomial Fraction Integer

axiom
q:=leftRankPolynomial()$G

(36)

Type: SparseUnivariatePolynomial? Fraction Polynomial Fraction Integer

axiom
map(factor,coefficients q)

(37)

Type: List Factored Fraction Polynomial Fraction Integer

axiomrightUnit()$G
   this algebra has no right unit

(38)

Type: Union("failed",...)

axiom
p1:=generic([x1,y1,z1,w1])$G

(39)

Type: GenericNonAssociativeAlgebra?(Fraction Integer,4,[e,i,j,k]?,[MATRIX,MATRIX,MATRIX,MATRIX]?)

axiom
p2:=generic([x2,y2,z2,w2])$G

(40)

Type: GenericNonAssociativeAlgebra?(Fraction Integer,4,[e,i,j,k]?,[MATRIX,MATRIX,MATRIX,MATRIX]?)

axiom
p3:=generic([x3,y3,z3,w3])$G

(41)

Type: GenericNonAssociativeAlgebra?(Fraction Integer,4,[e,i,j,k]?,[MATRIX,MATRIX,MATRIX,MATRIX]?)

axiom
leftRecip(p1)$G

(42)

Type: Union(GenericNonAssociativeAlgebra?(Fraction Integer,4,[e,i,j,k]?,[MATRIX,MATRIX,MATRIX,MATRIX]?),...)

axiomrightRecip(p1)$G
   this algebra has no right unit

(43)

Type: Union("failed",...)

axiom
leftRegularRepresentation(p1)

(44)

Type: Matrix Fraction Polynomial Fraction Integer

axiom
rightRegularRepresentation(p1)

(45)

Type: Matrix Fraction Polynomial Fraction Integer

axiom
associator(p1,p2,p3)$G

(46)

Type: GenericNonAssociativeAlgebra?(Fraction Integer,4,[e,i,j,k]?,[MATRIX,MATRIX,MATRIX,MATRIX]?)