Complete Reduction for Derivatives in a Primitive Tower
Du, Gao, Li et al.
A complete reduction $Ï$ for derivatives in a differential field is a linear operator on the field over its constant subfield. The reduction enables us to decompose an element $f$ as the sum of a derivative and the remainder $Ï(f)$. A direct application of $Ï$ is that $f$ is in-field integrable if and only if $Ï(f) = 0.$
In this paper, we present a complete reduction for derivatives in a primitive tower algorithmically. Typical examples for primitive towers are differential fields generated by (poly-)logarithmic functions and logarithmic integrals. Using remainders and residues, we provide a necessary and sufficient condition for an element from a primitive tower to have an elementary integral, and discuss how to construct telescopers for non-D-finite functions in some special primitive towers.
academic
Réduction Complète pour les Dérivées dans une Tour Primitive
Titre: Complete Reduction for Derivatives in a Primitive Tower
Auteurs: Hao Du (Université Postale de Pékin), Yiman Gao (Université Johannes Kepler), Wenqiao Li (Laboratoire Clé de Mécanisation Mathématique de l'Académie Chinoise des Sciences), Ziming Li (Laboratoire Clé de Mécanisation Mathématique de l'Académie Chinoise des Sciences)
Classification: cs.SC (Calcul Symbolique)
Conférence de Publication: ISSAC'25 (Symposium International sur le Calcul Symbolique et Algébrique)
La réduction complète ϕ des dérivées dans un corps différentiel est un opérateur linéaire du corps sur son sous-corps des constantes. Cette réduction nous permet de décomposer un élément f en une somme de dérivées et d'un reste ϕ(f). Une application directe de ϕ est que f est intégrable dans le corps si et seulement si ϕ(f)=0. Cet article propose algorithmiquement la réduction complète des dérivées dans une tour primitive. Les exemples typiques de tours primitives sont les corps différentiels engendrés par des fonctions (multi-)logarithmiques et l'intégrale logarithmique. En utilisant les restes et les résidus, nous fournissons des conditions nécessaires et suffisantes pour qu'un élément dans une tour primitive possède une intégrale élémentaire, et nous discutons de la construction de telescoper pour certaines fonctions non-D-finies dans des tours primitives particulières.
Problème central en intégration symbolique: En calcul symbolique, déterminer si une fonction possède une intégrale sous forme élémentaire est un problème fondamental. Pour les fonctions transcendantes de Liouville, ce problème est généralement décrit par des extensions monomiales.
Importance de la réduction complète: La réduction complète est un opérateur linéaire capable de décomposer tout élément d'un corps différentiel en une partie dérivée et un reste « minimal ». Cette décomposition est essentielle pour:
Juger l'intégrabilité des fonctions dans le corps
Le telescoping créatif basé sur la réduction
L'intégration (sommation) en termes finis
Limitations des méthodes existantes:
La décomposition additive n'est pas toujours une application linéaire, manquant de commodité théorique et pratique
Les réductions complètes existantes se concentrent principalement sur des types spécifiques: fonctions surexponentielles, fonctions algébriques, fonctions D-finies, etc.
Absence d'algorithme systématique de réduction complète pour la catégorie importante des tours primitives
Établissement d'un cadre algorithmique pour la réduction complète des dérivées dans les tours primitives: Proposition d'une méthode systématique en trois étapes pour construire la réduction complète
Développement d'algorithmes auxiliaires clés: Incluant les algorithmes de réduction auxiliaire (AuxiliaryReduction), de construction de base (Basis) et de projection (Projection)
Fourniture de conditions nécessaires et suffisantes pour les intégrales élémentaires: Critères de jugement basés sur les restes et les résidus pour les éléments dans les tours primitives
Extension de la méthode de construction de telescoper: Conditions suffisantes pour l'existence de telescoper pour certaines fonctions non-D-finies
Implémentation d'algorithmes efficaces: Les expériences montrent que cette méthode surpasse les méthodes existantes dans la plupart des cas
Étant donnée une tour primitive F0⊂F1⊂⋯⊂Fn, où Fi=Fi−1(ti) et ti est un monôme primitif sur Fi−1, l'objectif est de construire une réduction complète ϕ:Fn→Fn telle que:
Pour tout f∈Fn, il existe un unique g∈Fn et r∈im(ϕ) tels que f=g′+r
Pour l'extension par monôme primitif F(t), l'algorithme se divise en trois étapes:
Étape 1: Définition du sous-espace auxiliaire
Définir A=im(ϕ)⊗CC[t] comme le sous-espace auxiliaire de F[t]′ dans F[t], où ϕ:F→F est la réduction complète déjà existante sur F.
Étape 2: Détermination de la base de l'intersection
Construire une C-base {v0,v1,v2,…} de F[t]′∩A, où:
v0=ϕ(t′)
vi=ϕ(t′)ti−Mi,0(ti) (pour i≥1)
Étape 3: Fixation de l'espace supplémentaire
Déterminer l'espace supplémentaire Aθ de A dans F[t] par rapport à F[t]′ via des techniques de base efficaces.
Entrée: p ∈ F[t]
Sortie: (q,r) ∈ F[t] × A tel que p = q' + r
1. Initialiser p̃ ← p, q ← 0, r ← 0
2. tant que p̃ ≠ 0 faire
d ← deg(p̃), l ← lc(p̃)
Calculer la paire-R de l: (g, φ(l))
q ← q + gt^d, r ← r + φ(l)t^d
p̃ ← p̃ - lt^d - (dgt')t^(d-1)
3. retourner (q,r)
Algorithme 3.12 (Projection):
Projeter les éléments du sous-espace auxiliaire vers F[t]′ et l'espace supplémentaire θ.
Les expériences utilisent la tour primitive Q(x)(t1,t2,t3), où:
t1=log(x)
t2=log(x+1)
t3=log(t1)
Quatre groupes de fonctions intégrandes de complexité différente ont été testés, chaque groupe contenant des dérivées polynomiales de degrés différents.
La méthode CR surpasse globalement la méthode AD, avec de meilleures performances dans la plupart des cas de test
Par rapport à la fonction int de Maple, CR montre une performance supérieure pour les problèmes de complexité plus élevée, mais est légèrement plus lente dans les cas simples
Meilleure stabilité: CR et AD peuvent traiter certains problèmes d'intégration que la fonction int ne peut pas résoudre
Analyse des composants d'algorithme: HermiteReduce et AuxiliaryReduction sont les parties les plus coûteuses en temps, tandis que Projection est relativement efficace
Exemple 4.5: Pour la fonction
f=x2(x−1)t22((x−1)2t1+x)t23+x(x−1)t1
CR a trouvé avec succès son intégrale, tandis que Maple et Mathematica n'ont pas pu donner de résultat sous forme élémentaire.
Exemple 5.4: Démontre le processus complet de calcul d'intégrale élémentaire, incluant l'analyse des restes et le calcul des résidus.
Théorie de la réduction complète: Réductions complètes existantes pour les fonctions surexponentielles 5, les fonctions algébriques 12,15, les fonctions D-finies 6,13,35
Décomposition additive: Applications dans le telescoping créatif basé sur la réduction 2-4,24
Intégration symbolique: Algorithmes d'intégrale élémentaire pour les fonctions transcendantes de Liouville 8,17,26,28,34
Restriction de la portée d'application: Principalement orientée vers les tours primitives; recherche supplémentaire nécessaire pour d'autres types d'extensions transcendantes
Complexité de calcul: Pour les polynômes de haut degré, le temps de calcul reste relativement long
Espace d'optimisation de l'implémentation: Les algorithmes fondamentaux comme HermiteReduce offrent encore des possibilités d'optimisation
Densité de rédaction élevée: Contenu technique dense, exigeant un solide bagage mathématique du lecteur
Analyse de complexité insuffisante: Manque d'analyse théorique détaillée de la complexité
Portée expérimentale limitée: Tests principalement sur des tours primitives à trois niveaux; performances pour des cas de dimension supérieure inconnues
L'article cite 37 références pertinentes, couvrant les domaines importants de l'intégration symbolique, de la réduction complète, du telescoping créatif et autres, fournissant une base théorique solide pour cette recherche.