Commit 8921db95 authored by niklas.baumgarten's avatar niklas.baumgarten
Browse files

removing assemble transfer

parent 5ea922d3
......@@ -360,20 +360,3 @@ void MixedEllipticAssemble::SetFlux(const Vector &u, Vector &flux) {
}
}
}
void MixedEllipticAssemble::AssembleTransfer(TransferMatrix &TM) const {
// Only works for rt0_p0
TM = 0;
const matrixgraph &cg = TM.CoarseMatrixGraph();
for (cell c = cg.cells(); c != cg.cells_end(); ++c) {
for (int cf = 0; cf < c.Faces(); ++cf) {
for (int k = 0; k < 2; ++k)
TM(c.Face(cf), 0.5 * (c.Face(cf) + c.FaceCorner(cf, k)))[0] = 0.5;
if (c.Faces() == 3) {
for (int f = cf + 1; f < c.Faces(); ++f)
TM(c.Face(cf), 0.5 * (c.Face(cf) + c.Face(f)))[0] = 0.5;
} else
TM(c.Face(cf), 0.5 * (c.Face(cf) + c()))[0] = 0.5;
}
}
}
......@@ -57,8 +57,6 @@ public:
virtual void SetExactSolution(Vector &uEx) const override;
virtual void SetFlux(const Vector &u, Vector &flux) override;
void AssembleTransfer(TransferMatrix &TM) const override;
};
#endif
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