Parallelization of the SIP algorithm

This work is devoted to the development of an algorithm for the parallelization of the SIP-solver for pentadiagonal matrices. It bases on the so-called block Jacobi algorithm. The parallelization of the matrix solver is done using the OpenMP standard. The performance of the new algorithm is illustrated with the numerical solution of a 2D Laplace's equation problem and a 2D Navier-Stokes equation problem (lid driven cavity).