On the convergence of Newton-type methods under mild differentiability conditions

We introduce the new idea of recurrent functions to provide a new semilocal convergence analysis for Newton-type methods, under mild differentiability conditions. It turns out that our sufficient convergence conditions are weaker, and the error bounds are tighter than in earlier studies in some interesting cases (Chen, Ann Inst Stat Math 42:387–401, 1990; Chen, Numer Funct Anal Optim 10:37–48, 1989; Cianciaruso, Numer Funct Anal Optim 24:713–723, 2003; Cianciaruso, Nonlinear Funct Anal Appl 2009; Dennis 1971; Deuflhard 2004; Deuflhard, SIAM J Numer Anal 16:1–10, 1979; Gutiérrez, J Comput Appl Math 79:131–145, 1997; Hernández, J Optim Theory Appl 109:631–648, 2001; Hernández, J Comput Appl Math 115:245–254, 2000; Huang, J Comput Appl Math 47:211–217, 1993; Kantorovich 1982; Miel, Numer Math 33:391–396, 1979; Miel, Math Comput 34:185–202, 1980; Moret, Computing 33:65–73, 1984; Potra, Libertas Mathematica 5:71–84, 1985; Rheinboldt, SIAM J Numer Anal 5:42–63, 1968; Yamamoto, Numer Math 51: 545–557, 1987; Zabrejko, Numer Funct Anal Optim 9:671–684, 1987; Zinc̆ko 1963). Applications and numerical examples, involving a nonlinear integral equation of Chandrasekhar-type, and a differential equation are also provided in this study.     

Smoothing metrics on closed Riemannian manifolds through the Ricci flow

Under the assumption of the uniform local Sobolev inequality, it is proved that Riemannian metrics with an absolute Ricci curvature bound and a small Riemannian curvature integral bound can be smoothed to having a sectional curvature bound. This partly extends previous a priori estimates of Li (J Geom Anal 17:495–511, 2007; Adv Math 223:1924–1957, 2010).     

Several modifications of DPR estimator of the tail index

In this paper, we continue the investigation of an estimator proposed in [Yu. Davydov, V. Paulauskas, and A. Račkauskas, More on p-stable convex sets in Banach spaces, J. Theor. Probab., 13:39–64, 2000] and [V. Paulauskas, A new estimator for tail index, Acta Appl. Math., 79:55–67, 2003] and considered in [V. Paulauskas and M. Vaičiulis, Once more on comparison of tail index estimators, preprint, 2010]. We propose a class of modifications of the so-called DPR estimator and demonstrate that these modifications can have better asymptotic properties than the original DPR estimator.     

Parametric approximation of isotropic and anisotropic elastic flow for closed and open curves

Deckelnick and Dziuk (Math. Comput. 78(266):645–671, 2009) proved a stability bound for a continuous-in-time semidiscrete parametric finite element approximation of the elastic flow of closed curves in \mathbbR^d, d 3 2{\mathbb{R}^d, d\geq2} . We extend these ideas in considering an alternative finite element approximation of the same flow that retains some of the features of the formulations in Barrett et al. (J Comput Phys 222(1): 441–462, 2007; SIAM J Sci Comput 31(1):225–253, 2008; IMA J Numer Anal 30(1):4–60, 2010), in particular an equidistribution mesh property. For this new approximation, we obtain also a stability bound for a continuous-in-time semidiscrete scheme. Apart from the isotropic situation, we also consider the case of an anisotropic elastic energy. In addition to the evolution of closed curves, we also consider the isotropic and anisotropic elastic flow of a single open curve in the plane and in higher codimension that satisfies various boundary conditions.     

On a class of secant-like methods for solving nonlinear equations

We provide a semilocal convergence analysis for a certain class of secant-like methods considered also in Argyros (J Math Anal Appl 298:374–397, 2004, 2007), Potra (Libertas Mathematica 5:71–84, 1985), in order to approximate a locally unique solution of an equation in a Banach space. Using a combination of Lipschitz and center-Lipschitz conditions for the computation of the upper bounds on the inverses of the linear operators involved, instead of only Lipschitz conditions (Potra, Libertas Mathematica 5:71–84, 1985), we provide an analysis with the following advantages over the work in Potra (Libertas Mathematica 5:71–84, 1985) which improved the works in Bosarge and Falb (J Optim Theory Appl 4:156–166, 1969, Numer Math 14:264–286, 1970), Dennis (SIAM J Numer Anal 6(3):493–507, 1969, 1971), Kornstaedt (1975), Larsonen (Ann Acad Sci Fenn, A 450:1–10, 1969), Potra (L’Analyse Numérique et la Théorie de l’Approximation 8(2):203–214, 1979, Aplikace Mathematiky 26:111–120, 1981, 1982, Libertas Mathematica 5:71–84, 1985), Potra and Pták (Math Scand 46:236–250, 1980, Numer Func Anal Optim 2(1):107–120, 1980), Schmidt (Period Math Hung 9(3):241–247, 1978), Schmidt and Schwetlick (Computing 3:215–226, 1968), Traub (1964), Wolfe (Numer Math 31:153–174, 1978): larger convergence domain; weaker sufficient convergence conditions, finer error bounds on the distances involved, and a more precise information on the location of the solution. Numerical examples further validating the results are also provided.     

