2010-12-30 18:11:50 -04:00
|
|
|
#!/usr/bin/env python3
|
|
|
|
|
|
|
|
"""
|
|
|
|
A demonstration of classes and their special methods in Python.
|
|
|
|
"""
|
|
|
|
|
2010-07-05 18:44:05 -03:00
|
|
|
class Vec:
|
2010-12-30 18:11:50 -04:00
|
|
|
"""A simple vector class.
|
1992-08-13 09:14:11 -03:00
|
|
|
|
2010-12-30 18:11:50 -04:00
|
|
|
Instances of the Vec class can be constructed from numbers
|
1992-08-13 09:14:11 -03:00
|
|
|
|
2010-07-05 18:44:05 -03:00
|
|
|
>>> a = Vec(1, 2, 3)
|
|
|
|
>>> b = Vec(3, 2, 1)
|
1992-08-13 09:14:11 -03:00
|
|
|
|
2010-07-05 18:44:05 -03:00
|
|
|
added
|
|
|
|
>>> a + b
|
|
|
|
Vec(4, 4, 4)
|
|
|
|
|
|
|
|
subtracted
|
|
|
|
>>> a - b
|
|
|
|
Vec(-2, 0, 2)
|
|
|
|
|
|
|
|
and multiplied by a scalar on the left
|
|
|
|
>>> 3.0 * a
|
|
|
|
Vec(3.0, 6.0, 9.0)
|
1992-08-13 09:14:11 -03:00
|
|
|
|
2010-07-05 18:44:05 -03:00
|
|
|
or on the right
|
|
|
|
>>> a * 3.0
|
|
|
|
Vec(3.0, 6.0, 9.0)
|
2021-03-15 23:53:58 -03:00
|
|
|
|
|
|
|
and dot product
|
|
|
|
>>> a.dot(b)
|
|
|
|
10
|
|
|
|
|
|
|
|
and printed in vector notation
|
|
|
|
>>> print(a)
|
|
|
|
<1 2 3>
|
|
|
|
|
2010-07-05 18:44:05 -03:00
|
|
|
"""
|
2021-03-15 23:53:58 -03:00
|
|
|
|
2003-04-24 14:13:18 -03:00
|
|
|
def __init__(self, *v):
|
2003-09-10 18:12:59 -03:00
|
|
|
self.v = list(v)
|
1992-08-13 09:14:11 -03:00
|
|
|
|
2010-07-05 17:13:41 -03:00
|
|
|
@classmethod
|
|
|
|
def fromlist(cls, v):
|
2003-09-10 18:12:59 -03:00
|
|
|
if not isinstance(v, list):
|
2003-04-24 14:13:18 -03:00
|
|
|
raise TypeError
|
2010-07-05 17:13:41 -03:00
|
|
|
inst = cls()
|
|
|
|
inst.v = v
|
|
|
|
return inst
|
1992-08-13 09:14:11 -03:00
|
|
|
|
2003-04-24 14:13:18 -03:00
|
|
|
def __repr__(self):
|
2021-03-15 23:53:58 -03:00
|
|
|
args = ', '.join([repr(x) for x in self.v])
|
|
|
|
return f'{type(self).__name__}({args})'
|
|
|
|
|
|
|
|
def __str__(self):
|
|
|
|
components = ' '.join([str(x) for x in self.v])
|
|
|
|
return f'<{components}>'
|
1992-08-13 09:14:11 -03:00
|
|
|
|
2003-04-24 14:13:18 -03:00
|
|
|
def __len__(self):
|
|
|
|
return len(self.v)
|
1992-08-13 09:14:11 -03:00
|
|
|
|
2003-04-24 14:13:18 -03:00
|
|
|
def __getitem__(self, i):
|
|
|
|
return self.v[i]
|
1992-08-13 09:14:11 -03:00
|
|
|
|
2003-09-10 18:12:59 -03:00
|
|
|
def __add__(self, other):
|
2021-03-15 23:53:58 -03:00
|
|
|
"Element-wise addition"
|
2010-07-05 18:44:05 -03:00
|
|
|
v = [x + y for x, y in zip(self.v, other.v)]
|
2010-07-05 17:13:41 -03:00
|
|
|
return Vec.fromlist(v)
|
1992-08-13 09:14:11 -03:00
|
|
|
|
2003-09-10 18:12:59 -03:00
|
|
|
def __sub__(self, other):
|
2021-03-15 23:53:58 -03:00
|
|
|
"Element-wise subtraction"
|
2010-07-05 18:44:05 -03:00
|
|
|
v = [x - y for x, y in zip(self.v, other.v)]
|
2010-07-05 17:13:41 -03:00
|
|
|
return Vec.fromlist(v)
|
1992-08-13 09:14:11 -03:00
|
|
|
|
2003-04-24 14:13:18 -03:00
|
|
|
def __mul__(self, scalar):
|
2021-03-15 23:53:58 -03:00
|
|
|
"Multiply by scalar"
|
2010-07-05 18:44:05 -03:00
|
|
|
v = [x * scalar for x in self.v]
|
2010-07-05 17:13:41 -03:00
|
|
|
return Vec.fromlist(v)
|
1992-08-13 09:14:11 -03:00
|
|
|
|
2010-07-05 18:44:05 -03:00
|
|
|
__rmul__ = __mul__
|
1992-08-13 09:14:11 -03:00
|
|
|
|
2021-03-15 23:53:58 -03:00
|
|
|
def dot(self, other):
|
|
|
|
"Vector dot product"
|
|
|
|
if not isinstance(other, Vec):
|
|
|
|
raise TypeError
|
|
|
|
return sum(x_i * y_i for (x_i, y_i) in zip(self, other))
|
|
|
|
|
1992-08-13 09:14:11 -03:00
|
|
|
|
|
|
|
def test():
|
2010-07-05 18:44:05 -03:00
|
|
|
import doctest
|
|
|
|
doctest.testmod()
|
1992-08-13 09:14:11 -03:00
|
|
|
|
|
|
|
test()
|