546 research outputs found
On some approximation problems concerning sparse polynomials over finite fields
AbstractWe obtain new lower bounds on the number of non-zeros of sparse polynomials and give a fully polynomial time (ε, δ) approximation algorithm for the number of non-zeros of multivariate sparse polynomials over a finite field of q elements and degree less than q − 1. This partially answers an open problem of D. Grigoriev and M. Karpinski. Also, probabilistic and deterministic algorithms for testing identity to zero of a sparse polynomial given by a “black-box” are given. Finally, we propose an algorithm to estimate the size of the image of a univariate sparse polynomial
A local decision test for sparse polynomials
An ℓ-sparse (multivariate) polynomial is a polynomial containing at most ℓ-monomials in its explicit description. We assume that a polynomial is implicitly represented as a black-box: on an input query from the domain, the black-box replies with the evaluation of the polynomial at that input. We provide an efficient, randomized algorithm, that can decide whether a polynomial [MathML] given as a black-box is ℓ-sparse or not, provided that q is large compared to the polynomial's total degree. The algorithm makes only queries, which is independent of the domain size. The running time of our algorithm (in the bit-complexity model) is , where d is an upper bound on the degree of each variable. Existing interpolation algorithms for polynomials in the same model run in time . We provide a similar test for polynomials with integer coefficients
Development of symbolic algorithms for certain algebraic processes
This study investigates the problem of computing the exact greatest common divisor of two polynomials relative to an orthogonal basis, defined over the rational number field. The main objective of the study is to design and implement an effective and efficient symbolic algorithm for the general class of dense polynomials, given the rational number defining terms of their basis. From a general algorithm using the comrade matrix approach, the nonmodular and modular techniques are prescribed. If the coefficients of the generalized polynomials are multiprecision integers, multiprecision arithmetic will be required in the construction of the comrade matrix and the corresponding systems coefficient matrix. In addition, the application of the nonmodular elimination technique on this coefficient matrix extensively applies multiprecision rational number operations. The modular technique is employed to minimize the complexity involved in such computations. A divisor test algorithm that enables the detection of an unlucky reduction is a crucial device for an effective implementation of the modular technique. With the bound of the true solution not known a priori, the test is devised and carefully incorporated into the modular algorithm. The results illustrate that the modular algorithm illustrate its best performance for the class of relatively prime polynomials. The empirical computing time results show that the modular algorithm is markedly superior to the nonmodular algorithms in the case of sufficiently dense Legendre basis polynomials with a small GCD solution. In the case of dense Legendre basis polynomials with a big GCD solution, the modular algorithm is significantly superior to the nonmodular algorithms in higher degree polynomials. For more definitive conclusions, the computing time functions of the algorithms that are presented in this report have been worked out. Further investigations have also been suggested
Efficient spatial modelling using the SPDE approach with bivariate splines
Gaussian fields (GFs) are frequently used in spatial statistics for their
versatility. The associated computational cost can be a bottleneck, especially
in realistic applications. It has been shown that computational efficiency can
be gained by doing the computations using Gaussian Markov random fields (GMRFs)
as the GFs can be seen as weak solutions to corresponding stochastic partial
differential equations (SPDEs) using piecewise linear finite elements. We
introduce a new class of representations of GFs with bivariate splines instead
of finite elements. This allows an easier implementation of piecewise
polynomial representations of various degrees. It leads to GMRFs that can be
inferred efficiently and can be easily extended to non-stationary fields. The
solutions approximated with higher order bivariate splines converge faster,
hence the computational cost can be alleviated. Numerical simulations using
both real and simulated data also demonstrate that our framework increases the
flexibility and efficiency.Comment: 26 pages, 7 figures and 3 table
Robust Network Topology Inference and Processing of Graph Signals
The abundance of large and heterogeneous systems is rendering contemporary
data more pervasive, intricate, and with a non-regular structure. With
classical techniques facing troubles to deal with the irregular (non-Euclidean)
domain where the signals are defined, a popular approach at the heart of graph
signal processing (GSP) is to: (i) represent the underlying support via a graph
and (ii) exploit the topology of this graph to process the signals at hand. In
addition to the irregular structure of the signals, another critical limitation
is that the observed data is prone to the presence of perturbations, which, in
the context of GSP, may affect not only the observed signals but also the
topology of the supporting graph. Ignoring the presence of perturbations, along
with the couplings between the errors in the signal and the errors in their
support, can drastically hinder estimation performance. While many GSP works
have looked at the presence of perturbations in the signals, much fewer have
looked at the presence of perturbations in the graph, and almost none at their
joint effect. While this is not surprising (GSP is a relatively new field), we
expect this to change in the upcoming years. Motivated by the previous
discussion, the goal of this thesis is to advance toward a robust GSP paradigm
where the algorithms are carefully designed to incorporate the influence of
perturbations in the graph signals, the graph support, and both. To do so, we
consider different types of perturbations, evaluate their disruptive impact on
fundamental GSP tasks, and design robust algorithms to address them.Comment: Dissertatio
- …