LAPACK

La bibliothèque numérique LAPACK est disponible à l'adresse

http://www.netlib.org/lapack.

Cette bibliothèque est devenue un standard pour la résolution de grands systèmes linéaires par des méthodes directes (élimination de Gauss, décomposition LU voir chap. 5). Elle a l'avantage d'être gratuite et exempte de fautes.

Cette bibliothèque existe depuis de nombreuses années, et a été écrite en FORTRAN pour des raisons historiques et d'efficacité. Il est néanmoins possible d'appeler les routines de cette bibliothèque depuis des programmes C, C++, Pascal.

L'utilisation d'une telle bibliothèque a l'avantage de vous éviter de reprogrammer des routines standard. De plus, bon nombre de constructeurs d'ordinateurs proposent des versions optimisées de cette librairie.

Par exemple, voici la routine DGTSV permettant de résoudre un système linéaire tridiagonal grâce à l'algorithme d'élimination de Gauss.