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)
