A $\tau$ MATRIX BASED PRECONDITIONING METHOD FOR SOLVING THE SYSTEM OF LINEAR EQUATIONS FROM THE DISCRETIZATION OF THE RIESZ SPACE FRACTIONAL ADVECTION-DISPERSION EQUATIONS