#5109: [with patch, positive review] add support for Bell polynomials in Sage
---------------------------+------------------------------------------------
Reporter: mhansen | Owner: mhansen
Type: enhancement | Status: new
Priority: minor | Milestone: sage-3.3
Component: combinatorics | Resolution:
Keywords: |
---------------------------+------------------------------------------------
Changes (by wdj):
* summary: [with patch, needs review] add support for Bell polynomials
in Sage => [with patch, positive review] add
support for Bell polynomials in Sage
Comment:
This applies cleanly to 3.3.alpha2 and passes sage -t. The examples also
agree with some examples given on
http://en.wikipedia.org/wiki/Bell_polynomials, as well as this agreement:
{{{
sage: stirling_number2(6,2) == bell_polynomial(6,2)(1,1,1,1,1)
True
sage: stirling_number2(6,4) == bell_polynomial(6,4)(1,1,1)
True
sage: stirling_number2(7,4) == bell_polynomial(7,4)(1,1,1,1)
True
}}}
I ran sage -testall and got this failure:
{{{
sage -t "devel/sage/sage/rings/polynomial/toy_d_basis.py"
**********************************************************************
File "/Volumes/G-DRIVE-
MINI/sagestuff/sage-3.3.alpha2/devel/sage/sage/rings/polynomial/toy_d_basis.py",
line 91:
sage: d_basis(I)
Expected:
[x + 170269749119, y + 2149906854, z + 735710619426, 282687803443]
Got:
[x + 170269749119, y + 2149906854, z + 170335012540, 282687803443]
**********************************************************************
}}}
Though I don't know what a d_basis is, I think it is an unrelated failure
so I'm giving this a positive review.
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/5109#comment:2>
Sage <http://sagemath.org/>
Sage - Open Source Mathematical Software: Building the Car Instead of
Reinventing the Wheel
--~--~---------~--~----~------------~-------~--~----~
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
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-trac?hl=en
-~----------~----~----~----~------~----~------~--~---
