def __init__(self):
super(PyMUMPS, self).__init__()
self.utils = mu
# LAMBDAS!!!!
## For simplicity and to follow the documentation, the indices
## are Fotran-like : ICNTL(index) -> self.icntl(index-1)
self.ICNTL = lambda i, v : mu.iarr_setitem(self.icntl, i-1, np.int(v))
self.CNTL = lambda i, v : mu.darr_setitem(self.cntl, i-1, np.double(v))
self.INFO = lambda i : mu.iarr_getitem(self.info, i-1)
self.INFOG = lambda i : mu.iarr_getitem(self.infog, i-1)
self.RINFO = lambda i : mu.darr_getitem(self.rinfo, i-1)
self.RINFOG = lambda i : mu.darr_getitem(self.rinfog, i-1)
self.IRN = lambda v: self.set_int(self.a, v, offset=np.int(1))
self.A = lambda v: self.set_double(self.a, v)
self.RHS = lambda v: self.set_double(self.rhs, v)
def get_double(self, array_in, dim):
"""
Copies C @array_in into a Python array
"""
array_out = np.zeros(dim)
try:
for i in xrange(dim):
array_out[i] = mu.darr_getitem(array_in, i)
except Exception:
#@TODO:put here some stuffs!
pass
else:
return array_out
def get_rinfog(self, index):
mu.darr_getitem(self.rinfog, index-1)
