Rcpp port of a linear algebra function in R

Question

The following is an objective function for symmetrical non-negative matrix factorization that I\'m trying to port into Rcpp:

fit_H <- function(W,H, num.ite