spielwiese
|
Last change
on this file since 564136 was
c6963ca,
checked in by Hans Schoenemann <hannes@…>, 14 years ago
|
|
pagelength removed
git-svn-id: file:///usr/local/Singular/svn/trunk@13176 2c84dea3-7e68-4137-9b89-c4e89433aadc
|
-
Property mode set to
100644
|
|
File size:
1.0 KB
|
| Line | |
|---|
| 1 | LIB "tst.lib"; |
|---|
| 2 | tst_init(); |
|---|
| 3 | // |
|---|
| 4 | // test script for lift command |
|---|
| 5 | // |
|---|
| 6 | ring r1 = 32003,(x,y,z),(c,dp); |
|---|
| 7 | r1; |
|---|
| 8 | "-------------------------------"; |
|---|
| 9 | poly s1=y7+x3+xyz+z2; |
|---|
| 10 | ideal i1=jacob(s1); |
|---|
| 11 | ideal i2=std(i1); |
|---|
| 12 | ideal i=3x2+yz,7y6+2x2y+5xz; |
|---|
| 13 | i; |
|---|
| 14 | i1; |
|---|
| 15 | i2; |
|---|
| 16 | s1; |
|---|
| 17 | matrix T=lift(i2,i); |
|---|
| 18 | T; |
|---|
| 19 | matrix(i)-matrix(i2)*T; |
|---|
| 20 | "-----------------------------"; |
|---|
| 21 | poly s2=x4+y3+z6+2xyz; |
|---|
| 22 | ideal i3=jacob(s2); |
|---|
| 23 | ideal i4=std(i3); |
|---|
| 24 | s2; |
|---|
| 25 | i3; |
|---|
| 26 | i4; |
|---|
| 27 | matrix T1=lift(i4,i); |
|---|
| 28 | T1; |
|---|
| 29 | matrix(i)-matrix(i4)*T1; |
|---|
| 30 | "-------------------------------"; |
|---|
| 31 | listvar(all); |
|---|
| 32 | kill r1; |
|---|
| 33 | //From: Trond St|len Gustavsen <stolen@math.uio.no> |
|---|
| 34 | |
|---|
| 35 | ring r = 0, (x(1..5)), ds; |
|---|
| 36 | matrix m[2][4]= x(1), x(2), x(3), x(4), x(2), x(3), x(4), x(5); |
|---|
| 37 | |
|---|
| 38 | ideal j = minor(m, 2); |
|---|
| 39 | ideal i = std(j); |
|---|
| 40 | |
|---|
| 41 | list ires = mres(i, 2); |
|---|
| 42 | def ires(1..2) = ires[1..2]; |
|---|
| 43 | kill ires; |
|---|
| 44 | matrix jaco = jacob(ires(1)); |
|---|
| 45 | qring s = i; |
|---|
| 46 | s; |
|---|
| 47 | matrix mat = matrix(fetch(r,ires(2))); |
|---|
| 48 | matrix imat = transpose(mat); |
|---|
| 49 | matrix jac = fetch(r, jaco); |
|---|
| 50 | |
|---|
| 51 | list ll = nres(module(imat),3); |
|---|
| 52 | def imatres(1..3) = ll[1..3]; |
|---|
| 53 | |
|---|
| 54 | matrix T = lift(imatres(2), module(jac)); |
|---|
| 55 | T; |
|---|
| 56 | kill r; |
|---|
| 57 | tst_status(1);$; |
|---|
Note: See
TracBrowser
for help on using the repository browser.
