I have noticed that there are some functions related to model order reduction in matrix_comps.jl but they are not listed anywhere. Namely balreal() and baltrunc(). I only found them while going through the code. Is there any reason for hiding them (perhaps users are not encouraged to use them right now)? If not, it would be a pitty not to list them. I might even contribute some examples of usage.