Convergence of inexact Newton methods for generalized equations

For solving the generalized equation $f(x)+F(x) \ni 0$ , where $f$ is a smooth function and $F$ is a set-valued mapping acting between Banach spaces, we study the inexact Newton method described by $$\begin{aligned} \left( f(x_k)+ D f(x_k)(x_{k+1}-x_k) + F(x_{k+1})\right) \cap R_k(x_k, x_{k+1}) \ne \emptyset , \end{aligned}$$ where $Df$ is the derivative of $f$ and the sequence of mappings $R_k$ represents the inexactness. We show how regularity properties of the mappings $f+F$ and $R_k$ are able to guarantee that every sequence generated by the method is convergent either q-linearly, q-superlinearly, or q-quadratically, according to the particular assumptions. We also show there are circumstances in which at least one convergence sequence is sure to be generated. As a byproduct, we obtain convergence results about inexact Newton methods for solving equations, variational inequalities and nonlinear programming problems.     

Linear complementarity as absolute value equation solution

We consider the linear complementarity problem (LCP): $Mz+q\ge 0, z\ge 0, z^{\prime }(Mz+q)=0$ as an absolute value equation (AVE): $(M+I)z+q=|(M-I)z+q|$ , where $M$ is an $n\times n$ square matrix and $I$ is the identity matrix. We propose a concave minimization algorithm for solving (AVE) that consists of solving a few linear programs, typically two. The algorithm was tested on 500 consecutively generated random solvable instances of the LCP with $n=10, 50, 100, 500$ and 1,000. The algorithm solved $100\,\%$ of the test problems to an accuracy of $10^{-8}$ by solving 2 or less linear programs per LCP problem.     

On the compactness of the set of invariant Einstein metrics

Let $M = G/H$ be a connected simply connected homogeneous manifold of a compact, not necessarily connected Lie group $G$ . We will assume that the isotropy $H$ -module $\mathfrak{g/h }$ has a simple spectrum, i.e. irreducible submodules are mutually non-equivalent. There exists a convex Newton polytope $N=N(G,H)$ , which was used for the estimation of the number of isolated complex solutions of the algebraic Einstein equation for invariant metrics on $G/H$ (up to scaling). Using the moment map, we identify the space $\mathcal{M }_1$ of invariant Riemannian metrics of volume 1 on $G/H$ with the interior of this polytope $N$ . We associate with a point ${x \in \partial N}$ of the boundary a homogeneous Riemannian space (in general, only local) and we extend the Einstein equation to $\partial N$ . As an application of the Alekseevsksky–Kimel’fel’d theorem, we prove that all solutions of the Einstein equation associated with points of the boundary are locally Euclidean. We describe explicitly the set $T\subset \partial N$ of solutions at the boundary together with its natural triangulation. Investigating the compactification ${\overline{\mathcal{M }}}_{1} = N$ of $\mathcal{M }_1$ , we get an algebraic proof of the deep result by Böhm, Wang and Ziller about the compactness of the set $\mathcal{E }_1 \subset \mathcal{M }_1$ of Einstein metrics. The original proof by Böhm, Wang and Ziller was based on a different approach and did not use the simplicity of the spectrum. In Appendix, we consider the non-symmetric flag manifolds $G/H$ with the second Betti number $b_2=1$ . We calculate the normalized volumes $2,6,20,82,344$ of the corresponding Newton polytopes and discuss the number of complex solutions of the algebraic Einstein equation and the finiteness problem.     

A weak second-order split-step method for numerical simulations of stochastic differential equations

In an analogy from symmetric ordinary differential equation numerical integrators, we derive a three-stage, weak 2nd-order procedure for Monte-Carlo simulations of Itô stochastic differential equations. Our composite procedure splits each time step into three parts: an \(h/2\) -stage of trapezoidal rule, an \(h\) -stage martingale, followed by another \(h/2\) -stage of trapezoidal rule. In \(n\) time steps, an \(h/2\) -stage deterministic step follows another \(n-1\) times. Each of these adjacent pairs may be combined into a single \(h\) -stage, effectively producing a two-stage method with partial overlap between successive time steps.     

Discontinuous Galerkin method for hyperbolic equations involving \delta -singularities: negative-order norm error estimates and applications

