The Lambert W-functions and some of their integrals: a case study of high-precision computation

The real-valued Lambert W-functions considered here are w ₀(y) and w _− 1(y), solutions of we ^w  = y, − 1/e < y < 0, with values respectively in ( − 1,0) and ( − ∞ , − 1). A study is made of the numerical evaluation to high precision of these functions and of the integrals ò₁^￥ [-w₀(-xe^-x)]^a x^-bx\int_1^\infty [-w_0(-xe^{-x})]^\alpha x^{-\beta}\d x, α > 0, β ∈ ℝ, and ò₀¹ [-w_-1(-x e^-x)]^a x^-bx\int_0^1 [-w_{-1}(-x e^{-x})]^\alpha x^{-\beta}\d x, α > − 1, β < 1. For the latter we use known integral representations and their evaluation by nonstandard Gaussian quadrature, if α ≠ β, and explicit formulae involving the trigamma function, if α = β.     

Construction of multivariate compactly supported orthonormal wavelets

We propose a constructive method to find compactly supported orthonormal wavelets for any given compactly supported scaling function φ in the multivariate setting. For simplicity, we start with a standard dilation matrix 2I_2×2 in the bivariate setting and show how to construct compactly supported functions ψ₁,. . .,ψ_n with n>3 such that {2^kψ_j(2^kx−ℓ,2^ky−m), k,ℓ,m∈Z, j=1,. . .,n} is an orthonormal basis for L₂(ℝ²). Here, n is dependent on the size of the support of φ. With parallel processes in modern computer, it is possible to use these orthonormal wavelets for applications. Furthermore, the constructive method can be extended to construct compactly supported multi-wavelets for any given compactly supported orthonormal multi-scaling vector. Finally, we mention that the constructions can be generalized to the multivariate setting. Dedicated to Professor Charles A. Micchelli on the occasion of his 60th birthday Mathematics subject classifications (2000) 42C15, 42C30.     

Diagnostic tests for local coalescences of variables in Boolean functions

The article studies diagnostic tests for local k -fold coalescences of variables in Boolean functions f( [(x)\tilde]ⁿ )( 1 ￡ k ￡ n,  1 ￡ t ￡ 2^2k ) f\left( {{{\tilde{x}}^n}} \right)\left( {1 \leq k \leq n,\;1 \leq t \leq {2^{{2^k}}}} \right) . Upper and lower bounds are proved for the Shannon function of the length of the diagnostic test for local k -fold coalescences generated by the system of functions F_t^k \Phi_t^k . The Shannon function of the length of a complete diagnostic test for local k -fold coalescences behaves asymptotically as 2 ^k (n − k + 1) for n → ∞, k → ∞.     

Partitioning Posets

Given a poset P = (X, ≺ ), a partition X ₁, ..., X _k of X is called an ordered partition of P if, whenever x ∈ X _i and y ∈ X _j with x ≺ y, then i ≤ j. In this paper, we show that for every poset P = (X, ≺ ) and every integer k ≥ 2, there exists an ordered partition of P into k parts such that the total number of comparable pairs within the parts is at most (m − 1)/k, where m ≥ 1 is the total number of edges in the comparability graph of P. We show that this bound is best possible for k = 2, but we give an improved bound, , for k ≥ 3, where c(k) is a constant depending only on k. We also show that, given a poset P = (X, ≺ ) and an integer 2 ≤ k ≤ |X|, we can find an ordered partition of P into k parts that minimises the total number of comparable pairs within parts in time polynomial in the size of P. We prove more general, weighted versions of these results. Supported by an EPSRC doctoral training grant.     

Integral representation of even positive-definite functions of one variable

We obtain an integral representation of even positive-definite functions of one variable for which the kernel [k ₁(x + y) + k ₂ (x − y)] is positive definite.     

Two-Point Boundary Value Problems for Duffing Equations across Resonance

In this paper, we consider the equation y″+f(x,y)=0 with a nonresonance condition of the form A≤f _y (x,y)≤B, where (k−1)²<A≤k ²<⋅⋅⋅<m ²≤B<(m+1)², k,m∈ℤ⁺. With optimal control theory and the Schauder fixed-point theorem, by introducing a new cost functional, we obtain a new existence and uniqueness result for the above equation with two-point boundary-value conditions. This work was supported by NSFC Grant 10501017 and 985 Project of Jilin University.     

