**0**

# Is there a built-in Python vector class?

Started by Mar 09 2008 12:57 AM

,
6 replies to this topic

###
#2
Crossbones+ - Reputation: **1164**

Posted 09 March 2008 - 01:30 AM

If you want, here's the vector class that I wrote. No guarantees on anything, but it's working great for me. It's only for 2D vectors, but I'm guessing that that will be enough for your work. Note that I also define a class Point to be the same as Vector, which you might not want.

In addition, it seems that you can use NumPy for vector work, but I can't find much good documentation on how to do that. NumPy would require installation by end-users, but it's darn near standard for Python and if you're just playing around, that might not be an issue.

Hope that helps.

import math

class Vector:

'Represents a 2D vector.'

def __init__(self, x = 0, y = 0):

self.x = float(x)

self.y = float(y)

def __add__(self, val):

return Point( self[0] + val[0], self[1] + val[1] )

def __sub__(self,val):

return Point( self[0] - val[0], self[1] - val[1] )

def __iadd__(self, val):

self.x = val[0] + self.x

self.y = val[1] + self.y

return self

def __isub__(self, val):

self.x = self.x - val[0]

self.y = self.y - val[1]

return self

def __div__(self, val):

return Point( self[0] / val, self[1] / val )

def __mul__(self, val):

return Point( self[0] * val, self[1] * val )

def __idiv__(self, val):

self[0] = self[0] / val

self[1] = self[1] / val

return self

def __imul__(self, val):

self[0] = self[0] * val

self[1] = self[1] * val

return self

def __getitem__(self, key):

if( key == 0):

return self.x

elif( key == 1):

return self.y

else:

raise Exception("Invalid key to Point")

def __setitem__(self, key, value):

if( key == 0):

self.x = value

elif( key == 1):

self.y = value

else:

raise Exception("Invalid key to Point")

def __str__(self):

return "(" + str(self.x) + "," + str(self.y) + ")"

Point = Vector

def DistanceSqrd( point1, point2 ):

'Returns the distance between two points squared. Marginally faster than Distance()'

return ( (point1[0]-point2[0])**2 + (point1[1]-point2[1])**2)

def Distance( point1, point2 ):

'Returns the distance between two points'

return math.sqrt( DistanceSqrd(point1,point2) )

def LengthSqrd( vec ):

'Returns the length of a vector sqaured. Faster than Length(), but only marginally'

return vec[0]**2 + vec[1]**2

def Length( vec ):

'Returns the length of a vector'

return math.sqrt( LengthSqrd(vec) )

def Normalize( vec ):

'Returns a new vector that has the same direction as vec, but has a length of one.'

if( vec[0] == 0. and vec[1] == 0. ):

return Vector(0.,0.)

return vec / Length(vec)

def Dot( a,b ):

'Computes the dot product of a and b'

return a[0]*b[0] + a[1]*b[1]

def ProjectOnto( w,v ):

'Projects w onto v.'

return v * Dot(w,v) / LengthSqrd(v)

In addition, it seems that you can use NumPy for vector work, but I can't find much good documentation on how to do that. NumPy would require installation by end-users, but it's darn near standard for Python and if you're just playing around, that might not be an issue.

Hope that helps.

###
#4
Members - Reputation: **454**

Posted 09 March 2008 - 09:29 AM

Sorry -- that could have been a little more clear. I'm looking for a mathematical vector class, not a container. I think I would want to be familiar with pretty much one of the most basic features of the language before I go about writing programs in it :)

###
#5
Members - Reputation: **1682**

Posted 28 March 2008 - 05:16 AM

Quote:

Original post by silverphyre673

Sorry -- that could have been a little more clear. I'm looking for a mathematical vector class, not a container. I think I would want to be familiar with pretty much one of the most basic features of the language before I go about writing programs in it :)

One of the built-in Python types is 'complex', and complex numbers are (2D) points (or vectors).

###
#6
Members - Reputation: **340**

Posted 31 March 2008 - 10:27 PM

In almost all cases which make use of vector math, the array type in the numpy module is appropriate. All of the appropriate functions are implemented in optimized C, and the types are closely compatible with python lists.

Some 'common' operations on numpy arrays that you would expect on a custom 2 or 3 dimensional vector class aren't present (dot product), but are very easy to reproduce.

Some 'common' operations on numpy arrays that you would expect on a custom 2 or 3 dimensional vector class aren't present (dot product), but are very easy to reproduce.