Add bugfixes for special cases in sets.py

[jylambert]

In sets.py some lines are failling in some special cases.

For example in line 44 -> consumers = np.where(matrices["a_c"].sum(axis=1) == 1)[0] fails if there is only one heat sink, as the vector has only one dimension and not two.