Nico <nicolas.aunai@free.fr> wrote:

Je demande ici, car je n'ai pas la "culture" d'un numéricien, pourtant
je suis sûr que ce genre de probleme est "classique", et que poser
n=m=1024 ne me semble pas exagéré, il doit donc y avoir des biblio
etc... qui proposent des solutions.

Merci beaucoup

Nico

Tu peux regarder les méthodes pour les matrices "sparse" dans
matlab ou sinon scilab. Par exemple dans scilab, il y a la routine
"sparse", "lusolve" etc.
Dans la librairie "netlib" tu as des routines pour matrices sparse.
http://www.netlib.org/y12m/
notamment http://www.netlib.org/y12m/y12mfe.f

Dans lapack tu as des routines pour ça aussi:
http://www.netlib.org/lapack/lug/lapack_lug.html
notamment sgbsv dont la page man dit:

       SGBSV - compute the solution to a real system of linear equations
A * X
       = B, where A is a band matrix of order N with KL  subdiagonals
and  KU
       superdiagonals, and X and B are N-by-NRHS matrices

etc. etc.


--

Michel TALON