Ask Your Question

PeterL's profile - activity

2016-06-20 18:58:05 -0500 received badge  Nice Question (source)
2016-02-05 05:33:14 -0500 received badge  Notable Question (source)
2016-02-05 05:33:14 -0500 received badge  Popular Question (source)
2015-04-30 05:06:12 -0500 received badge  Student (source)
2015-04-30 00:51:50 -0500 received badge  Scholar (source)
2015-04-30 00:51:40 -0500 commented answer bell_polynomial(n,k).coefficients() raises an error

All polynomials with k=0 are affected, and the error is not even consistent: b = bell_polynomial(3,0); b.parent() leads to: AttributeError: 'int' object has no attribute 'parent'. tmontail, is there a convention how the coefficients of the null polynomial are represented in Sage? [] or 0?

2015-04-29 11:45:27 -0500 asked a question bell_polynomial(n,k).coefficients() raises an error

Consider the Bell polynomials:

for n in (0..4):
    print [bell_polynomial(n,k) for k in (0..n)]

[1]
[0, x_1]
[0, x_2, x_1^2]
[0, x_3, 3*x_1*x_2, x_1^3]
[0, x_4, 3*x_2^2 + 4*x_1*x_3, 6*x_1^2*x_2, x_1^4]

Extracting the coefficients I expect this triangle:

[[1]]
[[0], [1]]
[[0], [1], [1]]
[[0], [1], [3], [1]]
[[0], [1], [3, 4], [6], [1]]

However this call (which looks very natural) does not work:

for n in (0..4):
    print [bell_polynomial(n,k).coefficients() for k in (0..n)]

It gives the error message: AttributeError: 'sage.rings.rational.Rational' object has no attribute 'coefficients'. However I nowhere give a sage.rings.rational.Rational object as an input. My input is a polynomial if I can trust the name "bell_polynomial(n,k)". So may I consider this as an internal error/bug?

Ralf Stephan showed me a workaround:

def polynomial_coefficients(p):
    if (isinstance(p, (int, Integer, Rational))):
        return [p]
    else:
        return p.coefficients()

This gives indeed the expected output. However it might be even better to represent the null polynomial by the empty list (to comply with the understanding that only nonzero coefficients are shown).

if (isinstance(p, (int, Integer, Rational))):
    if p == 0 return [] else return [p]

This leads to:

[[1]]
[[], [1]]
[[], [1], [1]]
[[], [1], [3], [1]]
[[], [1], [3, 4], [6], [1]]

Is this a workaround just to fit my possibly ill-guided expectations or is this a reasonable patch for the general case? Or is it correct that an object which is called a polynomial raises an error if asked for the coefficients?