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) |--> 0Can 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: 4 years ago
Seen: 218 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.
