On Sat, Jul 20, 2013 at 06:32:20AM -0700, Yuancheng Peng wrote:
> In diffusion/diffusion_python.py:
> 
> > +# Reference: http://continuum.io/blog/the-python-and-the-complied-python
> > +
> > +import numpy as np
> > +
> > +
> > +def diffusePurePython(u, tempU, iterNum):
> > +    """
> > +    Apply nested iteration for the Forward-Euler Approximation
> > +    """
> > +    mu = .1
> > +    row = u.shape[0]
> > +    col = u.shape[1]
> > +
> > +    for n in range(iterNum):
> > +        for i in range(1, row - 1):
> > +            for j in range(1, col - 1):
> 
> I have tried, but it ran even slower on my machine, so I didn't add it.

It does not matter if it is slower, that's the way people write it. if
pytrhan cannot do anything good with it, that's pythran's problem, the
benchmark should not hide it :-$