Related papers: Functional norms, condition numbers and numerical …
In this article we give an implementation of the standard algorithm to segment a real algebraic plane curve defined implicitly. Our implementation is efficient and simpler than previous. We use global information to count the number of…
Gr{\"o}bner bases is one the most powerful tools in algorithmic non-linear algebra. Their computation is an intrinsically hard problem with a complexity at least single exponential in the number of variables. However, in most of the cases,…
This article presents a numerical illustration of a recently proposed strongly polynomial-time algorithm for the general linear programming (LP) problem. Each iteration of the proposed algorithm consists of two Gauss-Jordan pivoting…
We give a lower bound for the numerical index of two-dimensional real spaces with absolute and symmetric norm. This allows us to compute the numerical index of the two-dimensional real $L_p$-space for $3/2\leq p\leq 3$.
Enumerative Geometry is concerned with the number of solutions to a structured system of polynomial equations, when the structure comes from geometry. Enumerative real algebraic geometry studies real solutions to such systems, particularly…
Statistics and Optimization are foundational to modern Machine Learning. Here, we propose an alternative foundation based on Abstract Algebra, with mathematics that facilitates the analysis of learning. In this approach, the goal of the…
We compare algorithms for global optimization of polynomial functions in many variables. It is demonstrated that existing algebraic methods (Gr\"obner bases, resultants, homotopy methods) are dramatically outperformed by a relaxation…
We present a new family of quantum Weyl algebras where the polynomial part is the quantum analog of functions on homogeneous spaces corresponding to symmetric matrices, skew symmetric matrices, and the entire space of matrices of a given…
In the last decade major steps towards an algorithmic treatment of orthogonal polynomials and special functions (OP & SF) have been made, notably Zeilberger's brilliant extension of Gosper's algorithm on algorithmic definite hypergeometric…
We study the $l^p$ norms of a class of weighted mean matrices whose diagonal terms are given by $n^{\alpha}/\sum^{n}_{i=1}i^{\alpha}$ with $\alpha > -1$. The $l^p$ norms of such matrices are known for $p \geq 2, (\alpha+1)p >1$ and $1<p…
We consider $\ell_p$-direct sums ($1\leq p<\infty$) and $c_0$-direct sums of countably many normed spaces and find the duals of these spaces. We characterize the support functionals of arbitrary elements in these spaces to characterize…
Nonsmooth nonconvex optimization problems involving the $\ell^p$ quasi-norm, $p \in (0, 1]$, of a linear map are considered. A monotonically convergent scheme for a regularized version of the original problem is developed and necessary…
Partitioning and grouping of similar objects plays a fundamental role in image segmentation and in clustering problems. In such problems a typical goal is to group together similar objects, or pixels in the case of image processing. At the…
The simplex method in Linear Programming motivates several problems of asymptotic convex geometry. We discuss some conjectures and known results in two related directions -- computing the size of projections of high dimensional polytopes…
We study convolutional neural networks with monomial activation functions. Specifically, we prove that their parameterization map is regular and is an isomorphism almost everywhere, up to rescaling the filters. By leveraging on tools from…
Sparse matrix computations are ubiquitous in scientific computing. With the recent interest in scientific machine learning, it is natural to ask how sparse matrix computations can leverage neural networks (NN). Unfortunately, multi-layer…
A symmetric function of $N$ variables can be given in terms of symmetric polynomials of these variables. We determine those symmetric polynomials in which the dual differential operators take the neatest form when expressed in terms of our…
We present incomplete gamma kernels, a generalization of Locally Optimal Projection (LOP) operators. In particular, we reveal the relation of the classical localized $ L_1 $ estimator, used in the LOP operator for point cloud denoising, to…
Let $E$ be a Jordan rectifiable curve in the complex plane and let $G$ be the bounded component of $\mathbb{C}\backslash E$. Now let $n\in \mathbb{N}$, and let $m_{n,E}$ denote the extremal constants defined by \begin{equation*}m_{n,E}=\inf…
The analysis of complex nonlinear systems is often carried out using simpler piecewise linear representations of them. A principled and practical technique is proposed to linearize and evaluate arbitrary continuous nonlinear functions using…