LAUSR

This paper is devoted to the numerical computation of algebraic linear systems involving several matrix power functions; that is finding x solution to $$\sum _{\alpha \in \mathbb {R}}A^{\alpha }x=b$$ ∑ α ∈ R A α x = b . These systems will be referred to as Generalized Fractional Algebraic Linear Systems (GFALS). In this paper, we derive several gradient methods for solving these very computationally complex problems, which themselves require the solution to intermiediate standard Fractional Algebraic Linear Systems (FALS) $$A^{\alpha }x=b$$ A α x = b , with $$\alpha \in \mathbb {R_+}$$ α ∈ R + . The latter usually require the solution to many classical linear systems $$Ax=b$$ A x = b . We also show that in some cases, the solution to a GFALS problem can be obtained as the solution to a first-order hyperbolic system of conservation laws. We also discuss the connections between this PDE-approach and gradient-type methods. The convergence analysis is addressed and some numerical experiments are proposed to illustrate and compare the methods which are proposed in this paper.