     std::cout << "# of polynomial coefficients: "<< n <<std::endl;
     std::cout << "# of 2d cells                 "<< Nx*Ny <<std::endl;
     
+    std::cout << "... for a precision of "<< eps<<std::endl;
     t.tic();
     std::cout << "Number of pcg iterations "<< pcg( dA, dx, db, t2d_d, eps)<<std::endl;
     t.toc();
-    std::cout << "... for a precision of "<< eps<<std::endl;
     std::cout << "... on the device took "<< t.diff()<<"s\n";
     t.tic();
     dg::cg( dA, dx_, db_, t2d_d, eps, dx_.size());

    dg::DMatrix dA = dg::create::laplacianM( grid, dg::not_normed, dg::forward); 
    dg::DMatrix DX = dg::create::dx( grid);
    dg::HMatrix A = dA;
    t.toc();
    std::cout<< "Creation took "<<t.diff()<<"s\n";


#include "backend/cusp_thrust_backend.h"

#include "cg.h"
#include "laplace.h"

//leo3 can do 350 x 350 but not 375 x 375
const double ly = 2.*M_PI;
