Currently risk_model <= limit compiles to something like cp.sum_squares(...) <= limit**2. Instead we want cp.norm2(...) <= limit to improve execution performance and solver stability, especially with first order solvers. We'll need to add an alternative compilation route and use cp.hstack in the constraint object to concatenate.