File: associatedPrimes2-test.m2

package info (click to toggle)
macaulay2 1.21%2Bds-3
  • links: PTS, VCS
  • area: main
  • in suites: bookworm
  • size: 133,096 kB
  • sloc: cpp: 110,377; ansic: 16,306; javascript: 4,193; makefile: 3,821; sh: 3,580; lisp: 764; yacc: 590; xml: 177; python: 140; perl: 114; lex: 65; awk: 3
file content (27 lines) | stat: -rw-r--r-- 748 bytes parent folder | download | duplicates (4)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
-- Date: Thu, 27 Feb 2003 15:40:55 -0500 (EST)
-- From: Greg Smith <ggsmith@cpw.math.columbia.edu>
-- To: Mike Stillman <mike@polygon.math.cornell.edu>,
--    Dan Grayson <dan@math.uiuc.edu>
-- cc: jsidman@math.berkeley.edu
-- Subject: Bug report.
-- 
-- 
-- There appears to a bug in the Assassinator code.
-- In particular, "associatedPrimes P" doesn't return all the
-- associated primes in the following example.
-- 
S = QQ[a,b,c,d,e];
I1 = ideal(a,b,c);
I2 = ideal(a,b,d);
I3 = ideal(a,e);
P = I1*I2*I3;
L1 = associatedPrimes P;
L2 = apply(associatedPrimes monomialIdeal P, J -> ideal J);
assert( #L1 == #L2 )
-- 
-- Naturally, I'm happy to observe that the monomial
-- ideal code does return the correct answer.
-- 
-- Cheers,
--   Greg.
--