[sympy] Re: How to build a new printer

2009-07-04 Thread Luke
.py;h=40a8853b71e4642981096b463d46da7392fcae76;hp=587b44af846bbb76ef218559666a0f3635a88c1b;hb=0ffaf649242f4876c591158d19f33a6b0b17ee0b;hpb=162e518e9e28403ba78d27b523e946beefa4d47e Hopefully this will give you an idea of what you need to do to customize your printer. ~Luke On Jul 4, 7:49 am, Christophe

[sympy] Re: sphinx documentation for the rc1 release

2009-07-13 Thread Luke
I don't know if this will change when you push to the docs.sympy.org server, but the links for the books by Hestenes and Lasenby at the bottom of the GA Module currently don't work. Just thought I'd let you know. ~Luke On Jul 13, 7:59 pm, Alan Bromborsky wrote: > Ondrej Ce

[sympy] Automatic distribution of numbers versus symbols

2009-07-28 Thread Luke
2008a, this automatic distribution is not easily reversible by a simple call to factor. It seems to me that: 1) there should be a way to control whether or not distribution occurs automatically, and 2) that if distribution does occur automatically, it should be easily reversible

[sympy] Add combine question

2009-08-14 Thread Luke
cally'. It is interesting that Mul can auto combine things like: rf*(1 - (c5*s2 + c2*s4*s5))*(1 - (c5*s2 + c2*s4*s5))**(-0.5), but Add does not. Would this be hard to implement in Add? Thanks, ~Luke --~--~-~--~~~---~--~~ You received this message because you are

[sympy] Re: Add combine question

2009-08-14 Thread Luke
if: 1) Both Add terms (a1 and a2, say) are Mul instances and (set(a1.args) & set(a2.args)) != set([]) or 2) One Mul (a1) and one Symbol (a2) and (set(a1.args) & set(a2)) != set([]) then they have common factors and could be combined. Any thoughts? ~Luke On Aug 14, 9:49 pm, Luke wrote: &g

[sympy] Re: los alamos sprint

2009-08-15 Thread Luke
> > >>> we are having a little sprint in Los Alamos now, with me, Aaron and > >>> Luke. If you want to join us, we are on irc. > > >>> Me and Luke is merging the new trig branch, Aaron is working on   > >>> dsolve > >>> and related bug

[sympy] Issue 1616

2009-08-25 Thread Luke
ut The following traceback may be corrupted or invalid The error message is: ('EOF in multi-line statement', (96, 0)) --- ValueErrorTraceback (most recent call last) /home/luke/lib/pyt

[sympy] Re: Issue 1616

2009-08-25 Thread Luke
Bill has posted a patch to the Issue page on Google code. It fixes this problem for everything I could throw at it. ~Luke On Aug 25, 12:47 am, Luke wrote: > I don't know what changed, but this issue is now causing a lot of my > code in pydy to not be able to use .subs() on

[sympy] New code in bicycle branch

2009-09-24 Thread Luke
ems in the list of Relational objects. I have been using it to create functions that work well with scipy's odeint and fsolve and it has worked very well. To see it in action, run the rollingdisc.py script in the examples directory. It will print out the generated string exp

[sympy] Symbolic matrix inversion

2009-09-29 Thread Luke
titution dictionaries) for inverting a matrix full of sympy expressions. Does anybody know if this has been done by somebody somewhere, or have any other ideas on how it could be done better than the way I suggested? Thanks, ~Luke --~--~-~--~~~---~--~~ You receiv

[sympy] Re: Symbolic matrix inversion

2009-09-29 Thread Luke
On Sep 29, 1:09 pm, Ondrej Certik wrote: > On Tue, Sep 29, 2009 at 12:49 PM, Luke wrote: > > I'm using Sympy from within PyDy to generate the equations of motion for > > mechanical systems.  At the end of the day,  the equations can be most > > generally written a

[sympy] Re: Symbolic matrix inversion

2009-09-29 Thread Luke
as the ode formulation would be as generic as possible. ~Luke On Sep 29, 4:15 pm, Alan Bromborsky wrote: > Luke wrote: > > > On Sep 29, 1:09 pm, Ondrej Certik wrote: > > >> On Tue, Sep 29, 2009 at 12:49 PM, Luke wrote: > > >>> I'm using Sympy fr

