From 7f3b651c17b7f84232b6912f321448e8e3f908ef Mon Sep 17 00:00:00 2001 From: "hanouvit@vz" <vithanousek@seznam.cz> Date: Thu, 2 Mar 2017 18:12:09 +0100 Subject: [PATCH] =?UTF-8?q?Odstranen=C3=AD=20Debug=20Messag=C5=AF?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- .../Operators/diffusion/LinearDiffusion_impl.h | 17 +---------------- src/TNL/Solvers/ODE/Euler_impl.h | 5 +---- 2 files changed, 2 insertions(+), 20 deletions(-) diff --git a/src/TNL/Operators/diffusion/LinearDiffusion_impl.h b/src/TNL/Operators/diffusion/LinearDiffusion_impl.h index a458a2c8a2..17c39bde99 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 2637a0dc06..96e631bebc 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 ); } -- GitLab