ENH: Make irreducible parameters as local as possible
For pruning purposes and to keep the physical meaning of irreducible parameters it would be a good idea to make the irreducible parameters as "local" as possible.
First approach to this, as suggested by @freeriks , would be to run QR factorization on constraint_vectors. This works well for second order, however fails to fulfill our goal for higher orders. The desired outcome is:
The number of irreducbile parameters which only connects to one orbit should be maximized
Edited by Erik Fransson