Error bounds and parameter choice strategies for simplified regularization in Hilbert scales

Simplified regularization in the setting of Hilbert scales has been considered for obtaining stable approximate solutions for ill-posed operator equations. The derived error estimates using an a posteriori as well as an a priori parameter choice strategy are shown to be of optimal order with respect to certain natural assumptions on the ill-posedness of the equation.The work of M. Thamban Nair is partially supported by IC&SR, I.I.T., Madras     

Second order parallel algorithms for piecewise smooth displacement kernels

Second order parallel algorithms for Fredholm integral equations with piecewise smooth displacement kernels are derived. One is based on a difference scheme of Runge-Kutta type for an unusual partial differential equations for continuous functions of two variables. The other is based on the trapezoidal quadrature rule applied to a modified integral equations. It is found that the Runge-Kutta type algorithm exhibits certain advantages.The work of these authors was supported in part by the NSF Grant DMS-9007030The work of this author was supported in part by a grant from the National Science and Engineering Research Council of Canada     

Relations between convolution type operators on intervals and on the half-line

This paper is devoted to the question to obtain (algebraic and topologic) equivalence (after extension) relations between convolution type operators on unions of intervals and convolution type operators on the half-line. These operators are supposed to act between Bessel potential spaces,H ^s,p , which are the appropriate spaces in several applications. The present approach is based upon special properties of convenient projectors, decompositions and extension operators and the construction of certain homeomorphisms between the kernels of the projectors. The main advantage of the method is that it provides explicit operator matrix identities between the mentioned operators where the relations are constructed only by bounded invertible operators. So they are stronger than the (algebraic) Kuijper-Spitkovsky relation and the Bastos-dos Santos-Duduchava relation with respect to the transfer of properties on the prize that the relations depend on the orders of the spaces and hold only for non-critical orders:S – 1/p . For instance, (generalized) inverses of the operators are explicitly represented in terms of operator matrix factorization. Some applications are presented.This research was supported by Junta Nacional de Investigação Científica e Tecnológica (Portugal) and the Bundesminister für Forschung und Technologie (Germany) within the projectSingular Operators-new features and applications, and by a PRAXIS XXI project under the titleFactorization of Operators and Applications to Mathematical Physics.     

Circulant integral operators as preconditioners for Wiener-Hopf equations

In this paper, we study the solutions of finite-section Wiener-Hopf equations by the preconditioned conjugate gradient method. Our main aim is to give an easy and general scheme of constructing good circulant integral operators as preconditioners for such equations. The circulant integral operators are constructed from sequences of conjugate symmetric functions {C }. Letk(t) denote the kernel function of the Wiener-Hopf equation and be its Fourier transform. We prove that for sufficiently large if {C } is uniformly bounded on the real lineR and the convolution product of the Fourier transform ofC with uniformly onR, then the circulant preconditioned Wiener-Hopf operator will have a clustered spectrum. It follows that the conjugate gradient method, when applied to solving the preconditioned operator equation, converges superlinearly. Several circulant integral operators possessing the clustering and fast convergence properties are constructed explicitly. Numerical examples are also given to demonstrate the performance of different circulant integral operators as preconditioners for Wiener-Hopf operators.Research supported in part by HKRGC grant no. 221600070.     

Spline approximation methods for multidimensional periodic pseudodifferential equations

We investigate several numerical methods for solving the pseudodifferential equationAu=f on the n-dimensional torusT ⁿ . We examine collocation methods as well as Galerkin-Petrov methods using various periodical spline functions. The considered spline spaces are subordinated to a uniform rectangular or triangular grid. For given approximation method and invertible pseudodifferential operatorA we compute a numerical symbol ^C , resp. ^G , depending onA and on the approximation method. It turns out that the stability of the numerical method is equivalent to the ellipticity of the corresponding numerical symbol. The case of variable symbols is tackled by a local principle. Optimal error estimates are established.The second author has been supported by a grant of Deutsche Forschungsgemeinschaft under grant namber Ko 634/32-1.     

Blowing up behavior for a class of nonlinear VIEs connected with parabolic PDEs

We study the blowing-up behavior of solutions of a class of nonlinear integral equations of Volterra type that is connected with parabolic partial differential equations with concentrated nonlinearities. We present some analytic results and, in the case of the kernel of Abel-kind with power nonlinearity and fixed initial data, we give a numerical approximation by using one-point collocation methods.     

Schauder estimates for a class of second order elliptic operators on a cube

We consider a class of second order elliptic operators on a d-dimensional cube S_d. We prove that if the coefficients are of class C^k+δ(S_d), with k=0,1 and δ∈(0,1), then the corresponding elliptic problem admits a unique solution u belonging to C^k+2+δ(S_d) and satisfying non-standard boundary conditions involving only second order derivatives.     

Algorithms for inversion of diagonal plus semiseparable operator matrices

LetH be a Hilbert space andRHH be a bounded linear operator represented by an operator matrix which is a sum of a diagonal and of a semiseparable type of order one operator matrices. We consider three methods for solution of the operator equationRx=y. The obtained results yields new algorithms for solution of integral equations and for inversion of matrices.     

Galerkin method for Wiener-Hopf operators with piecewise continuous symbol

This paper is concerned with the applicability of maximum defect polynomial (Galerkin) spline approximation methods with graded meshes to Wiener-Hopf operators with matrix-valued piecewise continuous generating function defined on R. For this, an algebra of sequences is introduced, which contains the approximating sequences we are interested in. There is a direct relationship between the stability of the approximation method for a given operator and invertibility of the corresponding sequence in this algebra. Exploring this relationship, the methods of essentialization, localization and identification of the local algebras are used in order to derive stability criteria for the approximation sequences.Supported by grant Praxis XXI/BD/4501/94 from FCT.Partly supported by FCT/BMFT grant 423.     

