14 papers · ranked by Valyu relevance
Rolando D. Somma, Guang Hao Low, Dominic W. Berry, Ryan Babbush
We describe an efficient quantum algorithm for solving the linear matrix equation AX+XB = C, where A, B and C are given complex matrices and X is unknown. This is known as the Sylvester equation, a fundamental equation with applications in control theory and physics. Our approach constructs the solution matrix X/x in a…
Jason Rader, Terry Lyons, Patrick Kidger
We introduce Lineax, a library bringing linear solves and linear least-squares to the JAX+Equinox scientific computing ecosystem. Lineax uses general linear operators, and unifies linear solves and least-squares into a single, autodifferentiable API. Solvers and operators are user-extensible, without requiring the user…
Jing He, Qi-wei Kong, Ho-Chung Lui, Haitao Liu + 1 more
The definition of factor space and a unified optimization based classification model were developed for linear programming and supervised learning. Intelligent behaviour appeared in a decision process can be treated as a moving point y, the dynamic state observed and controlled by the agent, moving in a factor space…
Yair Censor, Yehuda Zur
Linear superiorization (abbreviated: LinSup) considers linear programming (LP) problems wherein the constraints as well as the objective function are linear. It allows to steer the iterates of a feasibilityseeking iterative process toward feasible points that have lower (not necessarily minimal) values of the objective…
Ho Yee Cheung, Tsz Chiu Kwok, Lap Lau
We consider the problem of computing the rank of an m × n matrix A over a field. We present a randomized algorithm to find a set of r = rank(A) linearly independent columns in O˜(|A| + r ω ) field operations, where |A| denotes the number of nonzero entries in A and ω < 2.38 is the matrix multiplication exponent.…
Nikica Hlupić, Ivo Beroš
An algorithm and associated strategy for solving polynomial systems within the optimization framework is presented. The algorithm and strategy are named, respectively, the penetrating gradient algorithm and the deepest descent strategy. The most prominent feature of penetrating gradient algorithm, after which it was…
Yair Censor
Linear superiorization considers linear programming problems but instead of attempting to solve them with linear optimization methods it employs perturbation resilient feasibility-seeking algorithms and steers them toward reduced (not necessarily minimal) target function values. The two questions that we set out to…
A. Ya. Rodionov
The purpose of this short article is to bring attention to unifying approach to software and hardware design suggested and developed by Grigory Litvinov, Viktor Maslov and coworkers[1]. The unifying approach is based on observation that many algorithms do not depend on particular models of a numerical domain and even…
Arya Chakraborty
— While time complexity and space complexity of an algorithm helps to determine its efficiency when time or space needs to be optimized respectively, they fail to determine the more efficient algorithm when time and space both need to be optimized simultaneously. This resulted in the development of the A1-Score Factor…
Minati De, Subhas C. Nandy, Sasanka Roy
Prune-and-search is an important paradigm for solving many important geometric problems. We show that the general prune-andsearch technique can be implemented where the objects are given in read-only memory. As examples we consider convex-hull in 2D, and linear programming in 2D and 3D. For the convex-hull problem…
Jianer Chen, Qin Huang, Iyad Kanj, Ge Xia
We study fundamental point-line covering problems in computational geometry, in which the input is a set S of points in the plane. The first is the Rich Lines problem, which asks for the set of all lines that each covers at least λ points from S, for a given integer parameter λ ≥ 2; this problem subsumes the…
Lily Li, Aleksandar Nikolov
Many problems in computer science and applied mathematics require rounding a vector w of fractional values lying in the interval [0, 1] to a binary vector x so that, for a given matrix A, Ax is as close to Aw as possible. For example, this problem arises in LP rounding algorithms used to approximate NP-hard…
Andrea Brilli, Morteza Kimiaei, Giampaolo Liuzzi, Stefano Lucidi
This paper is devoted to the analysis of worst case complexity bounds for linesearchtype derivative-free algorithms for the minimization of general non-convex smooth functions. We prove that two linesearch-type algorithms enjoy the same complexity properties which have been proved for pattern and direct search…
Myung Cho, Weiyu Xu
The null space condition of sensing matrices plays an important role in guaranteeing the success of compressed sensing. In this paper, we propose new efficient algorithms to verify the null space condition in compressed sensing (CS). Given an (n − m) × n (m > 0) CS matrix A and a positive k, we are interested in…