In this paper, we develop and analyze discontinuous Galerkin (DG) methods to solve hyperbolic equations involving $\delta $ -singularities. Negative-order norm error estimates for the accuracy of DG approximations to $\delta $ -singularities are investigated. We first consider linear hyperbolic conservation laws in one space dimension with singular initial data. We prove that, by using piecewise $k$ th degree polynomials, at time $t$ , the error in the $H^{-(k+2)}$ norm over the whole domain is $(k+1/2)$ th order, and the error in the $H^{-(k+1)}(\mathbb R \backslash \mathcal R _t)$ norm is $(2k+1)$ th order, where $\mathcal R _t$ is the pollution region due to the initial singularity with the width of order $\mathcal O (h^{1/2} \log (1/h))$ and $h$ is the maximum cell length. As an application of the negative-order norm error estimates, we convolve the numerical solution with a suitable kernel which is a linear combination of B-splines, to obtain $L^2$ error estimate of $(2k+1)$ th order for the post-processed solution. Moreover, we also obtain high order superconvergence error estimates for linear hyperbolic conservation laws with singular source terms by applying Duhamel’s principle. Numerical examples including an acoustic equation and the nonlinear rendez-vous algorithms are given to demonstrate the good performance of DG methods for solving hyperbolic equations involving $\delta $ -singularities.     

Minimal subgroups and the structure of finite groups

A subgroup $H$ of a finite group $G$ is weakly-supplemented in $G$ if there exists a proper subgroup $K$ of $G$ such that $G=HK$ . In this paper we prove that a finite group $G$ is $p$ -nilpotent if every minimal subgroup of $P\bigcap G^{N}$ is weakly-supplemented in $G$ , and when $p=2$ either every cyclic subgroup of $P\bigcap G^{N}$ with order 4 is weakly-supplemented in $G$ or $P$ is quaternion-free, where $p$ is the smallest prime number dividing the order of $G$ , $P$ a sylow $p$ -subgroup of $G$ .     

Fermat-type equations of signature (13,13,p) via Hilbert cuspforms

In this paper we prove that for $p > 13649$ equations of the form $x^{13} + y^{13} = Cz^{p}$ have no non-trivial primitive solutions $(a,b,c)$ such that $13 \not \mid c$ for an infinite family of values for $C$ . Our method consists on relating a solution $(a,b,c)$ to the previous equation to a solution $(a,b,c_1)$ of another Diophantine equation with coefficients in $\mathbb Q (\sqrt{13})$ . Then we attach to $(a,b,c_1)$ a Frey curve $E_{(a,b)}$ defined over $\mathbb Q (\sqrt{13})$ that is not a $\mathbb Q $ -curve. We prove a modularity result of independent interest for certain elliptic curves over totally real abelian number fields satisfying some local conditions at $3$ . This theorem, in particular, implies modularity of $E_{(a,b)}$ . This enables us to use level lowering results and apply the modular approach via Hilbert cuspforms over $\mathbb Q (\sqrt{13})$ to prove the non-existence of $(a,b,c_1)$ and, consequently, of $(a,b,c)$ .     

The Ubiquity of the Orders of Fractional Semifields of Even Characteristic

To each non-square integer \(2^{2N+1}\ge 2^5\) there correspond semifields \(D\) of order of \(2^{2N+1}\) that contain \(\text{ GF}(4)\) . Hence there exist affine planes for each non-square order \(2^{2N+1}\ge 2^{5}\) that contain subaffine planes of order \(2^2\) . Moreover, there also exists semifields \(D_1\) and \(D_2\) , with \(|D_1|= |D_2| =|D|\) such that \(D_1\) is commutative and \(D_2\) is non-commutative but neither \(D_1\) nor \(D_2\) contains \(\text{ GF}(4)\) .     

The {\mathcal {A}}-Truncated K-Moment Problem

