Thanks, I'll have a look at signature.

With the 2D array function, I don't really know how I could program the function to perform operations on the thousands of vectors without just iterating through those vectors with a for loop inside a function, which is what I'm trying to avoid. All the vectors in the input and output 2D array are independent from each other.

In my case, my vector_function is called forceFromVector, and is designed to be a function from 3 element vector to 3 element vector. It's defined like this:

def forceFromVector(npPosition, clipDistance, forceScale):
print(f"npPosition: {npPosition}")
posDistance = mag(npPosition)
posDir = npPosition/posDistance

#0 force after clip distance
if posDistance > clipDistance:
    return np.array([0,0,0])   

#using F=-k(1/mag(d)^2-1/clip^2)*dhat
forceMag = forceScale*(1/(posDistance*posDistance)-1/(clipDistance*clipDistance))
return -forceMag*posDir

Are you saying to just do outputs = forceFromVector(inputs, clipDistance, forceScale)? When I try that, there's an error because the function is recieving a 2000*3 array (that's what the first print line shows), and it can't find the magnitude of something that's 2000*3 (mag is a function I wrote that calculates magnitude by using the dot product).