[sympy] Re: Symbolic matrix inversion

2009-09-29 Thread Luke
s, you would end up only integrating 1 kinematic D.E. (the one associated with the lean of the disc), and 3 dynamic D.E.'s. So in these types of systems, it depends on what your needs are. ~Luke On Sep 29, 4:22 pm, Alan Bromborsky wrote: > Alan Bromborsky wrote: > > Luke wrote: >

[sympy] Re: Symbolic matrix inversion

2009-09-29 Thread Luke
I also check the GSL (GNU Scientific Library). They have a nice numerical integrator, but it doesn't allow for a mass matrix. ~Luke On Sep 29, 4:44 pm, Luke wrote: > Yes, this is something I should look into.  I am pretty sure that the > Netlib codes have this functionality, but it

[sympy] Re: Symbolic matrix inversion

2009-09-30 Thread Luke
also eliminates the iteration that is being done by the ODE solver during time integration. Thanks, ~Luke On Sep 29, 8:07 pm, Tim Lahey wrote: > On Sep 29, 2009, at 7:15 PM, Alan Bromborsky wrote: > > > Are there differential equation solvers where you don't have to invert >

[sympy] Re: Symbolic matrix inversion

2009-09-30 Thread Luke
immediately useful to those audiences with minimal hassle. ~Luke [0] -- J. P. Meijaard, Jim M. Papadopoulos, Andy Ruina, A. L. Schwab, 2007 ``Linearized dynamics equations for the balance and steer of a bicycle: a benchmark and review,'' Proceedings of the Royal Society A 463:1955-1982. On Se

[sympy] Re: Symbolic matrix inversion

2009-09-30 Thread Luke
On Sep 30, 8:39 am, Luke wrote: > The methods you suggest essentially takes care of the mass matrix > problem by solving a linear system numerically during numerical > integration.  I am familiar with tools out there that do this, but > this isn't what I'm looking to d

[sympy] Preventing expansion

2009-10-11 Thread Luke
Pow(x*y, 2, expand=False) but this didn't do it. Is there a way? Thanks, ~Luke --~--~-~--~~~---~--~~ 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 un

[sympy] GSoC project mentoring

2011-04-05 Thread Luke
development process should be lowered by the resources and experience of the people in my lab. This will help both students to hit the ground running a lot more quickly. Cheers, ~Luke [0] -- pydy.org -- You received this message because you are subscribed to the Google Groups "sympy&

[sympy] Re: GSoC project mentoring

2011-04-05 Thread Luke
reate NewtonianReferenceFrame and then call the declare_coords on it probably should be abstracted to a different class. This is something that we are going to think about more carefully than I did when I first work on PyDy. Cheers, ~Luke On Apr 5, 1:51 pm, Tim Lahey wrote: > On Tue, Apr

[sympy] Re: GSoC project mentoring

2011-04-05 Thread Luke
give it a look, but it hasn't been touched since June 26, 2009, so there may be quite a lot that breaks once I rebase from the current master. ~Luke On Apr 5, 3:02 pm, Aaron Meurer wrote: > I agree with Ondrej. > > Since PyDy is (still at least) separate form SymPy, another ide

[sympy] Re: GSoC project mentoring

2011-04-05 Thread Luke
would start working on. So are 3D symbolic vector operations and the framework for managing orientations and positions, something that users of SymPy find useful? If so, I think this would justify it's inclusion as a module, and I would like to help make it happen. Thoughts? Thanks,

[sympy] Re: GSoC project mentoring

2011-04-06 Thread Luke
ny comments, we would greatly appreciate it. ~Luke On Apr 5, 9:41 pm, Brian Granger wrote: > On Tue, Apr 5, 2011 at 8:11 PM, Luke wrote: > > Yeah, that is a good idea -- thanks! > > > Regardless of whether these two applicants end up being accepted, I > > would like to get back

[sympy] class design / modelling

2011-04-11 Thread Luke
n this approach to class design? ~Luke -- 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 opt

Re: [sympy] Question regarding vectors

2011-05-09 Thread Luke
nk all the other types of products have to make use of explicit method calls since there would be no way to know which type of product would be implied. ~Luke On Mon, May 9, 2011 at 4:23 PM, Gilbert Gede wrote: > In PyDy (which we plan to merge into SymPy.physics.classical this summer) > Ve

