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)