Divergence of inverse square
Hello,
I tried to calculate the divergence of 1/r^2e_r, but I got zero, where it should have been 4pidirac_delta(re_r)
E.<r,th,ph>=EuclideanSpace(coordinates="spherical")
v=E.vector_field([1/r^2,0,0])
v.div().display()
Results is
E^3 --> R
(r, th, ph) |--> 0
Can this be considered as a bug?