the function gradient from sympy.vector import gradient only works on scalar fields, it is missing the gradient of vector. Besides that, the divergence from sympy.vector import divergence only works on vectors, it is missing the divergence of dyadics.
I think we could solve it easily by implementing it manually for cylindrical and spherical coordinates. How can I check if the vector is in cylindrical or spherical coordinates?
the function gradient
from sympy.vector import gradient
only works on scalar fields, it is missing the gradient of vector. Besides that, the divergencefrom sympy.vector import divergence
only works on vectors, it is missing the divergence of dyadics.