I have a function f that compute the inner-product-like value of 2 vectors (the function is symmetrical, so f(x, y) = f(y, x)). It take 2 1d array
f
f(x, y) = f(y, x)