Alright, here's the pure pexpect based version for the Singular example:
sage: singular.lib('decodegb')
# as mentioned before x(1)...x(3) is really bad notation for Sage.
sage: r = singular.ring(3,'x(1..3)','dp')
---------------------------------------------------------------------------
Traceback (most recent call last)
...
RuntimeError: Singular error:
? error occurred in STDIN line 27: `if(defined((1..3)>0){kill (1..3;};`
? last reserved name was `defined`
skipping text from `;`
sage: r = singular.ring(3,'x(1),x(2),x(3)','dp')
---------------------------------------------------------------------------
RuntimeError Traceback (most recent call last)
# Thus we have use a 'trick'
sage: singular.eval('ring r=3,(x(1..3)),dp;')
'ring r=3,(x(1..3)),dp;'
sage: points = singular.pointsGen(3,1);
sage: points2 = singular.convPoints(points);
sage: p = singular.graspList(points2,1,11);
sage: p
[1]:
_[1,1]=0
_[2,1]=0
_[3,1]=0
[2]:
_[1,1]=0
_[2,1]=0
_[3,1]=1
[3]:
_[1,1]=0
_[2,1]=0
_[3,1]=-1
[4]:
_[1,1]=0
_[2,1]=1
_[3,1]=0
[5]:
_[1,1]=0
_[2,1]=1
_[3,1]=1
[6]:
_[1,1]=0
_[2,1]=1
_[3,1]=-1
[7]:
_[1,1]=0
_[2,1]=-1
_[3,1]=0
[8]:
_[1,1]=0
_[2,1]=-1
_[3,1]=1
[9]:
_[1,1]=0
_[2,1]=-1
_[3,1]=-1
[10]:
_[1,1]=1
_[2,1]=0
_[3,1]=0
[11]:
_[1,1]=1
_[2,1]=0
_[3,1]=1
sage: id = singular.vanishId(p);
sage: id
x(1)*x(2),
x(1)^2-x(1),
x(3)^3-x(3),
x(1)*x(3)^2-x(1)*x(3),
x(2)^3-x(2)
--
name: Martin Albrecht
_pgp: http://pgp.mit.edu:11371/pks/lookup?op=get&search=0x8EF0DC99
_otr: 47F43D1A 5D68C36F 468BAEBA 640E8856 D7951CCF
_www: http://www.informatik.uni-bremen.de/~malb
_jab: [email protected]
--
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org
