     for( unsigned i=0; i<multi; i++)
         arakawa( lhs, rhs, jac);
     t.toc();
-    if(rank==0)std::cout<<t.diff()/(double)multi<<" ";
+    if(rank==0)std::cout<<t.diff()/(double)multi<<" "<<std::flush;
     //The Elliptic scheme
     periods[0] = false, periods[1] = false;
     MPI_Comm commEll;

        if(rank==0)std::cout<<0.<<" ";
    if(rank==0)std::cout << res.d<< " "<<res.i<<" ";


    if(rank==0)std::cout <<std::endl;
    MPI_Finalize();


    dg::CartesianMPIGrid3d grid( 0, lx, 0, ly, 0,lz, n, Nx, Ny, Nz, bcx, bcy, dg::PER, comm);
    dg::Timer t;
    Vector w3d, lhs, rhs, jac, x, y, z;
    try{
        w3d = dg::transfer<Vector>( dg::create::weights( grid));
        lhs = dg::transfer<Vector>( dg::evaluate ( left, grid));
        rhs = dg::transfer<Vector>( dg::evaluate ( right,grid));
        jac = dg::transfer<Vector>( dg::evaluate ( jacobian,grid));
        x = y = z = lhs;
    }
    catch( std::exception& e)
    {
        if(rank==0)std::cout << std::endl;
        if(rank==0)std::cerr << "Caught std::exception: "<<e.what()<<std::endl;
        MPI_Finalize();
        return 0;
    }
    std::cout<< std::setprecision(6);
    unsigned multi=100;

