Un algorithme GMRES robuste au format tensor train

We consider the solution of linear systems with tensor product structure using a GMRES algorithm. In order to cope with the computational complexity in large dimension both in terms of floating point operations and memory requirement, our algorithm is based on low-rank tensor representation, namely the Tensor Train format. In a backward error analysis framework, we show how the tensor approximation affects the accuracy of the computed solution. With the bacwkward perspective, we investigate the situations where the $d$-dimensional problem to be solved results from the concatenation of a sequence of $(d-1)$-dimensional problems (like parametric linear operator or parametric right-hand side problems), we provide backward error bounds to relate the accuracy of the $d$-dimensional computed solution with the numerical quality of the sequence of $(d-1)$-dimensional solutions that can be extracted form it. This enables to prescribe convergence threshold when solving the $d$-dimensional problem that ensures the numerical quality of the $(d-1)$-dimensional solutions that will be extracted from the $d$-dimensional computed solution once the solver has converged. The above mentioned features are illustrated on a set of academic examples of varying dimensions and sizes.Nous considérons la résolution de systèmes linéaires avec une structure de produit tensoriel en utilisant un algorithme GMRES. Afin de faire face à la complexité de calcul en grande dimension, à la fois en termes d'opérations en virgule flottante et d'exigences de mémoire, notre algorithme est basé sur une représentation tensorielle à faible rang, à savoir le format Tensor Train. Dans un cadre d'analyse d'erreur inverse, nous montrons comment l'approximation tensorielle affecte la précision de la solution calculée. Dans une perspective d'erreur inverse, nous étudions les situations où le problème de dimension $d$ à résoudre résulte de la concaténation d'une séquence de problèmes de dimension $(d-1)$ (comme les problèmes d'opérateurs linéaires paramétriques ou de second membres paramétriques), nous fournissons des bornes d'erreur inverse pour relier la précision de la solution calculée de dimension $d$ à la qualité numérique de la séquence de solutions de dimension $(d-1)$ qui peut être extraite de celle-ci. Cela permet de prescrire un seuil de convergence lors de la résolution du problème à $d$ dimensions qui garantit la qualité numérique des solutions à $(d-1)$ dimensions qui seront extraites de la solution calculéeenà $d$ dimensions une fois que le solveur aura convergé. Les caractéristiques mentionnées ci-dessus sont illustrées sur un ensemble d'exemples académiques de dimensions et de tailles variables