After using BICGSTAB solver with ILUT preconditioner , I got bg.error()=-1.#IND, forn about matrix size of 20000, after 1 iteraiton, what does that mean? Can anyone help me? There must be something wrong I think, but it works well for small size matrix..