diff --git a/dune/tectonic/factories/strikeslipfactory.cc b/dune/tectonic/factories/strikeslipfactory.cc index 130e6dd25eeb7c49884131bf2520b51e86df29fa..d10478cd65c38474481db119c5b7f4bdf9cc9c51 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());