Commit 8f1dba97 authored by Dominic Kempf's avatar Dominic Kempf
Browse files

Remove debug statements

parent 89ae5390
......@@ -352,11 +352,6 @@ class EulerBernoulli2DLocalOperator
// The tangential vector for the curve
auto t = fibre->tangent(it->second);
auto n = ig.intersection().centerUnitOuterNormal();
if (( t*n < 0) != flipped) {
std::cout << "FLIPPED IS WRONG" << std::endl;
}
// Extract physical parameter of the fibre
auto E = fibre_modulus[fibindex];
auto d = fibre_radii[fibindex];
......
......@@ -58,13 +58,13 @@ class FibreReinforcedBulkOperator
virtual void alpha_volume(const EG& eg, const LFSU& lfsu, const X& x, const LFSV& lfsv, R& r) const override
{
//bulkoperator.alpha_volume(eg, lfsu, x, lfsv, r);
//fibreoperator.alpha_volume(eg, lfsu, x, lfsv, r);
bulkoperator.alpha_volume(eg, lfsu, x, lfsv, r);
fibreoperator.alpha_volume(eg, lfsu, x, lfsv, r);
}
virtual void alpha_boundary(const IG& ig, const LFSU& lfsu, const X& x, const LFSV& lfsv, R& r) const override
{
//bulkoperator.alpha_boundary(ig, lfsu, x, lfsv, r);
bulkoperator.alpha_boundary(ig, lfsu, x, lfsv, r);
}
virtual void alpha_skeleton(const IG& ig, const LFSU& lfsu_s, const X& x_s, const LFSV& lfsv_s, const LFSU& lfsu_n, const X& x_n, const LFSV& lfsv_n, R& r_s, R& r_n) const
......
......@@ -197,7 +197,7 @@ int main(int argc, char** argv)
ISTLV& istlv = Dune::PDELab::Backend::native(v);
v = 1.0;
set_constrained_dofs(cc,0.0,v); // dirichlet dofs are zero, others are 1, so we can use it as a mask!
// solver.apply(istlv,istlr,stat);
solver.apply(istlv,istlr,stat);
istldisplacement -= istlv;
//==================
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment