Using the hook formula, it is easy to compute the volume of order polytopes obtained from posets with partition shape, since this is the same as the number of linear extensions.
To my knowledge, computing the number of linear extensions of an arbitrary order polytope is #P-hard.
There is also the family of series-parallel posets, where computing the Ehrhart polynomial (and hence the volume) can be done in polynomial time, so there are natural families where computing the Ehrhart polynomial is easy.
My question is then: can one compute the Ehrhart polynomial of partition-shaped order polytopes quickly (and not just the volume)?