Error estimates for the numerical approximation of Neumann control problems

We continue the discussion of error estimates for the numerical analysis of Neumann boundary control problems we started in Casas et al. (Comput. Optim. Appl. 31:193–219, 2005). In that paper piecewise constant functions were used to approximate the control and a convergence of order O(h) was obtained. Here, we use continuous piecewise linear functions to discretize the control and obtain the rates of convergence in L ²(Γ). Error estimates in the uniform norm are also obtained. We also discuss the approach suggested by Hinze (Comput. Optim. Appl. 30:45–61, 2005) as well as the improvement of the error estimates by making an extra assumption over the set of points corresponding to the active control constraints. Finally, numerical evidence of our estimates is provided. The authors were supported by Ministerio de Ciencia y Tecnología (Spain).     

Error analysis of Jacobi derivative estimators for noisy signals

Recent algebraic parametric estimation techniques (see Fliess and Sira-Ramírez, ESAIM Control Optim Calc Variat 9:151–168, 2003, 2008) led to point-wise derivative estimates by using only the iterated integral of a noisy observation signal (see Mboup et al. 2007, Numer Algorithms 50(4):439–467, 2009). In this paper, we extend such differentiation methods by providing a larger choice of parameters in these integrals: they can be reals. For this, the extension is done via a truncated Jacobi orthogonal series expansion. Then, the noise error contribution of these derivative estimations is investigated: after proving the existence of such integral with a stochastic process noise, their statistical properties (mean value, variance and covariance) are analyzed. In particular, the following important results are obtained:

Non-Walker Self-Dual Neutral Einstein Four-Manifolds of Petrov Type III

The local structure of the manifolds named in the title is described. Although curvature homogeneous, they are not, in general, locally homogeneous. Not all of them are Ricci-flat, which answers an existence question about type III Jordan-Osserman metrics, raised by Díaz-Ramos, García-Río and Vázquez-Lorenzo (J. Geom. Anal. 16, 39–52, 2006). Work begun during the author’s visit to the University of Santiago de Compostela, supported by Grant MTM2006-01432 (Spain).     

On Constrained Nonlinear Hermite Subdivision

We determine shape-preserving regions and we describe a general setting to generate shape-preserving families for the 2-points Hermite subdivision scheme introduced by Merrien (Numer. Algorithms 2:187–200, [1992]). This general construction includes the shape-preserving families presented in Merrien and Sablonníere (Constr. Approx. 19:279–298, [2003]) and Pelosi and Sablonníere (C ¹ GP Hermite Interpolants Generated by a Subdivision Scheme, Prépublication IRMAR 06–23, Rennes, [2006]). New special families are presented as particular examples. Nonstationary and nonuniform versions of such schemes, which produce smoother limits, are discussed.      

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.     

Strong lower energy estimates for nonlinearly damped Timoshenko beams and Petrowsky equations

The purpose of this paper is to establish strong lower energy estimates for strong solutions of nonlinearly damped Timoshenko beams, Petrowsky equations in two and three dimensions and wave-like equations for bounded one-dimensional domains or annulus domains in two or three dimensions. We also establish weak lower velocity estimates for strong solutions of the nonlinearly damped Petrowsky equation in two and three dimensions. The feedbacks in consideration have arbitrary growth close to the origin. These results improve the strong lower energy decay rates obtained in our previous papers (Alabau-Boussouira in J Differ Equ 249:1145–1178, 2010; J Differ Equ 248:1473–1517, 2010) for strong solutions of the nonlinearly locally damped wave equation and extend to systems and to Petrowsky equation the method of Alabau-Boussouira (J Differ Equ 249:1145–1178, 2010; J Differ Equ 248:1473–1517, 2010). These results are the first ones for Timoshenko beams and Petrowsky equations.     

A derivative free iterative method for solving least squares problems

A derivative free iterative method for approximating a solution of nonlinear least squares problems is studied first in Shakhno and Gnatyshyn (Appl Math Comput 161:253–264, 2005). The radius of convergence is determined as well as usable error estimates. We show that this method is faster than its Secant analogue examined in Shakhno and Gnatyshyn (Appl Math Comput 161:253–264, 2005). Numerical example is also provided in this paper.     

