From 395bc7cb5ce0ffc49a3a316d847830b1dc61085f Mon Sep 17 00:00:00 2001 From: KatolaZ Date: Mon, 11 May 2015 23:25:18 +0100 Subject: Removed the last full matrix dangling around --- python/multired.py | 20 +++++++++++++------- 1 file changed, 13 insertions(+), 7 deletions(-) (limited to 'python') diff --git a/python/multired.py b/python/multired.py index 3952eec..0023d38 100644 --- a/python/multired.py +++ b/python/multired.py @@ -35,7 +35,7 @@ # -------------------------------------------- # # -- 2015/04/23 -- release 0.1 -# +# -- 2015/05/11 -- release 0.1.1 -- removed the last full matrices # @@ -123,8 +123,10 @@ class layer: elif matrix != None: self.adj_matr = copy.copy(matrix) self.N, _x = matrix.shape - K = np.multiply(self.adj_matr.sum(0), np.ones((self.N,self.N))) - D = np.diag(np.diag(K)) + #K = np.multiply(self.adj_matr.sum(0), np.ones((self.N,self.N))) + #D = np.diag(np.diag(K)) + K = self.adj_matr.sum(0) + D = K * eye(self.N) self.laplacian = csr_matrix(D - self.adj_matr) K = self.laplacian.diagonal().sum() self.resc_laplacian = csr_matrix(self.laplacian / K) @@ -135,8 +137,10 @@ class layer: self.N = N self.adj_matr = csr_matrix((self._ww, (self._ii, self._jj)), shape=(self.N, self.N)) self.adj_matr = self.adj_matr + self.adj_matr.transpose() - K = np.multiply(self.adj_matr.sum(0), np.ones((self.N,self.N))) - D = np.diag(np.diag(K)) + #K = np.multiply(self.adj_matr.sum(0), np.ones((self.N,self.N))) + #D = np.diag(np.diag(K)) + K = self.adj_matr.sum(0) + D = K * eye(self.N) self.laplacian = csr_matrix(D - self.adj_matr) K = self.laplacian.diagonal().sum() self.resc_laplacian = csr_matrix(self.laplacian / K) @@ -171,8 +175,10 @@ class layer: self.adj_matr = self.adj_matr + other_layer.adj_matr else: self.adj_matr = copy.copy(other_layer.adj_matr) - K = np.multiply(self.adj_matr.sum(0), np.ones((self.N,self.N))) - D = np.diag(np.diag(K)) + #K = np.multiply(self.adj_matr.sum(0), np.ones((self.N,self.N))) + #D = np.diag(np.diag(K)) + K = self.adj_matr.sum(0) + D = K * eye(self.N) self.laplacian = csr_matrix(D - self.adj_matr) K = self.laplacian.diagonal().sum() self.resc_laplacian = csr_matrix(self.laplacian / K) -- cgit v1.2.3