Phase-field models of microstructural pattern formation during alloy solidification are commonly solved numerically using the finite-difference method, which is ideally suited to carry out computationally efficient simulations on massively parallel computer architectures such as Graphic Processing Units. However, one known drawback this method that discretization differential terms involving sp...