Feedback for GSoC Proposal

80 views
Skip to first unread message

Sagar Bharadwaj

unread,
Apr 1, 2017, 5:40:38 AM4/1/17
to sympy
I am interested in Implementing multiple coordinate systems for Vectors in SymPy.
I have completed the timeline and here is my proposal. Please Review it and give me you valuable suggestions.


Alan Bromborsky

unread,
Apr 1, 2017, 8:05:58 AM4/1/17
to sy...@googlegroups.com
You proposal is very good.  My only suggestion is that you may want a class for vector differential operators (gradient operator).  If the coordinate system constructor returns the gradient operator you could (in future) overload * and ^ operators so that if f is a scalar field and V a vector field and grad the vector operator then -

gradient of scalar field = grad * f
divergence of vector field = grad * V
curl of vector field = grad ^ V

then the Laplacian would be grad * (grad * f)
--
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+un...@googlegroups.com.
To post to this group, send email to sy...@googlegroups.com.
Visit this group at https://groups.google.com/group/sympy.
To view this discussion on the web visit https://groups.google.com/d/msgid/sympy/b940f17a-d895-4c50-bfd4-70ddd1489c93%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.


Sagar Bharadwaj

unread,
Apr 1, 2017, 9:23:52 AM4/1/17
to sy...@googlegroups.com

Yes. There is a class called Del in sympy.vector.deloperator.py. But in current implementation, it does not return the del operator ( the grad operator). It only calculates gradient, curl and divergence.

Thank you for the suggestion.
Even I was thinking returning a del operator from del class would be better. I'll incorporate that in my proposal.

Thank you for reviewing my proposal.

On 1 Apr 2017 5:35 p.m., "Alan Bromborsky" <abrom...@gmail.com> wrote:
You proposal is very good.  My only suggestion is that you may want a class for vector differential operators (gradient operator).  If the coordinate system constructor returns the gradient operator you could (in future) overload * and ^ operators so that if f is a scalar field and V a vector field and grad the vector operator then -

gradient of scalar field = grad * f
divergence of vector field = grad * V
curl of vector field = grad ^ V

then the Laplacian would be grad * (grad * f)



On 04/01/2017 05:40 AM, Sagar Bharadwaj wrote:
I am interested in Implementing multiple coordinate systems for Vectors in SymPy.
I have completed the timeline and here is my proposal. Please Review it and give me you valuable suggestions.


--
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+unsubscribe@googlegroups.com.

To post to this group, send email to sy...@googlegroups.com.
Visit this group at https://groups.google.com/group/sympy.
To view this discussion on the web visit https://groups.google.com/d/msgid/sympy/b940f17a-d895-4c50-bfd4-70ddd1489c93%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

--
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+unsubscribe@googlegroups.com.

To post to this group, send email to sy...@googlegroups.com.
Visit this group at https://groups.google.com/group/sympy.

Sagar Bharadwaj

unread,
Apr 1, 2017, 11:55:33 AM4/1/17
to sy...@googlegroups.com

Any other suggestions? Anything more that needs to be done on the project other than what I mentioned?

Message has been deleted

Sagar Bharadwaj

unread,
Apr 1, 2017, 2:05:57 PM4/1/17
to sympy
I have made the required changes in both the wiki and the doc.
Please look into the 'Differential Geometry Functions' subsection in 'Implementation' section.
Also the 'June 15 - June 22' period in the Timeline.

szymon.m...@gmail.com

unread,
Apr 2, 2017, 3:49:18 AM4/2/17
to sympy
Hi Sagar,
I found mistake in your proposal. Lame coefficient are defined as you wrote, but you should also take a root from your formula.
See at wikipedia: https://en.wikipedia.org/wiki/Curvilinear_coordinates#Lam.C3.A9_coefficients
They are defined as a norm of vector Diff(r(q1,q2,q3), qi) 

Szymon

Sagar Bharadwaj

unread,
Apr 2, 2017, 4:34:11 AM4/2/17
to sympy
Hi Szymon,
Thank you for reviewing my proposal. 
In the Lame coefficients formula in my proposal, I have raised the whole expression to the power half, which is equivalent to taking it's root. I think you might have missed that because it's in small font.

Was that the mistake you were referring or have I committed a mistake elsewhere?

Sagar

szymon.m...@gmail.com

unread,
Apr 2, 2017, 5:46:17 AM4/2/17
to sympy

OK, I get it.



or have I committed a mistake elsewhere?
No, I just didn't see this power.

Szymon
Reply all
Reply to author
Forward
0 new messages