Hello, I discovered that the divergence() function from sympy.vector gives incorrect answers sometimes when using Cylindrical coordinates. For example:
divergence( 1 * rhat ) # returns 0 when it should give 1/R I found that this issue was raised last year, and a patch offered, but it was never integrated into the sympy release. https://github.com/sympy/sympy/pull/20047 I went into sympy/vector/operators.py and commented out lines 366-368 by hand. That seems to have fixed the problem. I also found this library, called "symfields" which does div, grad, curl properly in cylindrical, spherical, or any curvilinear coordinate system https://github.com/DocNan/SymFields it does not use CoordSys3D. -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving emails from it, send an email to sympy+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/sympy/fd4f32b4-8d1f-4873-9b25-3ff2177ace12n%40googlegroups.com.