Search

We complete a gap in the proof that exceptional polynomials are complete orthogonal systems in the associated Hilbert spaces.

International audience; In this paper, we consider a semi-linear fourth order parabolic equation in a bounded smooth domain Ω with homogeneous Dirichlet and Neumann boundary conditions. The main result of this paper is the null controllability and the exact controllability to the trajectories at any time T > 0 for the associated control system with a control function acting at the interior.; Dans ce papier, on considère uneéquation parabolique semi-linéaire de quatrième ordre dans un domaine borné régulier Ω avec des conditions aux limites de type Dirichlet et Neumann homogènes. Le résultat principal de ce papier concerne la contrôlabilitéà zéro et la contrôlabilité exacte pour tout T > 0 du système de contrôle associé avec un contrôle agissantà l'interieur.

The purpose of this paper is to show well-posedness results in L 2 -based Sobolev spaces for transmission, Dirichlet, Neumann, and mixed boundary value problems for the Brinkman system with L ∞ coefficients in Lipschitz domains on a compact Riemannian manifold of dimension m ≥ 2 . The Dirichlet, transmission, and mixed problems for the nonlinear Darcy-Forchheimer-Brinkman system with L ∞ coefficients are also analyzed. First, we focus on the well-posedness of linear transmission, Dirichlet and mixed boundary value problems for the Brinkman system with L ∞ coefficients in Lipschitz domains on compact Riemannian manifolds by using a variational approach that reduces such a boundary value problem to a mixed variational formulation defined in terms of two bilinear continuous forms, one of them satisfying a coercivity condition and another one the inf-sup condition. Further, we show the equivalence between each boundary value problem for the Brinkman system with L ∞ coefficients and its mixed variational counterpart, and then the well posedness in L 2 -based Sobolev spaces by using the Necas-Babuska-Brezzi technique. The second goal of this paper is the construction of the Newtonian and layer potential operators for the Brinkman system with L ∞ coefficients in Lipschitz domains on compact Riemannian manifolds by using the well-posedness results for the analyzed linear transmission problems. Various mapping properties of these operators are also obtained and used to describe the weak solutions of the Poisson problems with Dirichlet and Neumann conditions for the nonsmooth Brinkman system in terms of such potentials. Finally, we combine the well-posedness results of the Poisson problems of Dirichlet, transmission, and mixed type for the nonsmooth Brinkman system with a fixed point theorem in order to show the existence of a weak solution of the Poisson problem of Dirichlet, transmission, or mixed type for the (nonlinear) Darcy-Forchheimer-Brinkman system with L ∞ coefficients in L 2 -based Sobolev spaces in Lipschitz domains on compact Riemannian manifolds of dimension m ∈ { 2 , 3 } .

Though the boundary layer formation in the chemotactic process has been observed in experiment (cf. [63] ), the mathematical study on the boundary layer solutions of chemotaxis models is just in its infant stage. Apart from the sophisticated theoretical tools involved in the analysis, how to impose/derive physical boundary conditions is a state-of-the-art in studying the boundary layer problem of chemotaxis models. This paper will proceed with a previous work [24] in one dimension to establish the convergence of boundary layer solutions of the Keller–Segel model with singular sensitivity in a two-dimensional space (half-plane) with respect to the chemical diffusion rate denoted by e ≥ 0 . Compared to the one-dimensional boundary layer problem, there are many new issues arising from multi-dimensions such as possible Prandtl type degeneracy, curl-free preservation and well-posedness of large-data solutions. In this paper, we shall derive appropriate physical boundary conditions and gradually overcome these barriers and hence establish the convergence of boundary layer solutions of the singular Keller–Segel system in the half-plane as the chemical diffusion rate vanishes. Specially speaking, we justify that the boundary layer converges to the outer layer (solution with e = 0 ) plus the inner layer as e → 0 , where both outer and inner layer profiles are precisely derived and well understood. By doing this, the structure of boundary layer solutions is clearly characterized. We hope that our results and methods can shed lights on the understanding of underlying mechanisms of the boundary layer patterns observed in the experiment for chemotaxis such as the work by Tuval et al. [63] , and open a new window in the future theoretical study of chemotaxis models.

Diaconis and Griffiths (2014) study the multivariate Krawtchouk polynomials orthogonal on the multinomial distribution. In this paper we derive the reproducing kernel orthogonal polynomials Q n ( x , y ; N , p ) on the multinomial distribution which are sums of products of orthonormal polynomials in x and y of fixed total degree n = 0 , 1 , … , N . The Poisson kernel ∑ n = 0 N ρ n Q n ( x , y ; N , p ) arises naturally from a probabilistic argument. An application to a multinomial goodness of fit test is developed, where the chi-squared test statistic is decomposed into orthogonal components which test the order of fit. A new duplication formula for the reproducing kernel polynomials in terms of the 1-dimensional Krawtchouk polynomials is derived. The duplication formula allows a Lancaster characterization of all reversible Markov chains with a multinomial stationary distribution whose eigenvectors are multivariate Krawtchouk polynomials and where eigenvalues are repeated within the same total degree. The χ 2 cutoff time, and total variation cutoff time is investigated in such chains. Emphasis throughout the paper is on a probabilistic understanding of the polynomials and their applications, particularly to Markov chains.

