[sympy] Identifying repeated subexpressions in systems of equations

2008-06-13 Thread Luke
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@googlegroups.com

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

2008-06-16 Thread Luke
I'll get started on this over the summer. Thanks, ~Luke On Jun 15, 11:28 pm, Robert Kern [EMAIL PROTECTED] wrote: On Sun, Jun 15, 2008 at 17:46, Ondrej Certik [EMAIL PROTECTED] wrote: Anyway, I created an issue for this: http://code.google.com/p/sympy/issues/detail?id=891 I've attached

[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 [EMAIL PROTECTED] wrote

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

2009-01-16 Thread Luke
, 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 --~--~-~--~~~---~--~~ You received this message because you

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

2009-01-17 Thread Luke
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 abro...@verizon.net wrote: Luke wrote

[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
the ._args list for the above example then be: v1._args == [5*A[1], 6*A[2]] ? Thanks, ~Luke -- Forwarded message -- From: Ondrej Certik ond...@certik.cz Date: Thu, Apr 16, 2009 at 9:37 PM Subject: Re: Confused about something To: Luke hazelnu...@gmail.com On Thu, Apr 16, 2009 at 8:00

[sympy] Printing of functions

2009-04-17 Thread Luke
, ~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

[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 hazelnu...@gmail.com wrote: I'm writing some tests for some code that expresses a Vector expression in the coordinates of a different frame. I am

[sympy] Injecting symbols/variables into the workspace

2009-04-28 Thread Luke
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 --~--~-~--~~~---~--~~ You received this message because you are subscribed

[sympy] Re: Trigsimp error

2009-04-28 Thread Luke
, 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 akshaysriniva...@gmail.com wrote: I think trigsimp is too hack-ish. I'll try implementing the algorithm

[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 ond...@certik.cz wrote: On Tue, Apr 28, 2009 at 12:32 PM, Ondrej Certik ond...@certik.cz wrote

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

2009-05-11 Thread Luke
(most recent call last) /home/luke/Documents/PythonDynamics/ipython console in module() /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') 76

[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 ond...@certik.cz wrote: On Mon, May 11, 2009 at 9:23 PM, Luke hazelnu...@gmail.com wrote: Would there be any reason that the following should

[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 email to sympy

[sympy] Poll on usage of solve()

2009-05-21 Thread Luke
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 group

[sympy] Enhancing functionality of diff()

2009-05-22 Thread Luke
) 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 --~--~-~--~~~---~--~~ You received

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

2009-05-22 Thread Luke
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 same, so there must be some sort of bug in integrate. ~Luke On May 22, 2:27 am, Oyster lepto.pyt

[sympy] Re: Motion equation generation script sharing

2009-05-23 Thread Luke
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 control. ~Luke On May 23

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

2009-05-23 Thread Luke
}) Out[7]: 4*3**(1/2) + 14*I*3**(1/2) In [8]: 14.0*3.0**(1./2.) Out[8]: 24.248711305964282 [/sympy] Note the real part of result is nonzero, and that the imaginary part is different from the result that Matlab gives (really Maple under the hood, at least in 2008a). I'll file it under the issues. ~Luke

[sympy] Re: trigsimp

2009-05-25 Thread Luke
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 akshaysriniva...@gmail.com wrote: Luke wrote: I tried both of those options and had no luck :( Have you looked at the paper by Fu, Zhong, and Zeng: http://vv.cn/d

[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
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 because you are subscribed to the Google

[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 a.kito...@gmail.com wrote: On Mon, 25 May 2009 12:48:09 -0700 (PDT)Luke hazelnu...@gmail.com wrote: 1)  Why are there both 'is_Number' and a 'is_number' methods, and why do

[sympy] Re: trigsimp

2009-05-25 Thread Luke
? Thanks, ~Luke On May 25, 4:07 pm, Aaron S. Meurer asmeu...@gmail.com wrote: On May 25, 2009, at 11:56 AM, Luke wrote: Here is the link for the Maxima trigsimp() code.  It was written in 1981, according to the comments!!! http://maxima.cvs.sourceforge.net/viewvc/maxima/maxima/share/trigonom

[sympy] Re: trigsimp

2009-05-26 Thread Luke
.match(a*pi + b) {b_: -x, a_: 1}   e3.match(a*pi + b) {b_: x, a_: -1}   e4.match(a*pi + b) {b_: -x, a_: -1} just as you expected. Aaron Meurer On May 25, 2009, at 7:46 PM, Luke wrote: Aaron, Thanks for the clarification.  I think I get the idea, but I'm having trouble matching

[sympy] Re: Strange simplification

2009-05-26 Thread Luke
I get the same behavior on my machine. It seems like a bug to me. Not sure if the bug is in integrate or in simplify, but could you file this in the issues? Thanks, ~Luke On May 26, 12:16 pm, Neal Becker ndbeck...@gmail.com wrote: Out[80]: (-2*s**2*w**2 + w**4)/(s**4 + w**4) In [81

[sympy] Re: trigsimp

2009-05-26 Thread Luke
Symbol? ~Luke On May 26, 9:33 am, Luke hazelnu...@gmail.com wrote: Thanks. On May 25, 9:02 pm, Aaron S. Meurer asmeu...@gmail.com wrote: You need to use the exclude option I was telling you about.  Do,   a = Wild('a', exclude=[pi]) and   b = Wild('a', exclude=[pi]). and you get   e1

[sympy] Re: Strange simplification

2009-05-26 Thread Luke
))) --- NameError Traceback (most recent call last) /home/luke/lib/python/sympy/ipython console in module() NameError: name 'infty' is not defined In [6]: simplify(integrate(ia, (f, 0, oo))) Out[6]: 0 In [10]: oo

[sympy] Re: Strange simplification

2009-05-26 Thread Luke
[2] *Pi*w*x + 4*Pi^2*x^2]))/ (8*Sqrt[2]*Pi) They don't let you do definite integrals there. And the computation timed out on Wolfram Alpha. Anybody actually *know* what this integral should be? ~Luke On May 26, 5:25 pm, Robert Kern robert.k...@gmail.com wrote: On Tue, May 26, 2009 at 19:22

