Novel numerical method for studying the equilibrium crystal shape