A real algorithm for the Hermitian eigenvalue decomposition

Modifying complex plane rotations, we derive a new Jacobi-type algorithm for the Hermitian eigendecomposition, which uses only real arithmetic. When the fast-scaled rotations are incorporated, the new algorithm brings a substantial reduction in computational costs. The new method has the same convergence properties and parallelism as the symmetric Jacobi algorithm. Computational test results show that it produces accurate eigenvalues and eigenvectors and achieves great reduction in computational time.The work of this author was supported in part by the National Science Foundation grant CCR-8813493 and by the University of Minnesota Army High Performance Computing Research Center contract DAAL 03-89-C-0038.The work of this author was supported in part by the University of Minnesota Army High Performance Computing Research Center contract DAAL 03-89-C-0038.     

Diffusing particles with electrostatic repulsion

Summary. We study a diffusion model of an interacting particles system with general drift and diffusion coefficients, and electrostatic inter-particles repulsion. More precisely, the finite particle system is shown to be well defined thanks to recent results on multivalued stochastic differential equations (see [2]), and then we consider the behaviour of this system when the number of particles goes to infinity (through the empirical measure process). In the particular case of affine drift and constant diffusion coefficient, we prove that a limiting measure-valued process exists and is the unique solution of a deterministic PDE. Our treatment of the convergence problem (as ) is partly similar to that of T. Chan [3] and L.C.G. Rogers - Z. Shi [5], except we consider here a more general case allowing collisions between particles, which leads to a second-order limiting PDE. Received: 5 August 1996 / In revised form: 17 October 1996     

A note on minimization problems and multistep methods

Summary. We use the qualitative properties of the solution flow of the gradient equation to compute a local minimum of a real-valued function . Under the regularity assumption of all equilibria we show a convergence result for bounded trajectories of a consistent, strictly stable linear multistep method applied to the gradient equation. Moreover, we compare the asymptotic features of the numerical and the exact solutions as done by Humphries, Stuart (1994) and Schropp (1995) for one-step methods. In the case of -stable formulae this leads to an efficient solver for stiff minimization problems. Received July 10, 1995 / Revised version received June 27, 1996     

Thermal diffusivity measurements in porous ceramics by photothermal methods

Received: 16 September 1996/Accepted: 6 January 1997     

Poincaré’s inequalities and Talagrand’s concentration phenomenon for the exponential distribution

Summary. We present a simple proof, based on modified logarithmic Sobolev inequalities, of Talagrand’s concentration inequality for the exponential distribution. We actually observe that every measure satisfying a Poincaré inequality shares the same concentration phenomenon. We also discuss exponential integrability under Poincaré inequalities and its consequence to sharp diameter upper bounds on spectral gaps. Received: 10 June 1996 / In revised form: 9 August 1996     

The Jacobian and the Ginzburg-Landau energy

We study the Ginzburg-Landau functional for , where U is a bounded, open subset of . We show that if a sequence of functions satisfies , then their Jacobians are precompact in the dual of for every . Moreover, any limiting measure is a sum of point masses. We also characterize the -limit of the functionals , in terms of the function space B2V introduced by the authors in [16,17]: we show that I(u) is finite if and only if , and for is equal to the total variation of the Jacobian measure Ju. When the domain U has dimension greater than two, we prove if then the Jacobians are again precompact in for all , and moreover we show that any limiting measure must be integer multiplicity rectifiable. We also show that the total variation of the Jacobian measure is a lower bound for the limit of the Ginzburg-Landau functional. Received: 15 December 2000 / Accepted: 23 January 2001 / Published online: 25 June 2001     

Spline prewavelets for non-uniform knots

Summary In this paper, we develop a framework suitable for performing a multiresolution analysis using univariate spline spaces of arbitrary degree and with non-uniform knot-sequences. To this end, we show, among other things, the existence of compactly supported prewavelets and of prewavelets that are globally supported, but decay exponentially. In each case we obtain a decomposition of a fine spline space as a sum of a coarse spline space plus a spline space spanned by prewavelets.     

Spline collocation for strongly elliptic equations on the torus

Summary We prove convergence and error estimates in Sobolev spaces for the collocation method with tensor product splines for strongly elliptic pseudodifferential equations on the torus. Examples of applications include elliptic partial differential equations with periodic boundary conditions but also the classical boundary integral operators of potential theory on torus-shaped domains in three or more dimensions. For odd-degree splines, we prove convergence of nodal collocation for any strongly elliptic operator. For even-degree splines and midpoint collocation, we find an additional condition for the convergence which is satisfied for the classical boundary integral operators. Our analysis is a generalization to higher dimensions of the corresponding analysis of Arnold and Wendland [4].