 #include "parameters.h"
 // #include "geometry_circ.h"
 #include "dg/average.h"
+#include "dg/backend/linspace.cuh"
 
 #ifdef DG_BENCHMARK
 #include "dg/backend/timer.cuh"

    invert_invgammaN(invgammaNU,target,src); //=ne-1 = Gamma (ni-1)    
}




template<class Matrix, class container, class P>

    double mass_, energy_, diff_, ediff_,gammanex_,coupling_;
    std::vector<double> evec;
    //probe
     std::vector<container> probevec;
     const container probe_coord_X, probe_coord_Y;
     Matrix probeinterp;
//     container probevalue;    
//     dg::Grid1d<double> gy;
//     const container w1d;
