# [Scipy-svn] r4412 - trunk/scipy/cluster

scipy-svn@scip... scipy-svn@scip...
Wed Jun 4 03:08:32 CDT 2008

```Author: damian.eads
Date: 2008-06-04 03:08:30 -0500 (Wed, 04 Jun 2008)
New Revision: 4412

Modified:
trunk/scipy/cluster/hierarchy.py
Log:
Fixed issue with cosine.

Modified: trunk/scipy/cluster/hierarchy.py
===================================================================
--- trunk/scipy/cluster/hierarchy.py	2008-06-03 17:53:23 UTC (rev 4411)
+++ trunk/scipy/cluster/hierarchy.py	2008-06-04 08:08:30 UTC (rev 4412)
@@ -1463,18 +1463,18 @@
# Find out: Is there a dot subtraction operator so I can
# subtract matrices in a similar way to multiplying them?
# Need to get rid of as much unnecessary C code as possible.
-        elif mstr in set(['cosine_old', 'cos_old']):
+        elif mstr in set(['cosine', 'cos']):
norms = np.sqrt(np.sum(X * X, axis=1))
_hierarchy_wrap.pdist_cosine_wrap(_convert_to_double(X), dm, norms)
-        elif mstr in set(['cosine', 'cos']):
+        elif mstr in set(['old_cosine', 'old_cos']):
norms = np.sqrt(np.sum(X * X, axis=1))
nV = norms.reshape(m, 1)
# The numerator u * v
nm = np.dot(X, X.T)
# The denom. ||u||*||v||
de = np.dot(nV, nV.T);
-            dm = 1 - (nm / de)
-            dm[xrange(0,m),xrange(0,m)] = 0
+            dm = 1.0 - (nm / de)
+            dm[xrange(0,m),xrange(0,m)] = 0.0
dm = squareform(dm)
elif mstr in set(['correlation', 'co']):
X2 = X - X.mean(1)[:,np.newaxis]

```