Right now, a single reference category constraint is routed to the augmented lagrangian algorithm in fit_null(). However, this is a simpler case for estimation under the null hypothesis, and likely a faster algorithm could be constructed and implemented for this.