New design of pivoted banded linear equation solver versus CLAPACK library routines

A new method of pivoting applicable to banded unsymmetric linear equation systems has been introduced. It limits the fill--in, preserves basic band structure and makes full use of allocated memory. Two new features of this method include multiple pivoting and the threshold factor use for pivoting decision. The stability or rather the solution precision has been investigated basing on the examples of classic ill conditioned equation system (Hilbert) and quasi--random system with still more difficult conditioning characteristic. Comprehensive test results have been presented in the easily readable graphic form.They prove the overall efficiency of new solvers and compare them favourably with the software from CLAPACK library.

---

Przedstawiono nową metodę wyboru elementów głównych użyteczną w przypadku pasmowych niesymetrycznych układów równań liniowych. Metoda ta ogranicza wypełnianie, zachowuje podstawową strukturę pasmową i w pełni wykorzustuje zaalokowaną pamięć. Dwie nowe właściwości tej metody to wielokrotne wybieranie elementów głównych oraz stosowanie współczynników progowych przy podejmowaniu decyzji o wyborze elementów głównych. Stabilność, a właściwie dokładność rozwiązań zbadano w oparciu o przykład klasycznego źle uwarunkowanego układu równań (Hilbert) i pseudo--przypadkowego układu, jeszcze gorzej uwarunkowanego. Rezultaty ekstensywnych testów przedstawiono w dobrze czytelnej postaci graficznej. Potwierdzają one efektywność nowego oprogramowania, zwłaszcza w porównaniu z oprogramowaniem z biblioteki CLAPACK.