The code: ``` penalty <- 5 udr::ud_fit(fit0, control = list(lambda = penalty, penalty.type = "nu", unconstrained.update = "ed")) ``` should throw a warning / error, since the nuclear norm penalty is not implemented for the extreme deconvolution algorithm