         double gpa = gpR0*inv_aspect_ratio;
         double Rmin=gpR0-1.0*gpa;
         double Zmin=-1.0*gpa*1.00;
-        double Rmax=gpR0+1.0*gpa; 
+        double Rmax=gpR0+1.0*gpa;
         double Zmax=1.0*gpa*1.00;
         dg::CylindricalMPIGrid3d g3d( Rmin,Rmax, Zmin,Zmax, 0, 2.*M_PI, n, Nx ,Ny, Nz,dg::DIR, dg::DIR, dg::PER,commEll);
         dg::geo::TokamakMagneticField magfield = dg::geo::createGuentherField(gpR0, gpI0);

        dg::geo::guenther::FuncNeu funcNEU(gpR0,gpI0);
        Vector function = dg::evaluate( funcNEU, g3d) , dsTdsfb(function);

        t.tic(); 
        for( unsigned i=0; i<multi; i++)
            ds.symv(function,dsTdsfb);
        t.toc();
        if(rank==0)std::cout<<t.diff()/(double)multi<<" ";
    }
    else 
        if(rank==0)std::cout<<0.<<" ";
    if(rank==0)std::cout << res.d<< " "<<res.i<<" ";

    } catch( std::exception& e) {
        if(rank==0)std::cout << std::endl;

    MPI_Cart_get( comm, 3, dims, periods, coords);
    if(rank==0)
    {
        std::cout<< dims[0] <<" "<<dims[1]<<" "<<dims[2]<<" "<<dims[0]*dims[1]*dims[2];
        int num_threads = 1;
#ifdef _OPENMP
        num_threads = omp_get_max_threads( );
#endif //omp
        std::cout << " "<<num_threads;
        std::cout<<" "<< n <<" "<<Nx<<" "<<Ny<<" "<<Nz;
    }


