Poisson Solvers for Two-Dimensional Fluid Simulation in a Box

In this project, we implement and benchmark four numerical solvers for the voxelized Poisson problem in two dimension, namely Gauss--Seidel Relaxation, successive over-relaxation, and conjugate gradient with and without an incomplete Cholesky preconditioner. Algorithms are tested on simple scenes with heat sources placed inside a MAC-discretized box to generate velocity fields.

Implementation is done in Java with JBLAS.

Note: there is a slight delay in Scene III for the Gauss-Seidel implementation. This is a recording/syncing problem and not an implementation issue.