C afé Calcul Formel et Équations SYM Manuel Bronstein INRIA Chercheur DR (deceased in June 2005) Montserrat Argente INRIA Assistant TR, part-time Evelyne Hubert INRIA Chercheur CR. Acting team leader Alban Quadrat INRIA Chercheur CR Thomas Cluzeau INRIA PostDoc Postdoc, until August 2006 Nicolas Le Roux UnivFr PhD ATER, co-supervision with the University of Limoges, defended in November Overall Objectives Overall Objectives

We develop computer algebra methods and softwares for solving functional equations, i.e., equations where the unknowns represent functions rather than numerical values. We study in particular linear and nonlinear differential and ( q)-difference equations, partial and ordinary. Solving is meant in the broad sense of extracting qualitative features of the solutions. Our goal is to foster the use of those symbolic methods in science and engineering by producing the programs and tools necessary to apply them to academic and industrial problems. Because symbolic computation is based on algebraic structures while our problems are analytic in nature, we develop the necessary theory for their algebraization.

Scientific Foundations Differential ideals and D-modules differential systems differential algebra formal integrability differential elimination algebraic analysis D-modules formal integrability involution holonomic systems control theory

Algorithms based on algebraic theories are developed to investigate the structure of the solution set of general differential systems. Different algebraic and geometric theories are the sources of our algorithms.

Differential elimination and completion

Formal integrabilityis the first problem that our algorithms address. The idea is to completea system of partial differential equation so as to be in a position to compute the Hilbert differential dimension polynomialor equivalently, its coefficients, the Cartan characters. Those provide an accurate measure of the arbitrariness that comes in the solution set (how many arbitrary functions of so many variables). Closely related is the problem of determining the initial conditionsthat can be freely chosen for having a well-posed problem (i.e. that lead to the existence and uniqueness of solutions). This is possible if we can compute all the differential relations up to a given order, meaning that we cannot obtain equations of lower order by combining the existing equations in the system. Such a system is called formally integrableand numerous algorithms for making systems of partial differential equations formally integrable have been developed using different approaches by E. Cartan , C. Riquier , M. Janet and D. Spencer .

Differential eliminationis the second problem that our algorithms deal with. One typically wants to determine what are the lowest differential equations that vanish on the solution set of a given differential system. The sense in which lowesthas to be understood is to be specified. It can first be order-wise, as it is of use in the formal integrability problem. But one can also wish to find differential equations in a subset of the variables, allowing the model to be reduced.

The radical differential ideal generated by a set of differential polynomials , i.e. the left-hand side of differential equations where the right-hand side is zero, is the largest set of differential polynomials that vanish on the solution set of . This is the object that our algorithms manipulate and for which we compute adequate representations in order to answer the above questions.

In the nonlinear case the best we can hope for is to have information outside of some hypersurface. Actually, the radical differential ideal can be decomposed into components on which the answers to formal integrability and eliminations are different. For each component the characteristic setdelivers the information about the singular hypersurface together with the quasi-generating set and membership test.

Triangulation-decomposition algorithms perform the task of computing a characteristic set for all the components of the radical differential ideal generated by a finite set of differential polynomials. References for those algorithms are the book chapters written by E. Hubert  , . They are based on the differential algebra developed by Ritt and Kolchin .

The objectives for future research in the branch of triangulation-decomposition is the improvement of the algorithms, the development of alternative approaches to certain classes of differential systems and the study of the intrinsic complexity of differential systems.

Another problem, specific to the nonlinear case, is the understanding and algorithmic classification of the different behaviours of interference of the locus of one component on the locus of another. The problem becomes clear in the specific case of radical differential ideals generated by a single differential polynomial. One then wishes to understand the behaviour of non singular solutions in the vicinity of singular solutions. Only the case of the first order differential polynomial equations is clear. Singular solutions are either the envelope or the limit case of the non singular solutions and the classification is algorithmic , .

Algebraic analysis

Many systems coming from mathematical physics, applied mathematics and engineering sciences can be described by means of systems of ordinary or partial differential equations, difference equations, differential time-delay equations...  In the case of a linear system , the system can be defined by means of matrices with entries in non-commutative algebras of functional operators such as differential operators, shift operators, time-delay operators...  An important class of operators algebras is called Ore algebrasand was introduced in .

The methods of algebraic analysis(algebraic D-modules) , , give a way to intrinsically study a linear functional system by means of its associated finitely presented left module over an Ore algebra. Thanks to the works of B. Malgrange, V. Palamodov, J. Bernstein and M. Kashiwara, algebraic analysis yields new results and information about the algebraic and analytic properties of linear functional systems, their solutions and associated geometric invariants (e.g., characteristic varieties, dimensions).

Algebraic analysis is becoming algorithmic thanks to the recent development of Gröbner bases over Ore algebras (e.g., the Maple package Ore -algebradeveloped by F. Chyzak , Singular: Pluraldeveloped by G.-M. Greuel, V. Levandovskyy and H. Schönemann ) and involutive bases in Ore algebras (e.g., the Maple package JanetOredeveloped by D. Robertz ), enabling the implementation of efficient algorithms for making systems formally integrable, computing special closed-form solutions , or studying properties of linear functional systems by means of the algebraic properties of the underlying modules.

Within constructive algebraic analysis, it is now possible to algebraically study systems of linear functional equations, and our objectives in that field are: (i) to develop and implement efficient algorithms for computing the polynomial and rational solutions of such systems and, further, for factoring and decomposing their associated D-modules; (ii) to study the links between the algebraic and analytic properties of such systems (since the algorithmic determination of the algebraic properties yields information about the analytic properties); (iii) to apply the above algorithms to study problems coming from mathematical physics and control theory.

Groups of transformations symmetry nonlinear differential systems linear systems of partial differential equations differential invariants differential Galois group closed-form solutions variational equations dynamical systems Hamiltonian mechanics formal integrability Lie groups of transformations

Though not a major subject of expertise, the topic is at the crossroads of the algorithmic themes developed in the team.

The Lie group, or symmetry group, of a differential system is the (biggest) group of point transformations leaving the solution set invariant.

Besides the group structure, a Lie group has the structure of a differentiable manifold. This double structure allows to concentrate on studying the tangent space at the origin, the Lie algebra.

The Lie group and the Lie algebra thus capture the geometry of a differential system. This geometric knowledge is exploited to solvenonlinear differential systems.

The Lie algebra is described by the solution set of a system of linear partial differential equations, whose determination is algorithmic. The dimension of the solution space of that linear differential system is the dimension of the Lie group and can be determined by the tools described in Section . Explicit subalgebras can be determined thanks to the methods developed within the context of Section .

For a given group of transformations on a set of independent and dependent variables there exist invariant derivations and a finite set of differential invariants that generate all the differential invariants . This forms an intrinsic frame for expressing any differential system invariant under this group action.

This line of ideas took a pragmatic shape for computation with the general method of M. Fels and P. Olver for computing the generating set of invariants. The differential algebra they consider has features that go beyond classical differential algebra. We are engaged in investigating the algebraic and algorithmic aspects of the subject.

Galois groups of linear functional equations

Differential Galois theory, developed first by Picard and Vessiot, then algebraically by Kolchin, associates a linear algebraic group to a linear ordinary differential equation or system. Many properties of its solutions, in particular the existence of closed-form solutions, are then equivalent to group-theoretic properties of the associated Galois group . By developing algorithms that, given a differential equation, test such properties, Kovacic and Singer have made the existence of closed-form solutions decidable in the case of equations with polynomial coefficients. Furthermore, a generalization of differential Galois theory to linear ordinary difference equations has yielded an algorithm for computing their closed-form solutions . Those algorithms are however difficult to apply in practice (except for equations of order two) so that many algorithmic improvements have been published in the past 20 years. Our objectives in this field are to improve the efficiency of the basic algorithms and to produce complete implementations, as well as to generalize them and their building blocks to linear partial differential and difference equations.

An exciting application of differential Galois theory to dynamical systems is the Morales-Ramis theory, which arose as a development of the Kovalevskaya-Painlevé analysis and Ziglin's integrability theory , . By connecting the existence of first integrals with branching of solutions as functions of complex time to a property of the differential Galois group of a variational equation, it yields an effective method of proving non-integrability and detecting possible integrability of dynamical systems. Consider the system of holomorphic differential equations defined on a complex n-dimensional manifold M. If ( t)is a non-equilibrium solution of ( ), then the maximal analytic continuation of ( t)defines a Riemann surface with tas a local coordinate. Together with ( ) we consider its variational equations (VEs) restricted to , i.e. We can decrease the order of that system by considering the induced system on the normal bundle of : where is the projection. The system of s= n-1equations obtained in this way yields the so-called normal variational equations (NVEs). Their monodromy group is the image of the fundamental group 1( , t0)of obtained in the process of continuation of the local solutions defined in a neighbourhood of t0along closed paths with base point t0. A non-constant rational function f( z)of svariables z= ( z1, ..., zs)is called an integral (or invariant) of the monodromy group if f( g· z) = f( z)for all .

