make dimension part of the type

[tpapp]

This would allow allocation-free code (using SVector) for a lot of calculations. Of course optionally, as users can still go for Vectors etc.

Caveat: output of collocation_matrix may be too large. Opt-in to static matrices?