From spline theory it is well-known that univariate cubic spline interpolation, if carried out in its natural Hilbert space W 2 2 [ a , b ] and on point sets with fill distance h , converges only like O ( h 2 ) in L 2 [ a , b ] if no additional assumptions are made. But superconvergence up to order h 4 occurs if more smoothness is assumed and if certain additional boundary conditions are satisfied. This phenomenon was generalized in 1999 to multivariate interpolation in Reproducing Kernel Hilbert Spaces on domains Ω ⊂ R d for continuous positive definite Fourier-transformable shift-invariant kernels on R d . But the sufficient condition for superconvergence given in 1999 still needs further analysis, because the interplay between smoothness and boundary conditions is not clear at all. Furthermore, if only additional smoothness is assumed, superconvergence is numerically observed in the interior of the domain, but a theoretical foundation still is a challenging open problem. This paper first generalizes the “improved error bounds” of 1999 by an abstract theory that includes the Aubin–Nitsche trick and the known superconvergence results for univariate polynomial splines. Then the paper analyzes what is behind the sufficient conditions for superconvergence. They split into conditions on smoothness and localization, and these are investigated independently. If sufficient smoothness is present, but no additional localization conditions are assumed, it is numerically observed that superconvergence always occurs in the interior of the domain, and some supporting arguments are provided. If smoothness and localization interact in the kernel-based case on R d , weak and strong boundary conditions in terms of pseudodifferential operators occur. A special section on Mercer expansions is added, because Mercer eigenfunctions always satisfy the sufficient conditions for superconvergence. Numerical examples illustrate the theoretical findings.

It has been shown that natural interval extensions (NIE) can be used to calculate the largest positive Lyapunov exponent (LLE). However, the elaboration of NIE are not always possible for some dynamical systems, such as those modelled by simple equations or by Simulink-type blocks. In this paper, we use rounding mode of floating-point numbers to compute the LLE. We have exhibited how to produce two pseudo-orbits by means of different rounding modes; these pseudo-orbits are used to calculate the Lower Bound Error (LBE). The LLE is the slope of the line gotten from the logarithm of the LBE, which is estimated by means of a recursive least square algorithm (RLS). The main contribution of this paper is to develop a procedure to compute the LLE based on the LBE without using the NIE. Additionally, with the aid of RLS the number of required points has been decreased. Eight numerical examples are given to show the effectiveness of the proposed technique.

Motivated by the paper Gerard-Varet and Dormy (2010) [6] [JAMS, 2010] about the linear ill-posedness for the Prandtl equations around a shear flow with exponential decay in normal variable, and the recent study of well-posedness on the Prandtl equations in Sobolev spaces, this paper aims to extend the result in [6] to the case when the shear flow has general decay. The key observation is to construct an approximate solution that captures the initial layer to the linearized problem motivated by the precise formulation of solutions to the inviscid Prandtl equations.

Reduction theory has played a major role in the study of Hamiltonian systems. Whilst the Hamilton–Jacobi theory is one of the main tools to integrate the dynamics of certain Hamiltonian problems and a topic of research on its own. Moreover, the construction of several symplectic integrators relies on approximations of a complete solution of the Hamilton–Jacobi equation. The natural question that we address in this paper is how these two topics (reduction and Hamilton–Jacobi theory) fit together. We obtain a reduction and reconstruction procedure for the Hamilton–Jacobi equation with symmetries, even in a generalized sense to be clarified below. Several applications and relations to other reduction of the Hamilton–Jacobi theory are shown in the last section of the paper. It is remarkable that as by-product we obtain a generalization of the Ge–Marsden reduction procedure [18] and the results in [17] . Quite surprisingly, the classical ansatze available in the literature to solve the Hamilton–Jacobi equation (see [2] , [19] ) are also particular instances of our framework.

We study the approximation of compact linear operators defined over certain weighted tensor product Hilbert spaces. The information complexity is defined as the minimal number of arbitrary linear functionals needed to obtain an e -approximation for the d -variate problem which is fully determined in terms of the weights and univariate singular values. Exponential tractability means that the information complexity is bounded by a certain function that depends polynomially on d and logarithmically on e − 1 . The corresponding unweighted problem was studied in Hickernell et al. (2020) with many negative results for exponential tractability. The product weights studied in the present paper change the situation. Depending on the form of polynomial dependence on d and logarithmic dependence on e − 1 , we study exponential strong polynomial, exponential polynomial, exponential quasi-polynomial, and exponential ( s , t ) -weak tractability with max ( s , t ) ≥ 1 . For all these notions of exponential tractability, we establish necessary and sufficient conditions on weights and univariate singular values for which it is indeed possible to achieve the corresponding notion of exponential tractability. The case of exponential ( s , t ) -weak tractability with max ( s , t ) 1 is left for future study. The paper uses some general results obtained in Hickernell et al. (2020) and Kritzer and Woźniakowski (2019).

This paper studies the approximation of the Hilbert transform f ? = H f of continuous functions f with continuous conjugate f ? based on a finite number of samples. It is known that every sequence { H N f } N ? N which approximates f ? from samples of f diverges (weakly) with respect to the uniform norm. This paper conjectures that all of these approximation sequences even contain no convergent subsequence. A property which is termed strong divergence.The conjecture is supported by two results. First it is proven that the sequence of the sampled conjugate Fejer means diverges strongly. Second, it is shown that for every sample based approximation method { H N } N ? N there are functions f such that ? H N f ? ∞ exceeds any given bound for any given number of consecutive indices N .As an application, the later result is used to investigate a problem associated with a question of Ul'yanov on Fourier series which is related to the possibility to construct adaptive approximation methods to determine the Hilbert transform from sampled data. This paper shows that no such approximation method with a finite search horizon exists.