The asymptotic formula in Waring’s problem for one square and seventeen fifth powers

Let R _k,s (n) denote the number of solutions of the equation n = x² + y₁^k + y₂^k + ?+ y_s^k{n= x^2 + y_1^k + y_2^k + \cdots + y_s^k} in natural numbers x, y ₁, . . . , y _s . By a straightforward application of the circle method, an asymptotic formula for R _k,s (n) is obtained when k ≥ 3 and s ≥ 2 ^k–1 + 2. When k ≥ 6, work of Heath-Brown and Boklan is applied to establish the asymptotic formula under the milder constraint s ≥ 7 · 2 ^k–4 + 3. Although the principal conclusions provided by Heath-Brown and Boklan are not available for smaller values of k, some of the underlying ideas are still applicable for k = 5, and the main objective of this article is to establish an asymptotic formula for R _5,17(n) by this strategy.     

A bound for the representability of large numbers by ternary quadratic forms and nonhomogeneous waring equations

The solvability of the equation n = x ² + y ² + 6pz ² (p is a fixed large prime) is proved under some natural congruential conditions and the assumption nm ¹² > p ²¹. As an implication, the solvability of the equation n = x ² + y ² + u ³ + v ³ + z ⁴ + w ¹⁶ + t ^4k+1 for all sufficiently large n is established. Bibliography: 13 titles. Translated from Zapiski Nauchnykh Seminarov POMI, Vol. 357, 2008, pp. 5–21.     

A generalized mixed type of quartic–cubic–quadratic–additive functional equations

We determine the general solution of the functional equation f(x + ky) + f(x-ky) = g(x + y) + g(x-y) + h(x) + h(y) for fixed integers with k ≠ 0; ±1 without assuming any regularity conditions for the unknown functions f, g, h, and0020[(h)\tilde] \tilde{h} . The method used for solving these functional equations is elementary but it exploits an important result due to Hosszú. The solution of this functional equation can also be obtained in groups of certain type by using two important results due to Székelyhidi.     

Equations fonctionnelles et estimations de normes de matrices

We solve independently the equations 1/θ(x)θ(y)=ψ(x)−ψ(y)+φ(x−y)/θ(x−y) and 1/θ(x)θ(y)=σ(x)−σ(y)/θ(x−y)+τ(x)τ(y), τ(0)=0. In both cases we find θ²=aθ⁴+bθ²+c. We deduce estimates for the spectral radius of a matrix of type(1/θ(x _r −x _s )) (the accent meaning that the coefficients of the main diagonal are zero) and we study the case where thex _r are equidistant.

Dédié to à Monsieur le Professeur Otto Haupt à l'occasion de son cententiare avec les meilleurs voeux     

Numerical evaluation of Goursat’s infinite integral

The infinite integral ò₀^￥x dx/(1+x⁶sin²x)\int_0^{\infty}x\,dx/(1+x^6\sin^2x) converges but is hard to evaluate because the integrand f(x) = x/(1 + x ⁶sin² x) is a non-convergent and unbounded function, indeed f(kπ) = kπ→ ∞ (k→ ∞). We present an efficient method to evaluate the above integral in high accuracy and actually obtain an approximate value in up to 73 significant digits on an octuple precision system in C++.     

Accelerated hybrid conjugate gradient algorithm with modified secant condition for unconstrained optimization