Error estimates for abstract linear–quadratic optimal control problems using proper orthogonal decomposition

In this paper we investigate POD discretizations of abstract linear–quadratic optimal control problems with control constraints. We apply the discrete technique developed by Hinze (Comput. Optim. Appl. 30:45–61, 2005) and prove error estimates for the corresponding discrete controls, where we combine error estimates for the state and the adjoint system from Kunisch and Volkwein (Numer. Math. 90:117–148, 2001; SIAM J. Numer. Anal. 40:492–515, 2002). Finally, we present numerical examples that illustrate the theoretical results.     

A Fréchet derivative-free cubically convergent method for set-valued maps

We introduce a new iterative method in order to approximate a locally unique solution of variational inclusions in Banach spaces. The method uses only divided differences operators of order one. An existence–convergence theorem and a radius of convergence are given under some conditions on divided difference operator and Lipschitz-like continuity property of set-valued mappings. Our method extends the recent work related to the resolution of nonlinear equation in Argyros (J Math Anal Appl 332:97–108, 2007) and has the following advantages: faster convergence to the solution than all the previous known ones in Argyros and Hilout (Appl Math Comput, 2008 in press), Hilout (J Math Anal Appl 339:53–761, 2008, Positivity 10:673–700, 2006), and we do not need to evaluate any Fréchet derivative. We provide also an improvement of the ratio of our algorithm under some center-conditions and less computational cost. Numerical examples are also provided.      

A Note on the Subcritical Two Dimensional Keller-Segel System

The existence of solution for the 2D-Keller-Segel system in the subcritical case, i.e. when the initial mass is less than 8π, is reproved. Instead of using the entropy in the free energy and free energy dissipation, which was used in the proofs (Blanchet et al. in SIAM J. Numer. Anal. 46:691–721, 2008; Electron. J. Differ. Equ. Conf. 44:32, 2006 (electronic)), the potential energy term is fully utilized by adapting Delort’s theory on 2D incompressible Euler equation (Delort in J. Am. Math. Soc. 4:553–386, 1991).     

Fully computable robust a posteriori error bounds for singularly perturbed reaction–diffusion problems

A procedure for the construction of robust, upper bounds for the error in the finite element approximation of singularly perturbed reaction–diffusion problems was presented in Ainsworth and Babuška (SIAM J Numer Anal 36(2):331–353, 1999) which entailed the solution of an infinite dimensional local boundary value problem. It is not possible to solve this problem exactly and this fact was recognised in the above work where it was indicated that the limitation would be addressed in a subsequent article. We view the present work as fulfilling that promise and as completing the investigation begun in Ainsworth and Babuška (SIAM J Numer Anal 36(2):331–353, 1999) by removing the obligation to solve a local problem exactly. The resulting new estimator is indeed fully computable and the first to provide fully computable, robust upper bounds in the setting of singularly perturbed problems discretised by the finite element method.     

Asymptotic normality of the Parzen–Rosenblatt density estimator for strongly mixing random fields

We prove the asymptotic normality of the kernel density estimator (introduced by Rosenblatt, Proc Natl Acad Sci USA 42:43–47, 1956 and Parzen, Ann Math Stat 33:1965–1976, 1962) in the context of stationary strongly mixing random fields. Our approach is based on the Lindeberg’s method rather than on Bernstein’s small-block-large-block technique and coupling arguments widely used in previous works on nonparametric estimation for spatial processes. Our method allows us to consider only minimal conditions on the bandwidth parameter and provides a simple criterion on the strong mixing coefficients which do not depend on the bandwidth.     

One-dimensional chemotaxis kinetic model

