- added missing fillComplete() call when matrix has been instantiated. - added support for BiCGStab. - switched to right-preconditioning if using BiCGStab (left is not supported in Trilinos yet) - now calling matrix->apply for mvp which is the recommended way.