In his two fundamental papers , , Ziglin showed that if ( ) possesses a meromorphic first integral, then has a rational first integral. Ziglin found a necessary condition for the existence of a maximal number of first integrals (without involutivity property) for analytic Hamiltonian systems, when n= 2 m, in the language of the monodromy group. Namely, let us assume that there exists a non-resonant element . If the Hamiltonian system with mdegrees of freedom has mmeromorphic first integrals F1= H, ..., Fm, which are functionally independent in a connected neighbourhood of , then any other monodromy matrix transforms eigenvectors of ginto its eigenvectors.

There is a problem however in making that theory algorithmic: the monodromy group is known only for a few differential equations. To overcome that problem, Morales-Ruiz and Ramis recently generalized Ziglin's approach by replacing the monodromy group by the differential Galois group of the NVEs. They formulated a new criterion of non-integrability for Hamiltonian systems in terms of the properties of the connected identity component of : if a Hamiltonian system is meromorphically integrable in the Liouville sense in a neighbourhood of the analytic curve , then the identity component of the differential Galois group of NVEs associated with is Abelian.

Since always contains , the Morales-Ramis non-integrability theorem always yields stronger necessary conditions than the Ziglin criterion.

When applying the Morales-Ramis criterion, our first step is to find a non-equilibrium particular solution, which often lies on an invariant submanifold. Next, we calculate the corresponding VEs and NVEs. If we know that our Hamiltonian system possesses kfirst integrals in involution, then we can consider VEs on one of their common levels, and the order of NVEs is equal to s= 2( m- k) , . In the last step we have to check if the identity component of is Abelian, a task where the tools of algorithmic Galois theory (such has the Kovacic algorithm) become useful. In practice, we often check only whether that component is solvable (which is equivalent to check whether the NVEs have Liouvillian solutions), because the system is not integrable when that component is not solvable.

Our main objectives in that field are: (i) to apply the Morales-Ramis theory to various dynamical systems occurring in mechanics and astronomy; (ii) to develop algorithms that carry out effectively all the steps of that theory; (iii) to extend it by making use of non-homogeneous variational equations; (iv) to generalize it to various non-Hamiltonian systems, e.g. for systems with certain tensor invariants; (v) to formulate theorems about partial integrability of dynamical systems and about real integrability (for real dynamical systems) in the framework of the Morales-Ramis theory;

Mathematical web services Computer algebra formula databases deductive databases communication Web services OpenMath MathML

The general theme of this aspect of our work is to develop tools that make it possible to share mathematical knowledge or algorithms between different software systems running at arbitrary locations on the web.

Most computer algebra systems deal with a lot of non algorithmic knowledge, represented directly in their source code. Typical examples are the values of particular integrals or sums. A very natural idea is to group this knowledge into a database. Unfortunately, common database systems are not capable to support the kind of mathematical manipulations that are needed for an efficient retrieval (doing pattern-matching, taking into account commutativity, etc.). The design and implementation of a suitable database raise some interesting problems at the frontier of computer algebra. We are currently developing a prototype for such a database that is capable of doing some deductions. Part of our prototype could be applied to the general problem of searching through mathematical texts, a problem that we plan to address in the near future.

The computer algebra community recognized more than ten years ago that in order to share knowledge such as the above database on the web, it was first necessary to develop a standard for communicating mathematical objects (via interprocess communication, e-mail, archiving in databases). We actively participated in the definition of such a standard, OpenMath (partly in the course of a European project). We were also involved in the definition of MathML by the World Wide Web Consortium. The availability of these two standards is the first step needed to develop rich mathematical services and new architectures for computer algebra and scientific computation in general enabling a transparent and dynamic access to mathematical components. We are now working towards this goal by experimenting with our mathematical software and emerging technologies (Web Services) and participating in the further development of OpenMath.

Application Domains Panorama

We have applied our research on linear functional systems to mathematical physics (e.g., computation of quadratic first integral or conservation laws, computation of block-diagonal forms of systems appearing in mathematical physics) and to control theory (e.g., computation of parametrizations or autonomous elements, computation of block-diagonal forms of many systems appearing in control theory, stabilization of linear control systems).

We have also applied our algorithms and programs for computing differential Galois groups to determine necessary conditions for integrability in mechanical modeling and astronomy.

Other applications include biology, where our algorithms for solving recurrence equations have been applied to the steady-state equations of nonhomogeneous Markov chains modeling the evolution of microsatellites in genomes.

Software Maple package diffalg Evelyne Hubert Nonlinear differential systems differential elimination analysis of singular solutions differential algebra triangulation-decomposition algorithms

