Alexander,
I guess it should be:Number q = values[0] * p.values[0];
You are right, thanks for pointing this out. It looks like not many people are using dim>3 in tensor - otherwise a bug that fundamental would not have been undetected for half a year... :-) I fixed it in the development version and made sure we test it, too.
Best, Martin _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
