please include the computation of the objective function (for convergence
checks or at least print outs). Thanks.

while( iter < max_iterations ){
        iter = iter + 1;
        H = (H * (t(W) %*% (V/(W%*%H))))/t(colSums(W));
        W = (W * ((V/(W%*%H)) %*% t(H)))/t(rowSums(H));
        obj = as.scalar(colSums(W)%*%rowSums(H)) -sum(V * log(W%*%H));
        print("ITER=" + iter + " obj=" + obj);
}

Regards,
Matthias



From:   Frederick R Reiss/Almaden/IBM@IBMUS
To:     [email protected]
Date:   12/17/2015 01:41 AM
Subject:        Re: DML example on main SystemML website



We can use the Poisson nonnegative matrix factorization example from last
week's webcast:
      i = 0
      while(i < max_iterations) {
          H = (H * (t(W) %*% (V/(W%*%H + epsilon))))/ t(colSums(W))
          W = (W * ((V/(W%*%H) + epsilon) %*% t(H)))/ t(rowSums(H))
          i = i + 1;
      }



Sound ok to everyone?

Fred


Deron Eriksson ---12/16/2015 04:02:26 PM---Hi, I think the main SystemML
website at http://systemml.incubator.apache.org/

From: Deron Eriksson <[email protected]>
To: [email protected]
Date: 12/16/2015 04:02 PM
Subject: DML example on main SystemML website



Hi,

I think the main SystemML website at http://systemml.incubator.apache.org/
needs to be updated so that the DML example is an actual algorithm or at
least a fragment of an algorithm.

Does anyone have a recommendation for a short, concise example that shows
the power of DML?

Thanks!
Deron




Reply via email to