Revision history [back]

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?