     std::cout << "Number of pcg iterations "<< pcg( laplace, x, b, v3d, eps)<<std::endl;
     t.toc();
     std::cout << "... on the device took "<< t.diff()<<"s\n";
-    //x = dg::evaluate( initial, grid);
-    //t.tic();
-    //std::cout << "Number of mixed pcg iterations "<< pcg( laplace, x, b, inverse, v3d, eps)<<std::endl;
-    //t.toc();
-    //std::cout << "... on the device took "<< t.diff()<<"s\n";
+    x = dg::evaluate( initial, grid);
+    t.tic();
+    std::cout << "Number of mixed pcg iterations "<< pcg( laplace, x, b, inverse, v3d, eps)<<std::endl;
+    t.toc();
+    std::cout << "... on the device took "<< t.diff()<<"s\n";
     dg::DVec  error(  solution);
     dg::blas1::axpby( 1., x,-1., error);
 

//#define CUSP_DEVICE_BLAS_SYSTEM CUSP_DEVICE_BLAS_CUBLAS
#include <iostream>
#include <iomanip>


#include "backend/sparseblockmat.cuh"
#include "backend/derivatives.h"
#include "backend/typedefs.cuh"
#include "backend/cusp_thrust_backend.h"

#include "cg.h"
#include "elliptic.h"
