skip to main content


Title: On the stationarity for nonlinear optimization problems with polyhedral constraints
Abstract

For polyhedral constrained optimization problems and a feasible point$$\textbf{x}$$x, it is shown that the projection of the negative gradient on the tangent cone, denoted$$\nabla _\varOmega f(\textbf{x})$$Ωf(x), has an orthogonal decomposition of the form$$\varvec{\beta }(\textbf{x}) + \varvec{\varphi }(\textbf{x})$$β(x)+φ(x). At a stationary point,$$\nabla _\varOmega f(\textbf{x}) = \textbf{0}$$Ωf(x)=0so$$\Vert \nabla _\varOmega f(\textbf{x})\Vert $$Ωf(x)reflects the distance to a stationary point. Away from a stationary point,$$\Vert \varvec{\beta }(\textbf{x})\Vert $$β(x)and$$\Vert \varvec{\varphi }(\textbf{x})\Vert $$φ(x)measure different aspects of optimality since$$\varvec{\beta }(\textbf{x})$$β(x)only vanishes when the KKT multipliers at$$\textbf{x}$$xhave the correct sign, while$$\varvec{\varphi }(\textbf{x})$$φ(x)only vanishes when$$\textbf{x}$$xis a stationary point in the active manifold. As an application of the theory, an active set algorithm is developed for convex quadratic programs which adapts the flow of the algorithm based on a comparison between$$\Vert \varvec{\beta }(\textbf{x})\Vert $$β(x)and$$\Vert \varvec{\varphi }(\textbf{x})\Vert $$φ(x).

 
more » « less
Award ID(s):
2031213
NSF-PAR ID:
10415087
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
Mathematical Programming
Volume:
205
Issue:
1-2
ISSN:
0025-5610
Format(s):
Medium: X Size: p. 107-134
Size(s):
["p. 107-134"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Using extensive numerical simulation of the Navier–Stokes equations, we study the transition from the Darcy’s law for slow flow of fluids through a disordered porous medium to the nonlinear flow regime in which the effect of inertia cannot be neglected. The porous medium is represented by two-dimensional slices of a three-dimensional image of a sandstone. We study the problem over wide ranges of porosity and the Reynolds number, as well as two types of boundary conditions, and compute essential features of fluid flow, namely, the strength of the vorticity, the effective permeability of the pore space, the frictional drag, and the relationship between the macroscopic pressure gradient$${\varvec{\nabla }}P$$Pand the fluid velocityv. The results indicate that when the Reynolds number Re is low enough that the Darcy’s law holds, the magnitude$$\omega _z$$ωzof the vorticity is nearly zero. As Re increases, however, so also does$$\omega _z$$ωz, and its rise from nearly zero begins at the same Re at which the Darcy’s law breaks down. We also show that a nonlinear relation between the macroscopic pressure gradient and the fluid velocityv, given by,$$-{\varvec{\nabla }}P=(\mu /K_e)\textbf{v}+\beta _n\rho |\textbf{v}|^2\textbf{v}$$-P=(μ/Ke)v+βnρ|v|2v, provides accurate representation of the numerical data, where$$\mu$$μand$$\rho$$ρare the fluid’s viscosity and density,$$K_e$$Keis the effective Darcy permeability in the linear regime, and$$\beta _n$$βnis a generalized nonlinear resistance. Theoretical justification for the relation is presented, and its predictions are also compared with those of the Forchheimer’s equation.

     
    more » « less
  2. Abstract

    It has been recently established in David and Mayboroda (Approximation of green functions and domains with uniformly rectifiable boundaries of all dimensions.arXiv:2010.09793) that on uniformly rectifiable sets the Green function is almost affine in the weak sense, and moreover, in some scenarios such Green function estimates are equivalent to the uniform rectifiability of a set. The present paper tackles a strong analogue of these results, starting with the “flagship degenerate operators on sets with lower dimensional boundaries. We consider the elliptic operators$$L_{\beta ,\gamma } =- {\text {div}}D^{d+1+\gamma -n} \nabla $$Lβ,γ=-divDd+1+γ-nassociated to a domain$$\Omega \subset {\mathbb {R}}^n$$ΩRnwith a uniformly rectifiable boundary$$\Gamma $$Γof dimension$$d < n-1$$d<n-1, the now usual distance to the boundary$$D = D_\beta $$D=Dβgiven by$$D_\beta (X)^{-\beta } = \int _{\Gamma } |X-y|^{-d-\beta } d\sigma (y)$$Dβ(X)-β=Γ|X-y|-d-βdσ(y)for$$X \in \Omega $$XΩ, where$$\beta >0$$β>0and$$\gamma \in (-1,1)$$γ(-1,1). In this paper we show that the Green functionGfor$$L_{\beta ,\gamma }$$Lβ,γ, with pole at infinity, is well approximated by multiples of$$D^{1-\gamma }$$D1-γ, in the sense that the function$$\big | D\nabla \big (\ln \big ( \frac{G}{D^{1-\gamma }} \big )\big )\big |^2$$|D(ln(GD1-γ))|2satisfies a Carleson measure estimate on$$\Omega $$Ω. We underline that the strong and the weak results are different in nature and, of course, at the level of the proofs: the latter extensively used compactness arguments, while the present paper relies on some intricate integration by parts and the properties of the “magical distance function from David et al. (Duke Math J, to appear).

     
    more » « less
  3. Abstract

    In this article, we study the moduli of irregular surfaces of general type with at worst canonical singularities satisfying$$K^2 = 4p_g-8$$K2=4pg-8, for any even integer$$p_g\ge 4$$pg4. These surfaces also have unbounded irregularityq. We carry out our study by investigating the deformations of the canonical morphism$$\varphi :X\rightarrow {\mathbb {P}}^N$$φ:XPN, where$$\varphi $$φis a quadruple Galois cover of a smooth surface of minimal degree. These canonical covers are classified in Gallego and Purnaprajna (Trans Am Math Soc 360(10):5489-5507, 2008) into four distinct families, one of which is the easy case of a product of curves. The main objective of this article is to study the deformations of the other three, non trivial, unbounded families. We show that any deformation of$$\varphi $$φfactors through a double cover of a ruled surface and, hence, is never birational. More interestingly, we prove that, with two exceptions, a general deformation of$$\varphi $$φis two-to-one onto its image, whose normalization is a ruled surface of appropriate irregularity. We also show that, with the exception of one family, the deformations ofXare unobstructed even though$$H^2(T_X)$$H2(TX)does not vanish. Consequently,Xbelongs to a unique irreducible component of the Gieseker moduli space. These irreducible components are uniruled. As a result of all this, we show the existence of infinitely many moduli spaces, satisfying the strict Beauville inequality$$p_g > 2q-4$$pg>2q-4, with an irreducible component that has a proper quadruple sublocus where the degree of the canonical morphism jumps up. These components are above the Castelnuovo line, but nonetheless parametrize surfaces with non birational canonical morphisms. The existence of jumping subloci is a contrast with the moduli of surfaces with$$K^2 = 2p_g- 4$$K2=2pg-4, studied by Horikawa. Irreducible moduli components with a jumping sublocus also present a similarity and a difference to the moduli of curves of genus$$g\ge 3$$g3, for, like in the case of curves, the degree of the canonical morphism goes down outside a closed sublocus but, unlike in the case of curves, it is never birational. Finally, our study shows that there are infinitely many moduli spaces with an irreducible component whose general elements have non birational canonical morphism and another irreducible component whose general elements have birational canonical map.

     
    more » « less
  4. Abstract

    The azimuthal ($$\Delta \varphi $$Δφ) correlation distributions between heavy-flavor decay electrons and associated charged particles are measured in pp and p–Pb collisions at$$\sqrt{s_{\mathrm{{NN}}}} = 5.02$$sNN=5.02TeV. Results are reported for electrons with transverse momentum$$44<pT<16$$\textrm{GeV}/c$$GeV/c and pseudorapidity$$|\eta |<0.6$$|η|<0.6. The associated charged particles are selected with transverse momentum$$11<pT<7$$\textrm{GeV}/c$$GeV/c, and relative pseudorapidity separation with the leading electron$$|\Delta \eta | < 1$$|Δη|<1. The correlation measurements are performed to study and characterize the fragmentation and hadronization of heavy quarks. The correlation structures are fitted with a constant and two von Mises functions to obtain the baseline and the near- and away-side peaks, respectively. The results from p–Pb collisions are compared with those from pp collisions to study the effects of cold nuclear matter. In the measured trigger electron and associated particle kinematic regions, the two collision systems give consistent results. The$$\Delta \varphi $$Δφdistribution and the peak observables in pp and p–Pb collisions are compared with calculations from various Monte Carlo event generators.

     
    more » « less
  5. Abstract

    The double differential cross sections of the Drell–Yan lepton pair ($$\ell ^+\ell ^-$$+-, dielectron or dimuon) production are measured as functions of the invariant mass$$m_{\ell \ell }$$m, transverse momentum$$p_{\textrm{T}} (\ell \ell )$$pT(), and$$\varphi ^{*}_{\eta }$$φη. The$$\varphi ^{*}_{\eta }$$φηobservable, derived from angular measurements of the leptons and highly correlated with$$p_{\textrm{T}} (\ell \ell )$$pT(), is used to probe the low-$$p_{\textrm{T}} (\ell \ell )$$pT()region in a complementary way. Dilepton masses up to 1$$\,\text {Te\hspace{-.08em}V}$$TeVare investigated. Additionally, a measurement is performed requiring at least one jet in the final state. To benefit from partial cancellation of the systematic uncertainty, the ratios of the differential cross sections for various$$m_{\ell \ell }$$mranges to those in the Z mass peak interval are presented. The collected data correspond to an integrated luminosity of 36.3$$\,\text {fb}^{-1}$$fb-1of proton–proton collisions recorded with the CMS detector at the LHC at a centre-of-mass energy of 13$$\,\text {Te\hspace{-.08em}V}$$TeV. Measurements are compared with predictions based on perturbative quantum chromodynamics, including soft-gluon resummation.

     
    more » « less