spielwiese
|
Last change
on this file since 6fb723 was
894057,
checked in by Oleksandr Motsak <motsak@…>, 13 years ago
|
|
ADD: Tests from online manual (res+stat on mamawutz): short ones Tst/Manual/s.lst
|
-
Property mode set to
100644
|
|
File size:
676 bytes
|
| Line | |
|---|
| 1 | LIB "tst.lib"; tst_init(); |
|---|
| 2 | LIB "jacobson.lib"; |
|---|
| 3 | ring r = 0,(x,d),Dp; |
|---|
| 4 | def R = nc_algebra(1,1); setring R; // the 1st Weyl algebra |
|---|
| 5 | matrix m[2][2] = d,x,0,d; print(m); |
|---|
| 6 | list J = jacobson(m); // returns a list with 3 entries |
|---|
| 7 | print(J[2]); // a Jacobson Form D for m |
|---|
| 8 | print(J[1]*m*J[3] - J[2]); // check that U*M*V = D |
|---|
| 9 | /* now, let us do the same for the shift algebra */ |
|---|
| 10 | ring r2 = 0,(x,s),Dp; |
|---|
| 11 | def R2 = nc_algebra(1,s); setring R2; // the 1st shift algebra |
|---|
| 12 | matrix m[2][2] = s,x,0,s; print(m); // matrix of the same for as above |
|---|
| 13 | list J = jacobson(m); |
|---|
| 14 | print(J[2]); // a Jacobson Form D, quite different from above |
|---|
| 15 | print(J[1]*m*J[3] - J[2]); // check that U*M*V = D |
|---|
| 16 | tst_status(1);$ |
|---|
Note: See
TracBrowser
for help on using the repository browser.