An accelerated hybrid conjugate gradient algorithm represents the subject of this paper. The parameter β _k is computed as a convex combination of b_k^HS\beta_k^{HS} (Hestenes and Stiefel, J Res Nat Bur Stand 49:409–436, 1952) and b_k^DY\beta_k^{DY} (Dai and Yuan, SIAM J Optim 10:177–182, 1999), i.e. b_k^C = (1-q_k)b_k^HS + q_k b_k^DY\beta_k^C =\left({1-\theta_k}\right)\beta_k^{HS} + \theta_k \beta_k^{DY}. The parameter θ _k in the convex combinaztion is computed in such a way the direction corresponding to the conjugate gradient algorithm is the best direction we know, i.e. the Newton direction, while the pair (s _k , y _k ) satisfies the modified secant condition given by Li et al. (J Comput Appl Math 202:523–539, 2007) B _k + 1 s _k  = z _k , where z_k = y_k +(h_k / || s_k ||² )s_kz_k =y_k +\left({{\eta_k} / {\left\| {s_k} \right\|^2}} \right)s_k, h_k = 2( f_k -f_k+1 )+( g_k +g_k+1 )^Ts_k\eta_k =2\left( {f_k -f_{k+1}} \right)+\left( {g_k +g_{k+1}} \right)^Ts_k, s _k  = x _k + 1 − x _k and y _k  = g _k + 1 − g _k . It is shown that both for uniformly convex functions and for general nonlinear functions the algorithm with strong Wolfe line search is globally convergent. The algorithm uses an acceleration scheme modifying the steplength α _k for improving the reduction of the function values along the iterations. Numerical comparisons with conjugate gradient algorithms show that this hybrid computational scheme outperforms a variant of the hybrid conjugate gradient algorithm given by Andrei (Numer Algorithms 47:143–156, 2008), in which the pair (s _k , y _k ) satisfies the classical secant condition B _k + 1 s _k  = y _k , as well as some other conjugate gradient algorithms including Hestenes-Stiefel, Dai-Yuan, Polack-Ribière-Polyak, Liu-Storey, hybrid Dai-Yuan, Gilbert-Nocedal etc. A set of 75 unconstrained optimization problems with 10 different dimensions is being used (Andrei, Adv Model Optim 10:147–161, 2008).     

The primes contain arbitrarily long polynomial progressions

We establish the existence of infinitely many polynomial progressions in the primes; more precisely, given any integer-valued polynomials P ₁, …, P _k  ∈ Z[m] in one unknown m with P ₁(0) = … = P _k (0) = 0, and given any ε > 0, we show that there are infinitely many integers x and m, with 1 \leqslant m \leqslant x^e1 \leqslant m \leqslant x^\varepsilon, such that x + P ₁(m), …, x + P _k (m) are simultaneously prime. The arguments are based on those in [18], which treated the linear case P _j  = (j − 1)m and ε = 1; the main new features are a localization of the shift parameters (and the attendant Gowers norm objects) to both coarse and fine scales, the use of PET induction to linearize the polynomial averaging, and some elementary estimates for the number of points over finite fields in certain algebraic varieties.     

Geometric computation of the numerical radius of a matrix

We use methods of geometric computing combined with hermitean matrix eigenvalue/eigenvector evaluations to find the numerical radius w(A) of a real or complex square matrix A simply, quickly, and accurately. The numerical radius w(A) is defined as the maximal distance of points in the field of values F(A) = { x^* A x | ||x||₂ = 1 }F(A) = \{ x^* A x \mid \|x\|_2 = 1 \} from zero in ℂ. Its value is an indicator of the transient behavior of the discrete dynamical system f _k + 1 = Af _k . We describe and test a MATLAB code for solving this optimization problem that can have multiple and even infinitely many solutions with maximal distance.     

On the boundedness of one recurrent sequence in a banach space

We establish necessary and sufficient conditions under which a sequence x ₀ = y ₀ , x _n+1 = Ax _n  + y _n+1 , n ≥ 0, is bounded for each bounded sequence { y_n :n \geqslant 0 } ì { x ? è_{n = 1}^￥ D( Aⁿ ) |sup_{n \geqslant 0} || Aⁿ x || < ￥ }\left\{ {y_n :n \geqslant 0} \right\} \subset \left\{ {\left. {x \in \bigcup\nolimits_{n = 1}^\infty {D\left( {A^n } \right)} } \right|\sup _{n \geqslant 0} \left\| {A^n x} \right\| < \infty } \right\}, where A is a closed operator in a complex Banach space with domain of definition D(A) .     

Higher order Riesz transforms associated with Bessel operators

