Calculating the Galois group of \(Y'=AY+B\), \(Y'=AY\) completely reducible (Q1864902)

In his previous work with \textit{M. Singer} [J. Pure Appl. Algebra 139, 3-23 (1999; Zbl 0931.12009)] the author provided an algorithm for calculating the Galois group of a system of the form \(Y' = AY + B\). But that algorithm demanded the computation of a certain set of factorizations for a scalar linear differential polynomial associated with the system. This set is very large and difficult to compute in general. In this paper, the author gives a new and more efficient algorithm for calculating the Galois group of the system replacing the factorization tasks by a single block-diagonal decomposition of the coefficient matrix satisfying certain properties.