[sympy] Function vs. Symbol

2009-05-26 Thread Luke
the functionality I just described, but hopefully this would still allow it to retain the ability to work with all the builtin Sympy functions that work best with Symbol objects. Thoughts? ~Luke --~--~-~--~~~---~--~~ You received this message because you are subscribed

[sympy] Re: Function vs. Symbol

2009-05-27 Thread Luke
). But maybe it could still be done without subclassing? The other approach would be to make everything in Sympy work equally well with Function, Derivative, and Symbol. I don't know if this is a good idea or not, but it also seems like it would be a ton of work. What do you think? ~Luke W

[sympy] Re: Function vs. Symbol

2009-05-27 Thread Luke
extremely easy and for the stuff I'm doing, would let one focus on the problem, not on the syntax. Thoughts? ~Luke On May 27, 11:53 am, Luke hazelnu...@gmail.com wrote: Fabian,   I think the example you gave is good, but I think it would be better if you could imply that x == x(t) upon

[sympy] Re: Function vs. Symbol

2009-05-27 Thread Luke
the second approach for the auto substitution stuff. ~Luke On May 27, 4:49 pm, Ondrej Certik ond...@certik.cz wrote: We discuss this on IRC with Luke and Fabian. Now I understand -- Luke wants the result of differentiating not to be instances of Derivative() class, but rather some other symbols, e.g

[sympy] Re: Function vs. Symbol

2009-05-28 Thread Luke
I'm a little unclear about a few things with regards to how to properly subclass StrPrinter. On May 27, 5:08 pm, Ondrej Certik ond...@certik.cz wrote: On Wed, May 27, 2009 at 6:05 PM, Luke hazelnu...@gmail.com wrote: I like the first way for the fact that it just has 'x' instead of 'x (t

[sympy] Re: Function vs. Symbol

2009-05-28 Thread Luke
, as well as the code inside their cusstom class, and any other relevant code as a complete example of how to properly customize the printing? Thanks, ~Luke On May 28, 3:23 pm, Luke hazelnu...@gmail.com wrote: I'm a little unclear about a few things with regards to how to properly subclassStrPrinter

[sympy] Re: Function vs. Symbol

2009-05-28 Thread Luke
/printer.py, line 132, in _print res = getattr(expr, self.printmethod)() File /home/luke/lib/python/pydy/pydy.py, line 266, in _sympystr_ xsym('*') + k._sympystr_() UnicodeDecodeError: 'ascii' codec can't decode byte 0xe2 in position 11: ordinal not in range(128) when I go to ipython

[sympy] Re: Function vs. Symbol

2009-05-28 Thread Luke
to the point of understand how it works, but I'm not there yet. Thanks, ~Luke On May 28, 4:00 pm, Sebastian basti...@gmail.com wrote: Hi Luke, I think all your questions are answered in the docstring of printing/printer.py. There it tells you in which order it is tried to print an object. 1) Let

[sympy] Re: Function vs. Symbol

2009-05-28 Thread Luke
() It gives me an error because of the xsym('*') part: s += print_pydy(e.dict[k]) + xsym('*') + k.__str__() UnicodeDecodeError: 'ascii' codec can't decode byte 0xe2 in position 11: ordinal not in range(128) Any suggestions on how to fix this? Thanks, ~Luke On May 28, 4:35 pm, Luke hazelnu

[sympy] Re: Function vs. Symbol

2009-05-28 Thread Luke
the UnicodeDecodeError. I know my terminal is unicode enabled, the error is something else, I just don't know how to fix it. ~Luke On May 28, 6:34 pm, Sebastian basti...@gmail.com wrote: Luke wrote: Ok, thanks. Any ideas regarding the xsym(*) error? Not really, I don't even know where this xsym comes from

[sympy] UnicodeDecodeError when using xsym(*)

2009-05-28 Thread Luke
: And here is the error I'm getting: File t.py, line 56, in module print 'print A[1]', Vector((1+sin(q1))*A[1]) File /home/luke/lib/python/pydy/pydy.py, line 254, in __str__ return print_pydy(self) File /home/luke/lib/python/pydy/pydy.py, line 1157, in print_pydy return pp.doprint

[sympy] cos(x - y) == cos(y - x)

2009-05-29 Thread Luke
=True) In [5]: cos(x - y) == cos(y - x) Out[5]: True In [10]: x = Function('x')(z) In [11]: y = Function('y')(z) In [12]: cos(x - y) == cos(y - x) Out[12]: False In [13]: Fabian, are your assumptions going to have an is_even_function and is_odd_function attributes for Function? ~Luke