Polynomial Birth–Death Distribution Approximation in the Wasserstein Distance

The polynomial birth–death distribution (abbreviated, PBD) on ℐ={0,1,2,…} or ℐ={0,1,2,…,m} for some finite m introduced in Brown and Xia (Ann. Probab. 29:1373–1403, 2001) is the equilibrium distribution of the birth–death process with birth rates {α _i } and death rates {β _i }, where α _i ≥0 and β _i ≥0 are polynomial functions of i∈ℐ. The family includes Poisson, negative binomial, binomial, and hypergeometric distributions. In this paper, we give probabilistic proofs of various Stein’s factors for the PBD approximation with α _i =a and β _i =i+bi(i−1) in terms of the Wasserstein distance. The paper complements the work of Brown and Xia (Ann. Probab. 29:1373–1403, 2001) and generalizes the work of Barbour and Xia (Bernoulli 12:943–954, 2006) where Poisson approximation (b=0) in the Wasserstein distance is investigated. As an application, we establish an upper bound for the Wasserstein distance between the PBD and Poisson binomial distribution and show that the PBD approximation to the Poisson binomial distribution is much more precise than the approximation by the Poisson or shifted Poisson distributions.      

Semilinear Mixed Problems on Hilbert Complexes and Their Numerical Approximation

Arnold, Falk, and Winther recently showed (Bull. Am. Math. Soc. 47:281–354, 2010) that linear, mixed variational problems, and their numerical approximation by mixed finite element methods, can be studied using the powerful, abstract language of Hilbert complexes. In another recent article (arXiv:), we extended the Arnold–Falk–Winther framework by analyzing variational crimes (à la Strang) on Hilbert complexes. In particular, this gave a treatment of finite element exterior calculus on manifolds, generalizing techniques from surface finite element methods and recovering earlier a priori estimates for the Laplace–Beltrami operator on 2- and 3-surfaces, due to Dziuk (Lecture Notes in Math., vol. 1357:142–155, 1988) and later Demlow (SIAM J. Numer. Anal. 47:805–827, 2009), as special cases. In the present article, we extend the Hilbert complex framework in a second distinct direction: to the study of semilinear mixed problems. We do this, first, by introducing an operator-theoretic reformulation of the linear mixed problem, so that the semilinear problem can be expressed as an abstract Hammerstein equation. This allows us to obtain, for semilinear problems, a priori solution estimates and error estimates that reduce to the Arnold–Falk–Winther results in the linear case. We also consider the impact of variational crimes, extending the results of our previous article to these semilinear problems. As an immediate application, this new framework allows for mixed finite element methods to be applied to semilinear problems on surfaces.     

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.     

A priori error estimates for space–time finite element discretization of semilinear parabolic optimal control problems

In this paper, a priori error estimates for space–time finite element discretizations of optimal control problems governed by semilinear parabolic PDEs and subject to pointwise control constraints are derived. We extend the approach from Meidner and Vexler (SIAM Control Optim 47(3):1150–1177, 2008; SIAM Control Optim 47(3):1301–1329, 2008) where linear-quadratic problems have been considered, discretizing the state equation by usual conforming finite elements in space and a discontinuous Galerkin method in time. Error estimates for controls discretized by piecewise constant functions in time and cellwise constant functions in space are derived in detail and we explain how error estimate for further discretization approaches, e.g., cellwise linear discretization in space, the postprocessing approach from Meyer and R?sch (SIAM J Control Optim 43:970–985, 2004), and the variationally discrete approach from Hinze (J Comput Optim Appl 30:45–63, 2005) can be obtained. In addition, we derive an estimate for a setting with finitely many time-dependent controls.     

An extended doubly-adaptive quadrature method based on the combination of the Ninomiya and the FLR schemes

An improvement is made to an automatic quadrature due to Ninomiya (J. Inf. Process. 3:162–170, 1980) of adaptive type based on the Newton–Cotes rule by incorporating a doubly-adaptive algorithm due to Favati, Lotti and Romani (ACM Trans. Math. Softw. 17:207–217, 1991; ACM Trans. Math. Softw. 17:218–232, 1991). We compare the present method in performance with some others by using various test problems including Kahaner’s ones (Computation of numerical quadrature formulas. In: Rice, J.R. (ed.) Mathematical Software, 229–259. Academic, Orlando, FL, 1971).      

