Changed basis function __call__ to take xyz

This allows passing in a numpy array to save time.
For now, I just changed all of the doctests to pass
in a tuple(x,y,z) rather than x,y,z as *args.
Passing in an xyz object instead of the values to bf eval
Expanding test cases to try to get bf to eval on multiple points at once.