A fast preconditioned iterative method for two-dimensional options pricing under fractional differential models. (15th January 2020)