Numerical computation of convex risk measures