37 : gt(gt_), lb(lb_), Minv(0)
38 , qr(QuadratureRules<DF,dimD>::rule(gt, 2*lb.order()))
42 DUNE_THROW(Exception,
"size template parameter does not match size of "
45 const QRiterator qrend = qr.end();
46 for(QRiterator qrit = qr.begin(); qrit != qrend; ++qrit) {
48 lb.evaluateFunction(qrit->position(),base);
50 for(
unsigned int i = 0; i < size; ++i)
51 for(
unsigned int j = 0; j < size; ++j)
52 Minv[i][j] += qrit->weight() * base[i] * base[j];
70 const QRiterator qrend = qr.end();
71 for(QRiterator qrit = qr.begin(); qrit != qrend; ++qrit) {
73 R y = f(qrit->position());
76 lb.evaluateFunction(qrit->position(),base);
78 for(
unsigned int i = 0; i < size; ++i)
79 for(
unsigned int j = 0; j < size; ++j)
80 out[i] += Minv[i][j] * qrit->weight() * y * base[j];