In this paper we investigate Riesz transforms R _μ ^(k) of order k≥1 related to the Bessel operator Δ_μ f(x)=-f”(x)-((2μ+1)/x)f’(x) and extend the results of Muckenhoupt and Stein for the conjugate Hankel transform (a Riesz transform of order one). We obtain that for every k≥1, R _μ ^(k) is a principal value operator of strong type (p,p), p∈(1,∞), and weak type (1,1) with respect to the measure dλ(x)=x ^2μ+1  dx in (0,∞). We also characterize the class of weights ω on (0,∞) for which R _μ ^(k) maps L ^p (ω) into itself and L ¹(ω) into L ^1,∞(ω) boundedly. This class of weights is wider than the Muckenhoupt class of weights for the doubling measure dλ. These weighted results extend the ones obtained by Andersen and Kerman.     

Center conditions III: Parametric and model center problems

We consider an Abel equation (*)y’=p(x)y ² +q(x)y ³ withp(x), q(x) polynomials inx. A center condition for (*) (closely related to the classical center condition for polynomial vector fields on the plane) is thaty ₀=y(0)≡y(1) for any solutiony(x) of (*). Folowing [7], we consider a parametric version of this condition: an equation (**)y’=p(x)y ² +εq(x)y ³ p, q as above, ε ∈ ℂ, is said to have a parametric center, if for any ɛ and for any solutiony(ɛ,x) of (**)y(ɛ, 0)≡y(ɛ, 1).. We give another proof of the fact, shown in [6], that the parametric center condition implies vanishing of all the momentsm _k (1), wherem _k (x)=∫ ₀ ^{x pk} (t)q(t)(dt),P(x)=∫ ₀ ^x p(t)dt. We investigate the structure of zeroes ofm _k (x) and generalize a “canonical representation” ofm _k (x) given in [7]. On this base we prove in some additional cases a composition conjecture, stated in [6, 7] for a parametric center problem. The research of the first and the third author was supported by the Israel Science Foundation, Grant No. 101/95-1 and by the Minerva Foundation.     

Strict Betweennesses Induced by Posets as well as by Graphs

For a finite poset P = (V, ≤ ), let _s(P){\cal B}_s(P) consist of all triples (x,y,z) ∈ V ³ such that either x < y < z or z < y < x. Similarly, for every finite, simple, and undirected graph G = (V,E), let B_s(G){\cal B}_s(G) consist of all triples (x,y,z) ∈ V ³ such that y is an internal vertex on an induced path in G between x and z. The ternary relations B_s(P){\cal B}_s(P) and B_s(G){\cal B}_s(G) are well-known examples of so-called strict betweennesses. We characterize the pairs (P,G) of posets P and graphs G on the same ground set V which induce the same strict betweenness relation B_s(P)=B_s(G){\cal B}_s(P)={\cal B}_s(G).     

Center conditions II: Parametric and model center problems

We consider an Abel equation (*)y’=p(x)y ² +q(x)y ³ withp(x), q(x) polynomials inx. A center condition for (*) (closely related to the classical center condition for polynomial vector fields on the plane) is thaty ₀=y(0)≡y(1) for any solutiony(x) of (*). We introduce a parametric version of this condition: an equation (**)y’=p(x)y ² +εq(x)y ³ p, q as above, ℂ, is said to have a parametric center, if for any ε and for any solutiony(ε,x) of (**),y(ε,0)≡y(ε,1). We show that the parametric center condition implies vanishing of all the momentsm _k (1), wherem _k (x)=∫ ₀ ^{x pk} (t)q(t)(dt),P(x)=∫ ₀ ^x p(t)dt. We investigate the structure of zeroes ofm _k (x) and on this base prove in some special cases a composition conjecture, stated in [10], for a parametric center problem. The research of the first and the third author was supported by the Israel Science Foundation, Grant No. 101/95-1 and by the Minerva Foundation.     

Solution and stability of generalized mixed type cubic, quadratic and additive functional equation in quasi-Banach spaces

In this paper, we achieve the general solution and the generalized Hyers–Ulam–Rassias stability of the following functional equation

f(x+ky)+f(x−ky)=k²f(x+y)+k²f(x−y)+2(1−k²)f(x)