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}).