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?
In the current implementation, vector fields are assumed to be smooth functions on their domain. The vector field that you have defined is not (actually, it is not even defined at $r=0$), so I would say this is not a bug.
It would have been great, if it had given the dirac solution.
Asked: 2021-02-07 00:15:43 +0200
Seen: 147 times
Last updated: Feb 07 '21
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.