An FIO calculus for marine seismic imaging, II: Sobolev estimates

We establish sharp L ²-Sobolev estimates for classes of pseudodifferential operators with singular symbols [Guillemin and Uhlmann (Duke Math J 48:251–267, 1981), Melrose and Uhlmann (Commun Pure Appl Math 32:483–519, 1979)] whose non-pseudodifferential (Fourier integral operator) parts exhibit two-sided fold singularities. The operators considered include both singular integral operators along curves in \mathbb R²{\mathbb R^2} with simple inflection points and normal operators arising in linearized seismic imaging in the presence of fold caustics [Felea (Comm PDE 30:1717–1740, 2005), Felea and Greenleaf (Comm PDE 33:45–77, 2008), Nolan (SIAM J Appl Math 61:659–672, 2000)].     

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).     

Generalized Qualification and Qualification Levels for Spectral Regularization Methods

The concept of qualification for spectral regularization methods (SRM) for inverse ill-posed problems is strongly associated to the optimal order of convergence of the regularization error (Engl et al. in Regularization of inverse problems. Mathematics and its applications, vol. 375, Kluwer Academic, Dordrecht, 1996; Mathé in SIAM J. Numer. Anal. 42(3):968–973, 2004; Mathé and Pereverzev in Inverse Probl. 19(3):789–803, 2003; Vainikko in USSR Comput. Math. Math. Phys. 22(3): 1–19, 1982). In this article, the definition of qualification is extended and three different levels are introduced: weak, strong and optimal. It is shown that the weak qualification extends the definition introduced by Mathé and Pereverzev (Inverse Probl. 19(3):789–803, 2003), mainly in the sense that the functions associated with orders of convergence and source sets need not be the same. It is shown that certain methods possessing infinite classical qualification (e.g. truncated singular value decomposition (TSVD), Landweber’s method and Showalter’s method) also have generalized qualification leading to an optimal order of convergence of the regularization error. Sufficient conditions for a SRM to have weak qualification are provided and necessary and sufficient conditions for a given order of convergence to be strong or optimal qualification are found. Examples of all three qualification levels are provided and the relationships between them as well as with the classical concept of qualification and the qualification introduced in Mathé and Pereverzev (Inverse Probl. 19(3):789–803, 2003) are shown. In particular, SRMs having extended qualification in each one of the three levels and having zero or infinite classical qualification are presented. Finally, several implications of this theory in the context of orders of convergence, converse results and maximal source sets for inverse ill-posed problems, are shown. This work was supported by DARPA/SPO, NASA LaRC and the National Institute of Aerospace under Grant VT-03-1, 2535, by AFOSR Grants F49620-03-1-0243 and FA9550-07-1-0273, by Consejo Nacional de Investigaciones Científicas y Técnicas, CONICET, and by Universidad Nacional del Litoral, U.N.L., Argentina, through Project CAI+D 2006, P.E. 236.     

Determining <Emphasis Type="Italic">p</Emphasis>-values for systems cointegration tests with a prior adjustment for deterministic terms

Following Doornik (J Econ Surv 12:573–593, 1998) I present a procedure to approximate the asymptotic distributions of systems cointegration tests with a prior adjustment for deterministic terms suggested by Lütkepohl (Econometrica 72:647–662, 2004), Saikkonen and Lütkepohl (Econometric Theory 16:373–406, 2000a, J Business Econ Stat 18:451–464, 2000b, Time Series Anal 21:435–456, 2000c) and Saikkonen and Luukkonen (J Econ 81:93–126, 1997). These tests rely upon different assumptions as to the inclusion of deterministic components such as a constant, a linear trend or a level shift. The asymptotic distributions, which are functions of Brownian motions, are approximated by Gamma distributions. Only estimates of the mean and variance of the asymptotic test distributions are needed to fit the Gamma distributions. Such estimates are obtained from response surfaces. The required coefficients to compute the asymptotic moments are presented in this paper. Via the fitted Gamma distributions one can, then, easily derive p-values or arbitrary percentiles.     