Let \({\mathcal {A}}\subseteq {\mathbb {N}}^n\) be a finite set, and \(K\subseteq {\mathbb {R}}^n\) be a compact semialgebraic set. An \({\mathcal {A}}\) -truncated multisequence ( \({\mathcal {A}}\) -tms) is a vector \(y=(y_{\alpha })\) indexed by elements in \({\mathcal {A}}\) . The \({\mathcal {A}}\) -truncated \(K\) -moment problem ( \({\mathcal {A}}\) -TKMP) concerns whether or not a given \({\mathcal {A}}\) -tms \(y\) admits a \(K\) -measure \(\mu \) , i.e., \(\mu \) is a nonnegative Borel measure supported in \(K\) such that \(y_\alpha = \int _K x^\alpha \mathtt {d}\mu \) for all \(\alpha \in {\mathcal {A}}\) . This paper proposes a numerical algorithm for solving \({\mathcal {A}}\) -TKMPs. It aims at finding a flat extension of \(y\) by solving a hierarchy of semidefinite relaxations \(\{(\mathtt {SDR})_k\}_{k=1}^\infty \) for a moment optimization problem, whose objective \(R\) is generated in a certain randomized way. If \(y\) admits no \(K\) -measures and \({\mathbb {R}}[x]_{{\mathcal {A}}}\) is \(K\) -full (there exists \(p \in {\mathbb {R}}[x]_{{\mathcal {A}}}\) that is positive on \(K\) ), then \((\mathtt {SDR})_k\) is infeasible for all \(k\) big enough, which gives a certificate for the nonexistence of representing measures. If \(y\) admits a \(K\) -measure, then for almost all generated \(R\) , this algorithm has the following properties: i) we can asymptotically get a flat extension of \(y\) by solving the hierarchy \(\{(\mathtt {SDR})_k\}_{k=1}^\infty \) ; ii) under a general condition that is almost sufficient and necessary, we can get a flat extension of \(y\) by solving \((\mathtt {SDR})_k\) for some \(k\) ; iii) the obtained flat extensions admit a \(r\) -atomic \(K\) -measure with \(r\le |{\mathcal {A}}|\) . The decomposition problems for completely positive matrices and sums of even powers of real linear forms, and the standard truncated \(K\) -moment problems, are special cases of \({\mathcal {A}}\) -TKMPs. They can be solved numerically by this algorithm.     

On the geometry of von Neumann algebra preduals

Let \(M\) be a von Neumann algebra and let \(M_\star \) be its (unique) predual. We study when for every \(\varphi \in M_\star \) there exists \(\psi \in M_\star \) solving the equation \(\Vert \varphi \pm \psi \Vert =\Vert \varphi \Vert =\Vert \psi \Vert \) . This is the case when \(M\) does not contain type I nor type III \(_1\) factors as direct summands and it is false at least for the unique hyperfinite type III \(_1\) factor. We also characterize this property in terms of the existence of centrally symmetric curves in the unit sphere of \(M_\star \) of length \(4\) . An approximate result valid for all diffuse von Neumann algebras allows to show that the equation has solution for every element in the ultraproduct of preduals of diffuse von Neumann algebras and, in particular, the dual von Neumann algebra of such ultraproduct is diffuse. This shows that the Daugavet property and the uniform Daugavet property are equivalent for preduals of von Neumann algebras.     

Finite groups all of whose second maximal subgroups are {\mathcal{H }}_p-groups

Let $G$ be a finite group. A subgroup $H$ of $G$ is called an $\mathcal{H }$ -subgroup of $G$ if $N_G(H)\cap H^g\le H$ for all $g\in G$ . A group $G$ is said to be an ${\mathcal{H }}_p$ -group if every cyclic subgroup of $G$ of prime order or order 4 is an $\mathcal{H }$ -subgroup of $G$ . In this paper, the structure of a finite group all of whose second maximal subgroups are ${\mathcal{H }}_p$ -subgroups has been characterized.     

Shintani cocycles and the order of vanishing of p-adic Hecke L-series at s=0

Let $\chi $ be a Hecke character of finite order of a totally real number field $F$ . By using Hill’s Shintani cocycle we provide a cohomological construction of the $p$ -adic $L$ -series $L_p(\chi , s)$ associated to $\chi $ . This is used to show that $L_p(\chi , s)$ has a trivial zero at $s=0$ of order at least equal to the number of places of $F$ above $p$ where the local component of $\chi $ is trivial.     

On S-permutably embedded subgroups of finite groups

