## reduce-algebra-developers — Discussion of development, administration and support for Reduce

You can subscribe to this list here.

 2009 2010 2011 2012 2013 2014 2015 2016 2017 Jan (2) Feb (5) Mar Apr May (2) Jun (8) Jul (4) Aug Sep Oct (2) Nov (6) Dec Jan (1) Feb (1) Mar (3) Apr (2) May (2) Jun (2) Jul (18) Aug (13) Sep (7) Oct Nov Dec (2) Jan Feb (11) Mar Apr (4) May Jun (1) Jul (18) Aug (16) Sep (12) Oct (12) Nov (19) Dec (42) Jan (16) Feb (3) Mar (8) Apr (14) May (30) Jun (5) Jul (7) Aug (3) Sep (10) Oct (4) Nov (10) Dec (1) Jan (14) Feb (8) Mar (5) Apr (3) May (9) Jun (19) Jul Aug (27) Sep (5) Oct (18) Nov (12) Dec (8) Jan (5) Feb (8) Mar (20) Apr (22) May (28) Jun (9) Jul (6) Aug (46) Sep (40) Oct (15) Nov (8) Dec (34) Jan (20) Feb (15) Mar (18) Apr (20) May (3) Jun (13) Jul (10) Aug (19) Sep (8) Oct (31) Nov (26) Dec (13) Jan (13) Feb (4) Mar (14) Apr (28) May (19) Jun (7) Jul (1) Aug Sep (19) Oct (5) Nov (4) Dec (9) Jan (4) Feb (30) Mar Apr (5) May (1) Jun (1) Jul (3) Aug (2) Sep Oct Nov Dec
S M T W T F S

1

2

3

4
(1)
5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

Showing 1 results of 1

 [Reduce-algebra-developers] Matrices with noncommuting elements From: - 2010-01-04 14:02:40 ```Greetings, I need to write some reduce code to perform coalgebra/bialgebra calculations. I would like to be able to construct abstract (tensor) products such as Delta(A)_{n+1}=A_n*A+B_n*C in which * is a tensor product, A_n and B_n are unspecified non-commuting operators and A and C are 2x2 matrices. In practice A_n, B_n, ... depend on parameters. In particular I need to be able to do algebra with these, for example (A_n* A)(B_n*C)=(A_n B_n)*(AC) in which * is again a tensor procuct, and AC is the matrix product of A and C. I have tried various things, all of which fail or cause stack overflows. Do you have any suggestions? For example matrix An,Bn,Cn,Dn; A2:=mat((c,0),(0,b)); B2:=mat((0,0),(sqrt(ab),0)); C2:=mat((0,sqrt(ab)),(0,0)); D2:=mat((a,0),(0,0)); matrix up,down; up:=mat((1,0),(0,0)); down:=mat((0,0),(1,0)); operator An,Bn,Cn,Dn; noncom An,Bn,Cn,Dn; matrix M; load_package linalg; % coproducts in lowest dim'l rep DeltaA:=kronecker_product(A2,A2)+kronecker_product(B2,C2); DeltaB:=kronecker_product(A2,B2)+kronecker_product(B2,D2); % this fails operator dA; dA:=An*A2; for all i let for all j let dA(i)*dA(j)=(An(i)*An(j))*(A2*A2); dA(1)*dA(2)*up; Any suggestions would be welcome. Regards, Jeff Schmidt ```

Showing 1 results of 1