The diffalgis a library that has been part of the commercial release of Maple since Maple V.5 (with an initial version by F. Boulier) and has been developed up to Maple 7 by E. Hubert. The library implements a triangulation-decomposition algorithm for polynomially nonlinear systems and tools for the analysis of singular solutions. Newer releases are available on the web ( http://www.inria.fr/cafe/Evelyne.Hubert/diffalg) together with an extensive set of examples of applications and an efficiency test suite. Recent developments include optimized treatment of parameters, improvements for higher degree equations and generalization to non-commuting derivations.

Library OreModules Alban Quadrat correspondent Daniel Robertz Constructive algebraic analysis constructive homological algebra Ore algebras Gröbner bases linear functional systems control theory mathematical physics

The OreModulespackage of Ore -algebra (Ore -algebra is a part of the commercial release of Maple) is dedicated to the algebraic study of linear systems of functional equations defined over some Ore algebras (e.g., ordinary or partial differential equations, time-delay equations, difference equations) and their applications in mathematical physics (e.g., research of potentials, computations of field equations and conservation laws) and in linear control theory (e.g., parameterizability, flatness, autonomous elements, first integrals of motion, equivalences, controllability, observability, input-output behaviour).

The main novelty of OreModulesis to combine the recent developments of the Gröbner bases over some Ore algebras (non-commutative polynomial rings) with new algorithms of algebraic analysis in order to effectively check classical properties of module theory (e.g., existence of a non-trivial torsion submodule, torsion-freeness, reflexiveness, projectiveness, stably freeness, freeness), give their system-theoretical interpretations (existence of autonomous elements or successive parametrizations, existence of minimal/injective parametrizations, Bézout equations or generalized inverses) and compute important tools of homological algebra (e.g., (minimal) free resolutions, splitting, extension functor, projective dimension, Hilbert power series).

The main advantage of the language of homological algebra used in carries over the implementations in OreModules: up to the choice of the domain of functional operators which occurs in a given system, all algorithms are stated and implemented in sufficient generality such that linear systems defined over the Ore algebras developed in the Maple package of Ore -algebraare covered at the same time.

A library of more than 30 examples coming from engineering sciences and mathematical physics (e.g., a two pendulum mounted on a car, a wind tunnel model, stirred tank models, a two reflector antenna, an electric transmission line, linearized Einstein equations, Maxwell equations, linear elasticity, Lie-Poisson structures) illustrates the main functions and applications of OreModules.

Library Stafford Alban Quadrat correspondent Daniel Robertz Stafford theorems bases of free modules over the Weyl algebras generators of ideals over the Weyl algebras computation of flat outputs block-diagonalization and block-triangularization problems

The Staffordpackage of OreModulescontains an implementation of constructive versions of the J. T. Stafford's famous but tricky theorem stating that every ideal over the Weyl algebras and of differential operators with polynomial and rational coefficients can be generated by two generators. Based on this implementation and on algorithmic results recently obtained by the authors of this package, two algorithms have been implemented which compute bases of free modules over the Weyl algebras and .

The development of the Staffordpackage was mainly motivated by the wish to compute flat outputs of flat multidimensional linear systems with varying coefficients, a problem appearing in control theory. However, it has recently played an important role in the decomposition problem of linear systems of partial differential equations with polynomial or rational coefficients. For more details, see Section  .

Library Morphisms Thomas Cluzeau Alban Quadrat correspondent Factorization and decomposition of linear functional systems constructive homological algebra quadratic first integrals quadratic conservation laws algebraic analysis

The forthcoming Morphismspackage of OreModuleswas developed in order to constructively handle some homological tools such as computations of morphisms between two finitely presented modules over Ore algebras, compute kernel, coimage, image and cokernel of morphisms and projectors of the ring of endomorphisms. Using the packages Stafford(see Section  ) and Quillen-Suslin(see Section  ), these results allow us to effectively compute factorizations as well as finding decompositions of linear systems over Ore algebras. Hence, using the package Morphisms, we can highly simplify the form of linear functional systems appearing in mathematical physics and control theory.

In terms of module theory, the Morphismspackage gives some methods to test if two modules are isomorphic, if a given module contains a submodule (reducible modules) or if it can be written as the direct sum of two submodules. Applications of the Morphismspackage to mathematical physics and control theory are illustrated in a library of examples (e.g., computation of Galois symmetries of physical systems, quadratic first integral and conservation laws, decomposition of linear systems coming from mathematical physics and control theory, equivalence of systems).

Library QuillenSuslin Anna Fabiańska Alban Quadrat correspondent Quillen-Suslin theorem flat multidimensional systems Lin-Bose conjecture (weakly) doubly factorization block-diagonalization and block-triangularization problems

The forthcoming QuillenSuslinpackage of OreModules, developed by A. Fabiańska (University of Aachen) with the help of A. Quadrat, contains an implementation of the famous Quillen-Suslin theorem. In particular, this implementation allows us to compute bases of free modules over a commutative polynomial ring with coefficients in a field (mainly ) and in a principal ideal domain (mainly ). The development of this package was led by different applications of the Quillen-Suslin theorem in multidimensional systems theory (e.g., computation of flat outputs of flat multidimensional systems, equivalence problems, Lin-Bose conjecture, computation of (weakly) doubly coprime factorizations) and signal processing (e.g., parametrization of synthesis filters). Finally, it has recently played a crucial role in the decomposition problem of linear systems of partial differential equations with constant coefficients (see Section  ).

Library libaldor Manuel Bronstein Aldor standard data structures arithmetic

The libaldorlibrary, under development in the project for several years, became the standard Aldorlibrary, bundled with the compiler distribution since 2001.

Library Algebra Manuel Bronstein Marc Moreno-Maza University of Western Ontario Computer algebra linear algebra commutative algebra polynomials

The Algebralibrary, written in collaboration with Marc Moreno-Maza, is now also bundled with the Aldorcompiler, starting with version 1.0.2, which has been released in May 2004. It is intended to become a standard core for computer algebra applications written in Aldor.

Library i t Manuel Bronstein Computer algebra differential equations difference equations systems

The i tlibrarycontains our algorithms for solving functional equations. The stand-alone solver powers the the web service at http://www-sop.inria.fr/cafe/Manuel.Bronstein/sumit/bernina_demo.htmlthat has solved third order equations as well since the end of 2004.

New Results Fast algorithms for linear differential and difference equations Thomas Cluzeau Alin Bostan ALGO project Frédéric Chyzak ALGO project Bruno Salvy ALGO project Linear differential and difference equations polynomial and rational solutions complexity binary splitting baby step-giant step indefinite and definite summation

The search for rational (and thus polynomial) solutions of linear differential and difference equations lies at the heart of several important algorithms in computer algebra.

In , we investigate the problem of computing polynomial solutions of linear differential equations. Even for detecting the non-existence of non-zero polynomial solutions of a linear differential equation, there exists no general algorithm in polynomial complexity. However, high degree polynomial solutions are controlled by the differential equation. This allows to represent them with a compact data-structure: a linear recurrence and initial conditions. This compact representation enables us to compute polynomial solutions in quasi-optimal complexity (that is, optimal complexity modulo logarithm terms) using binary splitting and baby step/giant step techniques. The timings of our implementations in Maple and Magma confirm the complexity estimates in practice.

We then address the problem of computing rational solutions. The compact representation of polynomial solutions of linear differential equations also allows us to evaluate, in quasi-optimal complexity, their value and that of their derivatives at a rational or algebraic point. We can also perform divisions by high powers of a monomial which lead directly to an algorithm in quasi-optimal complexity to compute rational solutions of linear differential equations.

This year, we have continued our joint work by looking at the case of difference equations. The calculation of polynomial solutions works in the same way except that the coefficients of these polynomials satisfy a linear recurrence equation in the binomial basis instead of the monomial basis. The generalization of our techniques for computing rational solutions is a little bit more difficult because of the appearance of a new phenomenon called the dispersion. However, binary splitting and baby step/giant step techniques together with the compact representation of the numerator allow us to obtain an algorithm in quasi-optimal complexity. As a direct consequence, we obtained fast versions of Abramov's algorithm computing rational solutions of linear recurrence equations, Gosper's algorithm for indefinite summation, and Zeilberger's algorithm for definite summation. Our algorithms have been implemented in Maple. This work has been published in .

Formal solutions of partial differential systems at a singular point Nicolas Le Roux Moulay Barkatou University of Limoges Evelyne Hubert Pfaffian systems normal crossings singularities regular and irregular singularity

Completely integrable systems are well studied. They are encountered for instance in isomonodromic deformations or in differential geometry. We consider here completely integrable systems with normal crossings singularities. These systems can be rewritten locally in the form The Aiare m× mmatrices with entries in the ring and the piare non negative integers. Yis a m×1vector of unknowns with coordinates in an extension of .

We investigate algorithms for the computation of formal solutions of completely integrable systems with normal crossing singularities. We expect to extend the known algorithms for ordinary differential systems.

As in the ordinary case, there are three kinds of singularities at the origin. Singularities of the first kindare those for which pi= 0for all i. Every system of the first kind is equivalent, by a transformation Y= TZwith , to a system where the matrices have constant entries , . These constant matrices are the monodromy of the system. Singularities are regularif the system is equivalent to one of the first kind. All other singularities are called irregular. There are algorithms to decide whether a singularity is regular that is based on a criterion by P. Deligne and A. van den Essen . Yet this is not constructive as the algorithm does not produce the transformation. We obtained a dual version of the criterion by A. van den Essen. Our criterion is based on the duality between increasing and decreasing versions of Levelt algorithm and characterizes regularity of a singularity by means of stabilization of a decreasing sequence of lattices at the mth stage. Stabilization at the mth stage, not appearing in van den Essen criterion, is of main interest for algorithmic purposes.

For singularities of the first kind the solution space can be described in terms of the monodromy matrices. More precisely, one can find a basis of solutions , which form is where with det( T) 0and the C( i)s are m× mconstant matrices commuting with one another. We obtained from proofs in , methods for computing such a basis of solutions. We are going to implement these methods in Maple.

The description of the solution space at regular singularities is then obtained if we have the transformation that reduces the system to a system of the first kind.

For irregular singularities, H. Charrière was the first to describe the solution space in the case where the system only involves two independent variables . She showed that Turritin's form for linear differential system can be carried to this case. Her proof uses the lattice approach initiated by R. Gérard & A.H.M. Levelt , . Next the result was extended to the general case independently by A.H.M. Levelt & A. van den Essen and H. Charrière & R. Gérard .

Our first goal is to provide an algorithm for reducing the order of the singulatrity at the origin - that is the length of ( p 1, ..., p n)- by linear transformations. Such a reduction process would allow us to decide wether or not the singularity is regular. We have investigated the possibility to extend Moser's algorithm , and Levelt's algorithm that work in the ordinary case. We clarified the relationship between those two algorithms and gave a complexity study of both of them . The presentation with lattices by Levelt seems easier to extend. We have shown that it can be extended to the case n= 2 . It is also implemented in Maple. We have used parts of the implementation of Moser's algorithm by M. Barkatou . In the case n= 2, we are building an algorithm for computing a basis of formal solution for systems with irregular singularities, relying on the rank reduction algorithm previously proposed. This is to be written in a forthcoming paper.

For n>2, the difficulties we encounter are twofold. Firstly, the systems considered are not invariant under transformation Y= TZ. Secondly, there is no warranty that the lattices are free. The same problem arose in , .

Factorization and decomposition of linear functional systems Cluzeau Thomas Alban Quadrat Factorization and decomposition of linear functional systems constructive homological algebra quadratic first integrals quadratic conservation laws algebraic analysis

The factorization and decomposition problems of systems of linear ordinary differential or difference equations have long been studied in the mathematical and symbolic computation literatures , , , , , and in the CAFÉ project , . The main novelty of the approach developed in the CAFÉ project is to recast these problems within the algebraic analysis framework and to extend them to linear functional systems defined over Ore algebras. Except from the pioneering works of M. F. Singer and his collaborators, this approach is new in the case of finite-dimensional linear systems over Ore algebras and, to our knowledge, is the first one for general linear functional systems (determined/under-/over-determined). The second novelty is to use the expertise developed in the CAFÉ project on efficient implementations for linear systems over univariate Ore algebras (e.g., see Section  ). The last novelty is to apply our new results to different systems coming from mathematical physics and control theory.

Following ideas of M. F. Singer , , in the publications and , we have shown how to compute the morphisms from a left D-module M, finitely presented by a matrix Rwith entries in an Ore algebra D, to a left D-module Npresented by a matrix S. If Dis a commutative ring or Mand Nare two finite-dimensional k-vector spaces (e.g., case of ordinary differential systems or integrable connections), then we can effectively compute the set h o m D( M, N)of morphisms from Mto N. If the previous hypotheses are not fulfilled, then we can only compute the set of morphisms of a given order. A morphism from Mto Ndefines a transformation sending a solution of the system Sz= 0into a solution of Ry= 0. When S= R, the ring e n d D( M)of the endomorphisms of Mcorresponds to the Galois symmetries'' of the system Ry= 0. Moreover, we have explicitly characterized the kernel, coimage, image and cokernel of a morphism from Mto N, a fact that can be used to check the equivalence of the systems Ry= 0and Sz= 0.

We have proved that the existence of a non-injective endomorphism of a left D-module M, finitely presented by a matrix Rwith entries in an Ore algebra D, corresponded to a factorization of the form R= R1R2, where R1and R2are two matrices with entries in D. Using this result, we have shown how the integration of the system Ry= 0could be done by means of a cascade of integrations. Moreover, we have described how to compute the projectors of e n d D( M)and proved that they allowed us to decompose the system Ry= 0into two decoupled systems Q1y1= 0and Q2y2= 0, where Q1and Q2were two matrices with entries in D. Hence, the integration of the system Ry= 0becomes equivalent to the integrations of the two independent systems Q1y1= 0and Q2y2= 0. Then, under certain conditions on the projectors (freeness), we proved that the system Ry= 0was equivalent to a block diagonal system. In particular, these conditions always hold in the case of a univariate Ore algebra over a field of coefficients (i.e., ordinary differential/difference systems over the field of rational functions) and in the case of a multivariate commutative Ore algebras due to the Quillen-Suslin theorem (e.g., linear system of partial differential equations with constant coefficients) (see Section  ). Moreover, if some rank conditions on the projectors are fulfilled, then, using a result due to J. T. Stafford, we prove that a similar result also holds for the Weyl algebras An( k)and Bn( k)over a field kof characteristic 0 (i.e., linear system of partial differential equations with polynomial/rational coefficients over k) (see Section  ). Using the recent implementations of the Quillen-Suslin and Stafford theorems in the packages Staffordand QuillenSuslin, we give a constructive way to compute the block diagonal decomposition of Ry= 0when it exists. For all the previous results, no condition on the system Ry= 0is required such as D-finite, determined, underdetermined, overdetermined, i.e., this approach handles general linear systems over an Ore algebra. In the case of a finite-dimensional determined system over an Ore algebra with rational coefficients, we find again the results obtained in M. Wu's PhD thesis .

The different algorithms have been implemented in the package Morphisms(see Section  ) of OreModules(see Section  ). This package is available with a library of examples illustrating its main features and applications in mathematical physics (e.g., symmetries of the Maxwell/Euler/Saint-Venant/Einstein equations, quadratic first integrals of motion or conservation laws (e.g., linear elasticity, electromagnetism, hydrodynamics), decomposition of systems of partial differential equations (e.g., transmission lines, Cauchy-Riemann or Beltrami equations, acoustic or electromagnetic waves, two-dimensional rotational isentropic flow, incompressible fluid in rotation with a small fluid velocity), Monge parametrizations) and control theory (e.g., decomposition of differential time-delay systems appearing in applications such as a wind tunnel model, a flexible rod, a stirred tank, different tank models, a network model).

Linear systems over Ore algebras and applications Alban Quadrat Daniel Robertz Constructive algebraic analysis constructive homological algebra Ore algebras Gröbner bases linear functional systems control theory mathematical physics

In the publications , we study the structural properties of under-determined linear systems over Ore algebras. Using the recent development of Gröbner and Janet bases over Ore algebras, we show how to make effective some important concepts of homological algebra (e.g., free resolutions, split exact sequences, duality, extension functor, different dimensions, torsion-free degrees, grades). Using these results, we obtain some effective algorithms which check the different structural properties of under-determined linear systems over Ore algebras by means of the algebraic properties of the underlying modules (e.g., torsion/torsion-free/reflexive/projective/stably free/free modules). In particular, we explain why these properties are related to the possibility to successively parameterize all solutions of a system and its parametrizations. Moreover, we show how these properties generalize the well-known concepts of primeness, developed in the literature of multidimensional systems, to systems with varying coefficients. Then, using a dictionary between the structural properties of under-determined linear systems over Ore algebras and concepts of linear control theory, we show that the previous algorithms allow us to effectively check whether or not a multidimensional linear control system over certain Ore algebras is (weakly, strongly) controllable, observable, parameterizable, flat, -free...  The problem of parameterizing the solutions of multidimensional linear systems has been extensively studied in control theory by the school of M. Fliess and J. C. Willems. Our results allow us to effectively answer that problem for classes of functional systems. The different algorithms obtained in have been implemented by D. Robertz and A. Quadrat in package OreModules(see Section  and ) and have been illustrated on more than 30 different systems appearing in the literature. See for more details.

The purpose of , is to study the so-called Monge problemfor linear systems over some classes of Ore algebras. This problem was first studied by J. Hadamard and E. Goursat for systems of partial differential equations in the period 1900-1930. Using the algebraic analysis approach developed in , we obtain an algorithm to compute a Monge parametrizationof all the solutions of an underdetermined linear system over an Ore algebra obtained by gluing the autonomous elements of the system to the parameterizable subsystem. The autonomous elements and the parameterizable subsystem can be computed by means of the algorithms developed in . We point out that a Monge parametrization is more general than the parameterizations used in as it depends on arbitrary functions of a certain number of the independent variables whereas we only consider parameterizations depending on all the independent variables in . Effective algorithms checking the given condition are obtained and implemented in OreModules(see Section  and ). In particular, these results are useful if we want to parametrize the solutions of certain systems appearing in mathematical physics and control theory (e.g., tank models, flexible rod models, some linear systems of partial differential equations appearing in the theory of elasticity). See for different examples.

Finally, we have shown in how to use the Monge parametrization in the study of variational problems and optimal control. In particular, we proved that the existence of a Monge parametrization of the system allowed us to transform a variational problem with the system as a differential constraint into a variational problem without differential constraints (which can be solved by means of the Euler-Lagrange equations without Lagrangian multipliers).

Applications of the Stafford theorems and implementation Alban Quadrat Daniel Robertz Aachen University Stafford theorems bases of free modules over the Weyl algebras generators of ideals over the Weyl algebras computation of flat outputs block-diagonalization and block-triangularization problems

It is quite well-known that an analytic time-varying controllable ordinary differential system is flatoutside some singularities . In , we proved that every analytic time-varying controllable linear system was a projection of a flat system and we gave an explicit description of the flat system which projects onto a given controllable one. Moreover, we generalized this result to multidimensional linear systems over Ore algebras. Finally, we proved that every controllable multi-input ordinary differential linear system with polynomial (resp., analytic) coefficients is flat, answering a question asked by K. B. Datta. However, the problem of constructively finding the flat outputs of such a flat system was still open. We have recently studied this problem and solved it using techniques of algebraic analysis and symbolic computation.

A well-known and difficult result due to J. T. Stafford asserts that a stably free left module Mover the Weyl algebras D= An( k)or Bn( k)of differential operators with polynomial/rational coefficients -where kis a field of characteristic 0 -with r a n k D( M) 2is free. The purpose of is to present a constructive proof of this result as well as an effective algorithm for the computation of bases of M. This algorithm, based on the new constructive proofs , of J. T. Stafford's result on the number of generators of left ideals over D, performs Gaussian eliminations on the columns of the formal adjoint of the presentation matrix of M. Moreover, we show that J. T. Stafford's result is a particular case of a more general one asserting that a stably free left D-module Mwith r a n k D( M) s r ( R)is free, where s r ( R)denotes the stable range of a ring D. This result is constructive if the stability of unimodular vectors with entries in Dcan be tested. An algorithm which computes the left projective dimension of a general left D-module Mdefined by means of a finite free resolution is presented. In particular, it allows us to check constructively whether or not the left D-module Mis stably free. All these algorithms have recently been implemented in the package Stafford(see Section  ) based on OreModules(see Section  ). Hence, we can now effectively compute the two generators of any left ideal over Dand a basis of any free left D-module.

These results can be used in order to recognize whether or not a linear system of ordinary/partial differential equations with polynomial/rational coefficients admits an injective parametrization. This problem is interesting in the theory of flat systems , , . Finally, the solution of the decomposition problem obtained in highly uses the computation of bases of free modules. Hence, the decomposition problem can be constructively solved in the case of systems of partial differential equations with polynomial/rational coefficients using the package Stafford.

Applications of the Quillen-Suslin theorem and implementation Anna Fabiańska Aachen University Alban Quadrat Quillen-Suslin theorem flat multidimensional systems differential time-delay systems Lin-Bose conjecture (weakly) doubly factorization block-diagonalization and block-triangularization problems

In 1955, a famous conjecture due to Serre asserted that if we take a row vector vwith coefficients in a commutative polynomial ring A= k[ x1, ..., xn]( kis a field) which admits a right-inverse over A, then vcan be embedded into a square unimodular matrix, namely, a matrix with a non-zero constant determinant. This conjecture was proved independently by Quillen and Suslin in 1976 and is nowadays called the Quillen-Suslin theorem. In terms of module theory, Serre's conjecture means that a projective A-module is free and the computation of the unimodular matrix then gives an explicit basis for this module. Constructive proofs of the Quillen-Suslin theorem have been largely studied in the literature of constructive algebra (e.g., see ). The purpose of the internship of J. Evers (INRIA Sophia Antipolis, CAFÉ project, 2005) was to start implementing the difficult computations of bases for projective A-modules in OreModules. This project has recently been finished by A. Fabiańska (University of Aachen, Germany) within the framework of her PhD thesis supervised by Prof. W. Plesken with the help of A. Quadrat. A package called QuillenSuslinwill be soon available (see Section  ).

The Quillen-Suslin theorem has many applications in mathematical systems theory and signal processing and, in particular, it allows us to constructively compute flat outputs of flat or -flat shift-invariant multidimensional linear systems such as time-invariant differential time-delay control systems. Hence, the present work gives the first implementation for the computation of flat outputs of flat shift-invariant multidimensional linear systems . Applications of such a class of systems were largely demonstrated in the literature of control theory and, in particular, for the motion planning and tracking problems (e.g., see the library of examples of OreModulesand the references therein).

Moreover, using a constructive version of the Quillen-Suslin theorem, A. Fabiańska and A. Quadrat show in that every flat shift-invariant multidimensional linear system is equivalent to a controllable 1-D linear system obtained by setting to 0 all but one functional operators in the system matrix. It was also proved that a flat differential time-delay linear system is equivalent to the controllable ordinary differential linear system without delays, i.e., the system obtained by setting the time-delay amplitude to 0. The implementation of the Quillen-Suslin theorem gives a constructive way to get the precise invertible transforms between these systems. Applications of these results in the study of stabilization problems were started in but will be developed in forthcoming publications. Finally, other applications of the Quillen-Suslin theorem were demonstrated in and, in particular, we give constructive algorithms for proving the standard Lin-Bose's conjecture (generalization of the Serre conjecture), computing (weakly) doubly coprime factorizations of transfer matrices of multidimensional systems and computing bases for the decomposition problems (key point for this problem; see Section  for more details). These different algorithms have been implemented in QuillenSuslin.

Resolvent representation for systems of ordinary differential equations Thomas Cluzeau Evelyne Hubert Differential algebra resolvent representation triangulation-decomposition algorithms

We had previously shown the existence of resolvent representations for systems of (nonlinear) ordinary differential equations that define regular differential ideals. This representation is based on the fact that those systems are birationally equivalent to a single equation. The resolvent representation generalizes the representation used in geometric resolution which has led to probabilistic algorithms of bounded complexity for solving polynomial systems. In we provide practical algorithms for computing such representations. We propose two different approaches. The first one uses differential characteristic decompositions whereas the second one proceeds by prolongation and algebraic elimination. Both constructions depend on the choice of a tuple over the differential base field and their success relies on the chosen tuple to be separating. The probabilistic aspect of the algorithms comes from this choice. To control it, we exhibit a family of tuples for which we can bound the probability that one of its element is separating.

Algebraic and differential Invariants Evelyne Hubert Irina Kogan North Carolina State University Alexandre Sedoglavic University of Lille Eric Schost Ecole Polytechnique, University of Western Ontario Algebraic invariants differential invariants Lie groups algebraic groups cross-section

In we proposed an algorithm to compute a generating set of rational invariants for the rational action of an algebraic group. The algorithms are implemented in Maple ( http://www.inria.fr/cafe/Evelyne.Hubert/aida). In we showed the connection of this algebraic construction with the moving frame method of Fels & Olver . In particular we exhibit a family of algebraic invariants, i.e. algebraic functions of the rational invariants, that not only generate algebraically all algebraic invariants but also generate, functionally, the smooth invariants in the neighbourhood of a point whose orbit is of maximal dimension. Those algebraic invariants occur naturally in differential geometry, as for instance the (Euclidean) curvature.

An essential property of the algebraic invariants we consider is that the field they generate is isomorphic to the field of rational functions on a given cross-section to the orbits. Their algebraic structure is therefore simple and is the basis for determining the differential algebraic structure of the differential invariants. We have made some progress on the determination of this differential algebra.

A particularly favorable case is when the cross-section intersects generic orbits only once, in which case we say the cross-section is of degree one. Such cross-sections are known to exist for specialgroups , which include the general linear group and the connected solvable groups. The joint project with E. Schost has been to provide a construction for those cross-sections of degree 1.

Invariants are the natural objects to translate algebraically the geometry and symmetry of problems in diverse fields of science. They allow to reduce the dimension of problems by factoring out the symmetry. We applied those ideas for the nondimensionalisation of dynamical systems. Reducing the number of parameters in biological models retaining the qualitative properties of the dynamics has often been considered an art . One can observe though that the classical reductions performed correspond to the symmetry reduction for simple scalings or translations. We provided the algorithms to determine such symmetries, compute the invariants and the rewriting algorithms to perform the reductions. The whole problems boils down to linear algebra and the algorithms are shown to be polynomial time . An implementation in Maple applied to a library of examples from the mathematical biology litterature has been constructed.

Algebraic analysis approach to infinite-dimensional linear systems Alban Quadrat Infinite-dimensional systems analysis and synthesis problems internal/strong/simultaneous/robust/optimal stabilization parametrization of all the stabilizing controllers algebraic analysis

Until today, no one has seemed to be able to algebraically incorporate the initial and limit conditions of a system of partial differential equations within the differential module approach ( D-modules). From analytical and engineering points of view, the algebraic D-module approach is not yet totally satisfactory. Driven by applied questions arising in control theory, A. Quadrat has tried to develop a new algebraic analysis approach which could solve this negative point for a class of linear systems of partial differential equations depending on two or three independent variables (typically, one independent variable is the time variable and one or two independent variables are the space variables) usually encountered in engineering sciences. This class of systems is usually called infinite-dimensional linear systems as they have infinite-dimensional state-spaces or, equivalently, a certain number of functions need to be prescribed as initial or limit conditions in order to get a functional solution. Examples of infinite-dimensional systems are the systems of partial differential equations (e.g., wave, heat, Euler-Bernoulli, telegraphist equations), differential time-delay systems (e.g., transport equation, transmission lines, electrical lines) and fractional differential systems. Therefore, we have tried to develop a new algebraic analysis for infinite-dimensional linear systems and apply it to classical problems of control theory.

Using symbolic analysis, also called symbolic calculus, developed by O. Heaviside, P. Dirac and L. Schwartz, i.e., techniques based on Laplace transform, convolution products, distributions and symbolic integration (special functions), we can transform an infinite-dimensional linear system into a transfer matrix, i.e., into functional relations between the state and the outputs of the system and its inputs and initial conditions. The transfer matrix carries important information of the system and particularly its stability (e.g., unstable poles correspond to unstable modes). Combining natural functional algebras, which were introduced for studying the stability of the systems, and the fractional representation approach to systems given by the symbolic analysis approach, we have shown in the publications , , , , how to develop a new module-theoretic approach to infinite-dimensional linear systems.

Using the algebraic concepts of lattices, we have shown in the publications , , how to generalize the well-known Youla-Ku era parametrization of all stabilizing controllers for multi-input multi-ouput systems which do not necessarily admit coprime factorizations. The Youla-Ku era parametrization has been the cornerstone of the success of the -control in the recent years as this parametrization allows us to rewrite the problem of finding the optimal stabilizing controllers (for a certain norm) as affine, and thus, convex problems.

Finally, the algebraic analysis techniques also allowed us to positively solve in the publication the well-known conjecture of Z. Lin in the literature of multidimensional systems on the equivalence between internal stabilizability and the existence of doubly coprime factorizations for multidimensional systems.

Other Grants and Activities European initiatives PAI Amadeus

A collaboration with RISC (Research Institute in Symbolic Computation) in Hagenberg, Austria was funded by the PAI program AMADEUS. The coordinators are Ralf Hemmecke for the Austrian side and Evelyne Hubert replaced Manuel Bronstein in June 2005 for the French side. One objective of the collaboration is to develop and implement a very generic and efficient engine to compute Gröbner basis with ideals and modules over rings of linear operators, such as Ore algebras or Poincare-Birkhoff-de Witt algebras. This collaboration ends this year.

PAI Procope

The purpose of this collaboration of A. Quadrat with the Lehrstuhl B für Mathematik RWTH - Aachen, led by Prof. W. Plesken, is to implement new symbolic packages for the mathematical study of systems theory (e.g., the Stafford, Quillen-Suslin, Morphismspackages of OreModules) and develop their applications to control theory and mathematical physics. This collaboration will continue in 2007.

Visiting scientists France

Frédéric Nataf, Laboratoire J.L. Lions, University Paris 6, visited the CAFÉ project to give a talk on the use of the Smith forms in decomposition domains and perfectly matched layers and to work with Alban Quadrat on the Smith and Jacobson forms (16-17/01).

Eric Schost, Ecole Polytechnique, visited the CAFÉ project to collaborate with E. Hubert (18-20/04).

Jean-François Pommaret, Ecole Nationale des Ponts et Chaussées (CERMICS), visited the CAFÉ project for collaborating with A. Quadrat (25-29/09).

Europe

Within the PAI Procope, A. Fabiańska, University of Aachen, visited A. Quadrat from April 2nd to 7th and from November 20th to 24th. Moreover, D. Robertz, University of Aachen, will visit us in January (07-12/01/07) (see Section  ).

Ralf Hemmecke (RISC Linz, Austria) visited our project for one week in July (10-15) in order to collaborate with Alban Quadrat on the design of the libaldorand Algebralibraries (see Scetion  ) to include non-commutative polynomials such as the ones developed in OreModules. The visit is funded by the PAI Amadeus (see Section  ).

• Evelyne Hubert served as a referee for the PhD thesis of L. D'Alfonso, University of Buenos Aires, Argentina.

• Evelyne Hubert is member of the council of the association femmes & mathématiques. As such she participates to the monthly meetings of the council to propose, discuss and undertake the actions of the association. Those aim at promoting scientific education and careers for women.

• Evelyne Hubert is a member of the visiting scientist committee chaired by Denis Talay.

• Alban Quadrat has just been elected Associate editor of the international journal Multidimensional Systems and Signal Processing ( http://www.springerlink.com/content/1573-0824/) (Springer) for the period 2007-2010.

• Teaching
• Evelyne Hubert taught a course of practical computer algebra in ISIA(April 2006), a Master affiliated to the Ecole des Mines de Paris.

• Alban Quadrat gave four lectures at the Mathematics, Algorithms and Proofs (MAP) workshop(4 hours), Castro Urdiales, Spain (9-13/01) ( http://www.disi.unige.it/map/).

• Alban Quadrat gave two lectures (4 hours) on the factorization and decomposition problems of functional linear systems at the Institut de Recherche en Communications et en Cybernétique de Nantes ( http://www.irccyn.ec-nantes.fr/) (28-30/11).

• Alban Quadrat gave two lectures (4 hours) on the factorization and decomposition problems of functional linear systems at the Lehrstuhl B für Mathematik, Aachen University ( http://wwwb.math.rwth-aachen.de/index_en.html) (18-20/12).

• Dissertations and internships

Doctorates completed in 2006:

• Nicolas Le Roux, Université de Limoges : Solutions Formelles d'Equations aux dérivées partielles. Co-supervised by Moulay Barkatou (University of Limoges) and Evelyne Hubert.

• Doctorates in progress in the project:

• Anna Fabiańska is partially supervised by A. Quadrat in the framework with her PhD thesis at the University of Aachen (Germany) with Prof. W Plesken.

• Conferences and workshops, invited conferences

T. Cluzeau

• participated to Journées Matrice Structurées, University of Limoges, France (18-19/01) ( http://www.unilim.fr/pages_perso/jacques-arthur.weil/CF/matstruct.html).

• gave a talk at the conference International Symposium on Symbolic and Algebraic Computations - ISSAC 2006, ( http://issac2006.dima.unige.it/), Genova, Italy (July 09-12).

• presented his results at the seminar Laboratoire d'Arithmétique, de Calcul Formel et d'Optimisation, ( http://www.unilim.fr/laco/seminaires/calcul/index.html), University of Limoges, France, (07/04).

• gave a talk at the seminar at the Institut d'Informatique et Mathématiques Appliquées de Grenoble (IMAG), Laboratoire de Modélisation et Calcul, Grenoble, France, (13/04).

• E. Hubert

• talked at the seminar of the SALSA project-team in Paris, France (February 14th).

• was invited to present her work at the conference Gröbner bases in Symbolic Analysis http://www.ricam.oeaw.ac.at/srs/groebat the Radon Institute for Computational and Applied Mathematics in Linz, Austria (May 8-12).

• participated to the conference Théorie de l'Elimination http://www-sop.inria.fr/galaad/conf/jouanolou/at the Centre International de Recherche en Mathématiques, Luminy, France (May 15-19).

• gave a talk at the seminar Calcul formel et Complexitéat the Institut de Recherche Mathématique de Rennes, France (June 19th).

• contributed to a panel session at the conference Challenges in Computer Algebra Software( http://www.dagstuhl.de/06271) that took place in Dagstuhl, Germany (July 2-7).

• was invited to present her work at the occasion of the IMA summer program Overdetermined systems of Partial Differential Equations and Symmetries( http://www.ima.umn.edu/2005-2006/SP7.17-8.4.06) in Minneapolis, USA (July 17th - August 4th).

• presented a poster at the workshop Algorithms in Algebraic Geometry( http://www.ima.umn.edu/2006-2007) organized at the Institute For mathematics and its Applications, Minneopolis, USA (September 15th - 20th).

• was invited to give a talk at the workshop Software for Algebraic Geometry( http://www.ima.umn.edu/2006-2007) organized at the Institute For Mathematics and its Applications, Minneopolis, USA (October 22th - 27th).

• was an invited speaker at the workshop Global Integrability of Field Theories( http://www.gift2006.eu) that took place at the Cockcroft Institute in Daresbury, UK, (November 1st - 3rd).

• participated to the workshop on Mathematics in the Digital Age( http://www.ima.umn.edu/2006-2007) organized at the Institute For mathematics and its Applications, Minneopolis, USA (December 7th - 9th).

• participated to the Kolchin Seminar ( http://www.sci.ccny.cuny.edu/~ksda) in New York, USA (December 15th - 16th): she gave a tutorial talk at the Graduate Center series and presented her recent work at the Hunter series.

• was invited to give four lectures at the Mathematics, Algorithms and Proofs (MAP) workshop(4 hours), Castro Urdiales, Spain (9-13/01) ( http://www.disi.unige.it/map/).

• was invited to two talks at the conference Gröbner bases in Symbolic Analysis( http://www.ricam.oeaw.ac.at/srs/groeb) at the Radon Institute for Computational and Applied Mathematics in Linz, Austria (8-12/05) (workshops Gröbner Bases in Symbolic Analysis'' and Gröbner Bases, Control Theory and Signal Processing'').

• was invited to give a talk at the seminar of the Algorithms project, INRIA Rocquencourt, (26/05) ( http://algo.inria.fr/seminars/).

• organized a mini-symposium entitled Symbolic methods in multidimensional systems theory'' at the 17 thInternational Symposium on Mathematical Theory of Networks and Systems (MTNS2006), Kyoto Japan, (24-28/07) ( http://www-ics.acs.i.kyoto-u.ac.jp/mtns06/).

• The CAFÉ project

• organized the conference CAFE, Computer Algebra and Functional Equations. An international conference, in memory of Manuel Bronstein, Sophia Antipolis, France, (13/07) ( http://www-sop.inria.fr/cafe/Cafe2006eng.html).

• On solutions of linear ordinary difference equations in their coefficient field M. Bronstein M. Journal of Symbolic Computation 0747-7171 29 6 2000 841-877 Symbolic Integration I - Transcendental Functions M. Bronstein M. 2 n d Ed., 2004 Springer, Heidelberg 1997 Differential Algebra and Differential Equations (Special Issue of the Journal of Symbolic Computation) M. Bronstein M. W. Sit W. Academic Press, Londres 1999 Notes on triangular sets and triangulation-decomposition algorithms II: Differential Systems E. Hubert E. F. Winkler F. U Langer U. Symbolic and Numerical Scientific Computing LNCS 2630 Springer 2003 40–87 Notes on triangular sets and triangulation-decomposition algorithms I: Polynomial systems. E. Hubert E. F. Winkler F. U Langer U. Symbolic and Numerical Scientific Computing LNCS 2630 Springer 2003 1–39 Algebraic analysis of linear multidimensional control systems J.-F. Pommaret J.-F. A. Quadrat A. IMA Journal of Control and Information 0265-0754 16 3 1999 275-297 Solutions formelles d'équations aux dérivées partielles N. Le Roux N. Ph. D. Thesis Université de Limoges 2006 <span class="smallcap" align="left">OreModules</span>: A symbolic package for the study of multidimensional linear systems F. Chyzak F. A. Quadrat A. D. Robertz D. J. Chiasson J. J.-J. Loiseau J.-J. Applications of Time-Delay Systems LNCIS To appear Springer 2006 Rational Invariants of a Group Action. Construction and Rewriting E. Hubert E. I. Kogan I. Journal of Symbolic Computation 0747-7171 doi:10.1016/j.jsc.2006.03.005 2006 A lattice approach to analysis and synthesis problems A. Quadrat A. Mathematics of Control, Signals, and Systems 0932-4194 18 2 2006 147-186 On a generalization of the Youla-Ku <span class="math" align="left"><img align="bottom" width="5" height="9" src="math_image_25.png" xylemeAttach="25" border="0" alt="Im25 $\mover c\#711$"/></span>era parametrization. Part II: MIMO systems A. Quadrat A. Mathematics of Control, Signals, and Systems 0932-4194 18 3 2006 199-235 Low complexity algorithms for linear recurrences A. Bostan A. F. Chyzak F. T. Cluzeau T. B. Salvy B. ISSAC 2006 ACM press 2006 Using morphism computations for factoring and decomposing general linear functional systems T. Cluzeau T. A. Quadrat A. 17th International Symposium on Mathematical Theory of Networks and Systems (MTNS 2006) 2006 Flat multidimensional linear systems with constant coefficients are equivalent to controllable 1-D linear systems A. Fabiańska A. A. Quadrat A. 17th International Symposium on Mathematical Theory of Networks and Systems (MTNS 2006) 2006 Rank reduction of a class of pfaffian systems in two variables N. Le Roux N. M. Barkatou M. ISSAC '06: Proceedings of the 2006 international symposium on Symbolic and algebraic computation, New York, NY, USA ACM Press 2006 204–211 http://doi.acm.org/10.1145/1145768.1145804 Constructive computation of flat outputs of a class of multidimensional linear systems with variable coefficients A. Quadrat A. D. Robertz D. 17th International Symposium on Mathematical Theory of Networks and Systems (MTNS 2006) 2006 On the Monge problem and multidimensional optimal control A. Quadrat A. D. Robertz D. 17th International Symposium on Mathematical Theory of Networks and Systems (MTNS 2006) 2006 Les systèmes hamiltoniens et leur intégrabilité Cours Spécialisés 8, Collection SMF M. Audin M. SMF et EDP Sciences
Paris
2001
On the infinitesimal geometry of integrable systems A. Baider A. R. C. Churchill R. C. D. L. Rod D. L. M. F. Singer M. F. Mechanics day (Waterloo, ON, 1992), Providence, RI Fields Inst. Commun. 7 Amer. Math. Soc. 1996 5–56 A rational version of Moser's algorithm M. Barkatou M. Proceedings of the International Symposium on Symbolic and Algebraic Computation ACM Press July 1995 297-302 Die Irreduzibilität der homogenen linearen Differentialgleichungen E. Beke E. Math. Ann. 1894 Rings of Differential Operators J.-E. Bjork J.-E. 21 North-Holland Publishing Company 1979 Fast Algorithms for Polynomial Solutions of Linear Differential Equations A. Bostan A. T. Cluzeau T. B. Salvy B. M. Kauers M. Proceedings of ISSAC'2005 ACM Press 2005 45–52 Picard-Vessiot extensions for linear functional systems M. Bronstein M. Z. Li Z. M. Wu M. M. Kauer M. Proceedings of the 2005 International Symposium on Symbolic and Algebraic Computation ACM Press 2005 68-75 Les systèmes différentiels extérieurs et leurs applications géométriques E. Cartan E. Hermann 1945 Triangulation formelle de certains systèmes de Pfaff complètement intégrables et application à l' étude <span class="math" align="left"><img align="bottom" width="21" height="10" src="math_image_26.png" xylemeAttach="26" border="0" alt="Im26 $C^\#8734$"/></span>des systèmes non linéaires. H. Charrière H. Ann. Scuola Norm. Sup. Pisa CI. Sci. 7 4 1980 625-714 Formal Reduction of Integrable Linear Connexion having a certain kind of Irregular Singularities H. Charrière H. R. Gérard R. analysis 1 1981 85-115 F. Chyzak F. Gröbner Bases, Symbolic Summation and Symbolic Integration 251 LMS 1998 32-60 Effective algorithms for parametrizing linear control systems over Ore algebras F. Chyzak F. A. Quadrat A. D. Robertz D. Applicable Algebra in Engineering, Communications and Computing 0938-1279 16 5 2005 319-376 Algorithmique modulaire des équations différentielles linéaires T. Cluzeau T. Thèse de Mathématiques Université de Limoges 2004 Resolvent Representation for Regular Differential Ideals T. Cluzeau T. E. Hubert E. Applicable Algebra in Engineering, Communication and Computing 0938-1279 0938-1279 13 5 2003 395-425 Probabilistic Algorithms for Computing Resolvent Representations of Regular Differential Ideals T. Cluzeau T. E. Hubert E. Technical report HAL INRIA 2006 http://hal.inria.fr/inria-00089287 Using morphism computations for factoring and decomposing general linear functional systems T. Cluzeau T. A. Quadrat A. Technical report INRIA Report 5942 2006 https://hal.inria.fr/inria-00083224 An Introduction to Infinite-Dimensional Linear Systems Theory Texts in Applied Mathematics R. F. Curtain R. F. H. J. Zwart H. J. 21 Springer-Verlag 1991 Equations Différentielles à Points Singuliers Réguliers Lecture Notes In Mathematics P. Deligne P. 163 Springer-Verlag 1970 Regular Singularities along Normal Crossings A. Van Den Essen A. V. D. Gerard Ramis G. Systèmes de Pfaff et Equations Différentielles dans le Champ Complexe Lecture Notes 712 Springer-Verlag 1979 88-130 Irregular Singularities in Several Variables A. Van Den Essen A. V. D. A.H.M. Levelt A. Memoirs of AMS 40 270 1982 Applications of the Quillen-Suslin theorem in multidimensional systems theory A. Fabiańska A. A. Quadrat A. submitted 2006 Moving coframes. II. Regularization and theoretical foundations M. Fels M. P. J. Olver P. J. Acta Appl. Math. 55 2 1999 127–208 A Lie-Bäcklund approach to equivalence and flatness of nonlinear systems M. Fliess M. J. Lévine J. P. Martin P. P. Rouchon P. IEEE Trans. Automatic Control 0018-9286 49 1999 922-937 Controllability and observability of linear delay systems: an algebraic approach M. Fliess M. H. Mounier H. ESAIM: Control, Optimisation and Calculus of Variations 1292-8119 3 1998 301-314 Invariants mesurant l'irrégularité en un point singulier des systèmes d' équations différentielles linéaires R. Gérard R. A.H.M. Levelt A. Annales de l'Institut Fourier 0373-0956 23 1 1973 157-195 Sur les Connexions à Singularités Régulières dans le cas de Plusieurs Variables R. Gérard R. A.H.M. Levelt A. Funkcialaj Ekvacioj 19 2 1976 149-173 Solving Difference Equations in Finite Terms P.A. Hendriks P. M.F. Singer M. Journal of Symbolic Computation 0747-7171 27 3 March 1999 239-260 Towards an effective version of a theorem of Stafford A. Hillebrand A. W. Schmale W. J. Symbolic Computation 0747-7171 32 2001 699-716 Essential Components of an Algebraic Differential Equation E. Hubert E. Journal of Symbolic Computation 0747-7171 28 4-5 1999 657-680 Smooth and Algebraic Invariants of a Group Action. Local and Global Constructions E. Hubert E. I. Kogan I. submitted January 2006 Polynomial Time Nondimensionalisation of Ordinary Differential Equations via their Lie Point Symmetries E. Hubert E. A. Sedoglavic A. Technical report HAL-INRIA 2006 http://hal.inria.fr/inria-00001251 Leçons sur les systèmes d'équations aux dérivées partielles M. Janet M. Cahiers Scientifiques IV, Gauthier-Villars 1929 Algebraic Study of Systems of Partial Differential Equations Mémoires de la Société Mathématiques de France M. Kashiwara M. 63 1995 Differential Algebra and Algebraic Groups Pure and Applied Mathematics E. R. Kolchin E. R. 54 Academic Press, New York-London 1973 An Algorithm for Solving Second Order Linear Homogeneous Differential Equations J.J. Kovacic J. Journal of Symbolic Computation 0747-7171 2 1 March 1986 3-43 Non-commutative Computer Algebra for Polynomial Algebras: Gröbner Bases, Applications and Implementation V. Levandovskyy V. Ph. D. Thesis University of Kaiserslautern 2005 Algorithmic proofs of two theorems of Stafford A. Leykin A. J. Symbolic Computation 0747-7171 38 2004 1535-1550 Rational Solutions of Riccati-like Partial Differential Equations Ziming Li Z. Fritz Schwarz F. Journal of Symbolic Computation 0747-7171 31 6 June 2001 691–716 Factoring systems of linear PDEs with finite-dimensional solution spaces Ziming Li Z. Fritz Schwarz F. Serguei P. Tsarev S. P. Journal of Symbolic Computation 0747-7171 36 3–4 September/October 2003 443–471 A generalization of Serre's conjecture and related issues Z. Lin Z. N.K. Bose N. Linear Algebra and Its Applications 0024-3795 338 2001 125-138 Output feedback stabilizability and stabilization of linear <hi rend="italic">n</hi>-D systems Z. Lin Z. K. Galkowski K. J. Wood J. Multidimensional Signals, Circuits and Systems Taylor and Francis 1994 122-278 Algorithms for the Quillen-Suslin theorem A. Logar A. B. Sturmfels B. Journal of Algebra 145 1992 231-239 F. Malrait F. P. Martin P. P. Rouchon P. Dynamic feedback transformations of controllable linear time-varying systems 259 LNCIS 2001 55-62 Differential Galois theory and non-integrability of Hamiltonian systems Progress in Mathematics J. J. Morales Ruiz J. J. 179 Birkhäuser Verlag
Basel
1999
The Order of a Singularity in Fuchs' Theory J. Moser J. Mathematische Zeitschrift 72 1960 379-398 Mathematical Biology Interdisciplinary Applied Mathematics James D. Murray J. D. 17 Springer 2002 Algorithms for <hi rend="italic">D</hi>-modules <span class="math" align="left">-</span>restriction, tensor product, localization, and local cohomology groups T. Oaku T. N. Takayama N. J. Pure Appl. Algebra 0022-4049 156 2-3 2001 267-308 Polynomial and rational solutions of holonomic systems T. Oaku T. N. Takayama N. H. Tsai H. J. Pure Appl. Algebra 0022-4049 164 1-2 2001 199-220 Linear Differential Operators with Constant Coefficients V. Palamodov V. Springer-Verlag 1970 Invariant theory V. L. Popov V. L. E. B. Vinberg E. B. Algebraic geometry. IV 1994 122-278 On a generalization of the Youla-Ku <span class="math" align="left"><img align="bottom" width="6" height="9" src="math_image_23.png" xylemeAttach="23" border="0" alt="Im23 $\mover c\#711$"/></span>era parametrization. Part I: The fractional ideal approach to SISO systems A. Quadrat A. Systems and Control Letters 0167-6911 50 2 2003 135–148 The fractional representation approach to synthesis problems: an algebraic analysis viewpoint. Part I: (weakly) doubly coprime factorizations A. Quadrat A. SIAM J. Control and Optimization 0363-0129 42 1 2003 266-299 The fractional representation approach to synthesis problems: an algebraic analysis viewpoint. Part II: internal stabilization A. Quadrat A. SIAM J. Control and Optimization 0363-0129 42 1 2003 300-320 On a general structure of the stabilizing controllers based on stable range A. Quadrat A. SIAM J. Control and Optimization 0363-0129 42 6 2004 2264-2285 An algebraic interpretation to the operator-theoretic approach to stabilizability. Part I: SISO systems A. Quadrat A. Acta Applicandæ Mathematicæ 0167-8019 88 1 2005 1-45 Constructive computation of bases of free modules over the Weyl algebras A. Quadrat A. D. Robertz D. Technical report INRIA Report 5786 2005 http://hal.inria.fr/inria-00070235 On the blowing-up of stably free behaviours A. Quadrat A. D. Robertz D. Proceedings of CDC-ECC05 CDRom 2005 Parametrizing all solutions of uncontrollable multidimensional linear systems A. Quadrat A. D. Robertz D. Proceedings of the 16th IFAC World Congress CDRom 2005 On the Monge problem for multidimensional linear systems A. Quadrat A. D. Robertz D. submitted 2006 La Méthode des Fonctions Majorantes et les Systèmes d'Equations aux Dérivées Partielles C. Riquier C. Gauthier-Villars 1928 Differential Algebra Colloquium publications J. F. Ritt J. F. XXXIII American Mathematical Society 1950 http://www.ams.org/online_bks/coll33/ Formal Computational Methods for Control Theory D. Robertz D. Ph. D. Thesis University of Aachen 2006 Liouvillian Solutions of <span class="math" align="left"><hi rend="it">n</hi><sup>t h</sup></span>Order Homogeneous Linear Differential Equations M. Singer M. American Journal of Mathematics 103 1980 661-682 Testing Reducibility of Linear Differential Operators: a Group Theoretic Perspective M.F. Singer M. Applicable Algebra in Engineering, Communication and Computing 0938-1279 0938-1279 7 1996 77–104 Galois Theory of Linear Differential Equations Grundlehren der Mathematischen Wissenschaften M.F. Singer M. M. van der Put M. 328 Springer, Heidelberg 2003 Galois Theory of Difference Equations LNM 1666 Michael Singer M. Marius van der Put M. Springer 1997 Overdetermined systems of partial differential equations D. C. Spencer D. C. Bull. Amer. Math. Soc. 75 1965 1-114 On a Linear System of Pfaffian Equations with Regular Singular Points T. Takano. T. M. Yoshida M. Funkcialaj Ekvacioj 19 1976 175–189 Sur les invariants des groupes continus de transformations A. Tresse A. Acta Mathematica 18 1894 1-88 Control System Synthesis: A Factorization Approach M. Vidyasagar M. MIT Press 1985 Asymptotic Expansions of Ordinary Differential Equations W. Wasow W. Dover Phoenix Editions 2002 On Solutions of Linear Functional Systems and Factorization of Modules over Laurent-Ore Algebras M. Wu M. Ph. D. Thesis Université de Nice-Sophia Antipolis and Chinese Academy of Sciences 2005 Branching of Solutions and non-existence of first integrals in Hamiltonian mechanics. I S. L. Ziglin S. L. Functional Anal. Appl. 16 1982 181–189 Branching of Solutions and non-existence of first integrals in Hamiltonian mechanics. II S. L. Ziglin S. L. Functional Anal. Appl. 17 1983 6–17 Low complexity algorithms for linear recurrences A. Bostan A. F. Chyzak F. T. Cluzeau T. B. Salvy B. ISSAC 2006 ACM press 2006 Using morphism computations for factoring and decomposing general linear functional systems T. Cluzeau T. A. Quadrat A. 17th International Symposium on Mathematical Theory of Networks and Systems (MTNS 2006) 2006 Flat multidimensional linear systems with constant coefficients are equivalent to controllable 1-D linear systems A. Fabiańska A. A. Quadrat A. 17th International Symposium on Mathematical Theory of Networks and Systems (MTNS 2006) 2006 Solutions formelles d'équations aux dérivées partielles N. Le Roux N. Ph. D. Thesis Université de Limoges 2006 Rank reduction of a class of pfaffian systems in two variables N. Le Roux N. M. Barkatou M. ISSAC '06: Proceedings of the 2006 international symposium on Symbolic and algebraic computation, New York, NY, USA ACM Press 2006 204–211 http://doi.acm.org/10.1145/1145768.1145804 Constructive computation of flat outputs of a class of multidimensional linear systems with variable coefficients A. Quadrat A. D. Robertz D. 17th International Symposium on Mathematical Theory of Networks and Systems (MTNS 2006) 2006 On the Monge problem and multidimensional optimal control A. Quadrat A. D. Robertz D. 17th International Symposium on Mathematical Theory of Networks and Systems (MTNS 2006) 2006 A lattice approach to analysis and synthesis problems A. Quadrat A. Mathematics of Control, Signals, and Systems 0932-4194 18 2 2006 147-186 On a generalization of the Youla-Kuc ˇera parametrization. Part II: MIMO systems A. Quadrat A. Mathematics of Control, Signals, and Systems 0932-4194 18 3 2006 199-235