print regularity truncation info

dune/tectonic/utils/debugutils.hh

@@ -14,6 +14,34 @@
 
 //using namespace std;
 
+    template <class Friction, typename VectorType>
+    void printRegularityTruncation(const Friction& friction, const VectorType& x, double truncationTolerance = 1e8) {
+        using BitVector = Dune::BitSetVector<MY_DIM>;
+
+        BitVector truncationFlags(x.size());
+        truncationFlags.unsetAll();
+
+        size_t count = 0;
+        size_t vsmall = 0;
+
+
+        for (size_t i = 0; i < x.size(); ++i) {
+          //std::cout << f_.phi().regularity(i, x[i]) << " xnorm: " << x[i].two_norm() << std::endl;
+          auto tangential_x = x[i];
+          tangential_x[0] = 0.0;
+
+          if (tangential_x.two_norm()<1e-14) {
+              vsmall++;
+          }
+
+          if (friction.regularity(i, x[i]) > truncationTolerance) {
+            count++;
+          }
+        }
+
+         std::cout << "V<1e-14: " << vsmall << " regularityTruncation: " << count << std::endl;
+    }
+
     template <typename VectorType>
     auto average(const VectorType& vec) {
         using BlockType = typename VectorType::block_type;