= Laplace = Solves the Laplace equation `u_xx + u_yy = 0` in a rectangular domain, by iteratively applying a stencil function function. The first argument is the number of iterations to use. The second is a BMP file containing boundary conditions. {{{ repa-laplace 10000 pls-400x400.bmp pls-400x400-out.bmp }}} == Code == The main algorithm is at http://code.haskell.org/repa/repa-stable/repa-examples/Laplace/src/Solver.hs The wrapper is at http://code.haskell.org/repa/repa-stable/repa-examples/Laplace/src/Main.hs == Test Data == In the input file, greyscale gives the boundary conditions, blue is where to determine the solution. || pls-400x400.bmp || pls-400x400-out.bmp || || [[Image(Examples/Laplace:pls-400x400-thumb.png)]] || [[Image(Examples/Laplace:pls-400x400-out-thumb.jpg)]] || || [http://code.haskell.org/repa/wiki/images/pls-400x400.bmp full size] || [http://code.haskell.org/repa/wiki/images/pls-400x400-out.bmp full size] ||