Internal: fit a multinomial MBMM gate by direct marginal optimization

Usage

.mbmm_fit_multinom_direct(
  X,
  logL_mat,
  K,
  penalty,
  ridge,
  start = NULL,
  maxit = 4000
)