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
​full size	​full size