This paper considers the problem of approximating the spectral factor of continuous spectral densities with finite Dirichlet energy based on finitely many samples of these spectral densities. Although there exists a closed form expression for the spectral factor, this formula shows a very complicated behavior because of the non-linear dependency of the spectral factor from spectral density and because of a singular integral in this expression. Therefore approximation methods are usually applied to calculate the spectral factor. It is shown that there exists no sampling-based method which depends continuously on the samples and which is able to approximate the spectral factor for all densities in this set. Instead, to any sampling-based approximation method there exists a large set of spectral densities so that the approximation method does not converge to the spectral factor for every spectral density in this set as the number of available sampling points is increased. The paper will also show that the same results hold for sampling-based algorithms for the calculation of the outer function in the theory of Hardy spaces.

In this paper, we obtain an explicit formula for the discrepancy between the limit of the second inner variations of p -Laplace Allen–Cahn energies and the second inner variation of their Γ -limit which is the area functional. Our analysis explains the mysterious discrepancy term found in our previous paper [8] in the case p = 2 . The discrepancy term turns out to be related to the convergence of certain 4-tensors which are absent in the usual Allen–Cahn functional. These (hidden) 4-tensors suggest that, in the complex-valued Ginzburg–Landau setting, we should expect a different discrepancy term which we are able to identify. Along the way, we partially answer a question of Kohn and Sternberg [6] by giving a relation between the limit of second variations of the Allen–Cahn functional and the second inner variation of the area functional at local minimizers. Moreover, our analysis reveals an interesting identity connecting second inner variation and Poincare inequality for area-minimizing surfaces with volume constraint in the work of Sternberg and Zumbrun [16] .

In this paper we investigate quasi-Monte Carlo (QMC) integration using digital nets over Z b in reproducing kernel Hilbert spaces. The tent transformation (previously called baker’s transform) was originally used for lattice rules by Hickernell (2002) to achieve higher order convergence of the integration error for smooth non-periodic integrands, and later, has been successfully applied to digital nets over Z 2 by Cristea et al. (2007) and Goda (2015). The aim of this paper is to generalize the latter two results to digital nets over Z b for an arbitrary prime b . For this purpose, we introduce the b -adic tent transformation for an arbitrary positive integer b greater than 1, which is a generalization of the original (dyadic) tent transformation. Further, again for an arbitrary positive integer b greater than 1, we analyze the mean square worst-case error of QMC rules using digital nets over Z b which are randomly digitally shifted and then folded using the b -adic tent transformation in reproducing kernel Hilbert spaces. Using this result, for a prime b , we prove the existence of good higher order polynomial lattice rules over Z b among a smaller number of candidates as compared to the result by Dick and Pillichshammer (2007), which achieve almost the optimal convergence rate of the mean square worst-case error in unanchored Sobolev spaces of smoothness of arbitrary high order.

In this paper, we introduce and study an implicit iterative method to approximate a common solution of split equilibrium problem and fixed point problem for a nonexpansive semigroup in real Hilbert spaces. Further, we prove that the nets generated by the implicit iterative method converge strongly to the common solution of split equilibrium problem and fixed point problem for a nonexpansive semigroup. This common solution is the unique solution of a variational inequality problem and is the optimality condition for a minimization problem. Furthermore, we justify our main result through a numerical example. The results presented in this paper extend and generalize the corresponding results given by Plubtieng and Punpaeng [S. Plubtieng, R. Punpaeng, Fixed point solutions of variational inequalities for nonexpansive semigroups in Hilbert spaces, Math. Comput. Model. 48 (2008) 279–286] and Cianciaruso et al. [F. Cianciaruso, G. Marino, L. Muglia, Iterative methods for equilibrium and fixed point problems for nonexpansive semigroups in Hilbert space, J. Optim. Theory Appl. 146 (2010) 491–509].

We say that the polynomial sequence $(Q^{(\lambda)}_n)$ is a semiclassical Sobolev polynomial sequence when it is orthogonal with respect to the inner product $$ _S= +\lambda , $$ where ${\bf u}$ is a semiclassical linear functional, ${\mathscr D}$ is the differential, the difference or the $q$--difference operator, and $\lambda$ is a positive constant. In this paper we get algebraic and differential/difference properties for such polynomials as well as algebraic relations between them and the polynomial sequence orthogonal with respect to the semiclassical functional $\bf u$. The main goal of this article is to give a general approach to the study of the polynomials orthogonal with respect to the above nonstandard inner product regardless of the type of operator ${\mathscr D}$ considered. Finally, we illustrate our results by applying them to some known families of Sobolev orthogonal polynomials as well as to some new ones introduced in this paper for the first time., Comment: 23 pages, special issue dedicated to Professor Guillermo Lopez lagomasino on the occasion of his 60th birthday, accepted in Journal of Approximation Theory

In this paper our aim is to show some mean value inequalities for the modified Bessel functions of the first and second kinds. Our proofs are based on some bounds for the logarithmic derivatives of these functions, which are in fact equivalent to the corresponding Tur\'an type inequalities for these functions. As an application of the results concerning the modified Bessel function of the second kind we prove that the cumulative distribution function of the gamma-gamma distribution is log-concave. At the end of this paper several open problems are posed, which may be of interest for further research., Comment: 14 pages

In this paper, we rigorously prove the existence and stability of single-peaked patterns for the singularly perturbed Gierer–Meinhardt system on a compact two-dimensional Riemannian manifold without boundary which are far from spatial homogeneity. Throughout the paper we assume that the activator diffusivity ϵ 2 is small enough. We show that for the threshold ratio D ∼ 1 ϵ 2 of the activator diffusivity ϵ 2 and the inhibitor diffusivity D, the Gaussian curvature and the Green's function interact. A convex combination of the Gaussian curvature and the Green's function together with their derivatives are linked to the peak locations and the o ( 1 ) eigenvalues. A nonlocal eigenvalue problem (NLEP) determines the O ( 1 ) eigenvalues which all have negative part in this case.

