Could you please either show the code that gave that result or else add it 
as a test for #20047

/c

On Sunday, October 31, 2021 at 7:42:13 PM UTC-5 Sean Logan wrote:

> 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 [email protected].
To view this discussion on the web visit 
https://groups.google.com/d/msgid/sympy/3060a900-36ed-4a45-9514-e90a72853686n%40googlegroups.com.

Reply via email to