In this paper, we study a variation of the equations of a chemotaxis kinetic model and investigate it in one dimension. In fact, we use fractional diffusion for the chemoattractant in the Othmar–Dunbar–Alt system (Othmer in J Math Biol 26(3):263–298, 1988). This version was exhibited in Calvez in Amer Math Soc, pp 45–62, 2007 for the macroscopic well-known Keller–Segel model in all space dimensions. These two macroscopic and kinetic models are related as mentioned in Bournaveas, Ann Inst H Poincaré Anal Non Linéaire, 26(5):1871–1895, 2009, Chalub, Math Models Methods Appl Sci, 16(7 suppl):1173–1197, 2006, Chalub, Monatsh Math, 142(1–2):123–141, 2004, Chalub, Port Math (NS), 63(2):227–250, 2006. The model we study here behaves in a similar way to the original model in two dimensions with the spherical symmetry assumption on the initial data which is described in Bournaveas, Ann Inst H Poincaré Anal Non Linéaire, 26(5):1871–1895, 2009. We prove the existence and uniqueness of solutions for this model, as well as a convergence result for a family of numerical schemes. The advantage of this model is that numerical simulations can be easily done especially to track the blow-up phenomenon.     

Optimal Meshes for Finite Elements of Arbitrary Order

Given a function f defined on a bounded domain Ω⊂ℝ² and a number N>0, we study the properties of the triangulation T_N\mathcal{T}_{N} that minimizes the distance between f and its interpolation on the associated finite element space, over all triangulations of at most N elements. The error is studied in the norm X=L ^p for 1≤p≤∞, and we consider Lagrange finite elements of arbitrary polynomial degree m−1. We establish sharp asymptotic error estimates as N→+∞ when the optimal anisotropic triangulation is used, recovering the results on piecewise linear interpolation (Babenko et al. in East J. Approx. 12(1), 71–101, 2006; Babenko, submitted; Chen et al. in Math. Comput. 76, 179–204, 2007) and improving the results on higher degree interpolation (Cao in SIAM J. Numer. Anal. 45(6), 2368–2391, 2007, SIAM J. Sci. Comput. 29, 756–781, 2007, Math. Comput. 77, 265–286, 2008). These estimates involve invariant polynomials applied to the m-th order derivatives of f. In addition, our analysis also provides practical strategies for designing meshes such that the interpolation error satisfies the optimal estimate up to a fixed multiplicative constant. We partially extend our results to higher dimensions for finite elements on simplicial partitions of a domain Ω⊂ℝ ^d .     

Hedonic price function estimation in economics and marketing: revisiting Lancaster’s issue of “noncombinable” goods

Following Lancaster’s (J. Political Econ. 74(1):132–157, 1966; Variety, equity and efficiency, 1979) interpretation of his characteristics approach to consumer theory, this contribution focuses on theoretical and empirical arguments questioning the smoothness of traditional hedonic price estimation techniques. Lancaster argued strongly against “combinability”, i.e., that any efficient combination of characteristics is feasible and sensible. We therefore explicitly test the impact of convexity using a set of recent non-parametric estimators. The test is carried out on a sample of 114 digital cameras whose price evolution is followed over 6 months. The hypothesis of convexity is rejected using the Li (Econ. Rev. 15(3):261–274, 1996) test. The conclusions point out implications for economics and marketing.     

Asymptotic normality of location invariant heavy tail index estimator

Motivated by Fraga Alves (Extremes 4:199–217, 2001)’s work, a new class of location invariant Hill-type estimators for the tail index of a heavy tailed distribution is proposed in the paper. Its asymptotic behavior is derived, and the optimal choice of the sample fraction is discussed by mean squared error. Asymptotic comparisons and simulation studies are presented to show that the new estimator performs well compared to the known ones.     

On characters and formal degrees of discrete series of affine Hecke algebras of classical types

We address two fundamental questions in the representation theory of affine Hecke algebras of classical types. One is an inductive algorithm to compute characters of tempered modules, and the other is the determination of the constants in the formal degrees of discrete series (in the form conjectured by Reeder (J. Reine Angew. Math. 520:37–93, 2000)). The former is completely different from the Lusztig-Shoji algorithm (Shoji in Invent. Math. 74:239–267, 1983; Lusztig in Ann. Math. 131:355–408, 1990), and it is more effective in a number of cases. The main idea in our proof is to introduce a new family of representations which behave like tempered modules, but for which it is easier to analyze the effect of parameter specializations. Our proof also requires a comparison of the C ^∗-theoretic results of Opdam, Delorme, Slooten, Solleveld (J. Inst. Math. Jussieu 3:531–648, 2004; ; Int. Math. Res. Not., 2008; Adv. Math. 220:1549–1601, 2009; Acta Math. 205:105–187, 2010), and the geometric construction from Kato (Duke Math. J. 148:305–371, 2009; Am. J. Math. 133:518–553, 2011), Ciubotaru and Kato (Adv. Math. 226:1538–1590, 2011).     

