         dg::geo::guenther::FuncNeu funcNEU(gpR0,gpI0);
         Vector function = dg::evaluate( funcNEU, g3d) , dsTdsfb(function);
 
-        t.tic(); 
+        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<<" ";
+        if(rank==0)std::cout<<" "<<t.diff()/(double)multi;
+    }
+    else
+        if(rank==0)std::cout<<" 0.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<<" ";
    }



    t.tic();
    unsigned number = pcg(laplace, x, b, ellv3d, 1e-6);
    t.toc();
    if(rank==0)std::cout <<" "<< number << " "<<t.diff()/(double)number<<std::flush;
    dg::blas1::axpby( 1., solution, -1., x);
    exblas::udouble res;
    res.d = dg::blas2::dot( x, ellw3d, x);
