Loewner-based macromodeling with exact interpolation constraints