Introduced maxiter parameter to limit maximum number of iterations. Improved scaling with no. of processors.

Minor bug fix (courtesy of Sergey Koposov).