[ 
https://issues.apache.org/jira/browse/MAHOUT-846?page=com.atlassian.jira.plugin.system.issuetabpanels:comment-tabpanel&focusedCommentId=13175277#comment-13175277
 ] 

Hudson commented on MAHOUT-846:
-------------------------------

Integrated in Mahout-Quality #1266 (See 
[https://builds.apache.org/job/Mahout-Quality/1266/])
    MAHOUT-846: Cache pdf zProd2piR term constant over life of cluster

jeastman : 
http://svn.apache.org/viewcvs.cgi/?root=Apache-SVN&view=rev&rev=1222524
Files : 
* 
/mahout/trunk/core/src/main/java/org/apache/mahout/clustering/dirichlet/models/GaussianCluster.java

                
> Improve Scalability of Gaussian Cluster For Wide Vectors
> --------------------------------------------------------
>
>                 Key: MAHOUT-846
>                 URL: https://issues.apache.org/jira/browse/MAHOUT-846
>             Project: Mahout
>          Issue Type: Improvement
>    Affects Versions: 0.5
>            Reporter: Jeff Eastman
>            Assignee: Jeff Eastman
>             Fix For: 0.6
>
>
> The pdf() implementation in GaussianCluster is pretty lame. It is computing a 
> running product of the element pdfs which, for wide input vectors (Reuters is 
> 41,807), always underflows and returns 0. Here's the code:
> {noformat}
>   public double pdf(VectorWritable vw) {
>     Vector x = vw.get();
>     // return the product of the component pdfs
>     // TODO: is this reasonable? correct? It seems to work in some cases.
>     double pdf = 1;
>     for (int i = 0; i < x.size(); i++) {
>       // small prior on stdDev to avoid numeric instability when stdDev==0
>       pdf *= UncommonDistributions.dNorm(x.getQuick(i),
>           getCenter().getQuick(i), getRadius().getQuick(i) + 0.000001);
>     }
>     return pdf;
> {noformat}
>   }

--
This message is automatically generated by JIRA.
If you think it was sent incorrectly, please contact your JIRA administrators: 
https://issues.apache.org/jira/secure/ContactAdministrators!default.jspa
For more information on JIRA, see: http://www.atlassian.com/software/jira

        

Reply via email to