[sympy] Re: cos(x - y) == cos(y - x)

2009-05-29 Thread Luke
that is supposed to take care of normalizing cos(-x-y) to cos(x+y), but obviously the logic in there isn't catching this. ~Luke On May 29, 9:06 am, Aaron S. Meurer asmeu...@gmail.com wrote: On May 29, 2009, at 9:25 AM, Ondrej Certik wrote: On Fri, May 29, 2009 at 2:24 AM, Luke hazelnu

[sympy] Re: Motion equation generation script sharing

2009-06-01 Thread Luke
and Lasenby? Or is Doran and Lasenby available electronically? Thanks, ~Luke On Jun 1, 2:25 pm, Alan Bromborsky abro...@verizon.net wrote: Luke wrote: If you are familiar with Autolev, then you are probably familiar with some of the behavior that makes it very convenient and easy to derive

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

2009-06-18 Thread Luke
), won't get screened because currently (pi/2).is_Number==False, but (pi/2).is_number==True. Or should I keep both sets of if statements? I want to make all the trig functions have consistent behavior, so I just need some input on how this should be designed. Thanks, ~Luke On May 25, 9:15 pm

[sympy] Complex Infinity

2009-06-22 Thread Luke
functions all consistent. 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 sympy

[sympy] Re: Complex Infinity

2009-06-24 Thread Luke
Is there a reason it returns complex infinity versus just infinity? Does it have to do with the assumptions about the variables? Does anybody know an example where Mathematica returns just regular infinity? ~Luke On Jun 23, 10:23 am, Ondrej Certik ond...@certik.cz wrote: 2009/6/23 Roberto

[sympy] Re: Complex Infinity

2009-06-24 Thread Luke
ComplexInfinity, but I'm still not understanding why this behavior is more desirable than regular infinity. I'm fine with implementing it this way, but it would be nice to understand why this way is more correct or general, if indeed it is. Thanks, ~Luke On Jun 24, 8:27 am, Luke hazelnu...@gmail.com wrote

[sympy] How to properly subclass the pretty printer

2009-06-25 Thread Luke
it? 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 sympy+unsubscr...@googlegroups.com

[sympy] Re: How to properly subclass the pretty printer

2009-06-25 Thread Luke
) MyBreakfast = GreenEggsAndHam('I LOVE SYMPY') print_GreenEggsAndHam(MyBreakfast) pprint(MyBreakfast) ~Luke On Jun 25, 11:58 am, Ondrej Certik ond...@certik.cz wrote: On Thu, Jun 25, 2009 at 12:50 PM, Lukehazelnu...@gmail.com wrote: I'm trying to subclass the StrPrinter and PrettyPrinter

[sympy] Re: Solving equation sets in Mathematica

2009-06-26 Thread Luke
I would recommend looking at scipy.optimize.fsolve(), as Alan suggested. If you have a decent initial guess, it should work. ~Luke On Jun 26, 9:43 am, Alan Bromborsky abro...@verizon.net wrote: nandan jha wrote: Hello I am trying to solve a set of non-linear equations in Mathematica 7.0

