 
     ////////////////////////////////set up computations///////////////////////////
     Geometry grid( 0, p.lx, 0, p.ly, p.n, p.Nx, p.Ny, p.bc_x, p.bc_y
-        #ifdef TOEFL_MPI
+        #ifdef POET_MPI
         , comm
-        #endif //TOEFL_MPI
+        #endif //POET_MPI
     );
     Geometry grid_out( 0, p.lx, 0, p.ly, p.n_out, p.Nx_out, p.Ny_out, p.bc_x, p.bc_y
-        #ifdef TOEFL_MPI
+        #ifdef POET_MPI
         , comm
-        #endif //TOEFL_MPI
+        #endif //POET_MPI
     );
     //create RHS
     poet::Explicit< Geometry, DMatrix, DDiaMatrix, DCooMatrix, DVec > ex( grid, p);

using DDiaMatrix =  cusp::dia_matrix<int, dg::get_value_type<DVec>, cusp::device_memory>;
using DCooMatrix =  cusp::coo_matrix<int, dg::get_value_type<DVec>, cusp::device_memory>;
#define MPI_OUT
#endif //TOEFL_MPI



int main( int argc, char* argv[])
{
#ifdef TOEFL_MPI
    ////////////////////////////////setup MPI///////////////////////////////
    int provided;
    MPI_Init_thread( &argc, &argv, MPI_THREAD_FUNNELED, &provided);

using DDiaMatrix =  cusp::dia_matrix<int, dg::get_value_type<DVec>, cusp::device_memory>;
using DCooMatrix =  cusp::coo_matrix<int, dg::get_value_type<DVec>, cusp::device_memory>;
#define MPI_OUT if(rank==0)
#else //POET_MPI
using HVec = dg::HVec;
using DVec = dg::DVec;
using HMatrix = dg::HMatrix;
