diff --git a/src/sand-wedge.cc b/src/sand-wedge.cc index 33be5768821928696e823cabc602329ff82bb4ec..b1c6178d04c75e3038beb047f6bf392edcab35f3 100644 --- a/src/sand-wedge.cc +++ b/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);