Commit d05092df authored by niklas.baumgarten's avatar niklas.baumgarten
Browse files

bug fix outflow

parent d931a42e
......@@ -261,7 +261,7 @@ pair<double, double> EllipticAssemble::getInflowOutflow(const Vector &u) const {
BFParts bnd(u.GetMesh(), c);
if (!bnd.onBnd()) continue;
for (int face = 0; face < c.Faces(); ++face) {
if (bnd[face] < 0) continue;
if (bnd[face] != 1) continue;
auto *faceElem = getFaceElement(c, u, face);
for (int q = 0; q < faceElem->nQ(); ++q) {
double w = faceElem->QWeight(q);
......
Supports Markdown
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