Re: [sympy] Question regarding vectors

2011-05-09 Thread Luke
Can the vectors and multivectors in the GA module work with arbitrary sympy expressions? i.e, if v is a GA vector, and s is a sympy expression, does it make sense to do: s*v? Is the result of type Mul or of something else? ~Luke On Mon, May 9, 2011 at 5:42 PM, Alan Bromborsky wrote: > On

Re: [sympy] Question regarding vectors

2011-05-10 Thread Luke
asses that implement only the methods that make sense for them, including the __mul__ and __rmul__ methods, as Alan did in his MV class? ~Luke On Tue, May 10, 2011 at 5:44 AM, Alan Bromborsky wrote: > The assumption is that the expression multiplying the vector (multivector) > is a scalar and t

Re: [sympy] Re: Question regarding vectors

2011-05-10 Thread Luke
around the class methods, so the user can choose to use whichever approach is more natural to them. ~Luke On Tue, May 10, 2011 at 11:02 AM, Vinzent Steinberg wrote: > On 10 Mai, 03:30, Ronan Lamy wrote: >> Please, try to make the interface dot(v1, v2), and not v1.dot(v2). > > While

Re: [sympy] Re: Who uses SymPy?

2011-05-12 Thread Luke
but I think all the other benefits of FOSS will become clear to these people once they take the time to try it. But there is a whole generation of entrenchment that will eventually retire and make room for open source, so I see things getting better as time goes on. ~Luke On Thu, May 12, 2011 at 10:

Re: [sympy] Re: Question regarding vectors

2011-05-14 Thread Luke
ke sure we can accommodate many different reference frames which is common in multibody dynamics. ~Luke On Fri, May 13, 2011 at 6:29 PM, Gilbert gede wrote: > So I think we've decided to make a new Vector class, to replace the > previous UnitVector and Vector classes, and it will not ext

[sympy] Problem with user-defined functions

2011-06-02 Thread luke
Hi eveyone, I'm developing an web application which has to interact with "user- defined formulas" of some financial kpis. I decided to use sympy to have a more solid math engine. Basically the input I reiceve is very simple, it might be in the worst case something like: kpi -> "(log(sum('productio

[sympy] Re: Problem with user-defined functions

2011-06-02 Thread luke
would have been: >>> a = SUM('field') + SUM('field') + SUM('field') -> 3SUM('field') # one >>> single query >>> print a >>> 3*SUM('field') >>> print N(a) >>> 1234 just like other functions wo

[sympy] Re: Problem with user-defined functions

2011-06-02 Thread luke
e, "myresults") return result.find_one()['value'][unicode(cls.arg)] Because arg in _eval_evalf was strangely enough a float (more precisely 57, have no clue why :) ). Thank you all guys! On Jun 2, 6:07 pm, Ronan Lamy wrote: > Le jeudi 02 juin 2011 à 08:03 -0700, luke

[sympy] Re: Problem with user-defined functions

2011-06-02 Thread luke
a symbol but is a float (57 in my case, probably is an index of some internal variable.. I can't find any documentation about that). So it's pretty unusable, that's why I set cls.arg. However the code is working so thanks everybody. On Jun 2, 7:01 pm, Ronan Lamy wrote: > Le je

[sympy] Re: Problem with user-defined functions

2011-06-02 Thread luke
WAIT, I didn't see that you use self.arg[0]!! that's why my code was giving me 57, I just need to call cls.arg[0], great now I can remove the eval ;)! thanks again On Jun 2, 9:30 pm, luke wrote: > Actually you're wrong. Every instance of a class in python has its own > attri

[sympy] Re: Problem with user-defined functions

2011-06-02 Thread luke
Nope in spite of my enthusiasm cls.args wont work properly as it will give me a and not a :) I'll put back my eval since it worked. Hope that was of any help! On Jun 2, 9:32 pm, luke wrote: > WAIT, I didn't see that you use > self.arg[0]!! that's why my code was giving

[sympy] Re: Problem with user-defined functions

