CHAMP: A Cherednik Algebra Magma Package

We present a computer algebra package based on Magma for performing computations in rational Cherednik algebras at arbitrary parameters and in Verma modules for restricted rational Cherednik algebras. Part of this package is a new general Las Vegas algorithm for computing the head and the constituents of a module with simple head in characteristic zero which we develop here theoretically. This algorithm is very successful when applied to Verma modules for restricted rational Cherednik algebras and it allows us to answer several questions posed by Gordon in some specific cases. We could determine the decomposition matrices of the Verma modules, the graded G-module structure of the simple modules, and the Calogero-Moser families of the generic restricted rational Cherednik algebra for around half of the exceptional complex reflection groups. In this way we could also confirm Martino's conjecture for several exceptional complex reflection groups.Comment: Final version to appear in LMS J. Comput. Math. 41 pages, 3 ancillary files. CHAMP is available at http://thielul.github.io/CHAMP/. All results are listed explicitly in the ancillary PDF document (currently 935 pages). Please check the website for further update

Solving the "Isomorphism of Polynomials with Two Secrets" Problem for all Pairs of Quadratic Forms

We study the Isomorphism of Polynomial (IP2S) problem with m=2 homogeneous quadratic polynomials of n variables over a finite field of odd characteristic: given two quadratic polynomials (a, b) on n variables, we find two bijective linear maps (s,t) such that b=t . a . s. We give an algorithm computing s and t in time complexity O~(n^4) for all instances, and O~(n^3) in a dominant set of instances. The IP2S problem was introduced in cryptography by Patarin back in 1996. The special case of this problem when t is the identity is called the isomorphism with one secret (IP1S) problem. Generic algebraic equation solvers (for example using Gr\"obner bases) solve quite well random instances of the IP1S problem. For the particular cyclic instances of IP1S, a cubic-time algorithm was later given and explained in terms of pencils of quadratic forms over all finite fields; in particular, the cyclic IP1S problem in odd characteristic reduces to the computation of the square root of a matrix. We give here an algorithm solving all cases of the IP1S problem in odd characteristic using two new tools, the Kronecker form for a singular quadratic pencil, and the reduction of bilinear forms over a non-commutative algebra. Finally, we show that the second secret in the IP2S problem may be recovered in cubic time

On the computation of π\pi-flat outputs for differential-delay systems

We introduce a new definition of $\pi$-flatness for linear differential delay systems with time-varying coefficients. We characterize $\pi$- and $\pi$-0-flat outputs and provide an algorithm to efficiently compute such outputs. We present an academic example of motion planning to discuss the pertinence of the approach.Comment: Minor corrections to fit with the journal versio

Algebraic computations in derived categories

AbstractThis paper presents explicit algorithms for computations over a finite subspectroid of the bounded derived category of a finite spectroid. We will demonstrate methods for the construction of a projective resolution of a module and for finding the quiver of a finite spectroid given in terms of its radical spaces. This enables us to compute the endomorphism algebra of a tilting complex – or, in fact, any finite complex – in the derived category. In order to carry out these computations, we have to restrict to a finite base field or the field of rational numbers. We will show that it is possible to transfer the results to any extension of the base field, in particular to the algebraic closure

An Axiomatic Setup for Algorithmic Homological Algebra and an Alternative Approach to Localization

In this paper we develop an axiomatic setup for algorithmic homological algebra of Abelian categories. This is done by exhibiting all existential quantifiers entering the definition of an Abelian category, which for the sake of computability need to be turned into constructive ones. We do this explicitly for the often-studied example Abelian category of finitely presented modules over a so-called computable ring $R$, i.e., a ring with an explicit algorithm to solve one-sided (in)homogeneous linear systems over $R$. For a finitely generated maximal ideal $\mathfrak{m}$ in a commutative ring $R$ we show how solving (in)homogeneous linear systems over $R_{\mathfrak{m}}$ can be reduced to solving associated systems over $R$. Hence, the computability of $R$ implies that of $R_{\mathfrak{m}}$. As a corollary we obtain the computability of the category of finitely presented $R_{\mathfrak{m}}$-modules as an Abelian category, without the need of a Mora-like algorithm. The reduction also yields, as a by-product, a complexity estimation for the ideal membership problem over local polynomial rings. Finally, in the case of localized polynomial rings we demonstrate the computational advantage of our homologically motivated alternative approach in comparison to an existing implementation of Mora's algorithm.Comment: Fixed a typo in the proof of Lemma 4.3 spotted by Sebastian Posu

Computational Aspects of Retrieving a Representation of an Algebraic Geometry Code

Producción CientíficaCode-based cryptography is an interesting alternative to classic number-theoretic public key cryptosystem since it is conjectured to be secure against quantum computer attacks. Many families of codes have been proposed for these cryptosystems such as algebraic geometry codes. In [Designs, Codes and Cryptography, pages 1-16, 2012] -for so called very strong algebraic geometry codes $\mathcal C=C_L(\mathcal X, \mathcal P, E)$, where $\mathcal X$ is an algebraic curve over $\mathbb F_q$, $\mathcal P$ is an $n$-tuple of mutually distinct $\mathbb F_q$-rational points of $\mathcal X$ and $E$ is a divisor of $\mathcal X$ with disjoint support from $\mathcal P$ --- it was shown that an equivalent representation $\mathcal C=C_L(\mathcal Y, \mathcal Q, F)$ can be found. The $n$-tuple of points is obtained directly from a generator matrix of $\mathcal C$, where the columns are viewed as homogeneous coordinates of these points. The curve $\mathcal Y$ is given by $I_2(\mathcal Y)$, the homogeneous elements of degree $2$ of the vanishing ideal $I(\mathcal Y)$. Furthermore, it was shown that $I_2(\mathcal Y)$ can be computed efficiently as the kernel of certain linear map. What was not shown was how to get the divisor $F$ and how to obtain efficiently an adequate decoding algorithm for the new representation. The main result of this paper is an efficient computational approach to the first problem, that is getting $F$. The security status of the McEliece public key cryptosystem using algebraic geometry codes is still not completely settled and is left as an open problemThis research was partly supported by the Danish National Research Foundation and the National Science Foundation of China (Grant No.\ 11061130539) for the Danish-Chinese Center for Applications of Algebraic Geometry in Coding Theory and Cryptography and by Spanish grants MTM2007-64704, MTM2010-21580-C02-02 and MTM2012-36917-C03-03. Part of the research of the second author is also funded by the Vernon Wilson Endowed Chair at Eastern Kentucky University during his sabbatical leave

An introduction to constructive algebraic analysis and its applications

This text is an extension of lectures notes I prepared for les Journées Nationales de Calcul Formel held at the CIRM, Luminy (France) on May 3-7, 2010. The main purpose of these lectures was to introduce the French community of symbolic computation to the constructive approach to algebraic analysis and particularly to algebraic D-modules, its applications to mathematical systems theory and its implementations in computer algebra systems such as Maple or GAP4. Since algebraic analysis is a mathematical theory which uses different techniques coming from module theory, homological algebra, sheaf theory, algebraic geometry, and microlocal analysis, it can be difficult to enter this fascinating new field of mathematics. Indeed, there are very few introducing texts. We are quickly led to Björk's books which, at first glance, may look difficult for the members of the symbolic computation community and for applied mathematicians. I believe that the main issue is less the technical difficulty of the existing references than the lack of friendly introduction to the topic, which could have offered a general idea of it, shown which kind of results and applications we can expect and how to handle the different computations on explicit examples. To a very small extent, these lectures notes were planned to fill this gap, at least for the basic ideas of algebraic analysis. Since, we can only teach well what we have clearly understood, I have chosen to focus on my work on the constructive aspects of algebraic analysis and its applications.Ce texte est une extension des notes de cours que j'ai préparés pour les les Journées Nationales de Calcul Formel qui ont eu lieu au CIRM, Luminy (France) du 3 au 7 Mai 2010. Le but principal de ce cours était d'introduire la communauté française du calcul formel à l'analyse algébrique constructive, et particulièrement à la théorie des D-modules algébriques, à ses applications à la théorie mathématique des systèmes et à ses implantations dans des logiciels de calcul formel tels que Maple ou GAP4. Parce que l'analyse algébrique est une théorie mathématique qui utilise différentes techniques venant de la théorie des modules, de l'algèbre homologique, de la théorie des faisceaux, de la géométrie algébrique et de l'analyse microlocale, il peut être difficile d'entrer dans ce domaine, nouveau et fascinant, des mathématiques. En effet, il existe peu de textes introductifs. Nous sommes rapidement conduits aux livres de Björk qui, à première vue, peuvent sembler difficiles aux membres de la communauté de calcul formel et aux mathématiciens appliqués. Je pense que le problème vient moins de la difficulté technique de la littérature existante que du manque d'introductions pédagogiques qui donnent une idée globale du domaine, montrent quels types de résultats et d'applications on peut attendre et qui développent les différents calculs à mener sur des exemples explicites. A leur humble niveau, ces notes de cours ont pour but de combler ce manque, tout du moins en ce qui concerne les idées de base de l'analyse algébrique. Puisque l'on ne peut enseigner bien que les choses que l'on a bien comprises, j'ai choisi de restreindre cette introduction à mes travaux sur les aspects constructifs de l'analyse algébrique et sur ses applications

The degree-distance and transmission-adjacency matrices

Let $G$ be a connected graph with adjacency matrix $A(G)$. The distance matrix $D(G)$ of $G$ has rows and columns indexed by $V(G)$ with $uv$-entry equal to the distance $\mathrm{dist}(u,v)$ which is the number of edges in a shortest path between the vertices $u$ and $v$. The transmission $\mathrm{trs}(u)$ of $u$ is defined as $\sum_{v\in V(G)}\mathrm{dist}(u,v)$. Let $\mathrm{trs}(G)$ be the diagonal matrix with the transmissions of the vertices of $G$ in the diagonal, and $\mathrm{deg}(G)$ the diagonal matrix with the degrees of the vertices in the diagonal. In this paper we investigate the Smith normal form (SNF) and the spectrum of the matrices $D^{\mathrm{deg}}_+(G):=\mathrm{deg}(G)+D(G)$, $D^{\mathrm{deg}}(G):=\mathrm{deg}(G)-D(G)$, $A^{\mathrm{trs}}_+(G):=\mathrm{trs}(G)+A(G)$ and $A^{\mathrm{trs}}(G):=\mathrm{trs}(G)-A(G)$. In particular, we explore how good the spectrum and the SNF of these matrices are for determining graphs up to isomorphism. We found that the SNF of $A^{\mathrm{trs}}$ has an interesting behaviour when compared with other classical matrices. We note that the SNF of $A^{\mathrm{trs}}$ can be used to compute the structure of the sandpile group of certain graphs. We compute the SNF of $D^{\mathrm{deg}}_+$, $D^{\mathrm{deg}}$, $A^{\mathrm{trs}}_+$ and $A^{\mathrm{trs}}$ for several graph families. We prove that complete graphs are determined by the SNF of $D^{\mathrm{deg}}_+$, $D^{\mathrm{deg}}$, $A^{\mathrm{trs}}_+$ and $A^{\mathrm{trs}}$. Finally, we derive some results about the spectrum of $D^{\mathrm{deg}}$ and $A^{\mathrm{trs}}$.Comment: 19 page

Une approche par l’analyse algébrique effectivedes systèmes linéaires sur des algèbres de Ore

The purpose of this paper is to present a survey on the effective algebraic analysis approach to linear systems theory with applications to control theory and mathematical physics. In particular, we show how the combination of effective methods of computer algebra - based on Gröbner basis techniques over a class of noncommutative polynomial rings of functional operators called Ore algebras - and constructive aspects of module theory and homological algebra enables the characterization of structural properties of linear functional systems. Algorithms are given and a dedicated implementation, called OreAlgebraicAnalysis, based on the Mathematica package HolonomicFunctions, is demonstrated.Le but de ce papier est de présenter un état de l’art d’une approche par l’analyse algébrique effective de la théorie des systèmes linéaires avec des applications à la théorie du contrôle et à la physique mathématique.En particulier, nous montrons comment la combinaison des méthodes effectives de calcul formel - basées sur lestechniques de bases de Gröbner sur une classe d’algèbres polynomiales noncommutatives d’opérateurs fonctionnels appelée algèbres de Ore - et d’aspects constructifs de théorie des modules et d’algèbre homologique permet lacaractérisation de propriétés structurelles des systèmes linéaires fonctionnels. Des algorithmes sont donnés et uneimplémentation dédiée, appelée OREALGEBRAICANALYSIS, basée sur le package Mathematica HOLONOMIC-FUNCTIONS, est présenté