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 Frederick R Reiss---12/17/2015 01:41:00 AM---We can use the Poisson nonnegative matrix factorization example from last week's webcast:

From: Frederick R Reiss/Almaden/IBM@IBMUS
To: dev@systemml.incubator.apache.org
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 <deroneriksson@gmail.com>
To:
dev@systemml.incubator.apache.org
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