vectorFunctions.py

def add(vect, other):
    return tuple(a + b for a, b in zip(vect, other))

def sub(vect, other):
    return tuple(a - b for a, b in zip(vect, other))

def mul(vect, other):
    if hasattr(other, '__iter__'):
        return tuple(a * b for a, b in zip(vect, other))
    return tuple(a * b for a, b in zip(vect, [other] * len(vect)))

def div(vect, other):
    if hasattr(other, '__iter__'):
        return tuple(a / b for a, b in zip(vect, other))
    return tuple(a / b for a, b in zip(vect, [other] * len(vect)))
                 
def length(vect):
    total = sum(map(lambda a: operator.pow(a, 2), vect))
    divisor = math.sqrt(total)
    return divisor

def normalized(vect):
    divisor = [length(vect)] * len(vect)
    return div(vect, divisor)
   
def dot(left, right):
    return  sum(mul(left, right))