Showing 1–50 of 69 results for author: Koutschan, C

Representing Piecewise-Linear Functions by Functions with Minimal Arity

Authors: Christoph Koutschan, Anton Ponomarchuk, Josef Schicho

Abstract: Any continuous piecewise-linear function $F\colon \mathbb{R}^{n}\to \mathbb{R}$ can be represented as a linear combination of $\max$ functions of at most $n+1$ affine-linear functions. In our previous paper [``Representing piecewise linear functions by functions with small arity'', AAECC, 2023], we showed that this upper bound of $n+1$ arguments is tight. In the present paper, we extend this resul… ▽ More Any continuous piecewise-linear function $F\colon \mathbb{R}^{n}\to \mathbb{R}$ can be represented as a linear combination of $\max$ functions of at most $n+1$ affine-linear functions. In our previous paper [``Representing piecewise linear functions by functions with small arity'', AAECC, 2023], we showed that this upper bound of $n+1$ arguments is tight. In the present paper, we extend this result by establishing a correspondence between the function $F$ and the minimal number of arguments that are needed in any such decomposition. We show that the tessellation of the input space $\mathbb{R}^{n}$ induced by the function $F$ has a direct connection to the number of arguments in the $\max$ functions. △ Less

Submitted 4 June, 2024; originally announced June 2024.

How to generate all possible rational Wilf-Zeilberger forms?

Authors: Shaoshi Chen, Christoph Koutschan, Yisen Wang

Abstract: Wilf-Zeilberger pairs are fundamental in the algorithmic theory of Wilf and Zeilberger for computer-generated proofs of combinatorial identities. Wilf-Zeilberger forms are their high-dimensional generalizations, which can be used for proving and discovering convergence acceleration formulas. This paper presents a structural description of all possible rational such forms, which can be viewed as an… ▽ More Wilf-Zeilberger pairs are fundamental in the algorithmic theory of Wilf and Zeilberger for computer-generated proofs of combinatorial identities. Wilf-Zeilberger forms are their high-dimensional generalizations, which can be used for proving and discovering convergence acceleration formulas. This paper presents a structural description of all possible rational such forms, which can be viewed as an additive analog of the classical Ore-Sato theorem. Based on this analog, we show a structural decomposition of so-called multivariate hyperarithmetic terms, which extend multivariate hypergeometric terms to the additive setting. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Determinant evaluations inspired by Di Francesco's determinant for twenty-vertex configurations

Authors: Christoph Koutschan, Christian Krattenthaler, Michael Schlosser

Abstract: In his work on the twenty vertex model, Di Francesco [Electron. J. Combin. 28(4) (2021), Paper No. 4.38] found a determinant formula for the number of configurations in a specific such model, and he conjectured a closed form product formula for the evaluation of this determinant. We prove this conjecture here. Moreover, we actually generalize this determinant evaluation to a one-parameter family o… ▽ More In his work on the twenty vertex model, Di Francesco [Electron. J. Combin. 28(4) (2021), Paper No. 4.38] found a determinant formula for the number of configurations in a specific such model, and he conjectured a closed form product formula for the evaluation of this determinant. We prove this conjecture here. Moreover, we actually generalize this determinant evaluation to a one-parameter family of determinant evaluations, and we present many more determinant evaluations of similar type - some proved, some left open as conjectures. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: AmS-LaTeX, 37 pages

MSC Class: Primary 15A15; Secondary 05A15; 05A19 05B45 82B20

Journal ref: J. Symbol. Comput. 127 (2025), 102352

Diagonally symmetric alternating sign matrices

Authors: Roger E. Behrend, Ilse Fischer, Christoph Koutschan

Abstract: The enumeration of diagonally symmetric alternating sign matrices (DSASMs) is studied, and a Pfaffian formula is obtained for the number of DSASMs of any fixed size, where the entries for the Pfaffian are positive integers given by simple binomial coefficient expressions. This result provides the first known case of an exact enumeration formula for an alternating sign matrix symmetry class in whic… ▽ More The enumeration of diagonally symmetric alternating sign matrices (DSASMs) is studied, and a Pfaffian formula is obtained for the number of DSASMs of any fixed size, where the entries for the Pfaffian are positive integers given by simple binomial coefficient expressions. This result provides the first known case of an exact enumeration formula for an alternating sign matrix symmetry class in which a simple product formula does not seem to exist. Pfaffian formulae are also obtained for DSASM generating functions associated with several natural statistics, including the number of nonzero strictly upper triangular entries in a DSASM, the number of nonzero diagonal entries in a DSASM, and the column of the unique 1 in the first row of a DSASM. The proofs of these results involve introducing a version of the six-vertex model whose configurations are in bijection with DSASMs of fixed size, and obtaining a Pfaffian expression for its partition function. Various further results and conjectures are also obtained, including some related to the exact enumeration of off-diagonally symmetric alternating sign matrices, and some related to the asymptotic enumeration of DSASMs and other classes of alternating sign matrices. △ Less

Submitted 15 September, 2023; originally announced September 2023.

Comments: 69 pages

MSC Class: 05A05; 05A15; 05A16; 05A19; 05B20; 15B35; 82B20; 82B23

Collisionless shock region of the KdV equation and an entry in Gradshteyn and Ryzhik

Authors: Tewodros Amdeberhan, Victor Moll, John Lopez Santander, Ken McLaughlin, Christoph Koutschan

