Algorithm assumed that integration_element is the same for all

quadrature points! Fixed.

[[Imported from SVN: r1060]]

fem/l2projection.hh

@@ -44,7 +44,6 @@ namespace Dune
       for( ; it != endit ; ++it)
       {
         discFunc.localFunction( *it , lf );
-        double det = (*it).geometry().integration_element(quad.point(0));
 
         const typename FunctionSpaceType::BaseFunctionSetType & set =
           functionSpace_.getBaseFunctionSet(*it);
@@ -53,6 +52,7 @@ namespace Dune
         {
           for(int qP = 0; qP < quad.nop(); qP++)
           {
+            double det = (*it).geometry().integration_element(quad.point(qP));
             f.evaluate((*it).geometry().global( quad.point(qP) ), ret);
             set.eval(i,quad,qP,phi);
             lf[i] += det * quad.weight(qP) * (ret * phi);