A subgroup $A$ of a finite group $G$ is said to be $S$ -permutably embedded in $G$ if for each prime $p$ dividing the order of $A$ , every Sylow $p$ -subgroup of $A$ is a Sylow $p$ -subgroup of some $S$ -permutable subgroup of $G$ . In this paper we determine how the $S$ -permutable embedding of several families of subgroups of a finite group influences its structure.     

On some primary subgroups and the structure of finite groups

Suppose that $G$ is a finite group and $H$ is a subgroup of $G$ . $H$ is said to be an $s$ -quasinormally embedded in $G$ if for each prime $p$ dividing the order of $H$ , a Sylow $p$ -subgroup of $H$ is also a Sylow $p$ -subgroup of some $S$ -quasinormal subgroup of $G$ ; $H$ is said to be $c$ -normal in $G$ if $G$ has a normal subgroup $T$ such that $G=HT$ and $H\cap T\le H_{G}$ , where $H_{G}$ is the normal core of $H$ in $G$ . We fix in every non-cyclic Sylow subgroup $P$ of $G$ some subgroup $D$ satisfying $1<|D|<|P|$ and study the structure of $G$ under the assumption that every subgroup $H$ of $P$ with $|H|=|D|$ is either $s$ -quasinormally embedded or $c$ -normal in $G$ . Some recent results are generalized and unified.     

On an extension of the Frobenius^{\prime } theorem about p-nilpotency of a finite group

Suppose that \(G\) is a finite group and \(H\) is a subgroup of \(G\) . \(H\) is said to be \(s\) -quasinormally embedded in \(G\) if for each prime \(p\) dividing the order of \(H\) , a Sylow \(p\) -subgroup of \(H\) is also a Sylow \(p\) -subgroup of some \(s\) -quasinormal subgroup of \(G\) . We fix in every non-cyclic Sylow subgroup \(P\) of \(G\) some subgroup \(D\) satisfying \(1<|D|<|P|\) and study the \(p\) -nilpotency of \(G\) under the assumption that every subgroup \(H\) of \(P\) with \(|H|=|D|\) is \(s\) -quasinormally embedded in \(G\) . Some recent results and the Frobenius \(^{\prime }\) theorem are generalized.     

On matching cover of graphs

A k-matching cover of a graph \(G\) is a union of \(k\) matchings of \(G\) which covers \(V(G)\) . The matching cover number of \(G\) , denoted by \(mc(G)\) , is the minimum number \(k\) such that \(G\) has a \(k\) -matching cover. A matching cover of \(G\) is optimal if it consists of \(mc(G)\) matchings of \(G\) . In this paper, we present an algorithm for finding an optimal matching cover of a graph on \(n\) vertices in \(O(n^3)\) time (if use a faster maximum matching algorithm, the time complexity can be reduced to \(O(nm)\) , where \(m=|E(G)|\) ), and give an upper bound on matching cover number of graphs. In particular, for trees, a linear-time algorithm is given, and as a by-product, the matching cover number of trees is determined.     

Minimal generating and normally generating sets for the braid and mapping class groups of \mathbb{D }^{2}, \mathbb S ^{2} and \mathbb{R P^2}

We consider the (pure) braid groups $B_{n}(M)$ and $P_{n}(M)$ , where $M$ is the $2$ -sphere $\mathbb S ^{2}$ or the real projective plane $\mathbb R P^2$ . We determine the minimal cardinality of (normal) generating sets $X$ of these groups, first when there is no restriction on $X$ , and secondly when $X$ consists of elements of finite order. This improves on results of Berrick and Matthey in the case of $\mathbb S ^{2}$ , and extends them in the case of $\mathbb R P^2$ . We begin by recalling the situation for the Artin braid groups ( $M=\mathbb{D }^{2}$ ). As applications of our results, we answer the corresponding questions for the associated mapping class groups, and we show that for $M=\mathbb S ^{2}$ or $\mathbb R P^2$ , the induced action of $B_n(M)$ on $H_3(\widetilde{F_n(M)};\mathbb{Z })$ is trivial, $F_{n}(M)$ being the $n^\mathrm{th}$ configuration space of $M$ .     

Uniqueness of generalized p-area minimizers and integrability of a horizontal normal in the Heisenberg group

