Avoid c-instruction in basis gradient evaluation
Has to be handled with care. The attempt in d0a87348 failed because the jacobian inverse transposed is a DiagonalMatrix for YaspGrid and you don't get the correct entries with [] operator. If you wan't to do it the right way you have to detect the cases where dune uses a DiagonalMatrix and collapse your loop accordingly.