On Implicit Multifunction Theorems

We give implicit multifunction results generalizing to multifunctions the Robinson’s implicit function theorem (Robinson, Math Oper Res 16(2):292–309, 1991). To this end, we use parametric error bounds estimates for a suitable function refining the one given in Azé and Corvellec (ESAIM Control Optim Calc Var 10:409–425, 2004). Sharp approximations of the implicit multifunctions are given extending the results of Nachi and Penot (Control Cybernet 35:871–901, 2005). Dedicated to Boris Mordukhovich in honour of his 60th birthday.     

The robustness of mean and variance approximations in risk analysis

This paper examines further the problem of estimating the mean and variance of a continuous random variable from estimates of three points within the distribution, typically the median or mode and two extreme fractiles. The problem arises most commonly in PERT and risk analysis where it can usually be assumed that the distribution in question is bell-shaped and positively skewed, often typified by a Beta distribution. Over the years, a number of alternative approximations have been proposed, usually as modifications to the original PERT formulae. The accuracy of a number of these approximations is investigated based not only on a Beta distribution, but also for three other commonly used bell-shaped, positively skewed distributions, namely the Gamma, Lognormal and F distributions. It is shown that a balanced weighted average of the median and the 4% fractiles provides a consistent estimator of the distribution mean across all four distributions. Furthermore, reasonably accurate estimates of the variance can also be obtained by treating the three fractiles as defining an equivalent discrete distribution with the same probability weight as in the formula for the mean.     

On Rump’s characterization of <Emphasis Type="Italic">P</Emphasis>-matrices

The necessary and sufficient P-matrix condition by Rump (Linear Algebra Appl 363:237–250, 2003) is simplified by showing that one of its assumptions can be deleted without affecting validity of the result.     

Components of the Fundamental Category II

In this article we carry on the study of the fundamental category (Goubault and Raussen, Dihomotopy as a tool in state space analysis. In: Rajsbaum, S. (ed.) LATIN 2002: Theoretical Informatics. Lecture Notes in Computer Science, vol. 2286, Cancun, Mexico, pp. 16–37, Springer, Berlin Heidelberg New York, 2002; Goubault, Homology, Homotopy Appl., 5(2): 95–136, 2003) of a partially ordered topological space (Nachbin, Topology and Order, Van Nostrand, Princeton, 1965; Johnstone, Stone Spaces, Cambridge University Press, Cambridge, MA, 1982), as arising in e.g. concurrency theory (Fajstrup et al., Theor. Comp. Sci. 357: 241–278, 2006), initiated in (Fajstrup et al., APCS, 12(1): 81–108, 2004). The “algebra” of dipaths modulo dihomotopy (the fundamental category) of such a po-space is essentially finite in a number of situations. We give new definitions of the component category that are more tractable than the one of Fajstrup et al. (APCS, 12(1): 81–108, 2004), as well as give definitions of future and past component categories, related to the past and future models of Grandis (Theory Appl. Categ., 15(4): 95–146, 2005). The component category is defined as a category of fractions, but it can be shown to be equivalent to a quotient category, much easier to portray. A van Kampen theorem is known to be available on fundamental categories (Grandis, Cahiers Topologie Géom. Différentielle Catég., 44: 281–316, 2003; Goubault, Homology, Homotopy Appl., 5(2): 95–136, 2003), we show in this paper a similar theorem for component categories (conjectured in Fajstrup et al. (APCS, 12(1): 81–108, 2004). This proves useful for inductively computing the component category in some circumstances, for instance, in the case of simple PV mutual exclusion models (Goubault and Haucourt, A practical application of geometric semantics to static analysis of concurrent programs. In: Abadi, M., de Alfaro, L. (eds.) CONCUR 2005 – Concurrency Theory: 16th International Conference, San Francisco, USA, August 23–26. Lecture Notes in Computer Science, vol. 3653, pp. 503–517, Springer, Berlin Heidelberg New York, 2005), corresponding to partially ordered subspaces of R ⁿ minus isothetic hyperrectangles. In this last case again, we conjecture (and give some hints) that component categories enjoy some nice adjunction relations directly with the fundamental category.