diff --git a/disc/operators/p1mgtransfer.hh b/disc/operators/p1mgtransfer.hh index 2a5dbbe97600822ff2bfa356dd3beff98b18f344..3f0b11699b4b614770248565d0426da203bd74ae 100644 --- a/disc/operators/p1mgtransfer.hh +++ b/disc/operators/p1mgtransfer.hh @@ -286,28 +286,6 @@ namespace Dune treated[indexi] = true; } } - - - // eliminate interpolation in rows with essential boundary conditions - // rowiterator endi=(*A).end(); - // for (rowiterator i=(*A).begin(); i!=endi; ++i) - // { - // int indexi = i.index(); - - // // go through all components at this vertex - // for (int compi=0; compi<m; compi++) - // if ((*skipflag[compi])[indexi]) - // { - // // we have to set this row to zero - // coliterator colend = (*i).end(); - // for (coliterator colit=(*i).begin(); colit!=colend; ++colit) - // for (int compj=0; compj<m; compj++) - // (*colit)[compi][compj] = 0; - // } - // } - - // printmatrix(std::cout,*A,"prolongation matrix","row",9,1); - } // return the skip flag