1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
|
############################################################################
##
## Copyright (C) 2006-2007 University of Utah. All rights reserved.
##
## This file is part of VisTrails.
##
## This file may be used under the terms of the GNU General Public
## License version 2.0 as published by the Free Software Foundation
## and appearing in the file LICENSE.GPL included in the packaging of
## this file. Please review the following to ensure GNU General Public
## Licensing requirements will be met:
## http://www.opensource.org/licenses/gpl-license.php
##
## If you are unsure which license is appropriate for your use (for
## instance, you are interested in developing a commercial derivative
## of VisTrails), please contact us at contact@vistrails.org.
##
## This file is provided AS IS with NO WARRANTY OF ANY KIND, INCLUDING THE
## WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE.
##
############################################################################
import core.modules
import core.modules.module_registry
from core.modules.vistrails_module import Module, ModuleError
from SciPy import SciPy
from numpy import allclose, arange, eye, linalg, ones
from scipy import linsolve, sparse
#######################################################################
class Matrix(SciPy):
def setSize(self, size):
pass
def setMatrix(self, m):
self.matrix = m
def numElements(self):
return self.matrix.getnnz()
def maxNumElements(self):
return self.matrix.nzmax
def rows(self):
return self.matrix.shape[0]
def cols(self):
return self.matrix.shape[1]
def Reals(self):
out = SparseMatrix()
tmp = self.matrix.copy()
out.matrix = tmp._real()
return out
def Imaginaries(self):
out = SparseMatrix()
tmp = self.matrix.copy()
out.matrix = tmp._imag()
return out
def Conjugate(self):
out = SparseMatrix()
out.matrix = self.matrix.conjugate().copy()
return out
def GetRow(self, i):
return self.matrix.getrow(i)
def GetCol(self, i):
return self.matrix.getcol(i)
class SparseMatrix(Matrix):
def setSize(self, size):
self.matrix = sparse.csc_matrix((size, size))
self.matrix.setdiag(ones(size))
class DenseMatrix(Matrix):
def setSize(self, size):
self.matrix = sparse.csc_matrix((size, size))
self.matrix.setdiag(ones(size))
self.matrix.todense()
class DOKMatrix(Matrix):
def setSize(self, size):
self.matrix = sparse.dok_matrix((size, size))
self.matrix.setdiag(ones(size))
class COOMatrix(Matrix):
def __init__(self, mat):
self.matrix=mat
def setSize(self, size):
self.matrix = sparse.coo_matrix((size, size))
self.matrix.setdiag(ones(size))
class CSRMatrix(Matrix):
def __init__(self, mat):
self.matrix=mat
def setSize(self, size):
self.matrix = sparse.csr_matrix((size, size))
self.matrix.setdiag(ones(size))
class LILMatrix(Matrix):
def __init__(self, mat):
self.matrix=mat
def setSize(self, size):
self.matrix = sparse.lil_matrix((size, size))
self.matrix.setdiag(ones(size))
#######################################################################
|