Perturbed unstable variational inequalities with unbounded operators on approximately given sets

The stability and convergence of the solutions of perturbed and regularized variational inequality to the solutions of the primary (unstable a priori) variational inequality with proper monotone operator are investigated. All the objects of inequality: the operatorA, the right-hand partf and the set of constrains are to be perturbed. At the same time no assumptions of boundedness and smoothness of the operatorA are used. The connection between the parameters of perturbations, which guarantees strong convergence of approximate solutions, is established. It is proved that the existence of the solution to the unperturbed variational inequality is necessary and sufficient condition for convergence of the regularized perturbed inequality solutions.This research was supported in part by the Ministry of Science Grant 3481-1-91 and by the Ministry of Absorption Center for Absorption in Science.     

Finite sections for singular integral operators by weighted Chebyshev polynomials

A finite section method for the approximate solution of singular integral equations with piecewise continuous coefficients on intervals is considered. The problem is transformed in such a way that results which were previously obtained for singular integral equations on the unit circle using localization methods in Banach algebras are applicable to it. Thus, necessary and sufficient conditions for the stability of the approximation method can be proved.     

Elliptic regularization for the semi-linear telegraph system

The solution of the semi-linear telegraph system is compared with the solution of an elliptic regularization, to which one associates two-point boundary conditions. An asymptotic approximation for the solution of the elliptic regularization is constructed. The method employed here is the boundary function method due to Vishik and Lyusternik. The problem is singularly perturbed of elliptic-hyperbolic type. To conduct this analysis, high regularity with respect to t for the solutions of both problems is required. Finally, the order of this approximation is found in different spaces of functions.     

Singular integral operators with piecewise continuous coefficients in reflexive rearrangement-invariant spaces

The paper is devoted to some only recently uncovered phenomena emerging in the study of singular integral operators (SIO's) with piecewise continuous (PC) coefficients in reflexive rearrangement-invariant spaces over Carleson curves. We deal with several kinds of indices of submultiplicative functions which describe properties of spaces (Boyd and Zippin indices) and curves (spirality indices). We consider some disintegration condition which combines properties of spaces and curves, the Boyd and spirality indices.We show that the essential spectrum of SIO associated with the Riemann boundary value problem with PC coefficient arises from the essential range of the coefficient by filling in certain massive connected sets (so-called logarithmic leaves) between the endpoints of jumps.These results combined with the Allan-Douglas local principle and with the two projections theorem enable us to study the Banach algebra generated by SIO's with matrix-valued piecewise continuous coefficients. We construct a symbol calculus for this Banach algebra which provides a Fredholm criterion and gives a basis for an index formula for arbitrary SIO's from in terms of their symbols.     

Anti-selfdual Lagrangians II: unbounded non self-adjoint operators and evolution equations

New variational principles based on the concept of anti-selfdual (ASD) Lagrangians were recently introduced in “AIHP-Analyse non linéaire, 2006”. We continue here the program of using such Lagrangians to provide variational formulations and resolutions to various basic equations and evolutions which do not normally fit in the Euler-Lagrange framework. In particular, we consider stationary boundary value problems of the form as well ass dissipative initial value evolutions of the form where is a convex potential on an infinite dimensional space, A is a linear operator and is any scalar. The framework developed in the above mentioned paper reformulates these problems as and respectively, where is an “ASD” vector field derived from a suitable Lagrangian L. In this paper, we extend the domain of application of this approach by establishing existence and regularity results under much less restrictive boundedness conditions on the anti-selfdual Lagrangian L so as to cover equations involving unbounded operators. Our main applications deal with various nonlinear boundary value problems and parabolic initial value equations governed by transport operators with or without a diffusion term. Nassif Ghoussoub research was partially supported by a grant from the Natural Sciences and Engineering Research Council of Canada. The author gratefully acknowledges the hospitality and support of the Centre de Recherches Mathématiques in Montréal where this work was initiated. Leo Tzou’s research was partially supported by a doctoral postgraduate scholarship from the Natural Science and Engineering Research Council of Canada.     

On the spectral radius of positive operators on Banach sequence spaces

Let K₁,…,K_n be (infinite) non-negative matrices that define operators on a Banach sequence space. Given a function f:[0,∞)×…×[0,∞)→[0,∞) of n variables, we define a non-negative matrix and consider the inequality

     

A normalization problem for a class of singular integral operators with carleman shift and unbounded coefficients

The solution to a normalization problem for singular integral operators with Carleman shift and degenerate and unbounded coefficients inL _p () is obtained, where is either the unit circle or the real line. The approach followed consists mainly in two steps: the reduction to a singular integral operator with bounded coefficients and the use of the solution to an abstract normalization problem.This research was supported by JNICT under the grant PBIC/C/CEN/1040/92.     

Quadrature-difference methods for solving linear and nonlinear singular integro-differential equations

Here we propose and justify quadrature-difference methods for solving different kinds (linear, nonlinear and multidimensional) of periodic singular integro-differential equations.     

Fredholm differential operators with unbounded coefficients

We prove that a first-order linear differential operator G with unbounded operator coefficients is Fredholm on spaces of functions on with values in a reflexive Banach space if and only if the corresponding strongly continuous evolution family has exponential dichotomies on both and and a pair of the ranges of the dichotomy projections is Fredholm, and that the Fredholm index of G is equal to the Fredholm index of the pair. The operator G is the generator of the evolution semigroup associated with the evolution family. In the case when the evolution family is the propagator of a well-posed differential equation u′(t)=A(t)u(t) with, generally, unbounded operators , the operator G is a closure of the operator . Thus, this paper provides a complete infinite-dimensional generalization of well-known finite-dimensional results by Palmer, and by Ben-Artzi and Gohberg.