This paper illustrates an approach to generate multi-wing attractors in coupled Lorenz systems. In particular, novel four-wing (eight-wing) hyperchaotic attractors are generated by coupling two (three) identical Lorenz systems. The paper shows that the equilibria of the proposed systems have certain symmetries with respect to specific coordinate planes and the eigenvalues of the associated Jacobian matrices exhibit the property of similarity. In analogy with the original Lorenz system, where the two-wings of the butterfly attractor are located around the two equilibria with the unstable pair of complex-conjugate eigenvalues, this paper shows that the four-wings (eight-wings) of these attractors are located around the four (eight) equilibria with two (three) pairs of unstable complex-conjugate eigenvalues.

Comprehensive numerical simulations (reviewed in Dissipative Solitons, Akhmediev and Ankiewicz (Eds.), Springer, Berlin, 2005) of pulse solutions of the cubic–quintic Ginzburg–Landau Equation (CGLE), a canonical equation governing the weakly nonlinear behavior of dissipative systems in a wide variety of disciplines, reveal various intriguing and entirely novel classes of solutions. In particular, there are five new classes of pulse or solitary waves solutions, viz. pulsating, creeping, snake, erupting, and chaotic solitons. In contrast to the regular solitary waves investigated in numerous integrable and non-integrable systems over the last three decades, these dissipative solitons are not stationary in time. Rather, they are spatially confined pulse-type structures whose envelopes exhibit complicated temporal dynamics. The numerical simulations also reveal very interesting bifurcations sequences of these pulses as the parameters of the CGLE are varied. In this paper, we address the issues of central interest in the area, i.e., the conditions for the occurrence of the five categories of dissipative solitons, as well the dependence of both their shape and their stability on the various parameters of the CGLE, viz. the nonlinearity, dispersion, linear and nonlinear gain, loss and spectral filtering parameters. Our predictions on the variation of the soliton amplitudes, widths and periods with the CGLE parameters agree with simulation results. First, we elucidate the Hopf bifurcation mechanism responsible for the various pulsating solitary waves, as well as its absence in Hamiltonian and integrable systems where such structures are absent. Next, we develop and discuss a variational formalism within which to explore the various classes of dissipative solitons. Given the complex dynamics of the various dissipative solutions, this formulation is, of necessity, significantly generalized over all earlier approaches in several crucial ways. Firstly, the starting formulation for the Lagrangian is recent and not well explored. Also, the trial functions have been generalized considerably over conventional ones to keep the shape relatively simple (and the trial function integrable!) while allowing arbitrary temporal variation of the amplitude, width, position, speed and phase of the pulses. In addition, the resulting Euler–Lagrange equations are treated in a completely novel way. Rather than consider the stable fixed points which correspond to the well-known stationary solitons or plain pulses, we use dynamical systems theory to focus on more complex attractors viz. periodic, quasiperiodic, and chaotic ones. Periodic evolution of the trial function parameters on stable periodic attractors yield solitons whose amplitudes and widths are non-stationary or time dependent. In particular, pulsating and snake dissipative solitons may be treated in this manner. Detailed results are presented here for the pulsating solitary waves – their regimes of occurrence, bifurcations, and the parameter dependences of the amplitudes, widths, and periods agree with simulation results. Snakes and chaotic solitons will be addressed in subsequent papers. This overall approach fails only to address the fifth class of dissipative solitons, viz. the exploding or erupting solitons.

The cyclic projections algorithm is an important method for determining a point in the intersection of a finite number of closed convex sets in a Hilbert space. That is, for determining a solution to the ''convex feasibility'' problem. This is the third paper in a series on a study of the rate of convergence for the cyclic projections algorithm. In the first of these papers, we showed that the rate could be described in terms of the ''angles'' between the convex sets involved. In the second, we showed that these angles often had a more tractable formulation in terms of the ''norm'' of the product of the (nonlinear) metric projections onto related convex sets. In this paper, we show that the rate of convergence of the cyclic projections algorithm is also intimately related to the ''linear regularity property'' of Bauschke and Borwein, the ''normal property'' of Jameson (as well as Bakan, Deutsch, and Li's generalization of Jameson's normal property), the ''strong conical hull intersection property'' of Deutsch, Li, and Ward, and the rate of convergence of iterated parallel projections. Such properties have already been shown to be important in various other contexts as well.

A quasimode for a positive, symmetric and compact operator on a Hilbert space could be defined as a pair ( u , λ ), where u is a function approaching a certain linear combination of eigenfunctions associated with the eigenvalues of the operator in a “small interval” [ λ − r , λ + r ] . Its value in describing asymptotics for low and high frequency vibrations in certain singularly perturbed spectral problems, which depend on a small parameter e, has been made clear recently in many papers. In this paper, considering second order evolution problems, we provide estimates for the time t in which standing waves of the type e i λ t u approach their solutions u ( t ) when the initial data deal with quasimodes of the associated operators. We establish a general abstract framework and we extended it to the case where operators and spaces depend on the small parameter e: now λ and u can depend on e and also perform the estimates for t. We apply the results to vibrating systems with concentrated masses.

