diff --git a/dune/gfe/mixedriemanniantrsolver.cc b/dune/gfe/mixedriemanniantrsolver.cc index 166a063084e7633aac74d5558472883f66d0e747..2d4ed89a298f7eb1b0ee236e68072a6e1ce61f8b 100644 --- a/dune/gfe/mixedriemanniantrsolver.cc +++ b/dune/gfe/mixedriemanniantrsolver.cc @@ -129,7 +129,7 @@ setup(const GridType& grid, mmgStep0->basesolver_ = baseSolver0; mmgStep0->setSmoother(presmoother0, postsmoother0); mmgStep0->obstacleRestrictor_= new MandelObstacleRestrictor<CorrectionType0>(); - mmgStep0->verbosity_ = Solver::FULL; + mmgStep0->verbosity_ = Solver::QUIET; TrustRegionGSStep<MatrixType11, CorrectionType1>* presmoother1 = new TrustRegionGSStep<MatrixType11, CorrectionType1>; TrustRegionGSStep<MatrixType11, CorrectionType1>* postsmoother1 = new TrustRegionGSStep<MatrixType11, CorrectionType1>; @@ -141,7 +141,7 @@ setup(const GridType& grid, mmgStep1->basesolver_ = baseSolver1; mmgStep1->setSmoother(presmoother1, postsmoother1); mmgStep1->obstacleRestrictor_= new MandelObstacleRestrictor<CorrectionType1>(); - mmgStep1->verbosity_ = Solver::FULL; + mmgStep1->verbosity_ = Solver::QUIET; // ////////////////////////////////////////////////////////////////////////////////////// // Assemble a Laplace matrix to create a norm that's equivalent to the H1-norm @@ -428,7 +428,6 @@ void MixedRiemannianTrustRegionSolver<GridType,Basis0,TargetSpace0,Basis1,Target Dune::Timer solutionTimer; for (int ii=0; ii<200; ii++) { - std::cout << "Iteration " << ii << std::endl; residual0 = rhs_global0; stiffnessMatrix01.mmv(corr_global1, residual0); mmgStep0->setRhs(residual0); @@ -450,9 +449,6 @@ void MixedRiemannianTrustRegionSolver<GridType,Basis0,TargetSpace0,Basis1,Target double energy = 0.5 * (tmp0*corr_global0 + tmp1*corr_global1) - (rhs_global0*corr_global0 + rhs_global1*corr_global1); - - std::cout << "Energy: " << energy << std::endl; - if (energy > oldEnergy) //DUNE_THROW(Dune::Exception, "energy increase!"); std::cout << "Warning: energy increase!" << std::endl;