Add thrust for filling data

[jotachango]

To improve performance use thrust in that way:

    thrust::device_ptr<double> f_n_wrapped(f_n_ptr);
    thrust::fill(
            f_n_wrapped,
            f_n_wrapped + num_div_x*num_div_y,
            0.0);