We will hopefully soon get rid of such annoying #defines via our linear algebra interface, which will do all of this in a central place and also compute eigendecompositions using multiple CPUs or GPUs :)