mikaem/fenicstools — pull request #2 by MiroK (2014-03-25 21:35:16+00:00)
Title: Cr weights

This is a small modification to WeightedGradient.py that allows computation of gradient of CR function as matrix-vector product. There's also a demo showing that the computed gradient converges to the correct value as O(h^{1.6}) whereas the convergence rate of the projected gradient is only O(h^{1}).