Small Support Spline Riesz Wavelets in Low Dimensions

In Han and Shen (SIAM J. Math. Anal. 38:530–556, 2006), a family of univariate short support Riesz wavelets was constructed from uniform B-splines. A bivariate spline Riesz wavelet basis from the Loop scheme was derived in Han and Shen (J. Fourier Anal. Appl. 11:615–637, 2005). Motivated by these two papers, we develop in this article a general theory and a construction method to derive small support Riesz wavelets in low dimensions from refinable functions. In particular, we obtain small support spline Riesz wavelets from bivariate and trivariate box splines. Small support Riesz wavelets are desirable for developing efficient algorithms in various applications. For example, the short support Riesz wavelets from Han and Shen (SIAM J. Math. Anal. 38:530–556, 2006) were used in a surface fitting algorithm of Johnson et al. (J. Approx. Theory 159:197–223, 2009), and the Riesz wavelet basis from the Loop scheme was used in a very efficient geometric mesh compression algorithm in Khodakovsky et al. (Proceedings of SIGGRAPH, 2000).     

Asymptotic expansions for the conditional sojourn time distribution in the <Emphasis Type="Italic">M</Emphasis>/<Emphasis Type="Italic">M</Emphasis>/1-PS queue

We consider the M/M/1 queue with processor sharing. We study the conditional sojourn time distribution, conditioned on the customer’s service requirement, in various asymptotic limits. These include large time and/or large service request, and heavy traffic, where the arrival rate is only slightly less than the service rate. The asymptotic formulas relate to, and extend, some results of Morrison (SIAM J. Appl. Math. 45:152–167, [1985]) and Flatto (Ann. Appl. Probab. 7:382–409, [1997]). This work was partly supported by NSF grant DMS 05-03745.     

Asymptotic Approximations between the Hahn-Type Polynomials and Hermite, Laguerre and Charlier Polynomials

It has been shown in Ferreira et al. (Adv. Appl. Math 31:61–85, [2003]), López and Temme (Methods Appl. Anal. 6:131–196, [1999]; J. Cpmput. Appl. Math. 133:623–633, [2001]) that the three lower levels of the Askey table of hypergeometric orthogonal polynomials are connected by means of asymptotic expansions. In this paper we continue with that investigation and establish asymptotic connections between the fourth level and the two lower levels: we derive twelve asymptotic expansions of the Hahn, dual Hahn, continuous Hahn and continuous dual Hahn polynomials in terms of Hermite, Charlier and Laguerre polynomials. From these expansions, several limits between polynomials are derived. Some numerical experiments give an idea about the accuracy of the approximations and, in particular, about the accuracy in the approximation of the zeros of the Hahn, dual Hahn, continuous Hahn and continuous dual Hahn polynomials in terms of the zeros of the Hermite, Charlier and Laguerre polynomials.      

A Bilinear Pseudodifferential Calculus

In this paper we construct a new class of bilinear pseudodifferential operators which contains both the Coifman-Meyer class as well as the non-translation invariant class closely related both to the bilinear Hilbert transform and previously studied in Bényi et al. (J. Geom. Anal. 16(3):431–453, 2006), Bényi et al. (J. Anal. Math., 2009), Bernicot (Anal. PDE 1:1–27, 2008) as well as the bilinear Marcinkiewicz class studied in Grafakos and Kalton (Stud. Math. 146(2):115–156, 2001). We prove boundedness on Sobolev spaces for these operators as well as establish a symbolic calculus that exhibits the nice behavior of our new class under transposition and composition with linear operators.     

Stability Analysis for Minty Mixed Variational Inequality in Reflexive Banach Spaces

This paper is devoted to the stability analysis for a class of Minty mixed variational inequalities in reflexive Banach spaces, when both the mapping and the constraint set are perturbed. Several equivalent characterizations are given for the Minty mixed variational inequality to have nonempty and bounded solution set. A stability result is presented for the Minty mixed variational inequality with Φ-pseudomonotone mapping in reflexive Banach space, when both the mapping and the constraint set are perturbed by different parameters. As an application, a stability result for a generalized mixed variational inequality is also obtained. The results presented in this paper generalize and extend some known results in Fan and Zhong (Nonlinear Anal., Theory Methods Appl. 69:2566–2574, 2008) and He (J. Math. Anal. Appl. 330:352–363, 2007).