Abstract: The long-time behavior of solutions to the initial value problem for the Korteweg-de Vries equation on the whole line, with general initial conditions has been described uniformly using five different asymptotic forms. Four of these asymptotic forms were expected: the quiescent behavior (for |x| very large), a soliton region (in which the solution behaves as a collection of isolated solitary waves… ▽ More The long-time behavior of solutions to the initial value problem for the Korteweg-de Vries equation on the whole line, with general initial conditions has been described uniformly using five different asymptotic forms. Four of these asymptotic forms were expected: the quiescent behavior (for |x| very large), a soliton region (in which the solution behaves as a collection of isolated solitary waves), a self-similar region (in which the solution is described via a Painlevé transcendent), and a similarity region (where the solution behaves as a simple trigonometric function of the quantities t and x/t). A fifth asymptotic form, lying between the self-similar (Painlevé) and the similarity one, has been described in terms of classical elliptic functions. An integral of elliptic type, giving an explicit representation of the phase, has appeared in this context. The same integral has appeared in the table of integrals by Gradshteyn and Ryzhik. Our goal here is to confirm the validity of this entry. △ Less

Submitted 6 September, 2023; originally announced September 2023.

MSC Class: 33E05 (Primary) 35C20; 35P25 (Secondary)

Yay for Determinants!

Authors: Tewodros Amdeberhan, Christoph Koutschan, Doron Zeilberger

Abstract: In this {\it case study}, we hope to show why Sheldon Axler was not just wrong, but {\em wrong}, when he urged, in 1995: ``Down with Determinants!''. We first recall how determinants are useful in enumerative combinatorics, and then illustrate three versatile tools (Dodgson's condensation, the holonomic ansatz and constant term evaluations) to operate in tandem to prove a certain intriguing determ… ▽ More In this {\it case study}, we hope to show why Sheldon Axler was not just wrong, but {\em wrong}, when he urged, in 1995: ``Down with Determinants!''. We first recall how determinants are useful in enumerative combinatorics, and then illustrate three versatile tools (Dodgson's condensation, the holonomic ansatz and constant term evaluations) to operate in tandem to prove a certain intriguing determinantal formula conjectured by the first author. We conclude with a postscript describing yet another, much more efficient, method for evaluating determinants: `ask determinant-guru, Christian Krattenthaler', but advise people only to use it as a last resort, since if we would have used this last method right away, we would not have had the fun of doing it all by ourselves. △ Less

Submitted 8 September, 2023; v1 submitted 4 July, 2023; originally announced July 2023.

Comments: 15 pages; 2 new sections added beyond the first version; dedicated to Christian Krattenthaler on his 65th birthday; this paper appeared (in a revised form) in the journal, SLC (Séminaire Lotharingien de Combinatoire), under a new title "A case study in determinant evaluations"; {\tt https://www.mat.univie.ac.at/~slc/wpapers/s89amdkozeil.pdf}

Journal ref: A case study in determinant evaluations, Séminaire Lotharingien de Combinatoire, 89 (2023), Article B89a

Representing Piecewise Linear Functions by Functions with Small Arity

Authors: Christoph Koutschan, Bernhard Moser, Anton Ponomarchuk, Josef Schicho

Abstract: A piecewise linear function can be described in different forms: as an arbitrarily nested expression of $\min$- and $\max$-functions, as a difference of two convex piecewise linear functions, or as a linear combination of maxima of affine-linear functions. In this paper, we provide two main results: first, we show that for every piecewise linear function there exists a linear combination of… ▽ More A piecewise linear function can be described in different forms: as an arbitrarily nested expression of $\min$- and $\max$-functions, as a difference of two convex piecewise linear functions, or as a linear combination of maxima of affine-linear functions. In this paper, we provide two main results: first, we show that for every piecewise linear function there exists a linear combination of $\max$-functions with at most $n+1$ arguments, and give an algorithm for its computation. Moreover, these arguments are contained in the finite set of affine-linear functions that coincide with the given function in some open set. Second, we prove that the piecewise linear function $\max(0, x_{1}, \ldots, x_{n})$ cannot be represented as a linear combination of maxima of less than $n+1$ affine-linear arguments. This was conjectured by Wang and Sun in 2005 in a paper on representations of piecewise linear functions as linear combination of maxima. △ Less

Submitted 26 May, 2023; originally announced May 2023.

Some D-finite and Some Possibly D-finite Sequences in the OEIS

Authors: Manuel Kauers, Christoph Koutschan

Abstract: In an automatic search, we found conjectural recurrences for some sequences in the OEIS that were not previously recognized as being D-finite. In some cases, we are able to prove the conjectured recurrence. In some cases, we are not able to prove the conjectured recurrence, but we can prove that a recurrence exists. In some remaining cases, we do not know where the recurrence might come from. In an automatic search, we found conjectural recurrences for some sequences in the OEIS that were not previously recognized as being D-finite. In some cases, we are able to prove the conjectured recurrence. In some cases, we are not able to prove the conjectured recurrence, but we can prove that a recurrence exists. In some remaining cases, we do not know where the recurrence might come from. △ Less

Submitted 24 April, 2023; v1 submitted 5 March, 2023; originally announced March 2023.

MSC Class: 05A15 (Primary) 68W30; 33F10 (Secondary)

Journal ref: Journal of Integer Sequences, vol. 26, article 23.4.5, 2023

On Potentials Integrated by the Nikiforov-Uvarov Method

Authors: Lina Ellis, Ikumi Ellis, Christoph Koutschan, Sergei K. Suslov

Abstract: We discuss basic potentials of the nonrelativistic and relativistic quantum mechanics that can be integrated in the Nikiforov and Uvarov paradigm with the aid of a computer algebra system. This consideration may help the readers to study analytical methods of quantum physics. We discuss basic potentials of the nonrelativistic and relativistic quantum mechanics that can be integrated in the Nikiforov and Uvarov paradigm with the aid of a computer algebra system. This consideration may help the readers to study analytical methods of quantum physics. △ Less

Submitted 16 April, 2023; v1 submitted 4 March, 2023; originally announced March 2023.

Comments: 51 pages, 19 tables; 3 figures; complementary Mathematica notebook is available from the authors; It is also posted on Wolfram community https://community.wolfram.com/groups/-/m/t/2897057

Transcendence Certificates for D-finite Functions

Authors: Manuel Kauers, Christoph Koutschan, Thibaut Verron

Abstract: Although in theory we can decide whether a given D-finite function is transcendental, transcendence proofs remain a challenge in practice. Typically, transcendence is certified by checking certain incomplete sufficient conditions. In this paper we propose an additional such condition which catches some cases on which other tests fail. Although in theory we can decide whether a given D-finite function is transcendental, transcendence proofs remain a challenge in practice. Typically, transcendence is certified by checking certain incomplete sufficient conditions. In this paper we propose an additional such condition which catches some cases on which other tests fail. △ Less

Submitted 19 September, 2023; v1 submitted 13 February, 2023; originally announced February 2023.

Comments: 9 pages, 1 figure

Journal ref: Proceedings of International Symposium on Symbolic and Algebraic Computation 2023

A Unified Approach to Unimodality of Gaussian Polynomials

Authors: Christoph Koutschan, Ali K. Uncu, Elaine Wong

Abstract: In 2013, Pak and Panova proved the strict unimodality property of $q$-binomial coefficients $\binom{\ell+m}{m}_q$ (as polynomials in $q$) based on the combinatorics of Young tableaux and the semigroup property of Kronecker coefficients. They showed it to be true for all $\ell,m\geq 8$ and a few other cases. We propose a different approach to this problem based on computer algebra, where we establi… ▽ More In 2013, Pak and Panova proved the strict unimodality property of $q$-binomial coefficients $\binom{\ell+m}{m}_q$ (as polynomials in $q$) based on the combinatorics of Young tableaux and the semigroup property of Kronecker coefficients. They showed it to be true for all $\ell,m\geq 8$ and a few other cases. We propose a different approach to this problem based on computer algebra, where we establish a closed form for the coefficients of these polynomials and then use cylindrical algebraic decomposition to identify exactly the range of coefficients where strict unimodality holds. This strategy allows us to tackle generalizations of the problem, e.g., to show unimodality with larger gaps or unimodality of related sequences. In particular, we present proofs of two additional cases of a conjecture by Stanley and Zanello. △ Less

Submitted 31 August, 2023; v1 submitted 8 February, 2023; originally announced February 2023.

Comments: Supplementary material at https://wongey.github.io/unimodality

Journal ref: ISSAC 2023: Proceedings of the 2023 International Symposium on Symbolic and Algebraic Computation, July 2023, Pages 434-442

How does the Gerrymander Sequence Continue?

Authors: Manuel Kauers, Christoph Koutschan, George Spahn

Abstract: We compute the next few terms of the OEIS sequence A348456 and provide guessed equations for the generating functions of some sequences in its context. We compute the next few terms of the OEIS sequence A348456 and provide guessed equations for the generating functions of some sequences in its context. △ Less

Submitted 30 November, 2022; v1 submitted 5 September, 2022; originally announced September 2022.

Journal ref: Journal of Integer Sequences Vol 25 (2022) Article 22.9.7

Guessing with Little Data

Authors: Manuel Kauers, Christoph Koutschan

Abstract: Reconstructing a hypothetical recurrence equation from the first terms of an infinite sequence is a classical and well-known technique in experimental mathematics. We propose a variation of this technique which can succeed with fewer input terms. Reconstructing a hypothetical recurrence equation from the first terms of an infinite sequence is a classical and well-known technique in experimental mathematics. We propose a variation of this technique which can succeed with fewer input terms. △ Less

Submitted 8 July, 2022; v1 submitted 16 February, 2022; originally announced February 2022.

Journal ref: Proceedings of ISSAC 2022, pp 83--90

Realizations of Rigid Graphs

Authors: Christoph Koutschan

Abstract: A minimally rigid graph, also called Laman graph, models a planar framework which is rigid for a general choice of distances between its vertices. In other words, there are finitely many ways, up to isometries, to realize such a graph in the plane. Using ideas from algebraic and tropical geometry, we derive a recursive formula for the number of such realizations. Combining computational results wi… ▽ More A minimally rigid graph, also called Laman graph, models a planar framework which is rigid for a general choice of distances between its vertices. In other words, there are finitely many ways, up to isometries, to realize such a graph in the plane. Using ideas from algebraic and tropical geometry, we derive a recursive formula for the number of such realizations. Combining computational results with the construction of new rigid graphs via gluing techniques, we can give a new lower bound on the maximal possible number of realizations for graphs with a given number of vertices. △ Less

Submitted 3 January, 2022; originally announced January 2022.

Comments: In Proceedings ADG 2021, arXiv:2112.14770

Journal ref: EPTCS 352, 2021, pp. 4-13

Apéry limits for elliptic $L$-values

Authors: Christoph Koutschan, Wadim Zudilin

Abstract: For an (irreducible) recurrence equation with coefficients from $\mathbb Z[n]$ and its two linearly independent rational solutions $u_n,v_n$, the limit of $u_n/v_n$ as $n\to\infty$, when exists, is called the Apéry limit. We give a construction that realises certain quotients of $L$-values of elliptic curves as Apéry limits. For an (irreducible) recurrence equation with coefficients from $\mathbb Z[n]$ and its two linearly independent rational solutions $u_n,v_n$, the limit of $u_n/v_n$ as $n\to\infty$, when exists, is called the Apéry limit. We give a construction that realises certain quotients of $L$-values of elliptic curves as Apéry limits. △ Less

Submitted 16 November, 2021; originally announced November 2021.

Comments: 6 pages, 3 recurrence equations (including Apéry's for $ζぜーた(3)$)

Report number: RICAM Report 2021-34 MSC Class: Primary 11F67; Secondary 11G05; 11G40; 11J70; 11R06; 14K20; 33F10; 39A06

Journal ref: Bull. Austral. Math. Soc. 106 (2022), no. 2, 273--279

Legendre pairs of lengths $\ell\equiv0$ (mod 5)

Authors: Ilias Kotsireas, Christoph Koutschan, Dursun Bulutoglu, David Arquette, Jonathan Turner, Kenneth Ryan

Abstract: By assuming a type of balance for length $\ell=87$ and non-trivial subgroups of multiplier groups of Legendre pairs (LPs) for length $\ell=85$, we find LPs of these lengths. We then study the power spectral density (PSD) values of m-compressions of LPs of length 5m. We also formulate a conjecture for Legendre pairs of lengths $\ell \equiv 0$ (mod 5) and demonstrate how it can be used to decrease t… ▽ More By assuming a type of balance for length $\ell=87$ and non-trivial subgroups of multiplier groups of Legendre pairs (LPs) for length $\ell=85$, we find LPs of these lengths. We then study the power spectral density (PSD) values of m-compressions of LPs of length 5m. We also formulate a conjecture for Legendre pairs of lengths $\ell \equiv 0$ (mod 5) and demonstrate how it can be used to decrease the search space and storage requirements for finding such LPs. The newly found LPs decrease the number of integers in the range $\leq 200$ for which the existence question of LPs remains unsolved from 12 to 10. △ Less

Submitted 8 April, 2024; v1 submitted 3 November, 2021; originally announced November 2021.

MSC Class: 05B10; 05B20; 15B34

Journal ref: Special Matrices 2023; 11: 20230105

Binomial Determinants for Tiling Problems Yield to the Holonomic Ansatz

Authors: Hao Du, Christoph Koutschan, Thotsaporn Thanatipanonda, Elaine Wong

Abstract: We present and prove closed form expressions for some families of binomial determinants with signed Kronecker deltas that are located along an arbitrary diagonal in the corresponding matrix. They count cyclically symmetric rhombus tilings of hexagonal regions with triangular holes. We extend a previous systematic study of these families, where the locations of the Kronecker deltas depended on an a… ▽ More We present and prove closed form expressions for some families of binomial determinants with signed Kronecker deltas that are located along an arbitrary diagonal in the corresponding matrix. They count cyclically symmetric rhombus tilings of hexagonal regions with triangular holes. We extend a previous systematic study of these families, where the locations of the Kronecker deltas depended on an additional parameter, to families with negative Kronecker deltas. By adapting Zeilberger's holonomic ansatz to make it work for our problems, we can take full advantage of computer algebra tools for symbolic summation. This, together with the combinatorial interpretation, allows us to realize some new determinantal relationships. From there, we are able to resolve all remaining open conjectures related to these determinants, including one from 2005 due to Lascoux and Krattenthaler. △ Less

Submitted 21 September, 2021; v1 submitted 18 May, 2021; originally announced May 2021.

Comments: 45 pages; Supplementary material at https://wongey.github.io/binom-det

Journal ref: European Journal of Combinatorics, Volume 99, January 2022, 103437

There are EXACTLY 1493804444499093354916284290188948031229880469556 Ways to Derange a Standard Deck of Cards (ignoring suits) [and many other such useful facts]

Authors: Shalosh B. Ekhad, Christoph Koutschan, Doron Zeilberger

Abstract: In this memorial tribute to Joe Gillis, who taught us that Special Functions count, we show how the seminal Even-Gillis integral formula for the number of derangements of a multiset, in terms of Laguerre polynomials, can be used to efficiently compute not only the number of the title, but much harder ones, when it is interfaced with Wilf-Zeilberger algorithmic proof theory. In this memorial tribute to Joe Gillis, who taught us that Special Functions count, we show how the seminal Even-Gillis integral formula for the number of derangements of a multiset, in terms of Laguerre polynomials, can be used to efficiently compute not only the number of the title, but much harder ones, when it is interfaced with Wilf-Zeilberger algorithmic proof theory. △ Less

Submitted 25 January, 2021; originally announced January 2021.

Comments: Accompanied by a Maple package and output files that can be gotten from https://sites.math.rutgers.edu/~zeilberg/mamarim/mamarimhtml/multider.html

MSC Class: 33F10; 05A15; 33D45; 68W30

Tweaking the Beukers Integrals In Search of More Miraculous Irrationality Proofs A La Apery

Authors: Robert Dougherty-Bliss, Christoph Koutschan, Doron Zeilberger

Abstract: There are only aleph-zero rational numbers, while there are 2 to the power aleph-zero real numbers. Hence the probability that a randomly chosen real number would be rational is 0. Yet proving rigorously that any specific, natural, real constant, is irrational is usually very hard, witness that there are still no proofs of the irrationality of the Euler-Mascheroni constant, the Catalan constant, o… ▽ More There are only aleph-zero rational numbers, while there are 2 to the power aleph-zero real numbers. Hence the probability that a randomly chosen real number would be rational is 0. Yet proving rigorously that any specific, natural, real constant, is irrational is usually very hard, witness that there are still no proofs of the irrationality of the Euler-Mascheroni constant, the Catalan constant, or Zeta(5). Inspired by Frits Beukers' elegant rendition of Apery's seminal proofs of the irrationality of Zeta(2) and Zeta(3), and heavily using algorithmic proof theory, we systematically searched for other similar integrals, that lead to irrationality proofs. We found quite a few candidates for such proofs, including the square-root of Pi times Gamma(7/3)/Gamma(-1/6) and Gamma(19/6)/Gamma(8/3) divided by the square-root of Pi. △ Less

Submitted 20 January, 2021; originally announced January 2021.

Comments: 21 pages, accompanied by three Maple packages available from https://sites.math.rutgers.edu/~zeilberg/mamarim/mamarimhtml/beukers.html

Legendre pairs of lengths $\ell \equiv 0$ (mod $3$)

Authors: Ilias Kotsireas, Christoph Koutschan

Abstract: We prove a proposition that connects constant-PAF sequences and the corresponding Legendre pairs with integer PSD values. We show how to determine explicitly the complete spectrum of the $(\ell/3)$-rd value of the discrete Fourier transform for Legendre pairs of lengths $\ell \equiv 0 \, (\mod 3)$. This is accomplished by two new algorithms based on number-theoretic arguments. As an application, w… ▽ More We prove a proposition that connects constant-PAF sequences and the corresponding Legendre pairs with integer PSD values. We show how to determine explicitly the complete spectrum of the $(\ell/3)$-rd value of the discrete Fourier transform for Legendre pairs of lengths $\ell \equiv 0 \, (\mod 3)$. This is accomplished by two new algorithms based on number-theoretic arguments. As an application, we prove that Legendre pairs of the open lengths 117, 129, 133, and 147 exist by finding Legendre pairs of these lengths with a multiplier group of order at least 3. As a consequence, 85, 87, 115, 145, 159, 161, 169, 175, 177, 185, 187, 195 are the twelve integers in the range < 200 for which the question of existence of Legendre pairs remains unsolved. △ Less

Submitted 16 October, 2021; v1 submitted 8 January, 2021; originally announced January 2021.

Journal ref: Journal of Combinatorial Designs 29(12):870-887, 2021

Creative Telescoping on Multiple Sums

Authors: Christoph Koutschan, Elaine Wong

Abstract: We showcase a collection of practical strategies to deal with a problem arising from an analysis of integral estimators derived via quasi-Monte Carlo methods. The problem reduces to a triple binomial sum, thereby enabling us to open up the holonomic toolkit, which contains tools such as creative telescoping that can be used to deduce a recurrence satisfied by the sum. While applying these techniqu… ▽ More We showcase a collection of practical strategies to deal with a problem arising from an analysis of integral estimators derived via quasi-Monte Carlo methods. The problem reduces to a triple binomial sum, thereby enabling us to open up the holonomic toolkit, which contains tools such as creative telescoping that can be used to deduce a recurrence satisfied by the sum. While applying these techniques, a host of issues arose that partly needed to be resolved by hand. In other words, no creative telescoping implementation currently exists that can resolve all these issues automatically. Thus, we felt the need to compile the different strategies we tried and the difficulties that we encountered along the way. In particular, we highlight the necessity of the certificate in these computations and how its complexity can greatly influence the computation time. △ Less

Submitted 17 March, 2021; v1 submitted 17 October, 2020; originally announced October 2020.

Comments: 22 pages; Supplementary material at https://wongey.github.io/digital-nets-walsh/

Journal ref: Mathematics in Computer Science, Vol. 15(3), Pages 483-498 (2021)

Common Factors in Fraction-Free Matrix Decompositions

Authors: Johannes Middeke, David J. Jeffrey, Christoph Koutschan

Abstract: We consider LU and QR matrix decompositions using exact computations. We show that fraction-free Gauss--Bareiss reduction leads to triangular matrices having a non-trivial number of common row factors. We identify two types of common factors: systematic and statistical. Systematic factors depend on the reduction process, independent of the data, while statistical factors depend on the specific dat… ▽ More We consider LU and QR matrix decompositions using exact computations. We show that fraction-free Gauss--Bareiss reduction leads to triangular matrices having a non-trivial number of common row factors. We identify two types of common factors: systematic and statistical. Systematic factors depend on the reduction process, independent of the data, while statistical factors depend on the specific data. We relate the existence of row factors in the LU decomposition to factors appearing in the Smith--Jacobson normal form of the matrix. For statistical factors, we identify some of the mechanisms that create them and give estimates of the frequency of their occurrence. Similar observations apply to the common factors in a fraction-free QR decomposition. Our conclusions are tested experimentally. △ Less

Submitted 25 May, 2020; originally announced May 2020.

MSC Class: 00-01; 99-00

Journal ref: Mathematics in Computer Science 15(4), pp. 589-608, 2021

Diagonals of rational functions: from differential algebra to effective algebraic geometry (unabridged version)

Authors: Y. Abdelaziz, S. Boukraa, C. Koutschan, J-M. Maillard

Abstract: We show that the results we had obtained on diagonals of nine and ten parameters families of rational functions using creative telescoping, yielding modular forms expressed as pullbacked $ _2F_1$ hypergeometric functions, can be obtained, much more efficiently, calculating the $ j$-invariant of an elliptic curve canonically associated with the denominator of the rational functions. In the case whe… ▽ More We show that the results we had obtained on diagonals of nine and ten parameters families of rational functions using creative telescoping, yielding modular forms expressed as pullbacked $ _2F_1$ hypergeometric functions, can be obtained, much more efficiently, calculating the $ j$-invariant of an elliptic curve canonically associated with the denominator of the rational functions. In the case where creative telescoping yields pullbacked $ _2F_1$ hypergeometric functions, we generalize this result to other families of rational functions in three, and even more than three, variables. We also generalise this result to rational functions in more than three variables when the denominator can be associated to an algebraic variety corresponding to products of elliptic curves, foliation in elliptic curves. We also extend these results to rational functions in three variables when the denominator is associated with a {\em genus-two curve such that its Jacobian is a split Jacobian} corresponding to the product of two elliptic curves. We sketch the situation where the denominator of the rational function is associated with algebraic varieties that are not of the general type, having an infinite set of birational automorphisms. We finally provide some examples of rational functions in more than three variables, where the telescopers have pullbacked $ _2F_1$ hypergeometric solutions, the denominator corresponding to an algebraic variety being not simply foliated in elliptic curves, but having a selected elliptic curve in the variety explaining the pullbacked $ _2F_1$ hypergeometric solution. △ Less

Submitted 26 April, 2022; v1 submitted 3 February, 2020; originally announced February 2020.

Comments: 51 pages

MSC Class: 34M55; 47E05; 81Qxx; 32G34; 34Lxx; 34Mxx; 14Kxx

Calculation and Properties of Zonal Polynomials

Authors: Lin Jiu, Christoph Koutschan

Abstract: We investigate the zonal polynomials, a family of symmetric polynomials that appear in many mathematical contexts, such as multivariate statistics, differential geometry, representation theory, and combinatorics. We present two computer algebra packages, in SageMath and in Mathematica, for their computation. With the help of these software packages, we carry out an experimental mathematics study o… ▽ More We investigate the zonal polynomials, a family of symmetric polynomials that appear in many mathematical contexts, such as multivariate statistics, differential geometry, representation theory, and combinatorics. We present two computer algebra packages, in SageMath and in Mathematica, for their computation. With the help of these software packages, we carry out an experimental mathematics study of some properties of zonal polynomials. Moreover, we derive and prove closed forms for several infinite families of zonal polynomial coefficients. △ Less

Submitted 29 April, 2020; v1 submitted 30 January, 2020; originally announced January 2020.

MSC Class: Primary 05E05; Secondary 33C20 33C70 15B52 65C60

Journal ref: Mathematics in Computer Science 14:623-640, 2020

On Christol's conjecture

Authors: Y. Abdelaziz, C. Koutschan, J-M. Maillard

Abstract: We show that the unresolved examples of Christol's conjecture $ \, _3F_{2}\left([2/9,5/9,8/9],[2/3,1],x\right)$ and $_3F_{2}\left([1/9,4/9,7/9],[1/3,1],x\right)$, are indeed diagonals of rational functions. We also show that other $\, _3F_2$ and $\, _4F_3$ unresolved examples of Christol's conjecture are diagonals of rational functions. Finally we give two arguments that show that it is likely… ▽ More We show that the unresolved examples of Christol's conjecture $ \, _3F_{2}\left([2/9,5/9,8/9],[2/3,1],x\right)$ and $_3F_{2}\left([1/9,4/9,7/9],[1/3,1],x\right)$, are indeed diagonals of rational functions. We also show that other $\, _3F_2$ and $\, _4F_3$ unresolved examples of Christol's conjecture are diagonals of rational functions. Finally we give two arguments that show that it is likely that the $\, _3F_2([1/9, 4/9, 5/9], \, [1/3,1], \, 27 \cdot x)$ function is a diagonal of a rational function. △ Less

Submitted 1 January, 2020; v1 submitted 21 December, 2019; originally announced December 2019.

Comments: 13 pages

Heun functions and diagonals of rational functions (unabridged version)

Authors: Y. Abdelaziz, S. Boukraa, C. Koutschan, J-M. Maillard

Abstract: We provide a set of diagonals of simple rational functions of three and four variables that are squares of Heun functions. These Heun functions obtained through creative telescoping, turn out to be either pullbacked $_2F_1$ hypergeometric functions and in fact classical modular forms. We also obtain Heun functions that are Shimura curves as solutions of telescopers of rational functions. We provide a set of diagonals of simple rational functions of three and four variables that are squares of Heun functions. These Heun functions obtained through creative telescoping, turn out to be either pullbacked $_2F_1$ hypergeometric functions and in fact classical modular forms. We also obtain Heun functions that are Shimura curves as solutions of telescopers of rational functions. △ Less

Submitted 23 October, 2019; originally announced October 2019.

Comments: 64 pages

MSC Class: 34M55; 47E05; 81Qxx; 32G34; 34Lxx; 34Mxx; 14Kxx

Exact Lower Bounds for Monochromatic Schur Triples and Generalizations

Authors: Christoph Koutschan, Elaine Wong

Abstract: We derive exact and sharp lower bounds for the number of monochromatic generalized Schur triples $(x,y,x+ay)$ whose entries are from the set $\{1,\dots,n\}$, subject to a coloring with two different colors. Previously, only asymptotic formulas for such bounds were known, and only for $a\in\mathbb{N}$. Using symbolic computation techniques, these results are extended here to arbitrary… ▽ More We derive exact and sharp lower bounds for the number of monochromatic generalized Schur triples $(x,y,x+ay)$ whose entries are from the set $\{1,\dots,n\}$, subject to a coloring with two different colors. Previously, only asymptotic formulas for such bounds were known, and only for $a\in\mathbb{N}$. Using symbolic computation techniques, these results are extended here to arbitrary $a\in\mathbb{R}$. Furthermore, we give exact formulas for the minimum number of monochromatic Schur triples for $a=1,2,3,4$, and briefly discuss the case $0<a<1$. △ Less

Submitted 12 October, 2020; v1 submitted 3 April, 2019; originally announced April 2019.

Comments: 24 Pages

Journal ref: Algorithmic Combinatorics: Enumerative Combinatorics, Special Functions and Computer Algebra (Veronika Pillwein, Carsten Schneider, eds.), Texts & Monographs in Symbolic Computation, pages 223-248, Springer, 2020. ISBN 978-3-030-44558-4

On the singular value decomposition of n-fold integration operators

Authors: Ronny Ramlau, Christoph Koutschan, Bernd Hofmann

Abstract: In theory and practice of inverse problems, linear operator equations $Tx=y$ with compact linear forward operators $T$ having a non-closed range $\mathcal{R}(T)$ and mapping between infinite dimensional Hilbert spaces plays some prominent role. As a consequence of the ill-posedness of such problems, regularization approaches are required, and due to its unlimited qualification spectral cut-off is… ▽ More In theory and practice of inverse problems, linear operator equations $Tx=y$ with compact linear forward operators $T$ having a non-closed range $\mathcal{R}(T)$ and mapping between infinite dimensional Hilbert spaces plays some prominent role. As a consequence of the ill-posedness of such problems, regularization approaches are required, and due to its unlimited qualification spectral cut-off is an appropriate method for the stable approximate solution of corresponding inverse problems. For this method, however, the singular system $\{σしぐま_i(T),u_i(T),v_i(T)\}_{i=1}^\infty$ of the compact operator $T$ is needed, at least for $i=1,2,...,N$, up to some stopping index $N$. In this note we consider $n$-fold integration operators $T=J^n\;(n=1,2,...)$ in $L^2([0,1])$ occurring in numerous applications, where the solution of the associated operator equation is characterized by the $n$-th generalized derivative $x=y^{(n)}$ of the Sobolev space function $y \in H^n([0,1])$. Almost all textbooks on linear inverse problems present the whole singular system $\{σしぐま_i(J^1),u_i(J^1),v_i(J^1)\}_{i=1}^\infty$ in an explicit manner. However, they do not discuss the singular systems for $J^n,\;n \ge 2$. We will emphasize that this seems to be a consequence of the fact that for higher $n$ the eigenvalues $σしぐま^2_i(J^n)$ of the associated ODE boundary value problems obey transcendental equations, the complexity of which is growing with $n$. We present the transcendental equations for $n=2,3,...$ and discuss and illustrate the associated eigenfunctions and some of their properties. △ Less

Submitted 28 November, 2018; originally announced November 2018.

Report number: RICAM Reports Series 2018-29

Journal ref: In: Jin Cheng, Masahiro Yamamoto, Shuai Lu (eds.), Inverse Problems and Related Topics (2020), vol. 310 of Springer Proceedings in Mathematics & Statistics, pages 237-256

Diagonals of rational functions, pullbacked 2F1 hypergeometric functions and modular forms (unabrigded version)

Authors: Y. Abdelaziz, S. Boukraa, C. Koutschan, J-M. Maillard

Abstract: We recall that diagonals of rational functions naturally occur in lattice statistical mechanics and enumerative combinatorics. We find that a seven-parameter rational function of three variables with a numerator equal to one (reciprocal of a polynomial of degree two at most) can be expressed as a pullbacked 2F1 hypergeometric function. This result can be seen as the simplest non-trivial family of… ▽ More We recall that diagonals of rational functions naturally occur in lattice statistical mechanics and enumerative combinatorics. We find that a seven-parameter rational function of three variables with a numerator equal to one (reciprocal of a polynomial of degree two at most) can be expressed as a pullbacked 2F1 hypergeometric function. This result can be seen as the simplest non-trivial family of diagonals of rational functions. We focus on some subcases such that the diagonals of the corresponding rational functions can be written as a pullbacked 2F1 hypergeometric function with two possible rational functions pullbacks algebraically related by modular equations, thus showing explicitely that the diagonal is a modular form. We then generalise this result to eight, nine and ten parameters families adding some selected cubic terms at the denominator of the rational function defining the diagonal. We finally show that each of these previous rational functions yields an infinite number of rational functions whose diagonals are also pullbacked 2F1 hypergeometric functions and modular forms. △ Less

Submitted 12 May, 2018; originally announced May 2018.

Comments: 39 pages

MSC Class: 47E05; 81Qxx; 32G34; 34Lxx; 34Mxx; 14Kxx

Journal ref: Journal of Physics A: Mathematical and Theoretical, vol. 51, article 455201, 2018

Echelons of power series and Gabrielov's counterexample to nested linear Artin Approximation

Authors: M. E. Alonso, F. J. Castro-Jiménez, H. Hauser, C. Koutschan

Abstract: Gabrielov's famous example for the failure of analytic Artin approximation in the presence of nested subring conditions is shown to be due to a growth phenomenon in standard basis computations for echelons, a generalization of the concept of ideals in power series rings. Gabrielov's famous example for the failure of analytic Artin approximation in the presence of nested subring conditions is shown to be due to a growth phenomenon in standard basis computations for echelons, a generalization of the concept of ideals in power series rings. △ Less

Submitted 22 April, 2018; originally announced April 2018.

Comments: To appear in Bulletin of the London Mathematical Society

MSC Class: 32A05; 13F25; 13P10; 16W60; 14QXX; 14B12

Desingularization in the $q$-Weyl algebra

Authors: Christoph Koutschan, Yi Zhang

Abstract: In this paper, we study the desingularization problem in the first $q$-Weyl algebra. We give an order bound for desingularized operators, and thus derive an algorithm for computing desingularized operators in the first $q$-Weyl algebra. Moreover, an algorithm is presented for computing a generating set of the first $q$-Weyl closure of a given $q$-difference operator. As an application, we certify… ▽ More In this paper, we study the desingularization problem in the first $q$-Weyl algebra. We give an order bound for desingularized operators, and thus derive an algorithm for computing desingularized operators in the first $q$-Weyl algebra. Moreover, an algorithm is presented for computing a generating set of the first $q$-Weyl closure of a given $q$-difference operator. As an application, we certify that several instances of the colored Jones polynomial are Laurent polynomial sequences by computing the corresponding desingularized operator. △ Less

Submitted 25 February, 2018; v1 submitted 12 January, 2018; originally announced January 2018.

Journal ref: Advances in Applied Mathematics, 97:80-101, 2018

Lower bounds on the number of realizations of rigid graphs

Authors: Georg Grasegger, Christoph Koutschan, Elias Tsigaridas

Abstract: Computing the number of realizations of a minimally rigid graph is a notoriously difficult problem. Towards this goal, for graphs that are minimally rigid in the plane, we take advantage of a recently published algorithm, which is the fastest available method, although its complexity is still exponential. Combining computational results with the theory of constructing new rigid graphs by gluing, w… ▽ More Computing the number of realizations of a minimally rigid graph is a notoriously difficult problem. Towards this goal, for graphs that are minimally rigid in the plane, we take advantage of a recently published algorithm, which is the fastest available method, although its complexity is still exponential. Combining computational results with the theory of constructing new rigid graphs by gluing, we give a new lower bound on the maximal possible number of (complex) realizations for graphs with a given number of vertices. We extend these ideas to rigid graphs in three dimensions and we derive similar lower bounds, by exploiting data from extensive Gröbner basis computations. △ Less

Submitted 11 April, 2018; v1 submitted 23 October, 2017; originally announced October 2017.

MSC Class: 52C25; 05C04; 68W30

Journal ref: Experimental Mathematics, 2018

A Curious Family of Binomial Determinants That Count Rhombus Tilings of a Holey Hexagon

Authors: Christoph Koutschan, Thotsaporn Thanatipanonda

Abstract: We evaluate a curious determinant, first mentioned by George Andrews in 1980 in the context of descending plane partitions. Our strategy is to combine the famous Desnanot-Jacobi-Dodgson identity with automated proof techniques. More precisely, we follow the holonomic ansatz that was proposed by Doron Zeilberger in 2007. We derive a compact and nice formula for Andrews's determinant, and use it to… ▽ More We evaluate a curious determinant, first mentioned by George Andrews in 1980 in the context of descending plane partitions. Our strategy is to combine the famous Desnanot-Jacobi-Dodgson identity with automated proof techniques. More precisely, we follow the holonomic ansatz that was proposed by Doron Zeilberger in 2007. We derive a compact and nice formula for Andrews's determinant, and use it to solve a challenge problem that we posed in a previous paper. By noting that Andrews's determinant is a special case of a two-parameter family of determinants, we find closed forms for several one-parameter subfamilies. The interest in these determinants arises because they count cyclically symmetric rhombus tilings of a hexagon with several triangular holes inside. △ Less

Submitted 3 July, 2018; v1 submitted 8 September, 2017; originally announced September 2017.

Journal ref: Journal of Combinatorial Theory, Series A, vol. 166, pp. 352-381, 2019

Computing the number of realizations of a Laman graph

Authors: Jose Capco, Matteo Gallet, Georg Grasegger, Christoph Koutschan, Niels Lubbes, Josef Schicho

Abstract: Laman graphs model planar frameworks which are rigid for a general choice of distances between the vertices. There are finitely many ways, up to isometries, to realize a Laman graph in the plane. In a recent paper we provide a recursion formula for this number of realizations using ideas from algebraic and tropical geometry. Here, we present a concise summary of this result focusing on the main id… ▽ More Laman graphs model planar frameworks which are rigid for a general choice of distances between the vertices. There are finitely many ways, up to isometries, to realize a Laman graph in the plane. In a recent paper we provide a recursion formula for this number of realizations using ideas from algebraic and tropical geometry. Here, we present a concise summary of this result focusing on the main ideas and the combinatorial point of view. △ Less

Submitted 11 October, 2017; v1 submitted 12 July, 2017; originally announced July 2017.

Comments: Extended abstract; the long version is arxiv:1701.05500

MSC Class: 14T05; 14N99; 52C25; 05C99

Journal ref: Electronic Notes in Discrete Mathematics (Proceedings of Eurocomb 2017) 61, pp. 207-213, 2017

The number of realizations of a Laman graph

Authors: Jose Capco, Matteo Gallet, Georg Grasegger, Christoph Koutschan, Niels Lubbes, Josef Schicho

Abstract: Laman graphs model planar frameworks that are rigid for a general choice of distances between the vertices. There are finitely many ways, up to isometries, to realize a Laman graph in the plane. Such realizations can be seen as solutions of systems of quadratic equations prescribing the distances between pairs of points. Using ideas from algebraic and tropical geometry, we provide a recursive form… ▽ More Laman graphs model planar frameworks that are rigid for a general choice of distances between the vertices. There are finitely many ways, up to isometries, to realize a Laman graph in the plane. Such realizations can be seen as solutions of systems of quadratic equations prescribing the distances between pairs of points. Using ideas from algebraic and tropical geometry, we provide a recursive formula for the number of complex solutions of such systems. △ Less

Submitted 28 November, 2017; v1 submitted 19 January, 2017; originally announced January 2017.

Comments: 36 pages

MSC Class: 14T05; 14N99; 52C25; 05C99

Journal ref: SIAM Journal on Applied Algebra and Geometry, Volume 2, Issue 1, p. 94-125, 2018

Reduction-Based Creative Telescoping for Fuchsian D-finite Functions

Authors: Shaoshi Chen, Mark van Hoeij, Manuel Kauers, Christoph Koutschan

Abstract: Continuing a series of articles in the past few years on creative telescoping using reductions, we adapt Trager's Hermite reduction for algebraic functions to fuchsian D-finite functions and develop a reduction-based creative telescoping algorithm for this class of functions, thereby generalizing our recent reduction-based algorithm for algebraic functions, presented at ISSAC 2016. Continuing a series of articles in the past few years on creative telescoping using reductions, we adapt Trager's Hermite reduction for algebraic functions to fuchsian D-finite functions and develop a reduction-based creative telescoping algorithm for this class of functions, thereby generalizing our recent reduction-based algorithm for algebraic functions, presented at ISSAC 2016. △ Less

Submitted 22 November, 2016; originally announced November 2016.

Comments: arXiv admin note: text overlap with arXiv:1602.00424

Inverse Inequality Estimates with Symbolic Computation

Authors: Christoph Koutschan, Martin Neumüller, Cristian-Silviu Radu

Abstract: In the convergence analysis of numerical methods for solving partial differential equations (such as finite element methods) one arrives at certain generalized eigenvalue problems, whose maximal eigenvalues need to be estimated as accurately as possible. We apply symbolic computation methods to the situation of square elements and are able to improve the previously known upper bound, given in "p-… ▽ More In the convergence analysis of numerical methods for solving partial differential equations (such as finite element methods) one arrives at certain generalized eigenvalue problems, whose maximal eigenvalues need to be estimated as accurately as possible. We apply symbolic computation methods to the situation of square elements and are able to improve the previously known upper bound, given in "p- and hp-finite element methods" (Schwab, 1998), by a factor of 8. More precisely, we try to evaluate the corresponding determinant using the holonomic ansatz, which is a powerful tool for dealing with determinants, proposed by Zeilberger in 2007. However, it turns out that this method does not succeed on the problem at hand. As a solution we present a variation of the original holonomic ansatz that is applicable to a larger class of determinants, including the one we are dealing with here. We obtain an explicit closed form for the determinant, whose special form enables us to derive new and tight upper resp. lower bounds on the maximal eigenvalue, as well as its asymptotic behaviour. △ Less

Submitted 3 May, 2016; v1 submitted 3 February, 2016; originally announced February 2016.

MSC Class: 33F10; 65N12 (Primary) 65N30; 68W30; 65F15; 05A20; 15A15; 15A45 (Secondary)

Journal ref: Advances in Applied Mathematics 80:1-23, 2016

Holonomic Tools for Basic Hypergeometric Functions

Authors: Christoph Koutschan, Peter Paule

Abstract: With the exception of q-hypergeometric summation, the use of computer algebra packages implementing Zeilberger's "holonomic systems approach" in a broader mathematical sense is less common in the field of q-series and basic hypergeometric functions. A major objective of this article is to popularize the usage of such tools also in these domains. Concrete case studies showing software in action int… ▽ More With the exception of q-hypergeometric summation, the use of computer algebra packages implementing Zeilberger's "holonomic systems approach" in a broader mathematical sense is less common in the field of q-series and basic hypergeometric functions. A major objective of this article is to popularize the usage of such tools also in these domains. Concrete case studies showing software in action introduce to the basic techniques. An application highlight is a new computer-assisted proof of the celebrated Ismail-Zhang formula, an important q-analog of a classical expansion formula of plane waves in terms of Gegenbauer polynomials. △ Less

Submitted 1 February, 2016; originally announced February 2016.

Reduction-Based Creative Telescoping for Algebraic Functions

Authors: Shaoshi Chen, Manuel Kauers, Christoph Koutschan

Abstract: Continuing a series of articles in the past few years on creative telescoping using reductions, we develop a new algorithm to construct minimal telescopers for algebraic functions. This algorithm is based on Trager's Hermite reduction and on polynomial reduction, which was originally designed for hyperexponential functions and extended to the algebraic case in this paper. Continuing a series of articles in the past few years on creative telescoping using reductions, we develop a new algorithm to construct minimal telescopers for algebraic functions. This algorithm is based on Trager's Hermite reduction and on polynomial reduction, which was originally designed for hyperexponential functions and extended to the algebraic case in this paper. △ Less

Submitted 1 February, 2016; originally announced February 2016.

ACM Class: I.1.2

Lattice Green Functions: the d-dimensional face-centred cubic lattice, d=8, 9, 10, 11, 12

Authors: S. Hassani, C. Koutschan, J-M. Maillard, N. Zenine

Abstract: We previously reported on a recursive method to generate the expansion of the lattice Green function of the $d$-dimensional face-centred cubic lattice (fcc). The method was used to generate many coefficients for d=7 and the corresponding linear differential equation has been obtained. In this paper, we show the strength and the limit of the method by producing the series and the corresponding line… ▽ More We previously reported on a recursive method to generate the expansion of the lattice Green function of the $d$-dimensional face-centred cubic lattice (fcc). The method was used to generate many coefficients for d=7 and the corresponding linear differential equation has been obtained. In this paper, we show the strength and the limit of the method by producing the series and the corresponding linear differential equations for d=8, 9, 10, 11, 12. The differential Galois groups of these linear differential equations are shown to be symplectic for d=8, 10, 12 and orthogonal for d= 9, 11. The recursion relation naturally provides a 2-dimensional array $ T_d(n,j)$ where only the coefficients $ t_d(n,0)$ correspond to the coefficients of the lattice Green function of the d-dimensional fcc. The coefficients $ t_d(n,j)$ are associated to D-finite bivariate series annihilated by linear partial differential equations that we analyze. △ Less

Submitted 6 February, 2016; v1 submitted 21 January, 2016; originally announced January 2016.

Comments: 28 pages

MSC Class: 34M55; 47E05; 81Qxx; 32G34; 34Lxx; 34Mxx; 14Kxx

Exact ZF Analysis and Computer-Algebra-Aided Evaluation in Rank-1 LoS Rician Fading

Authors: Constantin Siriteanu, Akimichi Takemura, Christoph Koutschan, Satoshi Kuriki, Donald St. P. Richards, Hyundong Shin

Abstract: We study zero-forcing detection (ZF) for multiple-input/multiple-output (MIMO) spatial multiplexing under transmit-correlated Rician fading for an N_R X N_T channel matrix with rank-1 line-of-sight (LoS) component. By using matrix transformations and multivariate statistics, our exact analysis yields the signal-to-noise ratio moment generating function (m.g.f.) as an infinite series of gamma distr… ▽ More We study zero-forcing detection (ZF) for multiple-input/multiple-output (MIMO) spatial multiplexing under transmit-correlated Rician fading for an N_R X N_T channel matrix with rank-1 line-of-sight (LoS) component. By using matrix transformations and multivariate statistics, our exact analysis yields the signal-to-noise ratio moment generating function (m.g.f.) as an infinite series of gamma distribution m.g.f.'s and analogous series for ZF performance measures, e.g., outage probability and ergodic capacity. However, their numerical convergence is inherently problematic with increasing Rician K-factor, N_R , and N_T. We circumvent this limitation as follows. First, we derive differential equations satisfied by the performance measures with a novel automated approach employing a computer-algebra tool which implements Groebner basis computation and creative telescoping. These differential equations are then solved with the holonomic gradient method (HGM) from initial conditions computed with the infinite series. We demonstrate that HGM yields more reliable performance evaluation than by infinite series alone and more expeditious than by simulation, for realistic values of K , and even for N_R and N_T relevant to large MIMO systems. We envision extending the proposed approaches for exact analysis and reliable evaluation to more general Rician fading and other transceiver methods. △ Less

Submitted 19 May, 2016; v1 submitted 24 July, 2015; originally announced July 2015.

Comments: Accepted for publication by the IEEE Transactions on Wireless Communications, on April 7th, 2016; this is the final revision before publication

Proof of the Wilf-Zeilberger Conjecture for Mixed Hypergeometric Terms

Authors: Shaoshi Chen, Christoph Koutschan

Abstract: In 1992, Wilf and Zeilberger conjectured that a hypergeometric term in several discrete and continuous variables is holonomic if and only if it is proper. Strictly speaking the conjecture does not hold, but it is true when reformulated properly: Payne proved a piecewise interpretation in 1997, and independently, Abramov and Petkovsek in 2002 proved a conjugate interpretation. Both results address… ▽ More In 1992, Wilf and Zeilberger conjectured that a hypergeometric term in several discrete and continuous variables is holonomic if and only if it is proper. Strictly speaking the conjecture does not hold, but it is true when reformulated properly: Payne proved a piecewise interpretation in 1997, and independently, Abramov and Petkovsek in 2002 proved a conjugate interpretation. Both results address the pure discrete case of the conjecture. In this paper we extend their work to hypergeometric terms in several discrete and continuous variables and prove the conjugate interpretation of the Wilf-Zeilberger conjecture in this mixed setting. △ Less

Submitted 8 June, 2018; v1 submitted 17 July, 2015; originally announced July 2015.

Comments: 20 pages, 1 figure

MSC Class: 33C70; 33F10

Journal ref: Journal of Symbolic Computation, vol. 93, pages 133-147, 2019

Symbolic Derivation of Mean-Field PDEs from Lattice-Based Models

Authors: Christoph Koutschan, Helene Ranetbauer, Georg Regensburger, Marie-Therese Wolfram

Abstract: Transportation processes, which play a prominent role in the life and social sciences, are typically described by discrete models on lattices. For studying their dynamics a continuous formulation of the problem via partial differential equations (PDE) is employed. In this paper we propose a symbolic computation approach to derive mean-field PDEs from a lattice-based model. We start with the micros… ▽ More Transportation processes, which play a prominent role in the life and social sciences, are typically described by discrete models on lattices. For studying their dynamics a continuous formulation of the problem via partial differential equations (PDE) is employed. In this paper we propose a symbolic computation approach to derive mean-field PDEs from a lattice-based model. We start with the microscopic equations, which state the probability to find a particle at a given lattice site. Then the PDEs are formally derived by Taylor expansions of the probability densities and by passing to an appropriate limit as the time steps and the distances between lattice sites tend to zero. We present an implementation in a computer algebra system that performs this transition for a general class of models. In order to rewrite the mean-field PDEs in a conservative formulation, we adapt and implement symbolic integration methods that can handle unspecified functions in several variables. To illustrate our approach, we consider an application in crowd motion analysis where the dynamics of bidirectional flows are studied. However, the presented approach can be applied to various transportation processes of multiple species with variable size in any dimension, for example, to confirm several proposed mean-field models for cell motility. △ Less

Submitted 21 March, 2016; v1 submitted 29 June, 2015; originally announced June 2015.

Journal ref: Proceedings of the 17th International Symposium on Symbolic and Numeric Algorithms for Scientific Computing (SYNASC), pp. 27-33, 2015. ISBN 978-1-5090-0461-4

Planar Linkages Following a Prescribed Motion

Authors: Matteo Gallet, Christoph Koutschan, Zijia Li, Georg Regensburger, Josef Schicho, Nelly Villamizar

Abstract: Designing mechanical devices, called linkages, that draw a given plane curve has been a topic that interested engineers and mathematicians for hundreds of years, and recently also computer scientists. Already in 1876, Kempe proposed a procedure for solving the problem in full generality, but his constructions tend to be extremely complicated. We provide a novel algorithm that produces much simpler… ▽ More Designing mechanical devices, called linkages, that draw a given plane curve has been a topic that interested engineers and mathematicians for hundreds of years, and recently also computer scientists. Already in 1876, Kempe proposed a procedure for solving the problem in full generality, but his constructions tend to be extremely complicated. We provide a novel algorithm that produces much simpler linkages, but works only for parametric curves. Our approach is to transform the problem into a factorization task over some noncommutative algebra. We show how to compute such a factorization, and how to use it to construct a linkage tracing a given curve. △ Less

Submitted 5 April, 2016; v1 submitted 19 February, 2015; originally announced February 2015.

Comments: 33 pages, 12 figures

MSC Class: 70B15; 68W30; 70G55; 20G20; 16Z05; 14P05; 12Y05

Journal ref: Mathematics of Computation 86 (2017), pages 473-506

Integral D-Finite Functions

Authors: Manuel Kauers, Christoph Koutschan

Abstract: We propose a differential analog of the notion of integral closure of algebraic function fields. We present an algorithm for computing the integral closure of the algebra defined by a linear differential operator. Our algorithm is a direct analog of van Hoeij's algorithm for computing integral bases of algebraic function fields. We propose a differential analog of the notion of integral closure of algebraic function fields. We present an algorithm for computing the integral closure of the algebra defined by a linear differential operator. Our algorithm is a direct analog of van Hoeij's algorithm for computing integral bases of algebraic function fields. △ Less

Submitted 30 June, 2015; v1 submitted 15 January, 2015; originally announced January 2015.

ACM Class: I.1.2

Journal ref: Proceedings of the International Symposium on Symbolic and Algebraic Computation (ISSAC 2015), pages 251-258, 2015. ACM, New York, USA, ISBN 978-1-4503-3435-8

Fundamental Laser Modes in Paraxial Optics: From Computer Algebra and Simulations to Experimental Observation

Authors: Christoph Koutschan, Erwin Suazo, Sergei K. Suslov

Abstract: We study multi-parameter solutions of the inhomogeneous paraxial wave equation in a linear and quadratic approximation which include oscillating laser beams in a parabolic waveguide, spiral light beams, and other important families of propagation-invariant laser modes in weakly varying media. A "smart" lens design and a similar effect of superfocusing of particle beams in a thin monocrystal film a… ▽ More We study multi-parameter solutions of the inhomogeneous paraxial wave equation in a linear and quadratic approximation which include oscillating laser beams in a parabolic waveguide, spiral light beams, and other important families of propagation-invariant laser modes in weakly varying media. A "smart" lens design and a similar effect of superfocusing of particle beams in a thin monocrystal film are also discussed. In the supplementary electronic material, we provide a computer algebra verification of the results presented here, and of some related mathematical tools that were stated without proofs in the literature. We also demonstrate how computer algebra can be used to derive some of the presented formulas automatically, which is highly desirable as the corresponding hand calculations are very tedious. In numerical simulations, some of the new solutions reveal quite exotic properties which deserve further investigation including an experimental observation. △ Less

Submitted 18 October, 2015; v1 submitted 2 July, 2014; originally announced July 2014.

MSC Class: 35C05; 35Q55 (Primary); 68W30; 81Q05 (Secondary)

MIMO Zero-Forcing Performance Evaluation Using the Holonomic Gradient Method

Authors: Constantin Siriteanu, Akimichi Takemura, Satoshi Kuriki, Hyundong Shin, Christoph Koutschan

Abstract: For multiple-input multiple-output (MIMO) spatial-multiplexing transmission, zero-forcing detection (ZF) is appealing because of its low complexity. Our recent MIMO ZF performance analysis for Rician--Rayleigh fading, which is relevant in heterogeneous networks, has yielded for the ZF outage probability and ergodic capacity infinite-series expressions. Because they arose from expanding the conflue… ▽ More For multiple-input multiple-output (MIMO) spatial-multiplexing transmission, zero-forcing detection (ZF) is appealing because of its low complexity. Our recent MIMO ZF performance analysis for Rician--Rayleigh fading, which is relevant in heterogeneous networks, has yielded for the ZF outage probability and ergodic capacity infinite-series expressions. Because they arose from expanding the confluent hypergeometric function $ {_1\! F_1} (\cdot, \cdot, σしぐま) $ around 0, they do not converge numerically at realistically-high Rician $ K $-factor values. Therefore, herein, we seek to take advantage of the fact that $ {_1\! F_1} (\cdot, \cdot, σしぐま) $ satisfies a differential equation, i.e., it is a \textit{holonomic} function. Holonomic functions can be computed by the \textit{holonomic gradient method} (HGM), i.e., by numerically solving the satisfied differential equation. Thus, we first reveal that the moment generating function (m.g.f.) and probability density function (p.d.f.) of the ZF signal-to-noise ratio (SNR) are holonomic. Then, from the differential equation for $ {_1\! F_1} (\cdot, \cdot, σしぐま) $, we deduce those satisfied by the SNR m.g.f. and p.d.f., and demonstrate that the HGM helps compute the p.d.f. accurately at practically-relevant values of $ K $. Finally, numerical integration of the SNR p.d.f. produced by HGM yields accurate ZF outage probability and ergodic capacity results. △ Less

Submitted 15 April, 2015; v1 submitted 15 March, 2014; originally announced March 2014.

Comments: This manuscript was accepted in December 2014

Journal ref: IEEE Transactions on Wireless Communications, vol. 14, no. 4, April 2015, pp. 2322-2335

A Generalized Apagodu-Zeilberger Algorithm

Authors: Shaoshi Chen, Manuel Kauers, Christoph Koutschan

Abstract: The Apagodu-Zeilberger algorithm can be used for computing annihilating operators for definite sums over hypergeometric terms, or for definite integrals over hyperexponential functions. In this paper, we propose a generalization of this algorithm which is applicable to arbitrary $\partial$-finite functions. In analogy to the hypergeometric case, we introduce the notion of proper $\partial$-finite… ▽ More The Apagodu-Zeilberger algorithm can be used for computing annihilating operators for definite sums over hypergeometric terms, or for definite integrals over hyperexponential functions. In this paper, we propose a generalization of this algorithm which is applicable to arbitrary $\partial$-finite functions. In analogy to the hypergeometric case, we introduce the notion of proper $\partial$-finite functions. We show that the algorithm always succeeds for these functions, and we give a tight a priori bound for the order of the output operator. △ Less

Submitted 2 August, 2014; v1 submitted 11 February, 2014; originally announced February 2014.

ACM Class: I.1.2

Journal ref: Proceedings of the International Symposium on Symbolic and Algebraic Computation (ISSAC 2014), pages 107-114, 2014. ACM, New York, USA, ISBN 978-1-4503-2501-1

Creative Telescoping for Holonomic Functions

Authors: Christoph Koutschan

Abstract: The aim of this article is twofold: on the one hand it is intended to serve as a gentle introduction to the topic of creative telescoping, from a practical point of view; for this purpose its application to several problems is exemplified. On the other hand, this chapter has the flavour of a survey article: the developments in this area during the last two decades are sketched and a selection of r… ▽ More The aim of this article is twofold: on the one hand it is intended to serve as a gentle introduction to the topic of creative telescoping, from a practical point of view; for this purpose its application to several problems is exemplified. On the other hand, this chapter has the flavour of a survey article: the developments in this area during the last two decades are sketched and a selection of references is compiled in order to highlight the impact of creative telescoping in numerous contexts. △ Less

Submitted 17 July, 2013; originally announced July 2013.

Comments: Tutorial and survey article, 24 pages

Journal ref: In Carsten Schneider and Johannes Bluemlein (eds.): Computer Algebra in Quantum Field Theory: Integration, Summation and Special Functions. Texts & Monographs in Symbolic Computation, Springer-Verlag Wien 2013

Relativistic Coulomb Integrals and Zeilberger's Holonomic Systems Approach II

Authors: Christoph Koutschan, Peter Paule, Sergei K. Suslov

Abstract: We derive the recurrence relations for relativistic Coulomb integrals directly from the integral representations with the help of computer algebra methods. In order to manage the computational complexity of this problem, we employ holonomic closure properties in a sophisticated way. We derive the recurrence relations for relativistic Coulomb integrals directly from the integral representations with the help of computer algebra methods. In order to manage the computational complexity of this problem, we employ holonomic closure properties in a sophisticated way. △ Less

Submitted 27 February, 2014; v1 submitted 6 June, 2013; originally announced June 2013.

Journal ref: In: Algebraic and Algorithmic Aspects of Differential and Integral Operators, Lecture Notes in Computer Science 8372, pp. 135-145, Springer-Verlag Berlin Heidelberg, 2014. ISBN 978-3-642-54478-1