powered by
vp_cp-style PL objective (includes constants)
glmmfe_pl_objective( z_work, w_num, X, Z, beta, eta, tau2, inner, approx, vc_eps = 1e-12 )