On Thu, Mar 24, 2011 at 12:12 AM, Aaron S. Meurer <asmeu...@gmail.com> wrote:
> On Mar 23, 2011, at 10:20 PM, Tim Lahey wrote:
>
>> Hi,
>>
>> I've been thinking about GSoC and looking at the various ideas people have 
>> been posting. By the way, I'd like to see perturbation theory included in 
>> general (i.e., not just for quantum mechanics) since it's used in stability 
>> analysis.
>
> Yes definitely.  Another thing that should not be implemented only in physics.

My only hesitancy is that perturbation *theories* tend to be very
specific to the equations you are solving.  For example, even within
quantum mechanics, there are probably a dozen different perturbation
theories.  The reason for this is that:

* it is non-trivial to identify and fully solve the non-perturbed system.
* Writing a perturbation series that converges at all or quickly is
not an easy task.

I don't disagree that it would be nice to have this, but I think it is
not the place to *start* with perturbation theory for a GSoC project.
I am biased, but the various quantum perturbation theories are well
known and would be a good place for starting to identify the common
aspects of all perturbation theories (similar to how the operators,
commutators, etc. in physics.quantum will end up doing).

>>
>> I have a few linear algebra related ideas I'd like to see. One is "abstract" 
>> linear algebra. That is support for arbitrary matrix and vector 
>> calculations. So, one would just indicate that a symbol is a matrix/vector 
>> without stating the size. Then you could work with them. The non-commutative 
>> support is a start, but the transpose operation would need to be added along 
>> with derivatives. I have a Maple worksheet that has a crude version of this. 
>> This is used quite a bit in the Ritz approach to Finite Element Analysis.

Yes, this could be quite useful.

Cheers,

Brian

>>
>> The other linear algebra thing to add is support for block matrices. So, one 
>> could specify a matrix like
>>
>> M = [[A,B],[C,D]]
>>
>> where A, B, C, and D are arbitrary matrices. Then, you could take the 
>> inverse and other standard linear algebra. This fits in with the "abstract" 
>> linear algebra support. This kind of thing pops up a lot in control theory.
>>
>> Cheers,
>>
>> Tim.
>>
>> ---
>> Tim Lahey
>> PhD Candidate, Systems Design Engineering
>> University of Waterloo
>> http://about.me/tjlahey
>
> These sound like good ideas.
>
> By the way, Sherjil Ozair is also applying to do some work on the matrices, 
> but I think he is planning on doing different stuff, like sparse matrices and 
> integration with the polys ground types (gmpy).  Maybe you could collaborate 
> with him to apply to do different things.  There are enough things to do in 
> the matrices to cover two projects, I think.  That way (assuming you both had 
> high enough quality applications) you could both potentially be accepted at 
> the same time.
>
> Otherwise, if you apply to do the exact same thing, we will have to pick one 
> of you (again assuming both are good enough to be accepted).
>
> Aaron Meurer
>
> --
> You received this message because you are subscribed to the Google Groups 
> "sympy" group.
> To post to this group, send email to sympy@googlegroups.com.
> To unsubscribe from this group, send email to 
> sympy+unsubscr...@googlegroups.com.
> For more options, visit this group at 
> http://groups.google.com/group/sympy?hl=en.
>
>



-- 
Brian E. Granger, Ph.D.
Assistant Professor of Physics
Cal Poly State University, San Luis Obispo
bgran...@calpoly.edu
elliso...@gmail.com

-- 
You received this message because you are subscribed to the Google Groups 
"sympy" group.
To post to this group, send email to sympy@googlegroups.com.
To unsubscribe from this group, send email to 
sympy+unsubscr...@googlegroups.com.
For more options, visit this group at 
http://groups.google.com/group/sympy?hl=en.

Reply via email to