The fractional Laplacian and fractional gradient are operators which play fundamental role in modeling of anomalous diffusion in d-dimensional space with the fractional exponent \(\alpha \in (1,2)\) . The principal-value integrals are split into singular and regular parts where we avoid using any weight function for the approximation in the singularity neighborhood. The resulting approximation coefficients are calculated from optimal value of the singular domain radius which is only a function of the exponent \(\alpha \) and a given grid topology. Various difference schemes are presented for the regular rectangular grids with mesh size \(h>0\) , and also for the hexagonal and the dodecahedral ones. This technique enables to evaluate the fractional operators with the approximation error \(\textrm{O}(h^{4-\alpha })\) which is verified using testing functions with known analytical expression of their fractional Laplacian and fractional gradient. Resulting formulas can be also used for the numeric solution of the fractional partial differential equations.
