I have just returned from a one-day workshop at U. C. Santa Barbara honoring John Gilbert on his 70th birthday and his official retirement after 20 years on the UCSB faculty.... 더 읽어보기 >>

(The January 5 posting was premature and incomplete.)... 더 읽어보기 >>

Christian Reinsch and Roland Bulirsch both passed away recently, Reinsch on October 8 and Bulirsch on September 21. Reinsch was 88 years old and Bulirsch was 89. Both of them were retired professors of numerical analysis at the Technical University of Munich. Both of them were friends of mine. But in almost all other ways, they were very different people.... 더 읽어보기 >>

I do not recall having seen this matrix before, but I will not be surprised to learn that somebody else already knows all about it, especially if that person's name is Nick.... 더 읽어보기 >>

The answer: A\A is always I, except when it isn't.... 더 읽어보기 >>

In February and March I published three blog posts about an enhancement request for MATLAB concerning tie breakers with the round function, including the ability to round ties to even. Round One, Round Two and Round Three. Since then, a group of us at MathWorks, organized by my colleague (and academic great-great-grand-descendant) Heiko Weichelt, have been considering the request. We had a virtual design review recently. This post describes a prototype incorporating the features that were discussed at that review.... 더 읽어보기 >>

I am giving a five-minute talk today, May 26, at the virtual seminar on Complexity of Matrix Computations. Here are my slides. Two new MATLAB functions, tred and imtql, instrumented to count flops, are available in symeig.m.... 더 읽어보기 >>

Computing Eigenvalues of Symmetric MatricesSee revision.Get the MATLAB code (requires JavaScript) Published with MATLAB®... 더 읽어보기 >>

This is about linear systems with fewer equations than variables; A*x = b where the m -by- n matrix A has fewer rows that columns, so m < n . I have always called such systems wide or fat, but this is not respectful. So I consulted the Merriam-Webster Thesaurus and found commodious.... 더 읽어보기 >>

I recently published Round, With Ties to Even and followed that with Round Two. Then, in an email, Andy Bartlett pointed out that my new round function fails for some large values of x between flintmax/2 and flintmax.... 더 읽어보기 >>