On the Moser- and super-reduction algorithms of systems of linear differential equations and their complexity

The notion of irreducible forms of systems of linear differential equations with formal power series coefficients as defined by Moser [Moser, J., 1960. The order of a singularity in Fuchs' theory. Math. Z. 379-398] and its generalisation, the super-irreducible forms introduced in Hilali and Wazner [Hilali, A., Wazner, A., 1987. Formes super-irreductibles des systemes differentiels lineaires. Numer. Math. 50, 429-449], are important concepts in the context of the symbolic resolution of systems of linear differential equations [Barkatou, M., 1997. An algorithm to compute the exponential part of a formal fundamental matrix solution of a linear differential system. Journal of App. Alg. in Eng. Comm. and Comp. 8 (1), 1-23; Pflugel, E., 1998. Resolution symbolique des systemes differentiels lineaires. Ph.D. Thesis, LMC-IMAG; Pflugel, E., 2000. Effective formal reduction of linear differential systems. Appl. Alg. Eng. Comm. Comp., 10 (2) 153-187]. In this paper, we reduce the task of computing a super-irreducible form to that of computing one or several Moser-irreducible forms, using a block-reduction algorithm. This algorithm works on the system directly without converting it to more general types of systems as needed in our previous paper [Barkatou, M., Pflugel, E., 2007. Computing super-irreducible forms of systems of linear differential equations via Moser-reduction: A new approach. In: Proceedings of ISSAC'07. ACM Press, Waterloo, Canada, pp. 1-8]. We perform a cost analysis of our algorithm in order to give the complexity of the super-reduction in terms of the dimension and the Poincare-rank of the input system. We compare our method with previous algorithms and show that, for systems of big size, the direct block-reduction method is more efficient.

[1]  A. Storjohann Algorithms for matrix canonical forms , 2000 .

[2]  Mark Giesbrecht,et al.  Computing Rational Forms of Integer Matrices , 2002, J. Symb. Comput..

[3]  Donald G. Babbitt,et al.  Formal reduction theory of meromorphic differential equations: a group theoretic view. , 1983 .

[4]  Volker Dietrich Zur Reduktion von linearen Differentialgleichungssystemen , 1978 .

[5]  Moulay A. Barkatou,et al.  An algorithm for computing a companion block diagonal form for a system of linear differential equations , 1993, Applicable Algebra in Engineering, Communication and Computing.

[6]  W. Wasow Asymptotic expansions for ordinary differential equations , 1965 .

[7]  Moulay A. Barkatou,et al.  On Rational Solutions of Systems of Linear Differential Equations , 1999, J. Symb. Comput..

[8]  Moulay A. Barkatou,et al.  Computing super-irreducible forms of systems of linear differential equations via moser-reduction: a new approach , 2007, ISSAC '07.

[9]  Pamela B. Lawhead,et al.  Super-irreducible form of linear differential systems , 1986 .

[10]  J. Moser,et al.  The order of a singularity in Fuchs' theory , 1959 .

[11]  Moulay A. Barkatou,et al.  An algorithm to compute the exponential part of a formal fundamental matrix solution of a linear differential system , 2009, Applicable Algebra in Engineering, Communication and Computing.

[12]  Joachim von zur Gathen,et al.  Modern Computer Algebra , 1998 .

[13]  Eckhard Pflügel,et al.  Effective Formal Reduction of Linear Differential Systems , 2000, Applicable Algebra in Engineering, Communication and Computing.

[14]  Thomas Cluzeau,et al.  Factorization of differential systems in characteristic p , 2003, ISSAC '03.

[15]  A. Hilali,et al.  Formes super-irréductibles des systèmes différentiels linéaires , 1986 .

[16]  Marie-Pierre Quere Stuchlik Algorithmique des faisceaux lineaires de matrices, application a la theorie des systemes lineaires et a la resolution d'equations algebro-differentielles , 1997 .

[17]  K. Chu The solution of the matrix equations AXB−CXD=E AND (YA−DZ,YC−BZ)=(E,F) , 1987 .

[18]  Eckhard Pflügel An improved estimate for the maximal growth order of solutions of linear differential systems , 2004 .