Overloading for functions with different dimensions

[Scheiermann]

Background

Many functions are explicitly defined for 1D, 2D or 3D.
Thus more functions are created than needed, as they are basically the same (just changing dimension/parameters list)

Ideas how to solve

Use polymorphisms (overloading).