We aim to look at a faster approach for solving the typical finite difference system of elliptical two dimensional PDEs. The method uses the Fast Fourier Transform (FFT) algorithm to quickly solve the finite difference system. We will also look into if this method can be applied to non-uniform and/or non-square meshes. Following looking into this approach we will see how the results of the former method compare to a approach using integral equations. The two PDEs we will apply the two aforementioned methods to are the Laplace and Poisson equations.
Laplace's Equation:
Poisson's Equation:
We will look at a multitude of examples and compare the accuracy between the two methods.
This repo contains the report and presentation for our APPM 4660: Numerical Anaylsis II final project. The code for the project can be found under the matlab directory.