Re: [sympy] SciPy talk

2014-05-21 Thread Rathmann
Thanks -- this one makes sense. (It is shuffles the multiplications into a balanced tree so that the top-level multiplication is the product of all the odds (=n) times the product of the evens.) Sympy itself uses prime swing. On Friday, May 16, 2014 9:50:57 AM UTC-7, Richard Fateman wrote:

Re: [sympy] SciPy talk

2014-05-16 Thread Richard Fateman
You are right! So sorry! I copied the wrong program into my note. Here is a correct version. (Again, Peter's article is fun to read...) (defun ff(n m) (if (= n m) n (*(ff n (+ m m))(ff (- n m) (+ m m) For this one, (ff 10 1) computes 10! and does the following multiplications. (ff

Re: [sympy] SciPy talk

2014-05-15 Thread Richard Fateman
Gee, this is getting long. Here's the deal on the improved factorial. Perhaps you just haven't gotten large enough, or your multiplications are too slow, even for small numbers. Or there is something else odd about your program that I don't see. Certainly n=20 is way too small to be

Re: [sympy] SciPy talk

2014-05-15 Thread Rathmann
Right, balancing the multiplication tree can be a win, for the reasons you give. But, I don't think your code snippet is doing any of that. If I run a version (attached) that prints out the multiplies it evaluates, (fact 100) gives: (2 . 1) (3 . 2) (4 . 6) (5 . 24) (6 . 120) (7 . 720)

Re: [sympy] SciPy talk

2014-05-09 Thread Rathmann
Well, I wouldn't have used such colorful language, but I share Prof. Fateman's skepticism on the utility of lgamma for computing factorials. I tried the expression as given, (see attachment, if you want to try yourself or point out errors) and it begins to diverge from the correct value at n=17

Re: [sympy] SciPy talk

2014-05-03 Thread Richard Fateman
On Tuesday, April 1, 2014 2:58:34 PM UTC-7, Ondřej Čertík wrote: On Mon, Mar 31, 2014 at 9:14 PM, Matthew Rocklin mroc...@gmail.comjavascript: wrote: http://www.evanmiller.org/mathematical-hacker.html I reference that blog post pretty often. I fully intend to reference it

Re: [sympy] SciPy talk

2014-05-02 Thread Aaron Meurer
FYI, my SciPy talk for SymPy was not accepted (it was accepted for the poster session). My talk on conda was accepted, as was the SymPy tutorial. Aaron Meurer On Tue, Apr 1, 2014 at 4:58 PM, Ondřej Čertík ondrej.cer...@gmail.com wrote: On Mon, Mar 31, 2014 at 9:14 PM, Matthew Rocklin

Re: [sympy] SciPy talk

2014-05-02 Thread Richard Fateman
I think your arguments are weak, though given the audience, perhaps they would be appealing. Here's what I think constitute good arguments for people to know about CAS. Maybe even sympy. 1. Scientists, mathematicians and programmers all have a rich language and context for discussing the

Re: [sympy] SciPy talk

2014-05-02 Thread Ondřej Čertík
Hi Richard, On Fri, May 2, 2014 at 5:32 PM, Richard Fateman fate...@gmail.com wrote: I think your arguments are weak, though given the audience, perhaps they would be appealing. Here's what I think constitute good arguments for people to know about CAS. Maybe even sympy. 1. Scientists,

Re: [sympy] SciPy talk

2014-05-02 Thread Aaron Meurer
On May 2, 2014, at 11:23 PM, Ondřej Čertík ondrej.cer...@gmail.com wrote: Hi Richard, On Fri, May 2, 2014 at 5:32 PM, Richard Fateman fate...@gmail.com wrote: I think your arguments are weak, though given the audience, perhaps they would be appealing. Here's what I think constitute good

Re: [sympy] SciPy talk

2014-04-01 Thread Ondřej Čertík
On Mon, Mar 31, 2014 at 9:14 PM, Matthew Rocklin mrock...@gmail.com wrote: http://www.evanmiller.org/mathematical-hacker.html I reference that blog post pretty often. I fully intend to reference it again in my talk (if it is accepted). The interesting thing about the Factorial / Gamma /

Re: [sympy] SciPy talk

2014-03-31 Thread Matthew Rocklin
I like that you emphasized the utility for numerics, I think that this is likely to be a selling point for the SciPy crowd. I think that it's correct to inform the reader about what symbolics are but I think that the first couple of sentences (which do this) could be stronger/more direct. Right

Re: [sympy] SciPy talk

2014-03-31 Thread Matthew Rocklin
Last minute (this morning) I decided to submit a talk based around some of my more esoteric programming projects, multiple dispatch, pattern matching, strategies. I drew a lot of inspiration from SymPy when I was thinking about these things so if accepted then I'll probably use SymPy heavily for

Re: [sympy] SciPy talk

2014-03-31 Thread Aaron Meurer
On Mon, Mar 31, 2014 at 11:32 AM, Matthew Rocklin mrock...@gmail.comwrote: I like that you emphasized the utility for numerics, I think that this is likely to be a selling point for the SciPy crowd. Yes, this was very intentional. I may need some help gathering up some nice motivating

Re: [sympy] SciPy talk

2014-03-31 Thread Aaron Meurer
I like this talk. It sounds similar to the talks you gave last year. I've really been sold on this approach in the assumptions system, a la my branch. In the old assumptions and the current new assumptions, to write down a fact, you have to do a lot of logic. You maybe don't need to know the

Re: [sympy] SciPy talk

2014-03-31 Thread Tim Lahey
On 31 Mar 2014, at 20:29, Aaron Meurer wrote: On Mon, Mar 31, 2014 at 11:32 AM, Matthew Rocklin mrock...@gmail.comwrote: I like that you emphasized the utility for numerics, I think that this is likely to be a selling point for the SciPy crowd. Yes, this was very intentional. I may

Re: [sympy] SciPy talk

2014-03-31 Thread Aaron Meurer
That's a good point. One of the nicest things about symbolics, when you can get it, is that it can make things drastically more efficient by doing mathematical simplifications. Evaluating integrals symbolically is a nice example of this (especially for SymPy, which has some pretty nice algorithms

Re: [sympy] SciPy talk

2014-03-31 Thread Matthew Rocklin
http://www.evanmiller.org/mathematical-hacker.html I reference that blog post pretty often. I fully intend to reference it again in my talk (if it is accepted). The interesting thing about the Factorial / Gamma / loggamma example is that to find the solution you need to find someone who knows

Re: [sympy] SciPy talk

2014-03-31 Thread Matthew Rocklin
And regarding assumptions I agree with everything that you said. On Mon, Mar 31, 2014 at 8:14 PM, Matthew Rocklin mrock...@gmail.com wrote: http://www.evanmiller.org/mathematical-hacker.html I reference that blog post pretty often. I fully intend to reference it again in my talk (if it is

[sympy] SciPy talk

2014-03-29 Thread Aaron Meurer
Here is the detailed abstract I have so far for a talk for SciPy. Any suggestions are welcome. The deadline is April 1 (probably 5 PM central or thereabouts). I roughly based it on the matplotlib talk from last year http://conference.scipy.org/scipy2013/presentation_detail.php?id=211. Symbolic