Here is a more concise, but essentially equivalent, code snippet which
exhibits the problem. Just attach and run "prob()" Again, it takes
10-12 minutes to crash on my MacPro3,1 with 2 Quad-Core Xeon
Processors @ 2.8 Ghz, running 10.5.8, with Sage 4.3.2 self-compiled
with i686-apple-darwin9-gcc-4.0.1 (GCC) 4.0.1 (Apple Inc. build
5493).
Best,
Nathan
------- problem code begins -----
datum = ([0, 1, 0, 2, 0, 2], [([-1, -2, -1, 1, -1, 0], [1, 2, 1, 0, 1,
0], -1)])
coefficients, raw_natural_gluing_equations = datum
def eval_eqn( (a,b,c), z):
ans = c
for i, z in enumerate(z):
ans *= z**a[i] * (1 - z) ** b[i]
return ans
def first_order_equations():
n = len(coefficients)
a_vars = var(" ".join(["a%d" % i for i in range(n)]))
t = var('t')
base_eqns = raw_natural_gluing_equations
z = [ a * t ** e for a, e in zip(a_vars, coefficients) ]
eqns = [eval_eqn(eqn, z).subs(t=0) for eqn in base_eqns]
R = PolynomialRing(QQ, ['a%d' % i for i in range(n)])
poly_eqns = [ R(e.numerator()) for e in eqns]
def prob():
for i in xrange(1000000):
first_order_equations()
--
To post to this group, send an email to sage-devel@googlegroups.com
To unsubscribe from this group, send an email to
sage-devel+unsubscr...@googlegroups.com
For more options, visit this group at http://groups.google.com/group/sage-devel
URL: http://www.sagemath.org
To unsubscribe, reply using "remove me" as the subject.
