# source:git/Tst/Plural/doc-division2.tst@293722

spielwiese
Last change on this file since 293722 was 293722, checked in by Viktor Levandovskyy <levandov@…>, 14 years ago
*levandov: tests for new functions in doc, for plural git-svn-id: file:///usr/local/Singular/svn/trunk@12656 2c84dea3-7e68-4137-9b89-c4e89433aadc
• Property mode set to `100644`
File size: 1006 bytes
Line
1LIB "tst.lib";
2tst_init();
3LIB "dmod.lib";
4ring r = 0,(x,y),dp;
5poly f = x^3+xy;
6def S = Sannfs(f); setring S; // compute the annihilator of f^s
7LD; // is not a Groebner basis yet!
8poly f = imap(r,f);
9poly P = f*Dx-s*diff(f,x);
10ideal I = LD, f; // consider a bigger ideal
11module A = LD*gen(1), I*gen(2); module B = slimgb(A);
12vector V = [P,(s+1)^2]; NF(V,B); //0
13list LV = division(V,A);
14vector W = [1,s^2]; NF(W,B); //!=0
15list LW = division(W,A);
16matrix M = matrix(W); matrix N = matrix(A);
17matrix T = matrix(LW); matrix R = matrix(LW); matrix U  = matrix(LW);
18// test:
19transpose(U)*transpose(M) - transpose(T)*transpose(N) - transpose(R); // must be 0
20vector W2 = [s^2,s^2]; NF(W2,B); //!=0
21// matrix to matrix
22matrix M2 = 1,s^2,s^2,1;
23list LW2 = division(M2,A);
24matrix M = matrix(M2); matrix N = matrix(A);
25matrix T = matrix(LW2); matrix R = matrix(LW2); matrix U  = matrix(LW2);
26transpose(U)*transpose(M) - transpose(T)*transpose(N) - transpose(R); // must be 0
27tst_status(1);\$
Note: See TracBrowser for help on using the repository browser.