I construct a Lanczos process on a large and sparse matrix and use the resu
lts of this iteration to compute the inverse square root of the same matrix
. The algorithm is a stable version of an earlier proposal by the author. I
t can be used fur problems related to the matrix sign and polar decompositi
on. The application here comes from the theory of chiral fermions on the la
ttice. (C) 2000 Academic Press.