Skip to content
Snippets Groups Projects
Commit 91862454 authored by Sander, Oliver's avatar Sander, Oliver
Browse files

Disable some debugging output

parent 9175ef7d
No related branches found
No related tags found
No related merge requests found
...@@ -129,7 +129,7 @@ setup(const GridType& grid, ...@@ -129,7 +129,7 @@ setup(const GridType& grid,
mmgStep0->basesolver_ = baseSolver0; mmgStep0->basesolver_ = baseSolver0;
mmgStep0->setSmoother(presmoother0, postsmoother0); mmgStep0->setSmoother(presmoother0, postsmoother0);
mmgStep0->obstacleRestrictor_= new MandelObstacleRestrictor<CorrectionType0>(); mmgStep0->obstacleRestrictor_= new MandelObstacleRestrictor<CorrectionType0>();
mmgStep0->verbosity_ = Solver::FULL; mmgStep0->verbosity_ = Solver::QUIET;
TrustRegionGSStep<MatrixType11, CorrectionType1>* presmoother1 = new TrustRegionGSStep<MatrixType11, CorrectionType1>; TrustRegionGSStep<MatrixType11, CorrectionType1>* presmoother1 = new TrustRegionGSStep<MatrixType11, CorrectionType1>;
TrustRegionGSStep<MatrixType11, CorrectionType1>* postsmoother1 = new TrustRegionGSStep<MatrixType11, CorrectionType1>; TrustRegionGSStep<MatrixType11, CorrectionType1>* postsmoother1 = new TrustRegionGSStep<MatrixType11, CorrectionType1>;
...@@ -141,7 +141,7 @@ setup(const GridType& grid, ...@@ -141,7 +141,7 @@ setup(const GridType& grid,
mmgStep1->basesolver_ = baseSolver1; mmgStep1->basesolver_ = baseSolver1;
mmgStep1->setSmoother(presmoother1, postsmoother1); mmgStep1->setSmoother(presmoother1, postsmoother1);
mmgStep1->obstacleRestrictor_= new MandelObstacleRestrictor<CorrectionType1>(); 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 // 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 ...@@ -428,7 +428,6 @@ void MixedRiemannianTrustRegionSolver<GridType,Basis0,TargetSpace0,Basis1,Target
Dune::Timer solutionTimer; Dune::Timer solutionTimer;
for (int ii=0; ii<200; ii++) for (int ii=0; ii<200; ii++)
{ {
std::cout << "Iteration " << ii << std::endl;
residual0 = rhs_global0; residual0 = rhs_global0;
stiffnessMatrix01.mmv(corr_global1, residual0); stiffnessMatrix01.mmv(corr_global1, residual0);
mmgStep0->setRhs(residual0); mmgStep0->setRhs(residual0);
...@@ -450,9 +449,6 @@ void MixedRiemannianTrustRegionSolver<GridType,Basis0,TargetSpace0,Basis1,Target ...@@ -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); 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) if (energy > oldEnergy)
//DUNE_THROW(Dune::Exception, "energy increase!"); //DUNE_THROW(Dune::Exception, "energy increase!");
std::cout << "Warning: energy increase!" << std::endl; std::cout << "Warning: energy increase!" << std::endl;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment