How to speed up the construction of the random effects in joint models?

Also, check out the trick here.