The purpose of this paper is to derive product representations for generalizations of the Rogers–Ramanujan series. Special cases of the results presented here were first stated by Ramanujan in the “Lost Notebook” and proved by George Andrews. The analysis used in this paper is based upon the work of Andrews and the broad contributions made by Mourad Ismail and Walter Hayman. Each series considered is related to an extension of the Rogers–Ramanujan continued fraction and corresponds to an orthogonal polynomial sequence generalizing classical orthogonal sequences. Using Ramanujan's differential equations for Eisenstein series and corresponding analogues derived by V. Ramamani, the coefficients in the series representations of each zero are expressed in terms of certain Eisenstein series.

The aim of this paper is the investigation of the error which results from the method of approximate approximations applied to functions defined on compact intervals, only. This method, which is based on an approximate partition of unity, was introduced by Maz’ya in 1991 and has mainly been used for functions defined on the whole space up to now. For the treatment of differential equations and boundary integral equations, however, an efficient approximation procedure on compact intervals is needed.In the present paper we apply the method of approximate approximations to functions which are defined on compact intervals. In contrast to the whole space case here a truncation error has to be controlled in addition. For the resulting total error pointwise estimates and L1-estimates are given, where all the constants are determined explicitly.

This paper gives the conditions necessary for weighted convergence of Hermite–Fejér interpolation for a general class of even weights which are of exponential decay on the real line or at the end points of (-1,1). The results of this paper guarantee that the conditions of Theorem 2.3 in [11] are optimal.

This paper studies recent modifications of the limited memory BFGS (L-BFGS) method for solving large scale unconstrained optimization problems. Each modification technique attempts to improve the quality of the L-BFGS Hessian by employing (extra) updates in a certain sense. Because at some iterations these updates might be redundant or worsen the quality of this Hessian, this paper proposes an updates criterion to measure this quality. Hence, extra updates are employed only to improve the poor approximation of the L-BFGS Hessian. The presented numerical results illustrate the usefulness of this criterion and show that extra updates improve the performance of the L-BFGS method substantially.

We consider the Schrodinger equations with arbitrary (large) power non-linearity on the three-dimensional torus. We construct non-trivial probability measures supported on Sobolev spaces and show that the equations are globally well-posed on the supports of these measures, respectively. Moreover, these measures are invariant under the flows that are constructed. Therefore, the constructed solutions are recurrent in time. Also, we show slow growth control on the time evolution of the solutions. A generalization to any dimension is given. Our proof relies on a new approach combining the fluctuation-dissipation method and some features of the Gibbs measures theory for Hamiltonian PDEs. The strategy of the paper applies to other contexts.

In this paper we define a notion of uniform distribution and discrepancy of sequences in an abstract set E through reproducing kernel Hilbert spaces of functions on E. In the case of the finite-dimensional unit cube these discrepancies are very closely related to the worst case error obtained for numerical integration of functions in a reproducing kernel Hilbert space. In the compact case we show that the discrepancy tends to zero if and only if the sequence is uniformly distributed in our sense. Next we prove an existence theorem for such uniformly distributed sequences and investigate the relation to the classical notion of uniform distribution. Some examples conclude this paper.

This paper presents an analysis of the dynamical behaviour of a non-symmetric oscillator with piecewise-linearity. The Chen–Langford (C–L) method is used to obtain the averaged system of the oscillator. Using this method, the local bifurcation and the stability of the steady-state solutions are studied. A Runge–Kutta method, Poincaré map and the largest Lyapunov’s exponent are used to detect the complex dynamical phenomena of the system. It is found that the system with piecewise-linearity exhibits periodic oscillations, period-doubling, period-3 solution and then chaos. When chaos is found, it is detected by examining the phase plane, bifurcation diagram and the largest Lyapunov’s exponent. The results obtained in this paper show that the vibration system with piecewise-linearity do exhibit quite similar dynamical behaviour to the discrete system given by the logistic map.

We discuss and examine Weierstrass’ main contributions to approximation theory. §1. Weierstrass This is a story about Karl Wilhelm Theodor Weierstrass (Weierstras), what he contributed to approximation theory (and why), and some of the consequences thereof. We start this story by relating a little about the man and his life. Karl Wilhelm Theodor Weierstrass was born on October 31, 1815 at Ostenfelde near Munster into a liberal (in the political sense) Catholic family. He was the eldest of four children, none of whom married. Weierstrass was a very successful gymnasium student and was subsequently sent by his father to the University of Bonn to study commerce and law. His father seems to have had in mind a government post for his son. However neither commerce nor law was to his liking, and he “wasted” four years there, not graduating. Beer and fencing seem to have been fairly high on his priority list at the time. The young Weierstrass returned home, and after a period of “rest”, was sent to the Academy at Munster where he obtained a teacher’s certificate. At the Academy he fortuitously came under the tutelage and personal guidance of C. Gudermann who was professor of mathematics at Munster and whose basic mathematical love and interest was the subject of elliptic functions and power series. This interest he was successful in conveying to Weierstrass. In 1841 Weierstrass received his teacher’s certificate, and then spent the next 13 years as a teacher (for 6 years he was a teacher in a pregymnasium in the town of Deutsch-Krone (West Prussia), then for another 7 years in a gymnasium in Braunsberg (East Prussia)). During this period he continued learning mathematics, mainly by studying the work of Abel. He also published some mathematical papers. However these appeared in school journals and were quite naturally not discovered at that time by any who could understand or appreciate them. (Weierstrass’ collected works contain 7 papers from before 1854, the first of which On the development of modular functions (49 pp.) was written in 1840.) In 1854 Weierstrass published the paper On the theory of Abelian functions in Crelle’s Journal fur die Reine und Angewandte Mathematik (the first mathematical research journal, founded in 1826, and now referred to without Crelle’s name in the formal title). It created a sensation within the mathematical community. Here was a 39 year old school teacher whom no one within the mathematical community had heard of. And he had written a masterpiece, not only in its depth, but also in its mastery of an area. Recognition was

