Inexact block coordinate descent methods with application to the nonnegative matrix factorization