Version 5 (modified by benl, 5 years ago) |
---|
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 |
Attachments
- pls-400x400-thumb.png (2.6 kB) - added by benl 5 years ago.
- pls-400x400-out-thumb.jpg (6.0 kB) - added by benl 5 years ago.