[Output] Friction writing: Report (negative) horizontal instead of x component

src/sand-wedge.cc

@@ -100,11 +100,6 @@ template <class GridView> class SpecialWriter {
   using ElementPointer =
       typename GridView::Grid::template Codim<0>::EntityPointer;
 
-  void writeMagnitude(LocalVector const &v) {
-    writer_ << v[0] << " "; // Because of Dirichlet conditions, this is like a
-                            // signed norm
-  }
-
   void writeHorizontal(LocalVector const &v) {
     writer_ << MyGeometry::horizontalProjection(v) << " ";
   }
@@ -140,7 +135,7 @@ template <class GridView> class SpecialWriter {
         I(globalToLocal(MyGeometry::I)),
         U(globalToLocal(MyGeometry::U)),
         Z(globalToLocal(MyGeometry::Z)) {
-    writer_ << "|G| |H| |J| |I| Uv Uh Zv Zh" << std::endl;
+    writer_ << "Gh Hh Jh Ih Uv Uh Zv Zh" << std::endl;
   }
 
   void write(VirtualGridFunction<typename GridView::Grid, LocalVector> const &
@@ -148,16 +143,16 @@ template <class GridView> class SpecialWriter {
     LocalVector value;
 
     specialField.evaluateLocal(*G.first, G.second, value);
-    writeMagnitude(value);
+    writeHorizontal(value);
 
     specialField.evaluateLocal(*H.first, H.second, value);
-    writeMagnitude(value);
+    writeHorizontal(value);
 
     specialField.evaluateLocal(*J.first, J.second, value);
-    writeMagnitude(value);
+    writeHorizontal(value);
 
     specialField.evaluateLocal(*I.first, I.second, value);
-    writeMagnitude(value);
+    writeHorizontal(value);
 
     specialField.evaluateLocal(*U.first, U.second, value);
     writeVertical(value);
@@ -361,7 +356,7 @@ int main(int argc, char *argv[]) {
 
     FrictionWriter<ScalarVector, Vector> frictionWriter(
         vertexCoordinates, frictionalNodes, "friction",
-        [](LocalVector const &x) { return x[0]; });
+        MyGeometry::horizontalProjection);
     BoundaryWriter<ScalarVector, Vector> verticalSurfaceWriter(
         vertexCoordinates, surfaceNodes, "verticalSurface",
         MyGeometry::verticalProjection);