def testcreation2(self): a=g.symbol("a") b=g.symbol("b") x=g.symbol("x") y=g.symbol("y") m=g.matrix([[a,0],[0,b]]) self.assertEqual(repr(m),"[[a,0],[0,b]]") self.assertNotEqual(repr(m),"[[a,1],[0,b]]") m=g.matrix([[0,0],[0,0]]) self.assertEqual(repr(m),"[[0,0],[0,0]]") m[0,0]=a self.assertEqual(repr(m),"[[a,0],[0,0]]") m[1,1]=b self.assertEqual(repr(m),"[[a,0],[0,b]]") self.assertNotEqual(repr(m),"[[a,0],[b,0]]") self.assertEqual(m[0,0],a) self.assertNotEqual(m[0,0],b) self.assertEqual(m[0,1],0) self.assertNotEqual(m[0,1],b) self.assertEqual(m[1,1],b) self.assertNotEqual(m[1,1],0) M=g.matrix m[0,1]=a self.assertEqual(m,M([[a,a],[0,b]])) m[0,1]=b self.assertEqual(m,M([[a,b],[0,b]])) m[0,1]=b*a self.assertEqual(m,M([[a,a*b],[0,b]]))
def __init2__(self, l, symbs=None, time=None): self.i = len(l) self.j = len(l[0]) self.data = _g.matrix(l) self.spatial_symbs = symbs self.time = time if not symbs: self.symbs = range(self.j) self._lhs = []
def __init2__(self, l, symbs=None, time=None): self.i = len(l) self.j = len(l[0]) l = [[_toex(x) for x in sl] for sl in l] #Convert to nested list of swiginac types self.data = _g.matrix(l) self.spatial_symbs = symbs self.time = time if not symbs: self.symbs = range(self.j) self._lhs = []
def __init1__(self, i, j, symbs=None, time=None): self.i = i; self.j = j self.data = _g.matrix(i, j) self.spatial_symbs = symbs self.time = time if not symbs: self.symbs = range(self.j) for i in range(self.i): for j in range(self.j): self[i, j] = Expr() self._lhs = []