¿Existe alguna función que pueda utilizarse para calcular la divergencia del campo vectorial? (enMatlab) Espero que exista en numpy / scipy pero no puedo encontrarlo usando Google.

Necesito calculardiv[A * grad(F)], dónde

F = np.array([[1,2,3,4],[5,6,7,8]]) # (2D numpy ndarray)

A = np.array([[1,2,3,4],[1,2,3,4]]) # (2D numpy ndarray)

asi quegrad(F) es una lista de 2Dndarrays

Sé que puedo calcular divergencias comoesta Pero no quiero reinventar la rueda. (También esperaría algo más optimizado) ¿Alguien tiene sugerencias?