diff --git a/src/TNL/Operators/diffusion/LinearDiffusion_impl.h b/src/TNL/Operators/diffusion/LinearDiffusion_impl.h
index a458a2c8a2174fb0adcf335f1d64433d998c2b27..17c39bde99e163d1885d91e7efef5da6bc37f451 100644
--- a/src/TNL/Operators/diffusion/LinearDiffusion_impl.h
+++ b/src/TNL/Operators/diffusion/LinearDiffusion_impl.h
@@ -153,25 +153,10 @@ operator()( const PreimageFunction& u,
 {  
    static_assert( EntityType::entityDimensions == 2, "Wrong mesh entity dimensions." );
    static_assert( PreimageFunction::getEntitiesDimensions() == 2, "Wrong preimage function" );
-  // std::cerr << "Lin Diff op:" <<__LINE__ << std::endl<< std::flush;
    const typename EntityType::template NeighbourEntities< 2 >& neighbourEntities = entity.getNeighbourEntities();
-  // std::cerr << "Lin Diff op:" <<__LINE__ << std::endl<< std::flush;
    const RealType& hxSquareInverse = entity.getMesh().template getSpaceStepsProducts< -2, 0 >();
- //  std::cerr << "Lin Diff op:" <<__LINE__ << std::endl<< std::flush;
    const RealType& hySquareInverse = entity.getMesh().template getSpaceStepsProducts< 0, -2 >();
- //  std::cerr << "Lin Diff op:" <<__LINE__ << std::endl<< std::flush;
-   
-   //int stop;
-   //std::cin >> stop;
-   
- /*  std::cerr << "Lin Diff indexes:" <<__LINE__ << "\t"
-           << neighbourEntities.template getEntityIndex< -1,  0 >()<<"\t"
-           << neighbourEntities.template getEntityIndex< 1,  0 >()<<"\t"
-           <<neighbourEntities.template getEntityIndex< 0,  -1 >()<<"\t"
-           <<neighbourEntities.template getEntityIndex< 0,  1 >()<<"\t"
-           <<std::endl<< std::flush;*/
-   
-   //std::cin >> stop;
+
    return ( u[ neighbourEntities.template getEntityIndex< -1,  0 >() ]
           + u[ neighbourEntities.template getEntityIndex<  1,  0 >() ] ) * hxSquareInverse +
           ( u[ neighbourEntities.template getEntityIndex<  0, -1 >() ]
diff --git a/src/TNL/Solvers/ODE/Euler_impl.h b/src/TNL/Solvers/ODE/Euler_impl.h
index 2637a0dc06a0d3b330b9b0353d7a9a16cccd1972..96e631bebc4622b54a0b386fa4762a47fbf56f90 100644
--- a/src/TNL/Solvers/ODE/Euler_impl.h
+++ b/src/TNL/Solvers/ODE/Euler_impl.h
@@ -240,10 +240,7 @@ void Euler< Problem > :: computeNewTimeLevel( DofVectorPointer& u,
 
    
    
-   localResidue /= tau * ( RealType ) size;
-   
-   std::cout << time << "   " << tau << std:: endl;
-   
+   localResidue /= tau * ( RealType ) size;   
    MPIAllreduce( localResidue, currentResidue, 1, MPI_SUM, this->solver_comm );
 
 }