Solves the Laplace equation u_xx + u_yy = 0 in a rectangular domain, by iteratively applying a stencil 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
The main algorithm is at SolverStencil.hs
The wrapper is at Main.hs
In the input file, greyscale gives the boundary conditions, blue specifies where to determine the solution.
|full size||full size|
Compiled with GHC 6.13.20100309. 400x400 image. 10000 iterations.
Running on a Intel i7 iMac. 2.8Ghz, 4 cores x 2 threads/core. 256k L1, 8MB L2, 8GB main memory.
Times stated include IO.