Related papers: Condition Estimates for Pseudo-Arclength Continuat…
Singular equations with rank-deficient Jacobians arise frequently in algebraic computing applications. As shown in case studies in this paper, direct and intuitive modeling of algebraic problems often results in nonisolated singular…
We study a semismooth Newton-type method for the nearest doubly stochastic matrix problem where both differentiability and nonsingularity of the Jacobian can fail. The optimality conditions for this problem are formulated as a system of…
We describe inexact proximal Newton-like methods for solving degenerate regularized optimization problems and for the broader problem of finding a zero of a generalized equation that is the sum of a continuous map and a maximal monotone…
We introduce a new non-degeneracy condition at infinity for a real or a mixed polynomial mapping $F$ which allows us to approximate its bifurcation locus in terms of certain Newton polyhedra. We derive a sufficiency result for the Jacobian…
Many problems give rise to polynomial systems. These systems often have several parameters and we are interested to study how the solutions vary when we change the values for the parameters. Using predictor-corrector methods we track the…
We study the bifurcation of solutions of semilinear elliptic boundary value problems of the form \begin{align*} \begin{aligned} -\Delta u &= f_\lambda(|x|,u,|\nabla u|) &&\text{in }\Omega, u &= 0 &&\text{on }\partial\Omega, \end{aligned}…
By using the Hadamard matrix product concept, this paper introduces two generalized matrix formulation forms of numerical analogue of nonlinear differential operators. The SJT matrix-vector product approach is found to be a simple,…
This is a PhD thesis about generated Jacobian equations; our purpose is twofold. First, we provide an introduction to these equations, whilst, at the same time, collating some results scattered throughout the literature. The other goal is…
It is well-known since the work of Pardoux and Peng [12] that Backward Stochastic Differential Equations provide probabilistic formulae for the solution of (systems of) second order elliptic and parabolic equations, thus providing an…
We introduce a quadratically convergent semismooth Newton method for nonlinear semidefinite programming that eliminates the need for the generalized Jacobian regularity, a common yet stringent requirement in existing approaches. Our…
In this paper, we investigate condition numbers of eigenvalue problems of matrix polynomials with nonsingular leading coefficients, generalizing classical results of matrix perturbation theory. We provide a relation between the condition…
We proposed in this paper a new method, which we named the W4 method, to solve nonlinear equation systems. It may be regarded as an extension of the Newton-Raphson~(NR) method to be used when the method fails. Indeed our method can be…
We consider positive solutions, possibly unbounded, to the semilinear equation $-\Delta u=f(u)$ on continuous epigraphs bounded from below. Under the homogeneous Dirichlet boundary condition, we prove new monotonicity results for $u$, when…
Our main result asserts that a certain natural non-linear operator on Jacobi matrices built by a hyperbolic polynomial with real Julia set is a contraction in operator norm if the polynomial is sufficiently hyperbolic. This allows us to get…
The textbook Newton's iteration is practically inapplicable on solutions of nonlinear systems with singular Jacobians. By a simple modification, a novel extension of Newton's iteration regains its local quadratic convergence toward…
We describe a three precision variant of Newton's method for nonlinear equations. We evaluate the nonlinear residual in double precision, store the Jacobian matrix in single precision, and solve the equation for the Newton step with…
This article considers the semilinear boundary value problem given by the Poisson equation, -\Delta u=f(u) in a bounded domain \Omega\subset \R^{n} with smooth boundary. For the zero boundary value case, we approximate a solution using the…
We apply the methods of classical approximation theory (extreme properties of polynomials) to study the essential support $\Sigma_{ac}$ of the absolutely continuous spectrum of Jacobi matrices. First, we prove an upper bound on the measure…
Solving complex optimization problems in engineering and the physical sciences requires repetitive computation of multi-dimensional function derivatives. Commonly, this requires computationally-demanding numerical differentiation such as…
We construct and justify leading order weakly nonlinear geometric optics expansions for nonlinear hyperbolic initial value problems, including the compressible Euler equations. The technique of simultaneous Picard iteration is employed to…