Analysis of inhomogeneously filled waveguide devices by a hybrid mode matching-finite element technique