Commit 5ea922d3 authored by niklas.baumgarten's avatar niklas.baumgarten
Browse files

removing assemble transfer

parent 311793f7
......@@ -334,23 +334,3 @@ void LagrangeEllipticAssemble::SetFlux(const Vector &u, Vector &flux) {
}
}
}
void LagrangeEllipticAssemble::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.Corners(); ++cf) {
TM(c.Corner(cf), c.Corner(cf))[0] = 1;
}
for (int ce = 0; ce < c.Edges(); ++ce) {
for (int cf = 0; cf < 2; ++cf) {
TM(c.EdgeCorner(ce, cf), c.Edge(ce))[0] = 0.5;
}
}
if (c.Corners() == 4) {
for (int cf = 0; cf < c.Corners(); ++cf) {
TM(c.Corner(cf), c())[0] = 0.25;
}
}
}
}
......@@ -56,8 +56,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
......
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