[SciPy-Dev] Subversion scipy.stats irregular problem with source code example

[James Phillips]

On Fri, Oct 1, 2010 at 6:55 AM,  <josef.pktd@gmail.com> wrote:
> Which version of scipy are you using for testing?

Scipy 0.7.0 on Ubuntu Lucid Lynx with Python 2.6.5.


> I had to fix some python 2.5 incompatibilities (*args), but when I ran
> the script it seemed to get stuck after printing 4 or 5 iterations. I
> let it run for a few minutes but then killed the process.

At the top of diffev.py's solve() method is a loop that runs per
generation of the genetic algorithm, you might consider placing a
"print gen" statement there to see if it is running or stuck somehow.

The virtual server I'm using has 4 CPUs so I plan to run the
distribution fitting in parallel.  Here are current timing results
fitting both loc and scale, in parallel these would have taken 207
seconds (gamma time):

distribution: powerlaw took 69.5095369816 seconds to run
de parameters [ 1.31332767  0.55933617  2.45766383]
nnlf: 23.3610045864

distribution: beta took 47.2591450214 seconds to run
de parameters [ 0.70567183  0.6011756   0.607       2.41      ]
nnlf: 7.45846316363

distribution: gamma took 206.991358995 seconds to run
de parameters [  1.40727512e+02  -5.83784812e+00   5.55091421e-02]
nnlf: 26.7597491247

distribution: pareto took 86.9141609669 seconds to run
de parameters [  7.75987131e+13  -1.05369726e+14   1.05369726e+14]
nnlf: 35.2599606822


     James