In this paper we obtain the Plancherel formula for the spaces of L 2 -sections of line bundles over the complex projective hyperboloids G=H withGD U.p;qIC/ andHD U.1IC/ U.p 1;qIC/. The Plancherel formula is given in an explicit form by means of spherical distributions associated with a character of the subgroupH . We obtain the Plancherel formula by a special method which is also suitable for other problems, for example, for quantization in the spirit of Berezin. © 2000 Editions scientifiques et medicales Elsevier SAS Keywords: Plancherel formula, Spherical distributions, Complex hyperbolic spaces In this paper we obtain the Plancherel formula for the spaces of L 2 -sections of line bundles over the complex projective hyperboloids G=H with GD U.p;qIC/ and H D U.1IC/ U.p 1;qIC/, i.e. we present the decomposition of L 2 into irreducible represen- tations of the group G of class . In order to leave aside the well-known case of a hyperboloid with compact stabilizer subgroup, see (14), we assumep>1 ;q >0. The Plancherel formula is given in an explicit form by means of spherical distributions associated with a character of the subgroup H. We obtain the Plancherel formula by Molchanov's method, see (9). Namely, we follow the detailed scheme in (1), Sections 4, 7. This method deals with the spectral resolution of the radial part of the Laplace operator. The essential step is setting the boundary conditions at certain special points. Those conditions are prescribed by the behaviour of spherical distributions. Finally, it is necessary to perform various analytic continuations. This method is also suitable for other problems, for example, for quantization in the spirit of Berezin, namely, for the decomposition of the Berezin form. It is therefore why this method has to be preferred to the existing methods, described in (3). We use our results from (13). There we define -spherical distributions, study their asymptotic behaviour and express them by means of hypergeometric functions. We describe the irreducible unitary representations of the group G ,o f class associated with an isotropic cone. We give constructions for the Fourier and Poisson transform, define intertwining operators and diagonalize them. Some of those results are presented in Section 1.

The main result of this paper provides uniform decay rates obtained for the energy function associated with a three-dimensional structural acoustic model described by coupled system consisting of the wave equation and plate equation with the coupling on the interface between the acoustic chamber and the wall. The uniform stabilization is achieved by introducing a nonlinear dissipation acting via boundary forces applied at the edge of the plate and viscous or boundary damping applied to the wave equation. The results obtained in this paper extend, to the non-analytic, hyperbolic-like setting, the results obtained previously in the literature for acoustic problems modeled by structurally damped plates (governed by analytic semigroups). As a bypass product, we also obtain optimal uniform decay rates for the Euler Bernoulli plate equations with nonlinear boundary dissipation acting via shear forces only and without (i) any geometric conditions imposed on the domain ,(ii) any growth conditions at the origin imposed on the nonlinear function. This is in contrast with the results obtained previously in the literature ([22] and references therein).

We establish a connection between the theory of Lyapunov exponents and the properties of expansivity and sensitivity to initial conditions for a particular class of discrete time dynamical systems; cellular automata (CA). The main contribution of this paper is the proof that all expansive cellular automata have positive Lyapunov exponents for almost all the phase space configurations. In addition, we provide an elementary proof of the non-existence of expansive CA in any dimension greater than 1. In the second part of this paper we prove that expansivity in dimension greater than 1 can be recovered by restricting the phase space to asuitablesubset of the whole space. To this extent we describe a 2-dimensional CA which is expansive over adense uncountablesubset of the whole phase space. Finally, we highlight the different behavior of expansive and sensitive CA for what concerns the speed at which perturbations propagate.

