Numerical glacier and ice-sheet models have been widely used in the context of climate change and landform evolution. While simulations of ice flow were numerically expensive for a long time, their performance has recently been boosted to an unprecedented level by machine learning techniques. This paper aims at keeping classical numerics competitive by introducing a novel numerical scheme, which allows for simulations at spatial resolutions of 25 m or even finer on standard desktop PCs.