The Dynamical Systems Method for solving nonlinear equations with monotone operators

A review of the authors's results is given. Several methods are discussed for solving nonlinear equations $F(u)=f$, where $F$ is a monotone operator in a Hilbert space, and noisy data are given in place of the exact data. A discrepancy principle for solving the equation is formulated and justified. Various versions of the Dynamical Systems Method (DSM) for solving the equation are formulated. These methods consist of a regularized Newton-type method, a gradient-type method, and a simple iteration method. A priori and a posteriori choices of stopping rules for these methods are proposed and justified. Convergence of the solutions, obtained by these methods, to the minimal norm solution to the equation $F(u)=f$ is proved. Iterative schemes with a posteriori choices of stopping rule corresponding to the proposed DSM are formulated. Convergence of these iterative schemes to a solution to equation $F(u)=f$ is justified. New nonlinear differential inequalities are derived and applied to a study of large-time behavior of solutions to evolution equations. Discrete versions of these inequalities are established.Comment: 50p

A non-linear structure preserving matrix method for the low rank approximation of the Sylvester resultant matrix

A non-linear structure preserving matrix method for the computation of a structured low rank approximation S((f) over bar , (g) over bar) of the Sylvester resultant matrix S(f , g) of two inexact polynomials f = f(y) and g = g(y) is considered in this paper. It is shown that considerably improved results are obtained when f (y) and g(y) are processed prior to the computation of S((f) over bar , (g) over bar), and that these preprocessing operations introduce two parameters. These parameters can either be held constant during the computation of S((f) over bar , (g) over bar), which leads to a linear structure preserving matrix method, or they can be incremented during the computation of S((f) over bar, (g) over bar), which leads to a non-linear structure preserving matrix method. It is shown that the non-linear method yields a better structured low rank approximation of S((f) over bar , (g) over bar) and that the assignment of f (y) and g(y) is important because S((f) over bar , (g) over bar) may be a good structured low rank approximation of S(f, g), but S((f) over bar , (g) over bar) may be a poor structured low rank approximation of S (g f) because its numerical rank is not defined. Examples that illustrate the differences between the linear and non-linear structure preserving matrix methods, and the importance of the assignment off (y) and g(y), are shown. (C) 2010 Elsevier B.V. All rights reserved

Upper bounds on quantum query complexity inspired by the Elitzur-Vaidman bomb tester

Inspired by the Elitzur-Vaidman bomb testing problem [arXiv:hep-th/9305002], we introduce a new query complexity model, which we call bomb query complexity $B(f)$. We investigate its relationship with the usual quantum query complexity $Q(f)$, and show that $B(f)=\Theta(Q(f)^2)$. This result gives a new method to upper bound the quantum query complexity: we give a method of finding bomb query algorithms from classical algorithms, which then provide nonconstructive upper bounds on $Q(f)=\Theta(\sqrt{B(f)})$. We subsequently were able to give explicit quantum algorithms matching our upper bound method. We apply this method on the single-source shortest paths problem on unweighted graphs, obtaining an algorithm with $O(n^{1.5})$ quantum query complexity, improving the best known algorithm of $O(n^{1.5}\sqrt{\log n})$ [arXiv:quant-ph/0606127]. Applying this method to the maximum bipartite matching problem gives an $O(n^{1.75})$ algorithm, improving the best known trivial $O(n^2)$ upper bound.Comment: 32 pages. Minor revisions and corrections. Regev and Schiff's proof that P(OR) = \Omega(N) remove

Solving multivariate functional equations

This paper presents a new method to solve functional equations of multivariate generating functions, such as $$F(r,s)=e(r,s)+xf(r,s)F(1,1)+xg(r,s)F(qr,1)+xh(r,s)F(qr,qs),$$ giving a formula for $F(r,s)$ in terms of a sum over finite sequences. We use this method to show how one would calculate the coefficients of the generating function for parallelogram polyominoes, which is impractical using other methods. We also apply this method to answer a question from fully commutative affine permutations.Comment: 11 pages, 1 figure. v3: Main theorems and writing style revised for greater clarity. Updated to final version, to appear in Discrete Mathematic