Python QRMUMPSSolver.analyze Beispiele

Beispiel #1

Datei: test_qr_mumps_with_cysparse.py Projekt: PythonOptimizers/qr_mumps.py

 def test_analyze(self):
     solver = QRMUMPSSolver(self.A, verbose=False)
     solver.analyze()
     assert solver.analyzed is True

Beispiel #2

Datei: qr_mumps_cysparse_complex.py Projekt: PythonOptimizers/qr_mumps.py

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]

x = solver.solve(rhs)
print "x:"

Beispiel #3

Datei: test_qr_mumps_with_numpy.py Projekt: PythonOptimizers/qr_mumps.py

 def test_analyze(self):
     solver = QRMUMPSSolver((self.m, self.n, self.arow, self.acol, self.aval), verbose=False)
     solver.analyze()
     assert solver.analyzed==True