A penalty function for computing orthogonal non-negative matrix factorization