2011-06-02 Thread luke
s != undefined){ sum += doc.%(field)s; } }); return {%(field)s:sum}; };""" % {'field':cls.args[0]}) result = db.people.map_reduce(map, reduce, "myresults&qu

[sympy] Adding sympy.physics to documentation

2011-06-07 Thread Luke
n add some skeleton files that will get the basic sphinx framework started, and then the submodules can be more completely fleshed out and polished by the people who are most familiar with the various bits of code that sympy/physics includes. Thoughts? ~Luke [0] -- https://github.com/gilbertgede/

Re: [sympy] Re: Adding sympy.physics to documentation

2011-06-07 Thread Luke
cs.sympy.org in the past (other than lack of time)? > > No, it should be there, it is just that noone yet made sure it's > generated correctly when building the docs. > I'll get started on this. ~Luke -- You received this message because you are subscribed to the Google Gro

Re: [sympy] Re: Adding sympy.physics to documentation

2011-06-07 Thread Luke
ppy to go over the code and > document the parts that I use. Cool, thanks! I'll mail the list when I get this up and running on a github branch. ~Luke > > Stefan > > On 7 June 2011 20:34, Luke wrote: >> >> >> 1)  Should any of the code in sympy.physics b

[sympy] bugs in units

2011-06-08 Thread Luke
I was writing some documentation for physics/units.py and discovered some infinite recursion RuntimeError: >>> from sympy.physics.units import joule >>> joule .... File "/home/luke/repos/sympy/sympy/core/expr.py", line 140, in __lt__ return C.StrictInequality(sel

Re: [sympy] Re: Adding sympy.physics to documentation

2011-06-12 Thread Luke
+1 for sympy.physics.mechanics ~Luke On Thu, Jun 9, 2011 at 12:56 PM, Brian Granger wrote: > On Tue, Jun 7, 2011 at 10:28 PM, Luke wrote: >> On Tue, Jun 7, 2011 at 1:26 PM, krastanov.ste...@gmail.com >> wrote: >>> What would the classical module contain? Won'

[sympy] sympy/physics sphinx documentation

2011-06-26 Thread Luke
ished then push it in as a whole? ~Luke -- 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

Re: [sympy] Question on Interfaces

2011-08-11 Thread Luke
. If people have thoughts on any negative implications of implementing multiple ways to perform the same thing, especially related to long term maintenance of the code, it would be great to hear them. ~Luke On Thu, Aug 11, 2011 at 12:14 PM, Gilbert Gede wrote: > My Summer of Code project is writing

Re: [Sports Biomechanics] Re: [sympy] Question on Interfaces

2011-08-11 Thread Luke
your library (sympy.physics.mechanics). If you > rewrite the language definition, then it doesn't necessarily fit with the > rest of the code. > Jason > > On Thu, Aug 11, 2011 at 12:42 PM, Luke wrote: >> >> I ran this same question by my girlfriend who teaches under

[sympy] MathJax script tag

2011-08-16 Thread Luke
entered "make html" from the sympy/doc directory. Is this a bug, or am I doing something incorrectly? Thanks, ~Luke -- "Those who would give up essential liberty to purchase a little temporary safety deserve neither liberty nor safety." -- Benjamin Franklin, Historical Revie

Re: [sympy] MathJax script tag

2011-08-16 Thread Luke
Ok, I just upgraded to Sphinx 1.0.7 and it fixed the problem. Thanks, ~Luke On Tue, Aug 16, 2011 at 4:10 PM, Aaron Meurer wrote: > It works for me.  The only thing that I can think of is that your > version of Sphinx has a bug in it.  I am using 1.0.7, which generates >   src

[sympy] doctest testmod() globs keyword argument

2011-08-20 Thread Luke
ypes: python sympy/physics/doctest_globaldict.py[1] F [FAIL] __ sympy.physics.doctest_globaldict.foo File "/home/luke/repos/sympy/sympy/physics/doctest_globaldict.py", li

Re: [sympy] doctest testmod() globs keyword argument

2011-08-20 Thread Luke
Something in testmod() must be different from the default doctest module, I'll give it a look. > By the way, when you say "running Python on this file," do you mean > running "python -m doctest myfile.py"? > I mean "python myfile.py". ~Luke > Aar

Re: [sympy] doctest testmod() globs keyword argument

2011-08-21 Thread Luke
On Sun, Aug 21, 2011 at 10:42 AM, Aaron Meurer wrote: > On Sat, Aug 20, 2011 at 7:34 PM, Luke wrote: >> On Sat, Aug 20, 2011 at 3:37 PM, Aaron Meurer wrote: >>> I don't know of any doctests in SymPy that do this.  Why can't you put >>> imports in e

Re: [sympy] doctest testmod() globs keyword argument

2011-08-22 Thread Luke
hereas the Python standard library, doctest is a module with a large number of functions, testmod() being one of them. ~Luke On Mon, Aug 22, 2011 at 10:32 AM, Aaron Meurer wrote: > Well, I don't know why you are doing things this way, but clearly > python -m doctest file.py or ./bin/docte

Re: [sympy] doctest testmod() globs keyword argument

2011-08-25 Thread Luke
On Thu, Aug 25, 2011 at 7:00 AM, Ronan Lamy wrote: > > Le samedi 20 août 2011 à 18:34 -0700, Luke a écrit : > > On Sat, Aug 20, 2011 at 3:37 PM, Aaron Meurer wrote: > > > I don't know of any doctests in SymPy that do this.  Why can't you put > > > import

[sympy] Import errors with an "installed" version of Sympy

2011-10-16 Thread Luke
then try to import sympy, I get the following import errors: luke@ThinkPad-W510:~$ echo $PYTHONPATH luke@ThinkPad-W510:~$ python Python 2.7.2+ (default, Oct 4 2011, 20:06:09) [GCC 4.6.1] on linux2 Type "help", "copyright", "credits" or "license" for more i

Re: [sympy] Import errors with an "installed" version of Sympy

2011-10-16 Thread Luke
ted when I did the `python setup.py install` step. I deleted the build directory, and repeated my steps, and now everything works as expected. Thanks, ~Luke -- You received this message because you are subscribed to the Google Groups "sympy" group. To post to this group, send email to sympy@g

[sympy] Feynman Lectures Exercise Challenge

2011-11-14 Thread Luke
figure it out, apparently Feynman did. My solution is here: http://dlpeterson.com/FLP_Exercise_Challenge/solution.pdf ~Luke -- 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 u

Re: [sympy] Re: Feynman Lectures Exercise Challenge

2011-11-14 Thread Luke
Ondrej, Great comments! How are you!? We need to get a beer next time you are coming through -- do you have a regular schedule, maybe Google Calendar that we can share? > 1) How did you create the drawing in the pdf? I used TikZ. I looked into PStricks, which is very powerful, but it seems li

Re: [sympy] Re: Feynman Lectures Exercise Challenge

2011-11-14 Thread Luke
cently, I've used pydy to derive the equations of motion for a rattleback [0], a non-trivial nonholonomic system that has some very interesting dynamic behavior (and which can be predicted by Newtonian mechanics). I'll try to add the cart-pendulum system when I have time. ~Luke [0] -- h

Re: [sympy] Re: Feynman Lectures Exercise Challenge

2011-11-15 Thread Luke
I found an error in my calculations, I have corrected it in the most recent commit on the github page. ~Luke -- 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

Re: [sympy] Re: Feynman Lectures Exercise Challenge

2011-11-15 Thread Luke
rrect that the length of the pendulum is not needed if you want the displacement in the x or y directions (l cancels out). It is needed, however, if you want the amplitude of the angular displacements (see On Tue, Nov 15, 2011 at 4:32 PM, Alan Bromborsky wrote: > On 11/15/2011 02:18 PM, Luke w

Re: [sympy] Re: Feynman Lectures Exercise Challenge

2011-11-15 Thread Luke
On Tue, Nov 15, 2011 at 6:41 PM, Luke wrote: > Alan, >  Looks good, thanks.  Yeah, Lagrange's method is the way to go on > this problem, the reason I chose to do F=ma directly was purely > because the problem is targeted towards a freshman audience and they > won't be fa

Re: [sympy] Re: Feynman Lectures Exercise Challenge

2011-11-15 Thread Luke
have to use Lagrange multipliers, or is there another way? In the approach I outlined, you would simply dot the vector form of the motion equation into a direction aligned with the rod (-sin(theta)*\hat{i} - cos(theta)*\hat{j}), and then solve for F. ~Luke -- You received this message because you a

Re: [sympy] Re: Feynman Lectures Exercise Challenge

2011-11-15 Thread Luke
> Are you asking, for  example, how to calculate the tension in the pendulum > arm > using the Lagrangian method? Yes. -- 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 thi

Re: [sympy] Re: Feynman Lectures Exercise Challenge

2011-11-16 Thread Luke
ly obtain from a differential equation. But perhaps that relationship is derivable without appealing to differential equations. In any event, there must be some way to relate period to length, and if it isn't this relationship, I don't know what it is or how to rationalize it. ~Luke -- Y

Re: [sympy] Re: Feynman Lectures Exercise Challenge

2011-11-16 Thread Luke
Do you mean in my solution that involves differential equations? I think this dimensional analysis approach may have merit, I just need to see all the steps and make sure they can all be justified without referencing a differential equation. ~Luke On Wed, Nov 16, 2011 at 10:08 AM, Aaron Meurer

Re: [sympy] Re: Feynman Lectures Exercise Challenge

2011-11-16 Thread Luke
I don't feel I have any deeper understanding of pendulums, mechanics, or geometry, than when I started. If I sound like I'm criticizing you here, I apologize, I am not. I am primarily dissatisfied with this question and others like it, and I really don't see the point of them. ~Luke

Re: [sympy] Re: Feynman Lectures Exercise Challenge

2011-11-16 Thread Luke
On Wed, Nov 16, 2011 at 1:57 PM, Luke wrote: > A very similar solution was proposed on the forum and rejected because > he made use of the assumption of \omega_n^2 = sqrt(g/l). I meant \omega_n^2 = g/l ~Luke -- You received this message because you are subscribed to the Google Groups

Re: [sympy] Re: I am interested in applying to Google Summer of Code to work on a SymPy project

2012-03-27 Thread Luke
reliably, it makes sense to pursue the flexible body case. I don't think it make sense to pursue the flexible case until after that though, and code output will vary significantly between the two because you obtain PDE's in one case and ODE's in the other. Luke > > > On

Re: [sympy] Re: GSoC Idea

2012-04-02 Thread Luke
hink getting sympy.physics.mechanics into sage would solve this problem very well. Luke -- 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+unsubs

Re: [sympy] Re: I am interested in applying to Google Summer of Code to work on a SymPy project

2012-04-02 Thread Luke
asons to have a custom Lagrangian class. Luke -- 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 m

[sympy] Symbolic computation webinar

2012-06-07 Thread Luke
The Maple director of research, Dr. Jürgen Gerhard, is going to be interviewed June 28th in the following webinar: http://spectrum.ieee.org/webinar/2075921 I thought people here might find it interesting to hear what he has to say about symbolic computation. ~Luke -- "People call

[sympy] Identifying repeated subexpressions in systems of equations

2008-06-13 Thread Luke
hen this common subexpression replacement method is used. Any ideas or references on this subject? Thanks, ~Luke --~--~-~--~~~---~--~~ You received this message because you are subscribed to the Google Groups "sympy" group. To post to this group, send email to

[sympy] Re: Identifying repeated subexpressions in systems of equations

2008-06-15 Thread Luke
re yet how to approach. Thanks for the response and if you have any other ideas or references, I'd love to hear them! I'm guessing I'm going to need to study up on searching, sorting, and regular expressions. Thanks, ~Luke On Jun 14, 11:10 am, "Ondrej Certik" &

[sympy] Re: Identifying repeated subexpressions in systems of equations

2008-06-16 Thread Luke
ose subexpression and identify subsubexpressions), and introduce variables for any that are used more than once. If a subexpression is found, but doesn't occur more than once in any of the preceding equations, but then later shows up in another equation, it would then be introduced as

[sympy] Re: Identifying repeated subexpressions in systems of equations

2008-06-16 Thread Luke
Robert, The Mathematica link you provided is exactly what I'm trying to do. I haven't tried your python code yet but after reading it I think it should work great. I really appreciate your comments and your help! Thanks, ~Luke On Jun 16, 12:14 pm, "Robert Kern" <[

[sympy] Symbolic dot and cross products with implicit rotation matrices

2009-01-16 Thread Luke
and python fairly well, and I know the kinematics well, but I'm by no means an experience object oriented programmer, so I'm not sure about how the best way to structure things would be. Any suggestions on how I might start on something like this? Thanks, ~Luke --~--~-~--~~--

[sympy] Re: Symbolic dot and cross products with implicit rotation matrices

2009-01-17 Thread Luke
fine the position of one point q relative to another p: p_q_p = q2*a1 + q3*b2 and the be able to do: dot(p_q_p,a1) and get: q2 - sin(q1)*q3 (recall, b2 = -sin(q1)*a1+cos(q1)*a2) Hopefully this clarifies what I'm trying to do. Thanks, ~Luke On Jan 17, 5:12 am, Alan Bromborsky wrote: > L

[sympy] Identifying additive terms in an expression

2009-01-17 Thread Luke
If I have a symbolic expression like: expr = a*sin(x) + b*cos(y) - c*tan(z) I can get each individual additive term, sort of like: In[1]: terms(expr) Out[1]: [a*sin(x), b*cos(y), -c*tan(z)] ? Thanks, ~Luke --~--~-~--~~~---~--~~ You received this message

[sympy] Fwd: Confused about something

2009-04-16 Thread Luke
is how I implemented it. Should the ._args list for the above example then be: v1._args == [5*A[1], 6*A[2]] ? Thanks, ~Luke -- Forwarded message -- From: Ondrej Certik Date: Thu, Apr 16, 2009 at 9:37 PM Subject: Re: Confused about something To: Luke On Thu, Apr 16, 2009 at 8:

[sympy] Printing of functions

2009-04-17 Thread Luke
d I go about tackling this? Thanks, ~Luke --~--~-~--~~~---~--~~ 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 t

[sympy] Trigsimp error

2009-04-28 Thread Luke
Traceback (most recent call last) /home/luke/Documents/PythonDynamics/ in () /usr/lib/pymodules/python2.5/sympy/polys/factortools.pyc in factor(f, *symbols, **flags) 78 return f 79 ---> 80 coeff, factors = poly_factors(f, *symbols, **flags) 81 82 res

[sympy] Re: Trigsimp error

2009-04-28 Thread Luke
Incidentally, trying trigsimp with recursive=True, deep=True, or both True didn't give the desired result. On Tue, Apr 28, 2009 at 12:22 AM, Luke wrote: > I'm writing some tests for some code that expresses a Vector expression in > the coordinates of a different frame. I

[sympy] Injecting symbols/variables into the workspace

2009-04-28 Thread Luke
nctions of time. Even better would be something that allows for syntax such as: >>> gc('q', 3) But has the same behavior as above. Is there a way to do this but *not* inject it into the *global* namespace, only the local one? ~Luke --~--~-~--~~~-

[sympy] Re: Trigsimp error

2009-04-28 Thread Luke
file the bug report, and I'll try to look at the paper Akshay mentioned. Akshay, if you need more examples that should simplify, I can provide you with a handful more :) ~Luke On Apr 28, 8:28 am, Akshay Srinivasan wrote: > I think trigsimp is too hack-ish. I'll try implementing

[sympy] Re: Trigsimp error

2009-04-28 Thread Luke
be what we want to implement. I don't have much time to start coding on it for a few weeks, but I could probably look at a few papers and let it soak in.... ~Luke On Apr 28, 12:35 pm, Ondrej Certik wrote: > On Tue, Apr 28, 2009 at 12:32 PM, Ondrej Certik wrote: > > On Tue, Apr

[sympy] solve() and diff() with Symbol versus Function

2009-05-11 Thread Luke
Traceback (most recent call last) /home/luke/Documents/PythonDynamics/ in () /var/lib/python-support/python2.6/sympy/solvers/solvers.pyc in solve(f, *symbols, **flags) 73 74 if any(not s.is_Symbol for s in symbols): ---> 75 raise TypeError('not a Symbol&

[sympy] Re: solve() and diff() with Symbol versus Function

2009-05-13 Thread Luke
I'll work on it in the next week or so, I think should be able to get something that does the job. ~Luke On May 11, 11:06 pm, Ondrej Certik wrote: > On Mon, May 11, 2009 at 9:23 PM, Luke wrote: > > > Would there be any reason that the following should not be implemented:

[sympy] trigsimp

2009-05-20 Thread Luke
to tackle this together, let me know and we could figure out a reasonable approach. Thanks, ~Luke --~--~-~--~~~---~--~~ You received this message because you are subscribed to the Google Groups "sympy" group. To post to this group, send emai

[sympy] Re: trigsimp

2009-05-20 Thread Luke
worth considering. Before investing in a lot of coding time, it would nice to be sure that a good algorithm is being used, although I must say their comparisons with the other popular packages out there seem favorable. ~Luke On May 20, 7:50 am, Alan Bromborsky wrote: > Luke wrote: > > Las

[sympy] Simple question

2009-05-20 Thread Luke
within ipython? Thanks, ~Luke --~--~-~--~~~---~--~~ 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 symp

[sympy] Poll on usage of solve()

2009-05-21 Thread Luke
able to solve equation(s) for both Symbol objects and Function objects. It makes the most sense to me that it should only support lists and tuples. ~Luke --~--~-~--~~~---~--~~ You received this message because you are subscribed to the Google Groups "sympy&

[sympy] Re: Poll on usage of solve()

2009-05-22 Thread Luke
Chris, You're right, who knows when that functionality may be useful. And it isn't a problem to deal with all three easily, so we should leave it in there. ~Luke On May 22, 3:32 am, smichr wrote: > On May 22, 4:23 am, Luke wrote:> Does anybody use > solve() by passing

[sympy] Enhancing functionality of diff()

2009-05-22 Thread Luke
1) in diff() of sympy/core/basic.py 2) in Derivative of sympy/core/function.py It seems to me option 1 is the right choice but I'd like to hear input on this. Maybe it doesn't matter or maybe it could another place? Thanks, ~Luke --~--~-~--~~~---~--~~ Yo

[sympy] Re: sympy vs matlab integrate, get 'final' result?

2009-05-22 Thread Luke
it doesn't like it, as you pointed out. I'm not sure why this would happen, it seems like integrate should parse each additive term and try to integrate it, if this were the case, the above example should work like Matlab. Additionally, it doesn't seem that the results are the sa

[sympy] Re: Motion equation generation script sharing

2009-05-23 Thread Luke
ics automatically is in place, but there is definitely more to do in regards to how to deal with constraints, definitions of generalized speeds, etc... What sort of systems do you use Kane's method to study? My own area of research is two wheeled vehicle (bicycle/motorcycle) dynamics and co

[sympy] Re: sympy vs matlab integrate, get 'final' result?

2009-05-23 Thread Luke
least in 2008a). I'll file it under the issues. ~Luke In [9]: res.subs({X: 2.}) Out[9]: 4*3**(1/2) + 14*I*3**(1/2) On May 23, 7:31 am, Oyster wrote: > that is a bad news :( > I have thought to use py+sympy in my research work, but think I'd > better turn to matlab now. > >

[sympy] Re: trigsimp

2009-05-25 Thread Luke
to the existing matchers tuple? I'm kind of thinking out loud right now and trying to figure out the next step to take ~Luke On May 21, 7:42 am, Akshay Srinivasan wrote: > Luke wrote: > > I tried both of those options and had no luck :( > > > Have you looked

[sympy] Behavior of trig functions

2009-05-25 Thread Luke
enough to recognize arguments of the following form: x+/- n*pi and then return the correct result based upon some rule based lookup? Or would it better to require the user to call trigsimp(sin(n*pi +/- x)) to return the simplified result? Thoughts? ~Luke

[sympy] .is_Number .is_number and .is_NumberSymbol attributes

2009-05-25 Thread Luke
o me that if something is NumberSymbol, it should also be a Number (or number), but this isn't how pi currently is. I'm guessing this part is a work in progress with the assumptions system, but I'm not sure. ~Luke --~--~-~--~~~---~--~~ You received this message

[sympy] Re: .is_Number .is_number and .is_NumberSymbol attributes

2009-05-25 Thread Luke
Ok, that makes more sense now. Thanks for the clarification. ~Luke On May 25, 1:05 pm, Abderrahim Kitouni wrote: > On Mon, 25 May 2009 12:48:09 -0700 (PDT)Luke wrote: > > > 1)  Why are there both 'is_Number' and a 'is_number' methods, and why >

[sympy] Re: trigsimp

2009-05-25 Thread Luke
ere another simpler approach, or am I missing something fundamental here? Thanks, ~Luke On May 25, 4:07 pm, "Aaron S. Meurer" wrote: > On May 25, 2009, at 11:56 AM, Luke wrote: > > > > > > > Here is the link for the Maxima trigsimp() code.  It was written in &

  1   2   >