R: How to implement parallel processing when calculating polychoric correlation matrix?

Question

I am trying to implement parallel processing to speed up the calculation of a polychoric correlation matrix using the hetcor() function from the polycor