Usually my mixed models contain several categorical variables with a lot of unique levels, so X matrix is very sparse. I use glmmTMB packag
X
glmmTMB