Internal: fit a binary MBMM gate by direct marginal optimization

Usage

.mbmm_fit_binary_direct(
  X,
  logL_mat,
  penalty,
  ridge,
  start = NULL,
  maxit = 2000
)