Journal of Symbolic Computation最新文献

In a series of papers we develop symbolic computation algorithms to investigate finiteness of central configurations for the planar n-body problem. Our approach is based on Albouy-Kaloshin's work on finiteness of central configurations for the 5-body problems. In their paper, bicolored graphs called zw-diagrams were introduced for possible scenarios when the finiteness conjecture fails, and proving finiteness amounts to exclusions of central configurations associated to these diagrams. Following their method, the amount of computations becomes enormous when there are more than five bodies. In this paper we introduce matrix algebra for determination of possible diagrams and asymptotic orders, devise several criteria to reduce computational complexity, and determine possible zw-diagrams by automated deductions. For the planar six-body problem, we show that there are at most 86 zw-diagrams.

In recent years there has been increased use of machine learning (ML) techniques within mathematics, including symbolic computation where it may be applied safely to optimise or select algorithms. This paper explores whether using explainable AI (XAI) techniques on such ML models can offer new insight for symbolic computation, inspiring new implementations within computer algebra systems that do not directly call upon AI tools. We present a case study on the use of ML to select the variable ordering for cylindrical algebraic decomposition. It has already been demonstrated that ML can make the choice well, but here we show how the SHAP tool for explainability can be used to inform new heuristics of a size and complexity similar to those human-designed heuristics currently commonly used in symbolic computation.

We study linear PDE with constant coefficients. The constant rank condition on a system of linear PDEs with constant coefficients is often used in the theory of compensated compactness. While this is a purely linear algebraic condition, the nonlinear algebra concept of primary decomposition is another important tool for studying such system of PDEs. In this paper we investigate the connection between these two concepts. From the nonlinear analysis point of view, we make some progress in the study of weak lower semicontinuity of integral functionals defined on sequences of PDE constrained fields, when the PDEs do not have constant rank.

For any zero-dimensional polynomial ideal $I$ and any nonzero polynomial F, this paper shows that the union of the multi-set of zeros of the ideal sum $I+〈F〉$ and that of the ideal quotient $I:〈F〉$ is equal to the multi-set of zeros of $I$, where zeros are counted with multiplicities. Based on this zero relation and the computation of Gröbner bases, a complete multiplicity-preserved algorithm is proposed to decompose any zero-dimensional polynomial set into finitely many squarefree normal triangular sets, resulting in a squarefree normal representation for the zeros of the polynomial set. In the representation the multiplicities of the zeros of the triangular sets can be read out directly. Examples and experiments are presented to illustrate the algorithm and its performance.

For a wide family of formations $F$ it is proved that the $F$-residual of a permutation finite group can be computed in polynomial time. Moreover, if in the previous case $F$ is hereditary, then the $F$-subnormality of a subgroup can be checked in polynomial time.

In algebraic geometry, it is important to provide effective parametrizations for families of curves, both in theory and in practice. In this paper, we present such an effective parametrization for the moduli of genus-5 curves that are neither hyperelliptic nor trigonal. Subsequently, we construct an algorithm for a complete enumeration of non-special genus-5 curves having more rational points than a specified bound, where “non-special curve” means that the curve is non-hyperelliptic and non-trigonal with mild singularities of the associated sextic model that we propose. As a practical application, we implement this algorithm using the computer algebra system MAGMA, specifically for curves over the prime field of characteristic 3.

We describe an algorithm that constructs a list of all topological types of holomorphic actions of a finite group on a compact Riemann surface C of genus $g\ge 2$ with $C/G\cong P^1$.

This paper is concerned with Smith forms of bivariate polynomial matrices. For a bivariate polynomial square matrix with the determinant being the product of two distinct and irreducible univariate polynomials, we prove that it is equivalent to its Smith form. We design an algorithm to reduce this class of bivariate polynomial matrices to their Smith forms, and an example is given to illustrate the algorithm. Furthermore, we extend the above class of matrices to a more general case, and derive a necessary and sufficient condition for the equivalence of a matrix and one of its all possible existing Smith forms.

We lower bound the rank of a tensor by a linear combination of the ranks of three of its unfoldings, using Sylvester's rank inequality. In a similar way, we lower bound the symmetric rank by a linear combination of the symmetric ranks of three unfoldings. Lower bounds on the rank and symmetric rank of tensors are important for finding counterexamples to Comon's conjecture. A real counterexample to Comon's conjecture is a tensor whose real rank and real symmetric rank differ. Previously, only one real counterexample was known, constructed in a paper of Shitov. We divide the construction into three steps. The first step involves linear spaces of binary tensors. The second step considers a linear space of larger decomposable tensors. The third step is to verify a conjecture that lower bounds the symmetric rank, on a tensor of interest. We use the construction to build an order six real tensor whose real rank and real symmetric rank differ.

The MacWilliams' Extension Theorem is a classical result by Florence Jessie MacWilliams. It shows that every linear isometry between linear block-codes endowed with the Hamming distance can be extended to a linear isometry of the ambient space. Such an extension fails to exist in general for rank-metric codes, that is, one can easily find examples of linear isometries between rank-metric codes which cannot be extended to linear isometries of the ambient space. In this paper, we explore to what extent a MacWilliams' Extension Theorem may hold for rank-metric codes. We provide an extensive list of examples of obstructions to the existence of an extension, as well as a positive result.