We show that if the Generalized Riemann Hypothesis is true, the problem of deciding whether a system of polynomial equations in several complex variables has a solution is in the second level of the polynomial hierarchy. In fact, this problem is in AM, the ``Arthur-Merlin'''' class (recall that $\np \subseteq \am \subseteq \rp^{\tiny \np} \subseteq \Pi_2$). The best previous bound was PSPACE. An earlier version of this paper was distributed as NeuroCOLT Technical Report~96-44. The present paper includes in particular a new lower bound for unsatisfiable systems, and remarks on the Arthur-Merlin class.

In this paper we study sequences of linear operators which are "almost positive" outside sets of small Jordan measure. For them, we prove Korovkin-type theorems in terms of a modification of the R-convergence used previously by W. Dickmeis, K. Mevissen, R. J. Nessel, and E. Van Wickeren and the test families of functions which the author introduced in a previous paper.

In previous papers the first author extended the classical equiconvergence theorem of Walsh by the application of summability methods in order to enlarge the disk of equiconvergence to regions of equisummability. The aim of this paper is to study equisummability of sequences which arise from Hadamard products of a fixed power series with Lagrange polynomial interpolants.

We consider the following k-coupled nonlinear Schrodinger systems: { − Δ u j + λ j u j = μ j u j 3 + ∑ i = 1 , i ≠ j k β i , j u i 2 u j in R N , u j > 0 in R N , u j ( x ) → 0 as | x | → + ∞ , j = 1 , 2 , ⋯ , k , where N ≤ 3 , k ≥ 3 , λ j , μ j > 0 are constants and β i , j = β j , i ≠ 0 are parameters. There have been intensive studies for the above systems when k = 2 or the systems are purely attractive ( β i , j > 0 , ∀ i ≠ j ) or purely repulsive ( β i , j 0 , ∀ i ≠ j ); however very few results are available for k ≥ 3 when the systems admit mixed couplings and the components are organized into groups, i.e., there exist ( i 1 , j 1 ) and ( i 2 , j 2 ) such that β i 1 , j 1 > 0 and β i 2 , j 2 0 . In this paper we give the first systematic and an (almost) complete study on the existence of ground states when the systems admit mixed couplings and the components are organized into groups. We first divide these systems into repulsive-mixed and total-mixed cases. In the first case we prove nonexistence of ground states. In the second case we give a necessary condition for the existence of ground states and also provide estimates for Morse index. The key idea is the block decomposition of the systems (optimal block decompositions, eventual block decompositions), and the measure of total interaction forces between different blocks. Finally the assumptions on the existence of ground states are shown to be optimal in some special cases.

In this paper, we study large time behavior of complex-valued solutions to nonlinear Klein-Gordon equation with a gauge invariant quadratic nonlinearity in two spatial dimensions. To find a possible asymptotic behavior, we consider the final value problem. It turns out that one possible behavior is a linear solution with a logarithmic phase correction as in the real-valued case. However, the shape of the logarithmic correction term has one more parameter which is also given by the final data. In the real case the parameter is constant so one cannot see its effect. However, in the complex case it varies in general. The one dimensional case is also discussed., Comment: 25 papges, 2 figures

Bessel processes $(X_{t,k})_{t\ge0}$ in $N$ dimensions are classified via associated root systems and multiplicity constants $k\ge0$. They describe interacting Calogero-Moser-Suther\-land particle systems with $N$ particles and are related to $\beta$-Hermite and $\beta$-Laguerre ensembles. Recently, several central limit theorems were derived for fixed $t>0$, fixed starting points, and $k\to\infty$. In this paper we extend the CLT in the A-case from start in 0 to arbitrary starting distributions by using a limit result for the corresponding Bessel functions. We also determine the eigenvalues and eigenvectors of the covariance matrices of the Gaussian limits and study applications to CLTs for the intermediate particles for $k\to\infty$ and then $N\to\infty$., Comment: 20 pages

We make systematic developments on Lawson-Osserman constructions relating to the Dirichlet problem (over unit disks) for minimal surfaces of high codimension in their 1977 Acta paper. In particular, we show the existence of boundary functions for which infinitely many analytic solutions and at least one nonsmooth Lipschitz solution exist simultaneously. This newly-discovered amusing phenomenon enriches the understanding on the Lawson-Osserman philosophy., Supercedes arXiv:1610.08162. To appear in the Journal de Math\'ematiques Pures et Appliqu\'ees

We prove a partial fraction decomposition of a quotient of two functions E α ( i t x ) and I α ( i t ) which are defined in terms of the Bessel functions J α and J α + 1 of the first kind. This expansion leads naturally to the introduction of an orthonormal system with respect to the measure | x | 2 α + 1 d x 2 α + 1 Γ ( α + 1 ) in [ − 1 , 1 ] , which we call the Fourier–Dunkl system of the second kind. Euler–Dunkl polynomials E n , α ( x ) of degree n are defined by considering E α ( t x ) ∕ I α ( t ) as a generating function. It is shown that the sum ∑ m = 1 ∞ 1 ∕ j m , α 2 k , where j m , α are the positive zeros of J α , is equal (up to an explicit factor) to E 2 k − 1 , α ( 1 ) . For α = 1 ∕ 2 this leads to classical results of Euler since the function E 1 ∕ 2 ( x ) is the exponential function and E n , 1 ∕ 2 ( x ) are (essentially) the Euler polynomials. In the second part of the paper a sampling theorem of Whittaker–Shannon–Kotel’nikov type is established which is strongly related to the above-mentioned partial decomposition and which holds for all functions in the Payley–Wiener space defined by the Dunkl transform in [ − 1 , 1 ] .

The study of the restricted isometry property (RIP) of corrupted random matrices is particularly important in the field of compressed sensing (CS) with corruptions. If a matrix still satisfies the RIP after that a certain portion of rows are erased, then we say that this matrix has the strong restricted isometry property (SRIP). In the field of compressed sensing, random matrices which satisfy certain moment conditions are of particular interest. Among these matrices, those with entries generated from i.i.d. Gaussian or symmetric Bernoulli random variables are often typically considered. Recent studies have shown that matrices with entries generated from i.i.d. Gaussian random variables satisfy the SRIP under arbitrary erasure of rows with high probability. In this paper, we study the erasure robustness property of Bernoulli random matrices. Our main result shows that with overwhelming probability, the SRIP holds for Bernoulli random matrices. Moreover, our analysis leads to a robust version of the famous Johnson–Lindenstrauss lemma for Bernoulli random matrices.

The Boussinesq a b c d system is a 4-parameter set of equations posed in R t × R x , originally derived by Bona, Chen and Saut [11] , [12] as first order 2-wave approximations of the incompressible and irrotational, two dimensional water wave equations in the shallow water wave regime, in the spirit of the original Boussinesq derivation [17] . Among many particular regimes, depending each of them in terms of the value of the parameters ( a , b , c , d ) present in the equations, the generic regime is characterized by the setting b , d > 0 and a , c 0 . If additionally b = d , the a b c d system is Hamiltonian. The equations in this regime are globally well-posed in the energy space H 1 × H 1 , provided one works with small solutions [12] . In this paper, we investigate decay and the scattering problem in this regime, which is characterized as having (quadratic) long-range nonlinearities, very weak linear decay O ( t − 1 / 3 ) because of the one dimensional setting, and existence of non scattering solutions (solitary waves). We prove, among other results, that for a sufficiently dispersive a b c d systems (characterized only in terms of parameters a , b and c), all small solutions must decay to zero, locally strongly in the energy space, in proper subset of the light cone | x | ≤ | t | . We prove this result by constructing three suitable virial functionals in the spirit of works [27] , [28] , and more precisely [42] (valid for the simpler scalar “good Boussinesq” model), leading to global in time decay and control of all local H 1 × H 1 terms. No parity nor extra decay assumptions are needed to prove decay, only small solutions in the energy space.

Let $L$ be the generator of an analytic semigroup whose kernels satisfy Gaussian upper bounds and H\"older's continuity. Also assume that $L$ has a bounded holomorphic functional calculus on $L^2(\mathbb{R}^n)$. In this paper, we construct a frame decomposition for the functions belonging to the Hardy space $H_{L}^{1}(\mathbb{R}^n)$ associated to $L$, and for functions in the Lebesgue spaces $L^p$, $1, Comment: 37 pages, to appear in Journal of Approximation Theory

In this paper we elucidate the connection between various notions of differentiability in the Wasserstein space: some have been introduced intrinsically (in the Wasserstein space, by using typical objects from the theory of Optimal Transport) and used by various authors to study gradient flows, Hamiltonian flows, and Hamilton–Jacobi equations in this context. Another notion is extrinsic and arises from the identification of the Wasserstein space with the Hilbert space of square-integrable random variables on a non-atomic probability space. As a consequence, the classical theory of well-posedness for viscosity solutions for Hamilton–Jacobi equations in infinite-dimensional Hilbert spaces is brought to bear on well-posedness for Hamilton–Jacobi equations in the Wasserstein space.

Let ζ ( s ) be the Riemann zeta function. In 1911, Bohr showed that the set { ζ ( σ + i τ ) : σ > 1 , τ ∈ R } is dense in ℂ . By Voronin’s denseness theorems in 1972, the sets { ( ζ ( σ + i λ 1 + i τ ) , … , ζ ( σ + i λ n + i τ ) ) : σ ≥ 1 , τ ∈ R } with distinct λ 1 , … , λ n ∈ R and { ( ζ ( σ + i τ ) , ζ ′ ( σ + i τ ) , … , ζ ( n − 1 ) ( σ + i τ ) ) : σ ≥ 1 , τ ∈ R } are dense in ℂ n . By Voronin’s universality theorem, for any fixed 1 ∕ 2 σ 1 and any non-negative integer k , the set { ζ σ , τ ( k ) : τ ∈ R } is dense in C [ a , b ] , where ζ σ , τ ( k ) ( t ) ≔ ζ ( k ) ( σ + i t + i τ ) , t ∈ [ a , b ] . In the present paper, we prove that the set { ζ σ , τ ( k ) : σ ≥ 1 , τ ∈ R } ∩ C [ a , b ] is not dense in C [ a , b ] .

In this paper, we study the problem of whether the best approximation to any x in a real Hilbert space X from the closed convex set K ≔ C ∩ D can be characterized by the best approximation to a perturbation x − l of x from the set C for some l in a certain cone in X . The set C is a closed convex subset of X and D ≔ { x ∈ X : g j ( x ) ≤ 0 , ∀ j = 1 , 2 , … , m } , where the functions g j : X ⟶ R ( j = 1 , 2 , … , m ) are continuously Frechet differentiable that are not necessarily convex. We show under suitable conditions that this “perturbation property” is characterized by the strong conical hull intersection property of C and D at the point x 0 ∈ K . We prove this by first establishing a dual cone characterization of a nearly convex set. Our result shows that the convex geometry of K is critical for the characterization rather than the representation of D by convex inequalities, which is commonly assumed for the problems of best approximation from a convex set. In the special case where the set D is convex, we show that the Lagrange multiplier characterization of best approximation holds under the standard Slater’s constraint qualification together with a non-degeneracy condition. The lack of representation of D by convex inequalities is supplemented by the non-degeneracy condition, but the characterization, even in this special case, allows applications to problems with quasi-convex functions g j , j = 1 , 2 , … , m , as they guarantee the convexity of D . Simple numerical examples illustrate the nature of our assumptions.

In this paper we prove the existence of periodic solutions for gradient systems in finite and infinite dimensional spaces. The techniques of the proofs are based on the application of a global inverse functions theorem, the Schaefer fixed point theorem and the Faedou–Galerkin method. We apply our results in order to solve nonlinear reaction–diffusion equations with Dirichlet and Neumann boundary conditions.

In this paper we present a systematic way to describe exceptional Jacobi polynomials via two partitions. We give the construction of these polynomials and restate the known aspects of these polynomials in terms of their partitions. The aim is to show that the use of partitions is an elegant way to label these polynomials. Moreover, we prove asymptotic results according to the regular and exceptional zeros of these polynomials., 40 pages, 1 figure

A general theory of partial balayage on Riemannian manifolds is developed, with emphasis on compact manifolds. Partial balayage is an operation of sweeping measures, or charge distributions, to a prescribed density, and it is closely related to (construction of) quadrature domains for subharmonic functions, growth processes such as Laplacian growth and to weighted equilibrium distributions. Several examples are given in the paper, as well as some specific results. For instance, it is proved that, in two dimensions, harmonic and geodesic balls are the same if and only if the Gaussian curvature of the manifold is constant., Comment: 68 pages