From 846a30c81a953ea0b8530d8beb4f4de2c02f2371 Mon Sep 17 00:00:00 2001 From: podlesny <podlesny@zedat.fu-berlin.de> Date: Tue, 19 Jan 2021 22:43:43 +0100 Subject: [PATCH] 3D not implemented --- dune/tectonic/factories/strikeslipfactory.cc | 11 ++++++----- 1 file changed, 6 insertions(+), 5 deletions(-) diff --git a/dune/tectonic/factories/strikeslipfactory.cc b/dune/tectonic/factories/strikeslipfactory.cc index 130e6dd2..d10478cd 100644 --- a/dune/tectonic/factories/strikeslipfactory.cc +++ b/dune/tectonic/factories/strikeslipfactory.cc @@ -6,6 +6,7 @@ #include <dune/contact/projections/normalprojection.hh> +#include "../problem-data/grid/cuboidgeometry.hh" #include "../problem-data/bc.hh" #include "../problem-data/myglobalfrictiondata.hh" @@ -25,13 +26,13 @@ void StrikeSlipFactory<HostGridType, VectorTEMPLATE>::setBodies() { const auto& frictionParset = this->parset_.sub("boundary.friction"); #if MY_DIM == 3 // TODO: not implemented - std::array<double, 2> depths = {this->parset_.template get<double>("body0.depth"), this->parset_.template get<double>("body1.depth")}; + //std::array<double, 2> depths = {this->parset_.template get<double>("body0.depth"), this->parset_.template get<double>("body1.depth")}; - cuboidGeometries_[0] = std::make_shared<TriangleGeometry>(origins[0], lengths[0], heights[0], depths[0]); - cuboidGeometries_[0]->addWeakeningPatch(frictionParset, {origins[0][0], origins[0][1]+ heights[0], 0}, {origins[0][0] + lengths[0], origins[0][1]+ heights[0], 0}); + //triangleGeometries_[0] = std::make_shared<TriangleGeometry>(origins[0], lengths[0], heights[0], depths[0]); + //triangleGeometries_[0]->addWeakeningPatch(frictionParset, {origins[0][0], origins[0][1]+ heights[0], 0}, {origins[0][0] + lengths[0], origins[0][1]+ heights[0], 0}); - cuboidGeometries_[1] = std::make_shared<TriangleGeometry>(origins[1], lengths[1], heights[1], depths[1]); - cuboidGeometries_[1]->addWeakeningPatch(frictionParset, origins[1], {origins[1][0] + lengths[1], origins[1][1], 0}); + //triangleGeometries_[1] = std::make_shared<TriangleGeometry>(origins[1], lengths[1], heights[1], depths[1]); + //triangleGeometries_[1]->addWeakeningPatch(frictionParset, origins[1], {origins[1][0] + lengths[1], origins[1][1], 0}); #elif MY_DIM == 2 triangleGeometries_[0] = std::make_shared<TriangleGeometry>(origin, lengths[0], heights[0]); triangleGeometries_[0]->addWeakeningPatch(frictionParset, triangleGeometries_[0]->A(), triangleGeometries_[0]->C()); -- GitLab