We study the uniqueness of generalized \(p\) -minimal surfaces in the Heisenberg group. The generalized \(p\) -area of a graph defined by \(u\) reads \(\int |\nabla u+\vec {F}|+Hu.\) If \(u\) and \(v\) are two minimizers for the generalized \(p\) -area satisfying the same Dirichlet boundary condition, then we can only get \(N_{\vec {F}}(u) = N_{\vec {F}}(v)\) (on the nonsingular set) where \(N_{\vec {F}}(w) := \frac{\nabla w+\vec {F}}{|\nabla w+\vec {F}|}.\) To conclude \(u = v\) (or \(\nabla u = \nabla v)\) , it is not straightforward as in the Riemannian case, but requires some special argument in general. In this paper, we prove that \(N_{\vec {F}}(u) = N_{ \vec {F}}(v)\) implies \(\nabla u = \nabla v\) in dimension \(\ge \) 3 under some rank condition on derivatives of \(\vec {F}\) or the nonintegrability condition of contact form associated to \(u\) or \(v\) . Note that in dimension 2 ( \(n=1),\) the above statement is no longer true. Inspired by an equation for the horizontal normal \(N_{\vec {F}}(u),\) we study the integrability for a unit vector to be the horizontal normal of a graph. We find a Codazzi-like equation together with this equation to form an integrability condition.     

On generalizing transitivity, persistence, and sensitivity

A subgroup property $\alpha $ is transitive in a group $G$ if $U \alpha V$ and $V \alpha G$ imply that $U \alpha G$ whenever $U \le V \le G$ , and $\alpha $ is persistent in $G$ if $U \alpha G$ implies that $U \alpha V$ whenever $U \le V \le G$ . Even though a subgroup property $\alpha $ may be neither transitive nor persistent, a given subgroup $U$ may have the property that each $\alpha $ -subgroup of $U$ is an $\alpha $ -subgroup of $G$ , or that each $\alpha $ -subgroup of $G$ in $U$ is an $\alpha $ -subgroup of $U$ . We call these subgroup properties $\alpha $ -transitivity and $\alpha $ -persistence, respectively. We introduce and develop the notions of $\alpha $ -transitivity and $\alpha $ -persistence, and we establish how the former property is related to $\alpha $ -sensitivity. In order to demonstrate how these concepts can be used, we apply the results to the cases in which $\alpha $ is replaced with “normal” and the “cover-avoidance property.” We also suggest ways in which the theory can be developed further.     

On a Regularized Family of Models for Homogeneous Incompressible Two-Phase Flows

We consider a general family of regularized models for incompressible two-phase flows based on the Allen–Cahn formulation in \(n\) -dimensional compact Riemannian manifolds for \(n=2,3\) . The system we consider consists of a regularized family of Navier–Stokes equations (including the Navier–Stokes- \(\alpha \) -like model, the Leray- \(\alpha \) model, the modified Leray- \(\alpha \) model, the simplified Bardina model, the Navier–Stokes–Voight model, and the Navier–Stokes model) for the fluid velocity \(u\) suitably coupled with a convective Allen–Cahn equation for the order (phase) parameter \(\phi \) . We give a unified analysis of the entire three-parameter family of two-phase models using only abstract mapping properties of the principal dissipation and smoothing operators and then use assumptions about the specific form of the parameterizations, leading to specific models, only when necessary to obtain the sharpest results. We establish existence, stability, and regularity results and some results for singular perturbations, which as special cases include the inviscid limit of viscous models and the \(\alpha \rightarrow 0\) limit in \(\alpha \) models. Then we show the existence of a global attractor and exponential attractor for our general model and establish precise conditions under which each trajectory \(\left( u,\phi \right) \) converges to a single equilibrium by means of a Lojasiewicz–Simon inequality. We also derive new results on the existence of global and exponential attractors for the regularized family of Navier–Stokes equations and magnetohydrodynamics models that improve and complement the results of Holst et al. (J Nonlinear Sci 20(5):523–567, 2010). Finally, our analysis is applied to certain regularized Ericksen–Leslie models for the hydrodynamics of liquid crystals in \(n\) -dimensional compact Riemannian manifolds.