* Now this one is a little tricky. The problem is that we
  basically need a specialization for every combination
  of dimensions (since some functions may only be
  implemented for some dimensions).
* We now generate specialized declarations for all
  dimensions in the header.
* We also generate specialized definitions for all
  dimensions in the source and we require the
  implementer of a new functions to add the function in
  the source file.
* Adding a new function is a sub optimal task, since it
  requires some copying. This solution may not be the
  most elegant but it seemed to be *fastest*, since we
  can just copy/paste find/replace.
* Ideas and patches are welcome!