[Cleanup] Add toBool() convenience method

src/boundary_writer.cc

@@ -3,6 +3,7 @@
 #endif
 
 #include "boundary_writer.hh"
+#include "tobool.hh"
 
 template <class ScalarVector, class Vector>
 BoundaryWriter<ScalarVector, Vector>::BoundaryWriter(
@@ -16,7 +17,7 @@ BoundaryWriter<ScalarVector, Vector>::BoundaryWriter(
   std::fstream vertexCoordinateWriter(prefix + "Coordinates",
                                       std::fstream::out);
   for (size_t i = 0; i < boundaryNodes.size(); ++i)
-    if (boundaryNodes[i][0])
+    if (toBool(boundaryNodes[i]))
       vertexCoordinateWriter << vertexCoordinates[i] << std::endl;
   vertexCoordinateWriter.close();
 }
@@ -31,7 +32,7 @@ template <class ScalarVector, class Vector>
 void BoundaryWriter<ScalarVector, Vector>::writeKinetics(Vector const &u,
                                                          Vector const &v) {
   for (size_t i = 0; i < boundaryNodes.size(); ++i) {
-    if (!boundaryNodes[i][0])
+    if (!toBool(boundaryNodes[i]))
       continue;
 
     displacementWriter << projector_(u[i]) << " ";

src/friction_writer.cc

@@ -3,6 +3,7 @@
 #endif
 
 #include "friction_writer.hh"
+#include "tobool.hh"
 
 template <class ScalarVector, class Vector>
 FrictionWriter<ScalarVector, Vector>::FrictionWriter(
@@ -23,7 +24,7 @@ template <class ScalarVector, class Vector>
 void FrictionWriter<ScalarVector, Vector>::writeOther(
     ScalarVector const &coefficient, ScalarVector const &alpha) {
   for (size_t i = 0; i < boundaryNodes.size(); ++i) {
-    if (!boundaryNodes[i][0])
+    if (!toBool(boundaryNodes[i]))
       continue;
 
     coefficientWriter << coefficient[i] << " ";

src/sand-wedge.cc

@@ -68,6 +68,7 @@
 #include <dune/tectonic/globalnonlinearity.hh>
 
 #include "assemblers.hh"
+#include "tobool.hh"
 #include "enum_parser.cc"
 #include "enum_scheme.cc"
 #include "enum_state_model.cc"

src/state/dieterichstateupdater.hh

 #ifndef DIETERICH_STATE_UPDATER_HH
 #define DIETERICH_STATE_UPDATER_HH
 
+#include "tobool.hh"
 #include "stateupdater.hh"
 
 template <class ScalarVector, class Vector>
@@ -56,7 +57,7 @@ template <class ScalarVector, class Vector>
 void DieterichStateUpdater<ScalarVector, Vector>::solve(
     Vector const &velocity_field) {
   for (size_t i = 0; i < nodes.size(); ++i) {
-    if (not nodes[i][0])
+    if (not toBool(nodes[i]))
       continue;
 
     double const VoL = velocity_field[i].two_norm() / L;

src/state/ruinastateupdater.hh

 #ifndef RUINA_STATE_UPDATER_HH
 #define RUINA_STATE_UPDATER_HH
 
+#include "tobool.hh"
 #include "stateupdater.hh"
 
 template <class ScalarVector, class Vector>
@@ -42,7 +43,7 @@ template <class ScalarVector, class Vector>
 void RuinaStateUpdater<ScalarVector, Vector>::solve(
     Vector const &velocity_field) {
   for (size_t i = 0; i < nodes.size(); ++i) {
-    if (not nodes[i][0])
+    if (not toBool(nodes[i]))
       continue;
 
     double const VoL = velocity_field[i].two_norm() / L;

src/tobool.hh

+#ifndef TO_BOOL_HH
+#define TO_BOOL_HH
+
+template <class Alloc>
+bool toBool(Dune::BitSetVectorConstReference<1, Alloc> x) {
+  return x[0];
+}
+
+#endif