I like what you have written so far. You should think about how the DiracDelta and Heaviside functions that are already in SymPy should be used. Think about the different methods and attributes that the singularily function would need. I'd recommend showing an example of how you will use this to solve a specific beam problem. You can use an example to show off the types of methods the objects will have.
Keep up the good work. Jason moorepants.info +01 530-601-9791 On Tue, Mar 8, 2016 at 9:51 AM, SAMPAD SAHA <[email protected]> wrote: > Hello Jason, > > I have drafted some portions of my proposal here > <https://github.com/sampadsaha5/sympy/wiki/GSoC-2016-Application:-Singularity-Functions> > . > > I haven't completed it yet. I am on my first phase of preparing the > proposal. After completing "The Plan and Prototype" section I will Draft my > proposal on Sympy's wiki. > > Can you please review it ? > > I need your suggestion. > > > > -------------------- > With Regards, > > Sampad > > > > > > > > > Regards > Sampad Kumar Saha > Mathematics and Computing > I.I.T. Kharagpur > > On Fri, Mar 4, 2016 at 4:31 AM, SAMPAD SAHA <[email protected]> wrote: > >> Jason, >> >> Yah, That would be great. After drafting my proposal on wiki, I will let >> you know. >> >> >> >> Regards >> Sampad Kumar Saha >> Mathematics and Computing >> I.I.T. Kharagpur >> >> On Fri, Mar 4, 2016 at 4:31 AM, SAMPAD SAHA <[email protected]> >> wrote: >> >>> >>> Yah, That would be great. After drafting my proposal on wiki, I will let >>> you know. >>> >>> >>> >>> >>> Regards >>> Sampad Kumar Saha >>> Mathematics and Computing >>> I.I.T. Kharagpur >>> >>> On Fri, Mar 4, 2016 at 4:12 AM, Jason Moore <[email protected]> >>> wrote: >>> >>>> Seems reasonable. I recommend drafting up a python session like you >>>> have there showing the desired functionality. You can make this part of >>>> your proposal on the wiki. Once you have a reasonable amount of material we >>>> can discuss further. >>>> >>>> >>>> Jason >>>> moorepants.info >>>> +01 530-601-9791 >>>> >>>> On Thu, Mar 3, 2016 at 2:28 PM, SAMPAD SAHA <[email protected]> >>>> wrote: >>>> >>>>> Jason, >>>>> >>>>> Thanks. I will go through the printing packages. >>>>> >>>>> Will it be ok if we add the module for singularity function in >>>>> sympy/functions/special ? >>>>> >>>>> And I am thinking to do like this: >>>>> >>>>> There will be class which will be taking power and offset as argument >>>>> e.g. f = singularityFunc(power , offset , variable) and it would be >>>>> printed as <variable - offset> raised to "power" . >>>>> >>>>> Like this: >>>>> >>>>> In [1] : singularityFunc(n , a , x) >>>>> Out [1] : <x-a>^n >>>>> >>>>> >>>>> >>>>> Regards >>>>> Sampad Kumar Saha >>>>> Mathematics and Computing >>>>> I.I.T. Kharagpur >>>>> >>>>> On Thu, Mar 3, 2016 at 11:01 PM, Jason Moore <[email protected]> >>>>> wrote: >>>>> >>>>>> After you create some objects that represent the functions you will >>>>>> need to add functionality to the various printers that we have. See the >>>>>> sympy printing packages for examples. >>>>>> >>>>>> >>>>>> Jason >>>>>> moorepants.info >>>>>> +01 530-601-9791 >>>>>> >>>>>> On Thu, Mar 3, 2016 at 8:55 AM, SAMPAD SAHA <[email protected]> >>>>>> wrote: >>>>>> >>>>>>> Jason , >>>>>>> >>>>>>> >>>>>>> Suppose I want to print <x-2>**2 in isympy as <x-2> raised to 2 . >>>>>>> How we can do it? >>>>>>> >>>>>>> >>>>>>> >>>>>>> >>>>>>> >>>>>>> Regards >>>>>>> Sampad Kumar Saha >>>>>>> Mathematics and Computing >>>>>>> I.I.T. Kharagpur >>>>>>> >>>>>>> On Thu, Mar 3, 2016 at 1:29 AM, SAMPAD SAHA <[email protected]> >>>>>>> wrote: >>>>>>> >>>>>>>> Jason , >>>>>>>> >>>>>>>> Sorry, actually I didn't notice that I haven't added the list in >>>>>>>> cc. It was mistaken. >>>>>>>> >>>>>>>> I will let you know about my status of my proposal soon. >>>>>>>> >>>>>>>> Thank You. >>>>>>>> >>>>>>>> >>>>>>>> -------------------- >>>>>>>> Sampad >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> Regards >>>>>>>> Sampad Kumar Saha >>>>>>>> Mathematics and Computing >>>>>>>> I.I.T. Kharagpur >>>>>>>> >>>>>>>> On Thu, Mar 3, 2016 at 1:16 AM, Jason Moore <[email protected]> >>>>>>>> wrote: >>>>>>>> >>>>>>>>> Sampad, please try to keep conversations on list so that everyone >>>>>>>>> can benefit. >>>>>>>>> >>>>>>>>> Please draft your proposal on the SymPy wiki. >>>>>>>>> >>>>>>>>> >>>>>>>>> Jason >>>>>>>>> moorepants.info >>>>>>>>> +01 530-601-9791 >>>>>>>>> >>>>>>>>> On Wed, Mar 2, 2016 at 11:44 AM, SAMPAD SAHA < >>>>>>>>> [email protected]> wrote: >>>>>>>>> >>>>>>>>>> Jason, >>>>>>>>>> >>>>>>>>>> I have started working on my proposal. After adding some ideas on >>>>>>>>>> how to proceed with this project from myside , I will send you the >>>>>>>>>> link. >>>>>>>>>> >>>>>>>>>> Can you please suggest me , where should I write my proposal so >>>>>>>>>> that I can send it to you for reviewing? >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> Regards >>>>>>>>>> Sampad Kumar Saha >>>>>>>>>> Mathematics and Computing >>>>>>>>>> I.I.T. Kharagpur >>>>>>>>>> >>>>>>>>>> On Thu, Mar 3, 2016 at 1:05 AM, SAMPAD SAHA < >>>>>>>>>> [email protected]> wrote: >>>>>>>>>> >>>>>>>>>>> Jason, >>>>>>>>>>> >>>>>>>>>>> Thank You. I would try my level best to follow your suggestions. >>>>>>>>>>> I would try to send more patches to make my proposal more better. >>>>>>>>>>> >>>>>>>>>>> >>>>>>>>>>> >>>>>>>>>>> >>>>>>>>>>> Regards >>>>>>>>>>> Sampad Kumar Saha >>>>>>>>>>> Mathematics and Computing >>>>>>>>>>> I.I.T. Kharagpur >>>>>>>>>>> >>>>>>>>>>> On Thu, Mar 3, 2016 at 12:44 AM, Jason Moore < >>>>>>>>>>> [email protected]> wrote: >>>>>>>>>>> >>>>>>>>>>>> Sampad, >>>>>>>>>>>> >>>>>>>>>>>> Only one patch is required along with the proposal. But if you >>>>>>>>>>>> have a look at the winning proposals from previous years, I think >>>>>>>>>>>> you'll >>>>>>>>>>>> see that people with more patches and higher engagement with the >>>>>>>>>>>> community >>>>>>>>>>>> are much more likely to get accepted. I'd recommend doing as much >>>>>>>>>>>> as you >>>>>>>>>>>> can to make your proposal stand out from the rest. >>>>>>>>>>>> >>>>>>>>>>>> Jason >>>>>>>>>>>> >>>>>>>>>>>> >>>>>>>>>>>> Jason >>>>>>>>>>>> moorepants.info >>>>>>>>>>>> +01 530-601-9791 >>>>>>>>>>>> >>>>>>>>>>>> On Wed, Mar 2, 2016 at 11:06 AM, SAMPAD SAHA < >>>>>>>>>>>> [email protected]> wrote: >>>>>>>>>>>> >>>>>>>>>>>>> Jason, >>>>>>>>>>>>> >>>>>>>>>>>>> I had fixed an issue earlier . This >>>>>>>>>>>>> <https://github.com/sympy/sympy/pull/10225> was my PR. >>>>>>>>>>>>> Do I need to send any more patch file? >>>>>>>>>>>>> >>>>>>>>>>>>> ------------------- >>>>>>>>>>>>> With Regards >>>>>>>>>>>>> Sampad >>>>>>>>>>>>> >>>>>>>>>>>>> >>>>>>>>>>>>> Regards >>>>>>>>>>>>> Sampad Kumar Saha >>>>>>>>>>>>> Mathematics and Computing >>>>>>>>>>>>> I.I.T. Kharagpur >>>>>>>>>>>>> >>>>>>>>>>>>> On Sun, Feb 28, 2016 at 9:59 PM, Jason Moore < >>>>>>>>>>>>> [email protected]> wrote: >>>>>>>>>>>>> >>>>>>>>>>>>>> You, as the designer, will get to choose whatever names you >>>>>>>>>>>>>> want. We will discuss the naming of functions during deeper >>>>>>>>>>>>>> discussions. >>>>>>>>>>>>>> For now just pick whatever you think is appropriate. >>>>>>>>>>>>>> >>>>>>>>>>>>>> >>>>>>>>>>>>>> Jason >>>>>>>>>>>>>> moorepants.info >>>>>>>>>>>>>> +01 530-601-9791 >>>>>>>>>>>>>> >>>>>>>>>>>>>> On Sun, Feb 28, 2016 at 1:53 AM, SAMPAD SAHA < >>>>>>>>>>>>>> [email protected]> wrote: >>>>>>>>>>>>>> >>>>>>>>>>>>>>> Jason, >>>>>>>>>>>>>>> >>>>>>>>>>>>>>> Thank you for giving me a direction to think over this >>>>>>>>>>>>>>> project. >>>>>>>>>>>>>>> >>>>>>>>>>>>>>> Would it be fine if we name the module as Singularity? >>>>>>>>>>>>>>> >>>>>>>>>>>>>>> As you said, I will now focus on creating object which would >>>>>>>>>>>>>>> help us to use standard operators on different types of >>>>>>>>>>>>>>> singularity >>>>>>>>>>>>>>> functions. >>>>>>>>>>>>>>> I will also focus on creating a timeline to implement >>>>>>>>>>>>>>> different functionalities . I will let you know shortly about >>>>>>>>>>>>>>> what i am >>>>>>>>>>>>>>> upto. >>>>>>>>>>>>>>> >>>>>>>>>>>>>>> >>>>>>>>>>>>>>> >>>>>>>>>>>>>>> >>>>>>>>>>>>>>> >>>>>>>>>>>>>>> Regards >>>>>>>>>>>>>>> Sampad Kumar Saha >>>>>>>>>>>>>>> Mathematics and Computing >>>>>>>>>>>>>>> I.I.T. Kharagpur >>>>>>>>>>>>>>> >>>>>>>>>>>>>>> On Sun, Feb 28, 2016 at 12:11 AM, Jason Moore < >>>>>>>>>>>>>>> [email protected]> wrote: >>>>>>>>>>>>>>> >>>>>>>>>>>>>>>> Sampad, >>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>> You will likely need to create objects that represent >>>>>>>>>>>>>>>> different types of singularity functions. These objects should >>>>>>>>>>>>>>>> be able to >>>>>>>>>>>>>>>> be assembled into full expressions using standard operators: >>>>>>>>>>>>>>>> plus, minus, >>>>>>>>>>>>>>>> etc. Additionally you'll need some methods that can convert >>>>>>>>>>>>>>>> these >>>>>>>>>>>>>>>> expressions into the continuous piecewise forms and allow >>>>>>>>>>>>>>>> evaluation with >>>>>>>>>>>>>>>> .evalf(). The integrate and Integral objects should be able to >>>>>>>>>>>>>>>> integrate >>>>>>>>>>>>>>>> them. >>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>> I'd suggest writing up some psuedo code that would give an >>>>>>>>>>>>>>>> idea of what a sympy session might look like to solve some >>>>>>>>>>>>>>>> beam problems. >>>>>>>>>>>>>>>> Maybe something like: >>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>> >>> from sympy import SingularityFunction as f >>>>>>>>>>>>>>>> >>> w = -3 * f(x - 0)**(-1) + 6 * f(x - 2)**0 - 9 * f(x - >>>>>>>>>>>>>>>> 4)**-1 >>>>>>>>>>>>>>>> >>> pprint(w) >>>>>>>>>>>>>>>> -3 * f<x - 0>**(-1) + 6 * <x - 2>**0 - 9 * <x - 4>**-1 >>>>>>>>>>>>>>>> >>> integrate(w) >>>>>>>>>>>>>>>> -3 * <x - 0>**0 + 6 * <x - 2>**1 - 9 * <x - 4>**0 >>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>> Do some beam bending problems by hand and then think about >>>>>>>>>>>>>>>> how you'd like it to work if you could do it all with SymPy. >>>>>>>>>>>>>>>> For every >>>>>>>>>>>>>>>> example problem you do, write out pseudo code like above. >>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>> I could also imagine a higher level abstraction where you >>>>>>>>>>>>>>>> could do something like: >>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>> from sympy import Beam, DistributedLoad, PointLoad >>>>>>>>>>>>>>>> b = Beam(length=l) >>>>>>>>>>>>>>>> b.apply(DistrubutedLoad(start=l/2, end=l, value=w) >>>>>>>>>>>>>>>> b.apply(PointLoad(location=l/3, value=p) >>>>>>>>>>>>>>>> b.moment() # finds the moment vs x using singularity funcs >>>>>>>>>>>>>>>> in the background >>>>>>>>>>>>>>>> b.shear_moment_diagram() # uses the plot module to plot the >>>>>>>>>>>>>>>> piecewise functoin >>>>>>>>>>>>>>>> b.moment_sing_func() # returns the singularity funcs >>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>> But you'll need to implement the singularity function stuff >>>>>>>>>>>>>>>> first before the later. >>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>> Jason >>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>> Jason >>>>>>>>>>>>>>>> moorepants.info >>>>>>>>>>>>>>>> +01 530-601-9791 >>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>> On Fri, Feb 26, 2016 at 1:29 PM, SAMPAD SAHA < >>>>>>>>>>>>>>>> [email protected]> wrote: >>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>> Hi Jason, >>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>> Thank you for the explanation. It really helped me. >>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>> So, basically we want to start it, firstly, by creating a >>>>>>>>>>>>>>>>> module which would deal with the mathematical operations >>>>>>>>>>>>>>>>> performed on >>>>>>>>>>>>>>>>> Singularity Functions. After this whole module is prepared, >>>>>>>>>>>>>>>>> we would focus >>>>>>>>>>>>>>>>> on how to use this module for solving beam problems. Am I >>>>>>>>>>>>>>>>> correct? >>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>> Can you please explain me in brief that what are the >>>>>>>>>>>>>>>>> mathematical operations we wanted to implement on that module? >>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>> On Friday, February 26, 2016 at 4:54:59 PM UTC+5:30, >>>>>>>>>>>>>>>>> SAMPAD SAHA wrote: >>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>> Hi, >>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>> I am Sampad Kumar Saha , an Undergraduate Mathematics and >>>>>>>>>>>>>>>>>> Computing Student at I.I.T. Kharagpur. >>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>> I have gone through the idea page >>>>>>>>>>>>>>>>>> <https://github.com/sympy/sympy/wiki/GSoC-2016-Ideas> >>>>>>>>>>>>>>>>>> and I am interested in working on the project named >>>>>>>>>>>>>>>>>> Singularity >>>>>>>>>>>>>>>>>> Function >>>>>>>>>>>>>>>>>> <https://github.com/sympy/sympy/wiki/GSoC-2016-Ideas#singularity-functions> >>>>>>>>>>>>>>>>>> . >>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>> By going through the Idea, I understood that we want to >>>>>>>>>>>>>>>>>> add a package to Sympy which can be used for for solving >>>>>>>>>>>>>>>>>> beam bending >>>>>>>>>>>>>>>>>> stress and deflection problems using singularity function. >>>>>>>>>>>>>>>>>> Am I correct? >>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>> We can by this way:- >>>>>>>>>>>>>>>>>> While solving we will be having the moment function as an >>>>>>>>>>>>>>>>>> input which we can arrange in the form of singularity >>>>>>>>>>>>>>>>>> functions and then >>>>>>>>>>>>>>>>>> integrate it twice to get the deflection curve and we can >>>>>>>>>>>>>>>>>> give the plot or >>>>>>>>>>>>>>>>>> the equation obtained of deflection curve as an output. >>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>> I have gone through some documents available on internet >>>>>>>>>>>>>>>>>> which have brief studies on solving beam bending stress and >>>>>>>>>>>>>>>>>> deflection >>>>>>>>>>>>>>>>>> problems using singularity functions. >>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>> References:- >>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>> - Beam Deflection By Discontinuity Functions >>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>> <http://www.colorado.edu/engineering/CAS/courses.d/Structures.d/IAST.Lect12.d/IAST.Lect12.pdf> >>>>>>>>>>>>>>>>>> . >>>>>>>>>>>>>>>>>> - Beam Equation Using Singularity Functions >>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>> <https://www.civil.iitb.ac.in/~sghosh/CE327/singularity-1.pdf> >>>>>>>>>>>>>>>>>> . >>>>>>>>>>>>>>>>>> - Enhanced Student Learning in Engineering Courses >>>>>>>>>>>>>>>>>> with CAS Technology >>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>> <http://cei.ust.hk/files/public/sympo/2007/paper_r13-thomashu.pdf> >>>>>>>>>>>>>>>>>> . >>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>> Since there is just a brief idea given in the idea page, >>>>>>>>>>>>>>>>>> I have a doubt that what are the things other than solving >>>>>>>>>>>>>>>>>> beam bending >>>>>>>>>>>>>>>>>> stress and deflection problems to be implemented in the >>>>>>>>>>>>>>>>>> project? >>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>> Any type of suggestions are welcome. >>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>> ========================================================================================================================================== >>>>>>>>>>>>>>>>>> Regards >>>>>>>>>>>>>>>>>> Sampad Kumar Saha >>>>>>>>>>>>>>>>>> Mathematics and Computing >>>>>>>>>>>>>>>>>> I.I.T. Kharagpur >>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>> >>>>>>>>>>>>>>> >>>>>>>>>>>>>> >>>>>>>>>>>>> >>>>>>>>>>>> >>>>>>>>>>> >>>>>>>>>> >>>>>>>>> >>>>>>>> >>>>>>> >>>>>> >>>>> >>>> >>> >> > -- 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 post to this group, send email to [email protected]. 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/CAP7f1AiCh6ZVYU1zJaSOg0QQ_JtSBJG9hD7BBS4Kf0Dydp4MEw%40mail.gmail.com. For more options, visit https://groups.google.com/d/optout.
