scipy/speed — Issue #1: One-array, two-array Laplace ?
Author: denis-bz
Created: 2012-02-14 10:54:06+00:00
State: open

"One-array" Laplace, new = av( 2 new, 2 old ), is of course different from "two-array", new = av( 4 old );
they cache differently, converge differently.
NumPy u[1:-1,1:-1] = ((u[2:,1:-1]+u[:-2,1:-1]) ... is two-array,
Cython u[i,j] = ((u[i+1, j] + u[i-1, j]) ... is one-array,
Fortran ?
How about printing u.mean() at the end, as a check that all the codes are doing the same thing ?
cheers
-- denis