expose the functionality of automatic interface creation

[exaexa]

The current one is kinda hidden in flux_balance_constraints and not really reusable, despite potentially useful for custom models.

Also some of the semantics that people get from interfaces could get re-used in filtering/grouping the reactions sensibly, as requested in #78 .

Ideally if we expose the functionality, we should also somehow provide the universal reaction filter based on that.