def testNormSymmetricNotImplemented(self): def f(): return A.norm('1') def g(): return A.norm('inf') A = poisson.poisson2d_sym(self.n) self.failUnlessRaises(NotImplementedError, f) self.failUnlessRaises(NotImplementedError, g)
def setUp(self): import numpy self.n = 30 self.P = poisson.poisson1d(self.n) for i in range(self.n): self.P[i,i] = 4.0 self.A = poisson.poisson2d(self.n) self.S = poisson.poisson2d_sym(self.n) self.I = spmatrix.ll_mat_sym(self.n) for i in range(self.n): self.I[i,i] = -1.0 self.mask = numpy.zeros(self.n**2, 'l') self.mask[self.n/2*self.n:(self.n/2 + 1)*self.n] = 1 self.mask1 = numpy.zeros(self.n**2, 'l') self.mask1[(self.n/2 + 1)*self.n:(self.n/2 + 2)*self.n] = 1
def setUp(self): self.n = 20 self.A = poisson.poisson2d(self.n) self.S = poisson.poisson2d_sym(self.n) self.B = poisson.poisson2d_sym_blk(self.n)