Speeding up a hierarchical multinomial logit model

I posted my own question, with all my data and R code, here: Troubleshooting a "MaxDiff" hierarchical multinomial regression I would very much appreciate if you could take a look at my code and see what I may be doing incorrectly!