Done with a single list comprehension:
def cross3d(a, b): # 3d cross product return [a[i-2]*b[i-1]-a[i-1]*b[i-2] for i in range(3)]