本文整理汇总了Python中qr_mumps.solver.QRMUMPSSolver类的典型用法代码示例。如果您正苦于以下问题：Python QRMUMPSSolver类的具体用法？Python QRMUMPSSolver怎么用？Python QRMUMPSSolver使用的例子？那么恭喜您, 这里精选的类代码示例或许可以为您提供帮助。
在下文中一共展示了QRMUMPSSolver类的20个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于我们的系统推荐出更棒的Python代码示例。
示例1: test_dense_solve_single_rhs
def test_dense_solve_single_rhs(self):
solver = QRMUMPSSolver((self.m, self.n, self.arow, self.acol, self.aval), verbose=False)
solver.factorize()
e = np.ones(self.n, dtype=np.float64)
rhs = np.dot(self.A, e)
x = solver.solve(rhs)
assert np.allclose(np.dot(self.A,x), rhs, 1e-5)
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:7，代码来源:test_qr_mumps_with_numpy.py
示例2: test_dense_minimum_norm_multiple_rhs
def test_dense_minimum_norm_multiple_rhs(self):
solver = QRMUMPSSolver(self.A, verbose=False)
B = np.ones([self.m, 3], dtype=np.float64)
B[:, 1] = 2 * B[:, 1]
B[:, 2] = 3 * B[:, 2]
with pytest.raises(RuntimeError):
solver = solver.minimum_norm(B)
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:7，代码来源:test_qr_mumps_with_cysparse.py
示例3: test_dense_solve_single_rhs
def test_dense_solve_single_rhs(self):
solver = QRMUMPSSolver(self.A, verbose=False)
solver.factorize()
e = np.ones(self.n, dtype=np.complex64)
rhs = self.A * e
x = solver.solve(rhs)
assert np.allclose(x, e, 1e-5)
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:7，代码来源:test_qr_mumps_with_cysparse.py
示例4: test_dense_least_squares_multiple_rhs
def test_dense_least_squares_multiple_rhs(self):
solver = QRMUMPSSolver((self.m, self.n, self.arow, self.acol, self.aval), verbose=False)
B = np.ones([self.n, 3], dtype=np.complex64)
B[:, 1] = 2 * B[:, 1]
B[:, 2] = 3 * B[:, 2]
rhs = np.dot(self.A, B)
x = solver.least_squares(rhs)
assert np.allclose(x, B, 1e-5)
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:8，代码来源:test_qr_mumps_with_numpy.py
示例5: test_dense_minimum_norm_multiple_rhs
def test_dense_minimum_norm_multiple_rhs(self):
solver = QRMUMPSSolver((self.m, self.n, self.arow, self.acol, self.aval), verbose=False)
B = np.ones([self.n, 3], dtype=np.float64)
B[: ,1] = 2 * B[:,1]
B[: ,2] = 3 * B[:,2]
rhs = np.dot(self.A, B)
x = solver.minimum_norm(rhs)
assert np.allclose(np.dot(self.A,x), rhs, 1e-5)
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:8，代码来源:test_qr_mumps_with_numpy.py
示例6: test_dense_solve_multiple_rhs
def test_dense_solve_multiple_rhs(self):
solver = QRMUMPSSolver((self.m, self.n, self.arow, self.acol, self.aval), verbose=False)
solver.factorize()
B = np.ones([self.n, 3], dtype=np.float64)
B[: ,1] = 2 * B[:,1]
B[: ,2] = 3 * B[:,2]
rhs = np.dot(self.A,B)
x = solver.solve(rhs)
assert np.allclose(x, B, 1e-5)
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:9，代码来源:test_qr_mumps_with_numpy.py
示例7: test_dense_solve_multiple_rhs
def test_dense_solve_multiple_rhs(self):
solver = QRMUMPSSolver(self.A, verbose=False)
solver.factorize()
B = np.ones([self.m, 3], dtype=np.float32)
B[:, 1] = 2 * B[:, 1]
B[:, 2] = 3 * B[:, 2]
x = solver.solve(B)
assert np.allclose(self.A * x[:, 0], B[:, 0], 1e-5)
assert np.allclose(self.A * x[:, 1], B[:, 1], 1e-5)
assert np.allclose(self.A * x[:, 2], B[:, 2], 1e-5)
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:10，代码来源:test_qr_mumps_with_cysparse.py
示例8: test_dense_least_squares_multiple_rhs
def test_dense_least_squares_multiple_rhs(self):
solver = QRMUMPSSolver(self.A, verbose=False)
B = np.ones([self.n, 3], dtype=np.float64)
B[:, 1] = 2 * B[:, 1]
B[:, 2] = 3 * B[:, 2]
rhs = np.ones([self.m, 3], dtype=np.float64)
rhs[:, 0] = self.A * B[:, 0]
rhs[:, 1] = self.A * B[:, 1]
rhs[:, 2] = self.A * B[:, 2]
x = solver.least_squares(rhs)
assert np.allclose(x, B, 1e-5)
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:11，代码来源:test_qr_mumps_with_cysparse.py
示例9: QRMUMPSSolver
from qr_mumps.solver import QRMUMPSSolver
import numpy as np
m = 7
n = 5
arow = np.array([1,2,5,0,5,1,3,4,6,1,2,1,3], dtype=np.int32)
acol = np.array([0,0,0,1,1,2,2,2,2,3,3,4,4], dtype=np.int32)
aval = np.array([0.7,0.6,0.4,0.1,0.1,0.3,0.6,0.7,0.2,0.5,0.2,0.1,0.6], dtype=np.complex64)
solver = QRMUMPSSolver((m, n, arow, acol, aval), verbose=False)
solver.analyze()
solver.factorize()
rhs = np.ones(m, dtype=np.complex64)
print "rhs:"
print rhs
x = solver.solve(rhs)
print "x:"
print x
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:21，代码来源:qr_mumps_numpy_complex.py
示例10: test_dense_minimum_norm_wrong_size_rhs
def test_dense_minimum_norm_wrong_size_rhs(self):
solver = QRMUMPSSolver(self.A, verbose=False)
rhs = np.ones([self.m + 1, 1], dtype=np.float64)
with pytest.raises(ValueError):
x = solver.minimum_norm(rhs)
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:5，代码来源:test_qr_mumps_with_cysparse.py
示例11: test_least_squares
def test_least_squares(self):
solver = QRMUMPSSolver(self.A, verbose=False)
rhs = np.ones([self.m, 1], dtype=np.float64)
with pytest.raises(RuntimeError):
x = solver.least_squares(rhs)
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:5，代码来源:test_qr_mumps_with_cysparse.py
示例12: test_factorize
def test_factorize(self):
solver = QRMUMPSSolver((self.m, self.n, self.arow, self.acol, self.aval), verbose=False)
solver.factorize()
assert solver.analyzed==True
assert solver.factorized==True
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:5，代码来源:test_qr_mumps_with_numpy.py
示例13: test_dense_least_squares_wrong_size_rhs
def test_dense_least_squares_wrong_size_rhs(self):
solver = QRMUMPSSolver(self.A, verbose=False)
rhs = np.ones([self.m + 1, 1], dtype=np.float64)
with pytest.raises(ValueError):
solver = solver.least_squares(rhs)
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:5，代码来源:test_qr_mumps_with_cysparse.py
示例14: test_dense_minimum_norm_wrong_size_rhs
def test_dense_minimum_norm_wrong_size_rhs(self):
solver = QRMUMPSSolver((self.m, self.n, self.arow, self.acol, self.aval), verbose=False)
rhs = np.ones([self.m+1, 1], dtype=np.float32)
with pytest.raises(ValueError):
x = solver.minimum_norm(rhs)
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:5，代码来源:test_qr_mumps_with_numpy.py
示例15: test_factorize
def test_factorize(self):
solver = QRMUMPSSolver(self.A, verbose=False)
solver.factorize()
assert solver.analyzed is True
assert solver.factorized is True
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:5，代码来源:test_qr_mumps_with_cysparse.py
示例16: test_dense_least_squares_wrong_size_rhs
def test_dense_least_squares_wrong_size_rhs(self):
solver = QRMUMPSSolver((self.m, self.n, self.arow, self.acol, self.aval), verbose=False)
rhs = np.ones([self.m+1, 1], dtype=np.complex64)
with pytest.raises(ValueError):
x = solver.least_squares(rhs)
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:5，代码来源:test_qr_mumps_with_numpy.py
示例17: xrange
arow = np.array([0,0,0,1,1,2,2,2,2,3,3,4,4], dtype=np.int32)
acol = np.array([1,2,5,0,5,1,3,4,6,1,2,1,3], dtype=np.int32)
aval = np.array([0.7,0.6,0.4,0.1,0.1,0.3,0.6,0.7,0.2,0.5,0.2,0.1,0.6], dtype=np.float64)
print "A:"
print A
print "-"*80
print "Loop intantiating a solver and solving with one rhs"
print "-"*80
for i in xrange(0,10):
print " "+"-"*76
print " i = %d"%i
print " "+"-"*76
solver = QRMUMPSSolver((m, n, arow, acol, aval), verbose=False)
solver.factorize('scotch')
print solver.factorization_statistics
e = (i+1)*np.ones(n, dtype=np.float64)
rhs = np.dot(A, e)
print " rhs:"
print rhs
x = solver.solve(rhs)
print " x:"
print x
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:28，代码来源:qr_mumps_numpy_loop_multiple_rhs.py
示例18: QRMUMPSSolver
from qr_mumps.solver import QRMUMPSSolver
import numpy as np
m = 4
n = 4
A = np.array([[1, 2, 3, 4], [5, 0, 7, 8], [9, 10, 0, 12], [13, 14, 15, 0]], dtype=np.float64)
arow = np.array([0, 0, 0, 0, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3], dtype=np.int32)
acol = np.array([0, 1, 2, 3, 0, 1, 2, 3, 0, 1, 2, 3, 0, 1, 2, 3], dtype=np.int32)
aval = np.array([1, 2, 3, 4, 5, 0, 7, 8, 9, 10, 0, 12, 13, 14, 15, 0], dtype=np.float64)
print "A:"
print A
solver = QRMUMPSSolver((m, n, arow, acol, aval), verbose=False)
solver.analyze()
solver.factorize()
print "\n\nSolving Ax = b where b is a vector"
e = np.ones(n, dtype=np.float64)
b = np.dot(A, e)
print "b:"
print b
x = solver.solve(b)
np.testing.assert_almost_equal(x, e)
print "x: it should be a vector of ones"
print x
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:30，代码来源:qr_mumps_numpy_double.py
示例19: LLSparseMatrix
from qr_mumps.solver import QRMUMPSSolver
import numpy as np
import sys
A = LLSparseMatrix(mm_filename=sys.argv[1], itype=types.INT32_T, dtype=types.COMPLEX128_T)
(n, m) = A.shape
print "A:"
print A
e = np.ones(n, dtype=np.complex128)
rhs = A * e
solver = QRMUMPSSolver(A, verbose=True)
solver.analyze()
solver.factorize()
x = solver.solve(rhs)
print "x should be 1-column vector:"
print x
print "=" * 80
rhs = np.ones([n, 3], dtype=np.complex128)
rhs[:, 1] = 2 * rhs[:, 1]
rhs[:, 2] = 3 * rhs[:, 2]
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:30，代码来源:qr_mumps_cysparse_complex.py
示例20: test_least_squares
def test_least_squares(self):
solver = QRMUMPSSolver((self.m, self.n, self.arow, self.acol, self.aval), verbose=False)
rhs = np.ones([self.m, 1], dtype=np.float32)
with pytest.raises(RuntimeError):
x = solver.least_squares(rhs)
开发者ID:PythonOptimizers，项目名称:qr_mumps.py，代码行数:5，代码来源:test_qr_mumps_with_numpy.py
注：本文中的qr_mumps.solver.QRMUMPSSolver类示例由纯净天空整理自Github/MSDocs等源码及文档管理平台，相关代码片段筛选自各路编程大神贡献的开源项目，源码版权归原作者所有，传播和使用请参考对应项目的License；未经允许，请勿转载。 |
客服电话
APP下载
官方微信
请发表评论