Paper accepted in Journal of Computational Science. [Link]