[WIP] Strong rules for coordinate descent

please use

active_rule = np.abs(active_rule, active_rule)