Added a tridiagonal solver based on the distributed algorithm. Only periodic BCs are supported at the moment and works on a single rank only. In terms of performance, the distributed solver is around %5 slower than the thomas algorithm.
Next step is to introduce MPI communication to enable multiple ranks. As of now MPI is only mimicked by local data copies.
Then, support for Dirichlet and Neumann boundary conditions will be introduced.
Added a tridiagonal solver based on the distributed algorithm. Only periodic BCs are supported at the moment and works on a single rank only. In terms of performance, the distributed solver is around %5 slower than the thomas algorithm.
Next step is to introduce MPI communication to enable multiple ranks. As of now MPI is only mimicked by local data copies.
Then, support for Dirichlet and Neumann boundary conditions will be introduced.