function [v,lambda,rr]=InversePower(b,x0,tol,nmax);
% Inverse power method for the Hermitian Wilson operator
v=b/norm(b); rr=[];
while 1,
  u=bicg5(v,x0,1e-6,1000);
  u=mult_gamma5(u);
  lambda=v'*u;
  r=v-u/lambda;
  rnorm=norm(r); rr=[rr;rnorm];
  if rnorm<tol, break, end
  v=u/norm(u);
end
% Copyright (C) 2006-2007 Artan Borici.
% This program is a free software licensed under the terms of the GNU General Public License