From 473b833ba2c79d46497a539ebc6853688c65f791 Mon Sep 17 00:00:00 2001
From: podlesny <podlesny@zedat.fu-berlin.de>
Date: Mon, 15 Mar 2021 23:29:38 +0100
Subject: [PATCH] debug print
---
.../tectonic/spatial-solving/tnnmg/linearization.hh | 13 ++++++++++++-
1 file changed, 12 insertions(+), 1 deletion(-)
diff --git a/dune/tectonic/spatial-solving/tnnmg/linearization.hh b/dune/tectonic/spatial-solving/tnnmg/linearization.hh
index 3d986305..f266bf32 100644
--- a/dune/tectonic/spatial-solving/tnnmg/linearization.hh
+++ b/dune/tectonic/spatial-solving/tnnmg/linearization.hh
@@ -52,7 +52,7 @@ class Linearization {
}
}
- //std::cout << "regularityTruncation: " << count << std::endl;
+ std::cout << "regularityTruncation: " << count << std::endl;
}
template<class NV, class NBV, class T>
@@ -121,6 +121,7 @@ class Linearization {
BitVector obstacleTruncationFlags = ignore_;
BitVector regularityTruncationFlags = ignore_;
+ //regularityTruncationFlags.unsetAll();
//std::cout << "ignore truncation: " << truncationFlags_.count();
determineTruncation(x, f_.lowerObstacle(), f_.upperObstacle(), obstacleTruncationFlags, obstacleTruncationTolerance_); // obstacle truncation
@@ -177,11 +178,13 @@ class Linearization {
// nonlinearity part
phi.addHessian(x, hessian_);
+
//truncateMatrix(hessian_, regularityTruncationFlags, regularityTruncationFlags);
//std::cout << x[0] << " " << hessian_[0][0] << std::endl;
B -= hessian_;
+ //truncateMatrix(B, regularityTruncationFlags, regularityTruncationFlags);
//print(B, "phi hessian: ");
// compute gradient
@@ -198,6 +201,14 @@ class Linearization {
phi.addGradient(x, negativeGradient_);
C -= negativeGradient_;
+ //truncateVector(C, truncationFlags_);
+
+ /*for (size_t i=0; i<C.size(); i++) {
+ std::cout << "dof " << i << std::endl;
+ std::cout << "truncation: " << regularityTruncationFlags[i] << std::endl;
+ std::cout << "C[i]: " << C[i] << std::endl;
+ }*/
+ //print(truncationFlags_, "truncationFlags_: ");
//print(C, "phi gradient: ");
// -grad is needed for Newton step
--
GitLab