[sympy] Re: who still uses python2.4

2009-06-26 Thread Luke
I use 2.6 and have never used 2.4 since I started using python 2 years ago. Breaking sympy compatibility with python 2.4 is ok by me. It seems like the users of sympy are probably the type who can update more frequently than once every 3 years, so it doesn't seem like too much of an issue. On

[sympy] Re: Difierentiating With Respect to a Function

2009-06-28 Thread Luke
the most logical, we just need to ensure we can make the chain rule work, so this behavior is preserved: In[1]: sin(x).diff(t) Out[2]: x'*cos(x) ~Luke On Jun 28, 8:53 am, Aaron S. Meurer asmeu...@gmail.com wrote: The routines in solve() for solving for a function or a derivative   could

[sympy] Unicode combining diacritical characters

2009-07-02 Thread Luke
terminal, or the font I'm using in my terminal, doesn't support that particular character. I am using Konsole in Kubuntu 9.04, with the character encoding set to: Unicode-- UTF-8. Can anybody else get the 'COMBINING DOT ABOVE' character to work in their terminals? ~Luke

[sympy] Re: Working with sets

2009-07-03 Thread Luke
couldn't make an interval like: [2, 1]? ~Luke On Jul 3, 10:44 am, Christophe projet...@gmail.com wrote: Hello, I would like to know if it is possible to work with sets. A second question. Is the a way to do something like   intersection([1 ; 2];[1.5 ; 5])   and   union([1 ; 2];[1.5 ; 5])   where

[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 abro...@verizon.net wrote: Ondrej

[sympy] Automatic distribution of numbers versus symbols

2009-07-28 Thread Luke
or not distribution occurs automatically, and 2) that if distribution does occur automatically, it should be easily reversible, i.e. through something like factor Is there a way to keep 4*(x-y) as 4*(x-y) instead of -4*y+4*x? Thoughts? ~Luke

[sympy] Add combine question

2009-08-14 Thread Luke
'. 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 subscribed

[sympy] Re: Add combine question

2009-08-15 Thread Luke
) 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 hazelnu...@gmail.com wrote

[sympy] Issue 1616

2009-08-25 Thread Luke
be corrupted or invalid The error message is: ('EOF in multi-line statement', (96, 0)) --- ValueErrorTraceback (most recent call last) /home/luke/lib/python/pydy/examples/ipython console

