Add sorted result

[fbaeuerlein]

Add sort of result via tag
e.g. solve<>(matrix, sorted_t);