From b31e082d634bce8b63edfaa314f2d87a94655cb0 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Sun, 14 Jul 2013 15:18:42 +0200
Subject: [PATCH] [Cleanup] Kill an assertion (the code was buggy, too)
---
dune/tectonic/myblockproblem.hh | 20 --------------------
1 file changed, 20 deletions(-)
diff --git a/dune/tectonic/myblockproblem.hh b/dune/tectonic/myblockproblem.hh
index f1c8eae9..dc52328b 100644
--- a/dune/tectonic/myblockproblem.hh
+++ b/dune/tectonic/myblockproblem.hh
@@ -170,26 +170,6 @@ template <class MyConvexProblemTypeTEMPLATE> class MyBlockProblem {
// -grad is needed for Newton step
linearization.b *= -1.0;
- // b should be a descent direction
- {
- VectorType const direction = linearization.b;
- VectorType tmp = linearization.b; // b
- Arithmetic::addProduct(tmp, -1.0, linearization.A, u); // b-Au
- double const localA = tmp * direction; // <b-Au,v>
-
- linearization.A.mv(direction, tmp); // Av
- double const localb = tmp * direction; // <Av,v>
-
- MyDirectionalConvexFunction<
- Dune::GlobalNonlinearity<MatrixType, VectorType>> const
- psi(localA, localb, problem.phi, u, direction);
-
- Interval<double> D;
- psi.subDiff(0, D);
- if (!isnan(D[1]))
- assert(D[1] <= 0);
- }
-
// apply truncation to stiffness matrix and rhs
for (size_t row = 0; row < linearization.A.N(); ++row) {
auto const col_end = linearization.A[row].end();
--
GitLab