[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 hazelnu...@gmail.com 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

[sympy] New code in bicycle branch

2009-09-24 Thread Luke
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 expression. ~Luke

[sympy] Symbolic matrix inversion

2009-09-29 Thread Luke
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 received this message because you are subscribed

[sympy] Re: Symbolic matrix inversion

2009-09-29 Thread Luke
On Sep 29, 1:09 pm, Ondrej Certik ond...@certik.cz wrote: On Tue, Sep 29, 2009 at 12:49 PM, Luke hazelnu...@gmail.com 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

[sympy] Re: Symbolic matrix inversion

2009-09-29 Thread Luke
formulation would be as generic as possible. ~Luke On Sep 29, 4:15 pm, Alan Bromborsky abro...@verizon.net wrote: Luke wrote: On Sep 29, 1:09 pm, Ondrej Certik ond...@certik.cz wrote: On Tue, Sep 29, 2009 at 12:49 PM, Luke hazelnu...@gmail.com wrote: I'm using Sympy from within PyDy

[sympy] Re: Symbolic matrix inversion

2009-09-29 Thread Luke
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 abro...@verizon.net wrote: Alan Bromborsky wrote: Luke wrote: On Sep 29, 1:09 pm, Ondrej Certik ond...@certik.cz 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 hazelnu...@gmail.com wrote: Yes, this is something I should look into.  I am pretty sure that the Netlib codes have this functionality

[sympy] Re: Symbolic matrix inversion

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

[sympy] Re: Symbolic matrix inversion

2009-09-30 Thread Luke
be 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 Sep 29, 7:34

[sympy] Re: Symbolic matrix inversion

2009-09-30 Thread Luke
On Sep 30, 8:39 am, Luke hazelnu...@gmail.com 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 do.  I

[sympy] Preventing expansion

2009-10-11 Thread Luke
, 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 unsubscribe from this group

[sympy] GSoC project mentoring

2011-04-05 Thread Luke
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 group. To post to this group, send email

[sympy] Re: GSoC project mentoring

2011-04-05 Thread Luke
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 asmeu...@gmail.com wrote: I agree with Ondrej. Since PyDy is (still at least) separate form SymPy, another idea regarding GSoC would

[sympy] class design / modelling

2011-04-12 Thread Luke
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 options, visit this group

Re: [sympy] Question regarding vectors

2011-05-09 Thread Luke
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 gilbertg...@gmail.com wrote: In PyDy (which we plan to merge into SymPy.physics.classical this summer) Vector is one of the classes already implemented (along

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 abro

Re: [sympy] Question regarding vectors

2011-05-10 Thread Luke
be standalone classes 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 abro...@verizon.net wrote: The assumption is that the expression multiplying the vector

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 vinzent.steinb...@googlemail.com wrote: On 10 Mai, 03:30, Ronan Lamy ronan.l...@gmail.com wrote: Please, try to make the interface dot(v1

Re: [sympy] Re: Who uses SymPy?

2011-05-12 Thread Luke
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:32 AM, Sherjil Ozair

Re: [sympy] Re: Question regarding vectors

2011-05-14 Thread Luke
dynamics. ~Luke On Fri, May 13, 2011 at 6:29 PM, Gilbert gede gilbertg...@gmail.com 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 extend any other sympy classes.  It will use sympy's matrix class.  The current

[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 -

[sympy] Re: Problem with user-defined functions

2011-06-02 Thread luke
') # one single query print a 3*SUM('field') print N(a) 1234 just like other functions work. e.g. log(10)+log(10)+log(10) 3*log(10) On Jun 2, 9:38 am, Mateusz Paprocki matt...@gmail.com wrote: Hi, On 2 June 2011 09:07, luke luca.giaco...@gmail.com wrote: Hi eveyone, I'm

[sympy] Re: Problem with user-defined functions

2011-06-02 Thread luke
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 ronan.l...@gmail.com wrote: Le jeudi 02 juin 2011 à 09:35 -0700, luke a écrit : Yes, that worked. But I had to restructure

[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 luca.giaco...@gmail.com wrote: Actually you're wrong. Every instance of a class in python has its own

[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 property and not a symbol :) I'll put back my eval since it worked. Hope that was of any help! On Jun 2, 9:32 pm, luke luca.giaco...@gmail.com wrote: WAIT, I didn't see that you use self.arg[0]!! that's why my code

[sympy] Re: Problem with user-defined functions

2011-06-02 Thread luke
; } }); return {%(field)s:sum}; }; % {'field':cls.args[0]}) result = db.people.map_reduce(map, reduce, myresults) return result.find_one()['value'][unicode(cls.args[0])] On Jun 2, 9:36 pm, luke luca.giaco...@gmail.com

[sympy] Adding sympy.physics to documentation

2011-06-07 Thread Luke
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/sympy/tree

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

2011-06-07 Thread Luke
when building the docs. I'll get started on this. ~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

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

2011-06-07 Thread Luke
when I get this up and running on a github branch. ~Luke Stefan On 7 June 2011 20:34, Luke hazelnu...@gmail.com wrote: 1)  Should any of the code in sympy.physics be bundled into it's own submodule (similar to quantum and/or classical)? You could have either a 'physics' submodule

[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(self, other) File /home/luke/repos

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 elliso...@gmail.com wrote: On Tue, Jun 7, 2011 at 10:28 PM, Luke hazelnu...@gmail.com wrote: On Tue, Jun 7, 2011 at 1:26 PM, krastanov.ste...@gmail.com krastanov.ste...@gmail.com wrote: What would

[sympy] sympy/physics sphinx documentation

2011-06-27 Thread Luke
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 options, visit this group

Re: [sympy] Question on Interfaces

2011-08-11 Thread Luke
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 gilbertg...@gmail.com wrote: My Summer of Code project is writing a submodule

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

2011-08-11 Thread Luke
with the rest of the code. Jason On Thu, Aug 11, 2011 at 12:42 PM, Luke hazelnu...@gmail.com wrote: I ran this same question by my girlfriend who teaches undergraduate physics classes.  She isn't an experienced programmer, so concepts of operator overloading and object oriented vs functional

[sympy] MathJax script tag

2011-08-16 Thread Luke
doesn't make sense with the result that all the LaTeX code appears undrendered. I am using Sphinx 1.0.1, and Python 2.7.1, and I simply 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

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 asmeu...@gmail.com 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   script

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

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

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

2011-08-22 Thread Luke
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 asmeu...@gmail.com wrote: Well, I don't know why you are doing things this way, but clearly python -m doctest file.py or ./bin/doctest

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

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

[sympy] Import errors with an installed version of Sympy

2011-10-16 Thread Luke
', and 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 information. from sympy import * Traceback (most

  1   2   >