diff --git a/src/lbm/SuViscoelasticity.h b/src/lbm/SuViscoelasticity.h index 9783b3981a85c1aa3995af9d35e385979981a3f1..493a3fa4987fb27d911317e93ee7686d5d1142e3 100644 --- a/src/lbm/SuViscoelasticity.h +++ b/src/lbm/SuViscoelasticity.h @@ -150,8 +150,6 @@ public: gradu(b, a) += velocity->get(cell + *d)[a] * real_c(d.c(b)) * real_c(0.5); } else if(boundaryHandling->isDomain(cell - *d)){ gradu(b, a) += velocity->get(cell - *d)[a] * real_c(d.c(b)) * real_c(-0.5) + velocity->get(cell)[a] * real_c(d.c(b)); - } else { - gradu(b, a) += real_t(0); } } } @@ -224,8 +222,6 @@ public: f[i] += stress->get(cell + *d)(k, i) * real_c(d.c(k)); } else if(boundaryHandling->isDomain(cell - *d)){ f[i] += -stress->get(cell - *d)(k, i) * real_c(d.c(k)) + stress->get(cell)(k, i) * real_c(d.c(k)) * real_c(2.0); - } else { - f[i] += real_t(0); } } }