The doubling algorithms are very efficient iterative methods for computing the unique minimal nonnegative solution to an M-matrix algebraic Riccati equation (MARE). They are globally and quadratically convergent, except for MARE in the critical case where convergence is linear with the linear rate 1/2. However, the initialization phase and the doubling iteration kernel of any doubling algorithm involve inverting nonsingular M-matrices. In particular, for MARE in the critical case, the M-matrices in the doubling iteration kernel, although nonsingular, move towards singular M-matrices at convergence. These inversions are causes of concerns on entrywise relative accuracy of the eventually computed minimal nonnegative solution. Fortunately, a nonsingular M-matrix can be inverted by the GTH-like algorithm of Alfa, Xue, and Ye [Math. Comp., 71:217--236, 2002] to almost full entrywise relative accuracy, provided a triplet representation of the matrix is known. Recently, Nguyen and Poloni [Numer. Math., 130(4):763--792, 2015] discovered a way to construct triplet representations in a cancellation-free manner for all involved M-matrices in the doubling iteration kernel, for a special class of MAREs arising from Markov-modulated fluid queues. In this paper, we extend Nguyen's and Poloni's work to all MAREs by also devising a way to construct the triplet representations cancellation-free. Our construction, however, is not a straightforward extension of theirs. It is made possible by an introduction of novel recursively computable auxiliary nonnegative vectors. As the second contribution, we propose an entrywise relative residual for an approximate solution. The residual has an appealing feature of being able to reveal the entrywise relative accuracies of all entries, large and small, of the approximation. This is in marked contrast to the usual legacy normalized residual which reflects relative accuracies of large entries well but not so much those of very tiny entries. Numerical examples are presented to demonstrate and confirm our claims.

This volume presents selections of Peter J. Bickels major papers, along with comments on their novelty and impact on the subsequent development of statistics as a discipline. Each of the eight parts concerns a particular area of research and provides new commentary by experts in the area. The parts range from Rank-Based Nonparametrics to Function Estimation and Bootstrap Resampling. Peters amazing career encompasses the majority of statistical developments in the last half-century or about about half of the entire history of the systematic development of statistics. This volume shares insights on these exciting statistical developments with future generations of statisticians. The compilation of supporting material about Peters life and work help readers understand the environment under which his research was conducted. The material will also inspire readers in their own research-based pursuits. This volume includes new photos of Peter Bickel, his biography, publication list, and a list of his students. These give the reader a more complete picture of Peter Bickel as a teacher, a friend, a colleague, and a family man.

No abstract uploaded!

No abstract uploaded!

We consider reflector imaging in a weakly random waveguide. We address the situation in which the source is farther from the reflector to be imaged than the energy equipartition distance, but the receiver array is closer to the reflector to be imaged than the energy equipartition distance. As a consequence, the reflector is illuminated by a partially coherent field and the signals recorded by the receiver array are noisy. This paper shows that migration of the recorded signals cannot give a good image, but an appropriate migration of the cross correlations of the recorded signals can give a very good image. The resolution and stability analysis of this original functional shows that the reflector can be localized with an accuracy of the order of the wavelength even when the receiver array has small aperture, and that broadband sources are necessary to ensure statistical stability, whatever the aperture of the array.