On 2/10/11 9:53 AM, Dox wrote:
Nicalas... Your suggestion almost work, and in fact it is exactly what I'm talking about!Specifically, my idea is to work with connections with values in a non- Abelian Lie algebra, SU(2), so there are 3 generators. Therefore, the first entry of my function is a form, and the second is a Lie algebra generator!. I tried this class nAform(object): def __init__(self, a, b): self._form = a self._matrix = b def __add__(self, other): if isinstance(other, nAform): if (self._matrix == other._matrix): return nAform(self._form + other._form, self._matrix) else: raise NotImplemented raise NotImplemented def __mul__(self, other): if isinstance(other, nAform): return nAform(self._form.wedge(other._form), self._matrix.commutator(other._matrix)) raise NotImplemented def __repr__(self): return str((self._form, self._matrix)) def __str__(self): return self.__repr__() And for monomials form it can multiply. But I got stacked in the definition of an addition, if the generators of the algebra are different! Which is, of course, one of the most important part of the class. Any additional help is welcome!
I made a mistake in my original code when I did "raise NotImplemented". I believe it should be "return NotImplemented" instead (i.e., NotImplemented is a constant, not an error to be raised). See the last line of the documentation for these functions:
http://docs.python.org/reference/datamodel.html#object.__add__ Thanks, Jason -- To post to this group, send an email to sage-devel@googlegroups.com To unsubscribe from this group, send an email to sage-devel+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-devel URL: http://www.sagemath.org
