Ask Your Question
0

Divergence of inverse square

asked 2021-02-07 00:15:43 +0100

curios_mind gravatar image

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?

edit retag flag offensive close merge delete

1 Answer

Sort by ยป oldest newest most voted
0

answered 2021-02-07 14:24:51 +0100

eric_g gravatar image

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.

edit flag offensive delete link more

Comments

It would have been great, if it had given the dirac solution.

curios_mind gravatar imagecurios_mind ( 2021-02-08 14:54:22 +0100 )edit

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

Stats

Asked: 2021-02-07 00:15:43 +0100

Seen: 187 times

Last updated: Feb 07 '21