|
|
|
last edited 11 years ago by test1 |
Edit detail for SandBoxTensorProduct revision 7 of 9
| 1 2 3 4 5 6 7 8 9 | ||
|
Editor: Bill Page
Time: 2009/10/12 00:23:03 GMT-7 |
||
| Note: update | ||
changed: -FreeModule is defined in poly.spad with a local category. -FreeModuleCat is defined in xpoly.spad. -FreeModuleCat does not require OrderedSet, but FreeModule does. -FreeModule1 is of category FreeModuleCat, its Rep however is FreeModule. FreeModule is defined in poly.spad. FreeModuleCategory is defined in xpoly.spad. changed: - M1 : FreeModuleCat(R, B1), M2 : FreeModuleCat(R, B2)): TPcat == TPimp where - TPcat == Join(TensorProductCategory(R,M1,M2),FreeModuleCat(R,Product(B1,B2))) with M1 : FreeModuleCategory(R, B1), M2 : FreeModuleCategory(R, B2)): TPcat == TPimp where TPcat == Join(TensorProductCategory(R,M1,M2),FreeModuleCategory(R,Product(B1,B2))) with changed: - TPimp == FreeModule1(R, Product(B1, B2)) add TPimp == FreeModule(R, Product(B1, B2)) add changed: - ltx1:List TERM1 := ListOfTerms x1 - ltx2:List TERM2 := ListOfTerms x2 ltx1:List TERM1 := listOfTerms x1 ltx2:List TERM2 := listOfTerms x2 changed: - for t1 in ListOfTerms x1 repeat - for t2 in ListOfTerms x2 repeat for t1 in listOfTerms x1 repeat for t2 in listOfTerms x2 repeat changed: -M:=FreeModule1(Integer,Symbol) -N:=FreeModule1(Integer,Symbol) M:=FreeModule(Integer,Symbol) N:=FreeModule(Integer,Symbol)
Date: Fri, 15 May 2009 21:03:00 +0200 Franz Lehner wrote:
Attached is a prototype for tensor products. It is free modules over commutative rings.
FreeModule? is defined in poly.spad. FreeModuleCategory? is defined in xpoly.spad.
spad
)abbrev category TENSORC TensorProductCategory
TensorProductCategory(R:CommutativeRing, M : Module(R), N : Module(R)):Category == Module(R) with
product: (M, N) -> %
)abbrev category TENSORP TensorProductProperty
TensorProductProperty(R:CommutativeRing, M : Module(R), N : Module(R), _
MxN : TensorProductCategory(R, M, N), S : Module(R)): Category == with
eval: (MxN, (M, N) -> S) -> S
)abbrev package TENSORD TensorProduct
TensorProduct(R : CommutativeRing, B1 : OrderedSet, B2 : OrderedSet, _
M1 : FreeModuleCategory(R, B1), M2 : FreeModuleCategory(R, B2)): TPcat == TPimp where
TPcat == Join(TensorProductCategory(R,M1,M2),FreeModuleCategory(R,Product(B1,B2))) with
if M1 has Algebra(R) and M2 has Algebra(R) then Algebra(R)
TERM1 == Record(k: B1, c: R)
TERM2 == Record(k: B2, c: R)
B1xB2 == Product(B1,B2)
TERM == Record(k: B1xB2, c: R)
TPimp == FreeModule(R, Product(B1, B2)) add
import Rep, TERM1, TERM2, TERM, B1xB2
product(x1:M1,x2:M2):% ==
zero? x1 or zero? x2 => return 0
ltx1:List TERM1 := listOfTerms x1
ltx2:List TERM2 := listOfTerms x2
res : List TERM := []
for s1 in ltx1 repeat
for s2 in ltx2 repeat
res := concat!(res,[makeprod(s1.k, s2.k), s1.c*s2.c]$TERM)
res pretend %
if M1 has Algebra(R) and M2 has Algebra(R) then
(x1:% * x2:%):% ==
res : % := 0
for t1 in listOfTerms x1 repeat
for t2 in listOfTerms x2 repeat
-- the coefficients
t1c:R := t1.c
t2c:R := t2.c
-- the basis elements
t1k:B1xB2 := t1.k
t2k:B1xB2 := t2.k
t1a: M1 := monom(selectfirst t1k,1)
t1b: M2 := monom(selectsecond t1k,1)
t2a: M1 := monom(selectfirst t2k,1)
t2b: M2 := monom(selectsecond t2k,1)
res:= res + t1.c*t2.c *product(t1a*t2a,t1b*t2b)
res
spad
Compiling FriCAS source code from file
/var/zope2/var/LatexWiki/6037869360649691585-25px001.spad using
old system compiler.
TENSORC abbreviates category TensorProductCategory
------------------------------------------------------------------------
initializing NRLIB TENSORC for TensorProductCategory
compiling into NRLIB TENSORC
;;; *** |TensorProductCategory| REDEFINED
Time: 0.04 SEC.
finalizing NRLIB TENSORC
Processing TensorProductCategory for Browser database:
--->-->TensorProductCategory((product (% M N))): Not documented!!!!
--->-->TensorProductCategory(constructor): Not documented!!!!
--->-->TensorProductCategory(): Missing Description
; compiling file "/var/zope2/var/LatexWiki/TENSORC.NRLIB/TENSORC.lsp" (written 30 OCT 2009 02:28:19 PM):
; compiling (/VERSIONCHECK 2)
; compiling (DEFPARAMETER |TensorProductCategory;CAT| ...)
; compiling (DEFPARAMETER |TensorProductCategory;AL| ...)
; compiling (DEFUN |TensorProductCategory| ...)
; compiling (DEFUN |TensorProductCategory;| ...)
; /var/zope2/var/LatexWiki/TENSORC.NRLIB/TENSORC.fasl written
; compilation finished in 0:00:00.031
------------------------------------------------------------------------
TensorProductCategory is now explicitly exposed in frame initial
TensorProductCategory will be automatically loaded when needed from
/var/zope2/var/LatexWiki/TENSORC.NRLIB/TENSORC
TENSORP abbreviates category TensorProductProperty
------------------------------------------------------------------------
initializing NRLIB TENSORP for TensorProductProperty
compiling into NRLIB TENSORP
;;; *** |TensorProductProperty| REDEFINED
Time: 0.02 SEC.
finalizing NRLIB TENSORP
Processing TensorProductProperty for Browser database:
--->-->TensorProductProperty((eval (S MxN (Mapping S M N)))): Not documented!!!!
--->-->TensorProductProperty(constructor): Not documented!!!!
--->-->TensorProductProperty(): Missing Description
; compiling file "/var/zope2/var/LatexWiki/TENSORP.NRLIB/TENSORP.lsp" (written 30 OCT 2009 02:28:20 PM):
; compiling (/VERSIONCHECK 2)
; compiling (DEFPARAMETER |TensorProductProperty;CAT| ...)
; compiling (DEFPARAMETER |TensorProductProperty;AL| ...)
; compiling (DEFUN |TensorProductProperty| ...)
; compiling (DEFUN |TensorProductProperty;| ...)
; /var/zope2/var/LatexWiki/TENSORP.NRLIB/TENSORP.fasl written
; compilation finished in 0:00:00.017
------------------------------------------------------------------------
TensorProductProperty is now explicitly exposed in frame initial
TensorProductProperty will be automatically loaded when needed from
/var/zope2/var/LatexWiki/TENSORP.NRLIB/TENSORP
TENSORD abbreviates package TensorProduct
------------------------------------------------------------------------
initializing NRLIB TENSORD for TensorProduct
compiling into NRLIB TENSORD
importing Rep
importing Record(k: B1,c: R)
importing Record(k: B2,c: R)
importing Record(k: Product(B1,B2),c: R)
importing Product(B1,B2)
compiling exported product : (M1,M2) -> $
Time: 0.11 SEC.
****** Domain: M1 already in scope
augmenting M1: (Algebra R)
****** Domain: M2 already in scope
augmenting M2: (Algebra R)
compiling exported * : ($,$) -> $
Time: 0.12 SEC.
****** Domain: M1 already in scope
augmenting M1: (Algebra R)
****** Domain: M2 already in scope
augmenting M2: (Algebra R)
(time taken in buildFunctor: 10)
;;; *** |TensorProduct| REDEFINED
;;; *** |TensorProduct| REDEFINED
Time: 0.02 SEC.
Warnings:
[1] not known that (OrderedSet) is of mode (CATEGORY domain (IF (has B1 (Finite)) (IF (has B2 (Finite)) (ATTRIBUTE (Finite)) noBranch) noBranch) (IF (has B1 (Monoid)) (IF (has B2 (Monoid)) (ATTRIBUTE (Monoid)) noBranch) noBranch) (IF (has B1 (AbelianMonoid)) (IF (has B2 (AbelianMonoid)) (ATTRIBUTE (AbelianMonoid)) noBranch) noBranch) (IF (has B1 (CancellationAbelianMonoid)) (IF (has B2 (CancellationAbelianMonoid)) (ATTRIBUTE (CancellationAbelianMonoid)) noBranch) noBranch) (IF (has B1 (Group)) (IF (has B2 (Group)) (ATTRIBUTE (Group)) noBranch) noBranch) (IF (has B1 (AbelianGroup)) (IF (has B2 (AbelianGroup)) (ATTRIBUTE (AbelianGroup)) noBranch) noBranch) (IF (has B1 (OrderedAbelianMonoidSup)) (IF (has B2 (OrderedAbelianMonoidSup)) (ATTRIBUTE (OrderedAbelianMonoidSup)) noBranch) noBranch) (IF (has B1 (OrderedSet)) (IF (has B2 (OrderedSet)) (ATTRIBUTE (OrderedSet)) noBranch) noBranch) (SIGNATURE makeprod ($ B1 B2)) (SIGNATURE selectfirst (B1 $)) (SIGNATURE selectsecond (B2 $)))
Cumulative Statistics for Constructor TensorProduct
Time: 0.25 seconds
--------------non extending category----------------------
.. TensorProduct(#1,#2,#3,#4,#5) of cat
(|Join| (|TensorProductCategory| |#1| |#4| |#5|)
(|FreeModuleCategory| |#1| (|Product| |#2| |#3|))
(CATEGORY |package|
(IF (|has| |#4| (|Algebra| |#1|))
(IF (|has| |#5| (|Algebra| |#1|))
(ATTRIBUTE (|Algebra| |#1|))
|noBranch|)
|noBranch|))) has no
(IF (|has| |#1| (|CommutativeRing|))
(ATTRIBUTE (|Module| |#1|))
|noBranch|) finalizing NRLIB TENSORD
Processing TensorProduct for Browser database:
--->-->TensorProduct(): Missing Description
; compiling file "/var/zope2/var/LatexWiki/TENSORD.NRLIB/TENSORD.lsp" (written 30 OCT 2009 02:28:21 PM):
; compiling (/VERSIONCHECK 2)
; compiling (DEFUN |TENSORD;product;M1M2$;1| ...)
; compiling (DEFUN |TENSORD;*;3$;2| ...)
; compiling (DEFUN |TensorProduct| ...)
; compiling (DEFUN |TensorProduct;| ...)
; compiling (MAKEPROP (QUOTE |TensorProduct|) ...)
; /var/zope2/var/LatexWiki/TENSORD.NRLIB/TENSORD.fasl written
; compilation finished in 0:00:00.132
------------------------------------------------------------------------
TensorProduct is now explicitly exposed in frame initial
TensorProduct will be automatically loaded when needed from
/var/zope2/var/LatexWiki/TENSORD.NRLIB/TENSORDaxiom
M:=FreeModule(Integer,Symbol)
 | (1) |
Type: Domain
axiom
N:=FreeModule(Integer,Symbol)
 | (2) |
Type: Domain
axiom
a1:='a1::M
 | (3) |
Type: FreeModule?(Integer,Symbol)
axiom
a2:='a2::M
 | (4) |
Type: FreeModule?(Integer,Symbol)
axiom
b1:='b1::N
 | (5) |
Type: FreeModule?(Integer,Symbol)
axiom
b2:='b2::N
 | (6) |
Type: FreeModule?(Integer,Symbol)
axiom
MxN:=TensorProduct(Integer,Symbol,Symbol,M,N);
Type: Domain
axiom
t:=product(a1+a2,b1+b2)$MxN;
axiom
t
 | (7) |
axiom
leadingMonomial t
 | (8) |
Type: Product(Symbol,Symbol)
axiom
numberOfMonomials t
 | (9) |
Type: PositiveInteger?
Demonstrating the axioms of the tensor product:
axiom
x:M
Type: Void
axiom
y:M
Type: Void
axiom
u:M
Type: Void
axiom
p:=2*x+3*u
 | (10) |
Type: FreeModule?(Integer,Symbol)
axiom
q:=5*x+7*y+11*u
 | (11) |
Type: FreeModule?(Integer,Symbol)
axiom
MxM:=TensorProduct(Integer,Symbol,Symbol,M,M);
Type: Domain
axiom
r:=product(p,q)$MxM
 | (12) |
axiom
w:= 13*y+17*y+19*u
 | (13) |
Type: FreeModule?(Integer,Symbol)
axiom
test( product(p+q,w)$MxM = product(p,w)$MxM + product(q,w)$MxM )
 | (14) |
Type: Boolean
axiom
test( product(p,q+w)$MxM = product(p,q)$MxM + product(p,w)$MxM )
 | (15) |
Type: Boolean
axiom
test( product(p,23*w)$MxM = 23*product(p,w)$MxM )
 | (16) |
Type: Boolean
axiom
test( product(23*p,w)$MxM = 23*product(p,w)$MxM )
 | (17) |
Type: Boolean
first example of a bialgebra --Bill Page, Wed, 13 May 2009 08:06:57 -0700 reply
SandBoxHopfAlgebra?