= 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] ||