pourquoi l algorithme d'Arnoldi s arrête exacetment en n iterations supposant A inversible.
pour ceux qui ont une idée svp au secours
l algorithme peut etre generalisé pour la resolution des systemes lineaires a second memebre multiple AX=B . X et B sont nxs
le produit scalaire est au sens de Frobenius ( <X,Y>_F=trace(transpose(X).Y)
cette generalisation s arrête aussi en n iterations !!! pourquoi ?!!