On the augmented system approach to sparse least-squares problems

Summary We study the augmented system approach for the solution of sparse linear least-squares problems. It is well known that this method has better numerical properties than the method based on the normal equations. We use recent work by Arioli et al. (1988) to introduce error bounds and estimates for the components of the solution of the augmented system. In particular, we find that, using iterative refinement, we obtain a very robust algorithm and our estimates of the error are accurate and cheap to compute. The final error and all our error estimates are much better than the classical or Skeel's error analysis (1979) indicates. Moreover, we prove that our error estimates are independent of the row scaling of the augmented system and we analyze the influence of the Björck scaling (1967) on these estimates. We illustrate this with runs both on large-scale practical problems and contrived examples, comparing the numerical behaviour of the augmented systems approach with a code using the normal equations. These experiments show that while the augmented system approach with iterative refinement can sometimes be less efficient than the normal equations approach, it is comparable or better when the least-squares matrix has a full row, and is, in any case, much more stable and robust.This author was visiting Harwell and was funded by a grant from the Italian National Council of Research (CNR), Istituto di Elaborazione dell'Informazione-CNR, via S. Maria 46, I-56100 Pisa, ItalyThis author was visiting Harwell from Faculty of Mathematics and Computer Science of the University of Amsterdam     

Inclusive\bar K^{0*} (892) andK 0*(892) production on silicon by 200 GeVK − andπ −

The production of the neutralK ^? (892) resonances by 200 GeVK ^? andπ ^? has been studied over the kinematic range 0.0<x _f<1.0 andp _t ² <5.0 GeV². Longitudinal and transverse momentum distributions are presented. In addition the decay angular distributions inK ^? fragmentation to \(\bar K^{0*} \) have been investigated.     

Lifetimes of charged and neutralD mesons

We have measured the lifetimes of hadronically produced charged and neutralD mesons using silicon microstrip detectors and an active silicon target in the NA32 spectrometer at the CERN SPS. We obtainτ _D± = (10.9± _1.5 ^1.9 )·10^?13s andτ _D ^(?)10 = (4.2±0.5)·10^?13s based on 59 and 90 fully reconstructed decays respectively, giving a ratioτ _D±/τ _D ^(?)10 of 2.6 ±0.5.     

Joint production of ϕ mesons andπ ±,π 0,p,\bar p,K s 0 andK ± in hadronic interactions

The joint production of ? mesons andπ ^±,π ⁰,p, \(\bar p\) ,K _s ⁰ andK ^± is investigated using a sample of 600,000 inclusive ? meson events obtained in hadron Be interactions with incidentπ ^±,p, \(\bar p\) andK ^± beams. Evidence is presented for the joint production of ? mesons and strange particles produced with non-strange incident beams. With incidentK ^± beam the number of additional strange particles is suppressed. The results are found to be in agreement with the qualitative predictions of a parton fusion model. The comparison with the Lund model for lowp _T processes is fair.     

InclusiveK *0(890) and\bar K*^0 (890) production on beryllium byK − and π− at 175 GeV

The reactions \(K^ - Be \to {}^(\bar K^) *^0 (890)X,\pi ^ - Be \to {}^(\bar K^) *^0 (890)X\) , have been studied in a 175 GeV unseparated hadron beam in the kinematic range 0<x _F <1.0 andp _T ² <5 GeV². Integrated cross-sections and the dependence of the cross-sections on the longitudinal and transverse momentum are presented, together with quark counting rules predictions. The nuclear dependence ofK ^? fragmentation intoK ^*0(890) with respect to Feynmanx is investigated from hydrogen to beryllium.