Commit 311793f7 authored by niklas.baumgarten's avatar niklas.baumgarten
Browse files

removing assemble transfer

parent cc216235
......@@ -554,19 +554,3 @@ double HybridEllipticAssemble::GoalFunctional(const Vector &u) const {
}
return PPM->Sum(goal);
}
void HybridEllipticAssemble::AssembleTransfer(TransferMatrix &TM) const {
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;
}
}
}
\ No newline at end of file
......@@ -60,8 +60,6 @@ public:
virtual void SetFlux(const Vector &u, Vector &flux) override;
void AssembleTransfer(TransferMatrix &TM) const override;
void LocalProblem(const cell &c, const Vector &u,
RMatrix &A, RVector &B, RVector &R) const;
......
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