itive = Symbol('x', positive=True)
> > >>> diff(x**2, x_positive)
> > 0
> > >>> x == x_positive
> > False
> >
> > If you use assumptions, it is good practice to always use the same
> > assumptions for any given symbol name. It also help
so I can't reproduce this half-remembered claimed behaviour.)
On Wednesday, June 3, 2020 at 11:30:31 PM UTC+1, Oscar wrote:
>
> >> >> On Wed, Jun 3, 2020 at 12:53 PM James Bateman
> wrote:
> >> >>>
> >> >>> I've just discovered a bug in my code
>>
>> Out[4]: int
>>
>> In [5]: a == b
>>
>> Out[5]: True
>>
>> In [6]: a is b
>>
>
;>
> >> In [1]: a = 1
> >>
> >> In [2]: b = 1
> >>
> >> In [3]: type(a)
> >> Out[3]: int
> >>
> >> In [4]: type(b)
> >> Out[4]: int
> >>
> >> In [5]: a == b
> >> Out[5]: True
> &g
I've just discovered a bug in my code which boiled down to the following,
where a symbol "y" was given the same SymPy name as an existing symbol.
import sympy as sp
x = sp.Symbol('x')
y = sp.Symbol('y')
x == y # True
x is y # True; expected False
x + y # 2*x; expected x + x (which would have
ows).
>
> Probably there is a tutorial somewhere that explains how to use
> google-cloud-sdk more generally. It might be worth going through
> something like that before trying to get sympy-live to work.
>
>
> On Mon, 20 Jan 2020 at 05:29, James . wrote:
> >
> > if yo
if you dont mind can u pls tell me in specific
while running this command(../../google-cloud-sdk/bin/dev_appserver.py .)
its showing like this "'..' is not recognized as an internal or external
command,
operable program or batch file.""
my cloud sdk in this path "C:\Users
github.com/sympy/sympy-live/pull/139
> After merging that PR sympy live will not work locally but will work
> when pushed to google.
>
>
> On Sun, 19 Jan 2020 at 12:59, James . wrote:
> >
> > Hi everyone,
> > If anyone using sympy live shell on their local ma
from 2 days but still
couldn't solve .
On Sat, Jan 18, 2020, 01:21 James . wrote:
> working now above things, but at this link
> https://github.com/sympy/sympy-live this command "
> ../google_appengine/dev_appserver.py"
> giving error
>
--
You received this message b
working now above things, but at this link
https://github.com/sympy/sympy-live this command "
../google_appengine/dev_appserver.py"
giving error
--
You received this message because you are subscribed to the Google Groups
"sympy" group.
To unsubscribe from this group and stop receiving emails
where is this "
dev_appserver.py
"
i cant find it anywhere
On Fri, Jan 17, 2020 at 11:08 PM James . wrote:
> just now added ssh key,still same error coming..
>
>
> On Fri, Jan 17, 2020 at 10:38 PM Aaron Meurer wrote:
>
>> I think you will need to add yo
>
> On Fri, Jan 17, 2020 at 6:41 AM James . wrote:
> >
> > https://github.com/sympy/sympy-live/blob/master/README.rst
> > from this i have followed every step,but im not able to see the gui of
> sympy live shell in my local server,i think its because of this command "
""
git submodule update
Cloning into 'C:/Users/James/AppData/Local/Google/Cloud
SDK/sympy-live/static'...
Host key verification failed.
fatal: Could not read from remote repository.
Please make sure you have the correct access rights
and the repository exists.
fatal: clone of 'g...@gith
rom the commit history of the SymPy Live
>> repo.
>>
>> The advantage of the App Engine is that it does autoscaling and
>> automatic sandboxing.
>>
>> Aaron Meurer
>>
>> On Tue, Jan 14, 2020 at 11:32 AM James .
>> wrote:
>>
Hi,My name is James, i am a 3rd year undergraduate student from India. I am
using python from 3 years, from a couple of months ago I started using
sympy and I am very much interested in refactoring sympy live.can someone
help me what are the prerequisites I need to do this?
On Tue, Jan 14, 2020
Is there a way to write a replacement rule for a function f with an
arbitrary number of arguments that makes it linear in all its arguments?
An example for when f has three arguments:
1) f( x1+x4 , x2 , x3 ) = f(x4,x2,x3) + f(x1,x2,x3)
2) f( x1 , x2+x4 , x3 ) = f(x1,x2,x3) + f(x1,x4,x3)
3) f( x1
6 for sol2s in sol2r: 7897
sol1 = solve(Integral(1/F.subs(v, sol2s.rhs), u).doit() - t - C2, u) 7898
sol = []
TypeError: 'Equality' object is not iterable
As far as I can tell I setup the problem correctly, does anyone have any
advice?
James
--
You received this message because
To kind of expand on what Jason's saying a 3D pendulum can be completely
defined using just (x, y, z) and you can deduce the angles from these
coordinates. In your case the pendulum only has two degrees of freedom (x
and y for instance and z be calculated because the pendulum has a fixed
Hi all,
I am not seeing an up to date version of the docs when I visit
docs.sympy.org. I have tried the /latest and /dev versions an neither seems
to work. The docs appear to not have been updated in awhile either as the
change I'm looking at was merged on April 2nd in PR #10878
Hi all,
I will be at the conference next week and I need a place to stay the night
of July 17th. I was wondering if any of the other attendees would mind me
crashing of the floor where ever they are staying?
Thanks,
Brandon
--
You received this message because you are subscribed to the
Theoretically you should be able to transform the results of Lagrange's
>> method to those of Kane's. I'm not sure if that procedure is laid out
>> anywhere in the literature. But it may have some complications.
>>
>> The f's are simply vector equations of those variable
These statements are found in the Kane's method and Lagrange's method docs
and are seemingly contradictory
"In mechanics we are assuming there are 5 basic sets of equations needed to
describe a system."
"In mechanics we are assuming there are 3 basic sets of equations needed to
describe a
Sorry, here is the url https://github.com/sympy/planet-sympy/pull/38.
Brandon
On Friday, May 6, 2016 at 11:20:04 AM UTC-5, James Milam wrote:
>
> Here is the link to PR url adding my blog to Planet Sympy
>
> Thanks,
> Brandon
>
> On Monday, May 2, 2016 at 11:47:12 AM UTC-5
Here is the link to PR url adding my blog to Planet Sympy
Thanks,
Brandon
On Monday, May 2, 2016 at 11:47:12 AM UTC-5, Ondřej Čertík wrote:
>
> Hi,
>
> For GSoC students:
>
> Please send a PR against https://github.com/sympy/planet-sympy, adding
> your blog.
>
> Please reply to this thread
530-601-9791
>
> On Tue, Apr 26, 2016 at 6:56 PM, James Goppert <james@gmail.com
> > wrote:
>
>> Hi Jason,
>>
>> Thanks for you interest:
>>
>> I think Modelica is a nicer format for the core mathematical model of
>> dynamic systems since it
ple will find it useful.
On Tuesday, April 26, 2016 at 6:02:55 PM UTC-4, James Goppert wrote:
>
> I'm working on a Modelica compiler written in pure python that outputs a
> sympy model. Please let me know if anyone is interested in contributing.
>
> Please see the jupyter notebook
eneral-Ideas
Thanks,
Brandon
On Sunday, April 10, 2016 at 8:14:29 PM UTC-4, James Milam wrote:
>
> Hi all,
>
> I was looking through the wiki for an ideas page where I could post some
> thoughts on a contribution idea that I do not have the time to implement.
> While l
Hi all,
I was looking through the wiki for an ideas page where I could post some
thoughts on a contribution idea that I do not have the time to implement.
While looking, I found a couple of idea pages but they were very old and
did not seem to be relevant. I was wondering if a general ideas
Hi Jason,
I've updated the application and am planning on doing final revisions and
submission tomorrow and was wondering if you or anyone else would be able
to look at it and provide feedback one last time.
Thank you for all of your assistance
--
You received this message because you are
I believe it is just code to demonstrate your work process, however, if you
can relate it to your project you'll have a better understanding of the
relevant code that you will need to talk about in your proposal.
--
You received this message because you are subscribed to the Google Groups
Hi Jason,
I was looking at also producing a Newton-Euler method for equations of
motion generation but I saw an edit you made on another proposal saying to
focus on the base class and speed ups OR base class and Newton-Euler method
so that the scope of the project is not too big. I was
Last I've heard mention of that project is in this message thread
https://groups.google.com/forum/?fromgroups=#!topic/sympy/KFyfC4gQUSQ
On Monday, March 14, 2016 at 5:21:01 PM UTC-4, Prince Kumar wrote:
>
> Can anyone please give a detailed information on the current status of
> project
I apologize for the previous post needing a scroll bar to the left. I
honestly have no idea how I turned off word wrap. I just posted a rough
draft for this project to wiki. Any general comments/revisions would be
greatly appreciated, especially regarding the scope/timeline for the
project.
I have gone through the example in the linked thread and have been able to
get snakeviz to create a profile when running
pydy-code-gen/pydy_code_gen/model's
generate_n_link_pendulum_on_cart_equations_of_motion()
function for a 10 link system. I do not have any experience working in C++
so I
Hi all,
I am currently pursuing a master's degree in mechanical engineering at the
University of Florida and am interested in participating in GSoC with the
sympy project. I have recently completed a course on analytical dynamics
which lead to my interest in working with the equations of
Hello GSoC Students (and mentors!),
In the next few days/weeks, GSoC pull requests will start coming in. This
is really exciting! After the pull requests come in, it's important to
start reviewing these in a timely fashion. While each student has a
mentor/mentors who will work with them
Dear all,
The matrix in sympy save as dictionary. I want to save a big matrix with
size 3x3. Is there any method to save it to hard-disk driver? I
tried with numpy but
it cannot work.
Thank you,
--
You received this message because you are subscribed to the Google Groups
sympy group.
On Friday, April 17, 2015 at 8:17:56 AM UTC-5, Harsh Gupta wrote:
Great initiative, I suggest adding `needs decision` to the classifier tag.
It is for the issues for which we are not sure if they are valid or not.
`Needs Decision` already exists, and is a special (pink) tag. The reason
for
wrote:
Glad to see you taking this on. Quite a few issues are either
duplicate or already fixed, so there is definitely some cleanup
possible.
On Fri, Apr 17, 2015 at 12:40 AM, James Crist crist...@umn.edu
wrote:
As of this writing, SymPy has 1648 issues open. That's more than
numpy
On Friday, April 17, 2015 at 1:34:55 PM UTC-5, Joachim Durchholz wrote:
- Help labeling.
I haven't found a way to enable tagging for people that do not have full
write access to the repository. I suspect GitHub does not offer any.
This is also the reason why tagging does not happen very
making the PR. I'd
like to remove it unless someone makes a strong argument to the contrary.
Perhaps, we can reword it. How about Needs changes or Waiting for
changes?
Sudhanshu Mishra
On Sat, Apr 18, 2015 at 12:35 AM, James Crist cris...@umn.edu
javascript: wrote:
On Friday, April 17
...@durchholz.org wrote:
One thing I saw proposed elsewhere is to group labels, by naming them
group:label.
Am 17.04.2015 um 07:40 schrieb James Crist:
*Submodule tags (html #FF, blue):*
Everything after `sympy.` for the specific submodule. Keep the naming and
casing consistent with the sympy
Do people still use this?
Many recent beginner PRs have been tagged with it. If no one is attached,
I'd like to remove it.
- Jim
On Fri, Apr 17, 2015 at 6:28 PM, Aaron Meurer asmeu...@gmail.com wrote:
On Fri, Apr 17, 2015 at 2:05 PM, James Crist crist...@umn.edu wrote:
On Friday, April 17
As of this writing, SymPy has 1648 issues open. That's more than numpy,
scipy, or pandas (more than numpy and scipy combined!). Further, our issue
tagging system is a mess. We can do better than this!
Many of these issues are imported from google code, and may be already
fixed in master (some
I expect this to turn into a flamewar, so I'm going to ask everyone to only
give one opinion (voting style). Please please please don't fight about
this.
I just got back from PyCon, and there was a presentation there on pep 8
formatting, and one of the main points was that pure formatting PRs
, Joachim Durchholz j...@durchholz.org wrote:
Am 14.04.2015 um 19:54 schrieb James Crist:
For example, say I make a tiny bug fix in function foo - I could also
clean
up some of the code in foo. That way the last person to touch foo is not
someone who added a space between an operator, but someone
to be
the main symbolics engine for sage, as the increased user base would
probably result in an increased development team.
On Tue, Apr 14, 2015 at 4:31 PM, James Crist crist...@umn.edu wrote:
Have you tried git blame -CCC -M? According to the manpage it is
supposed to be smarter about stuff.
I had
.
On Thu, Apr 2, 2015 at 11:54 AM, Joachim Durchholz j...@durchholz.org wrote:
Am 02.04.2015 um 18:46 schrieb James Crist:
Performance:
Using symbols() in all contexts might have performance ramifications,
creating new Symbol() objects means more memory pressure than reusing
precreated symbols
Performance:
Using symbols() in all contexts might have performance ramifications,
creating new Symbol() objects means more memory pressure than reusing
precreated symbols from sympy.abc (which happen 521 times in SymPy
itself, hopefully just in test code).
We cache symbol creation,
One way that's pretty common is to do something like:
result = whatever your testing
solution = something you type in
assert simplify(result - solution) == 0
This is better than constructing a special case form, because that form may
change in the future.
On Saturday, March 28, 2015 at 9:28:00
Clearly Jim needs to move out to the West Coast.
I'm working on it :)
I would suggest deciding on a technology beforehand
If we're hosting our docs at readthedocs (as Jason said on the wiki), then
website will be a gallery, blog, some prose pages, and that's it. Any of
the popular static
`lambdify` is intended for numeric evaluation (but can be made to evaluate
symbolically). By default, the functions `sin`, `cos`, etc... are pulled
from `math` or `numpy`, which expect floats (and if not given a float, they
attempt to convert to float).
What's happening here, is you're passing
Isn't this similar to SymPy Gamma? http://gamma.sympy.org/
On Wednesday, January 7, 2015 5:36:57 PM UTC-6, Anas Belkhadir wrote:
did you understand me what i want to say ??
2015-01-07 21:28 GMT+00:00 Francesco Bonazzi franz@gmail.com
javascript::
On Wednesday, January 7, 2015
Bleagh, I still need to fix this, completely forgot about it. Somewhere in
all the changes made in my GSoC project, this test stopped working. Due to
the large size of the expression involved (just printing it takes forever),
debugging the failure is tricky. Currently it's left with the
Hello,
Read this: https://github.com/sympy/sympy/wiki/introduction-to-contributing
On Friday, December 19, 2014 10:39:56 AM UTC-6, GAURAV MISHRA wrote:
Hi everyone
I am Gaurav and i am pursuing my B.Tech in Computer Science (3 year) from
IIIT Hyderabad (India).I am a Mathematics and
rules in a ruleset is small (~none for non-AC patterns), but
it does exist.
- Jim Crist
On Sat, Nov 29, 2014 at 10:29 AM, Aaron Meurer asmeu...@gmail.com wrote:
On Fri, Nov 28, 2014 at 6:40 PM, Richard Fateman fate...@gmail.com
wrote:
On Thursday, November 27, 2014 7:49:30 PM UTC-8, James
predicates for each path will get increasingly expensive.
On Sat, Nov 29, 2014 at 11:39 AM, James Crist crist...@umn.edu wrote:
@Richard,
Thanks for the Jenks paper, that was a good read. I also read your paper
on semantic matching
http://dl.acm.org.ezp2.lib.umn.edu/citation.cfm?id=806300, which
Oh boy, this is going to be a big post. Responding to everyone in turn:
*@Aaron:*
Nonlinear, AC pattern matching is NP complete. Linear AC pattern matches
can be found in polynomial time.
Interesting. Why is that?
Joachim got it right, having each match constrained by other matches,
All,
In my spare time, I've been working on implementing a fast pattern matcher
that accounts for Associative and Commutative symbols. It's going to be a
while before I'm ready to release the code (it needs some serious cleanup),
but as of now it is partly functional. Some notation:
T = set
You need to include the assumptions in the call to refine:
refine(sqrt((x - y)**2), Q.positive(x - y))
x - y
If I remember correctly, the `assuming` context manager also works for
this, but I can't be certain.
On Monday, November 10, 2014 6:30:25 PM UTC-6, Andrew Spielberg wrote:
Hi Aaron,
I really want to get this one in before the release:
https://github.com/sympy/sympy/pull/7824.
On Tuesday, October 28, 2014 12:24:02 PM UTC-5, Aaron Meurer wrote:
We are just about ready to create a release candidate for SymPy 0.7.6.
The only remaining blocking PR to my knowledge is
.
On Thursday, October 23, 2014 1:23:07 PM UTC-5, James Crist wrote:
For the large expressions that we see in `mechanics`, calls to `simplify`
can take an extremely long time. However, for simple expressions,
simplification is desirable. Currently we don't simplify by default inside
any
For the large expressions that we see in `mechanics`, calls to `simplify`
can take an extremely long time. However, for simple expressions,
simplification is desirable. Currently we don't simplify by default inside
any of our library code, as it's impossible to tell whether the expression
can
still be a problem though. This was just a proposal - I'm not adamant
that sympy needs such a feature.
On Thursday, October 23, 2014 2:39:01 PM UTC-5, Joachim Durchholz wrote:
Am 23.10.2014 um 20:23 schrieb James Crist:
However, this isn't the easiest thing to do in Python. The best
like doing that until I get some validation on the concept.
On Thursday, October 23, 2014 3:06:22 PM UTC-5, James Crist wrote:
What's composable in this context?
Easy to write without intruding too much into the actual function.
That would not affect SymPy itself
A second proof of concept PR, this time using a context manager. I actually
like this more, but it has its own issues as well.
https://github.com/sympy/sympy/pull/8297
On Thursday, October 23, 2014 4:14:27 PM UTC-5, James Crist wrote:
Proof of concept PR here https://github.com/sympy/sympy
to confirm, the latest official release is 0.7.5 and does not work
with your code. Is that right? Do you know when the next release will be?
On Friday, September 26, 2014 10:12:23 AM UTC-5, James Crist wrote:
@david:
This is all in sympy master. The codegen stuff has had major work done
@david:
This is all in sympy master. The codegen stuff has had major work done to
it in the last development cycle.
On Fri, Sep 26, 2014 at 10:10 AM, David Shin shin.da...@gmail.com wrote:
James,
Your code doesn't seem to work as expected:
$ python demo.py 10
Traceback (most recent call
The `test_kane3` test I think is not a good test in general. It's going to
run slow, and personally I've never had it complete (although my laptop is
slow). It's a fine example of the capabilities of our library, but as a
test - not so good. I doubt the rest of the mechanics team would agree
instead. That's what I get for trying to do things the easy way :/
Anyway, problem is *tentatively* solved. Still need to do the same for
`MatrixExpr` types, this only works for `Expr`. Can't imagine it will be
that different though.
- Jim
On Thursday, September 25, 2014 12:07:51 PM UTC-5, James
asme...@gmail.com
javascript: wrote:
On Fri, Sep 19, 2014 at 9:20 AM, James Crist cris...@umn.edu
javascript: wrote:
it looks like in the first example that the expression is returned
directly while in the second case it is not(?)
Ideally, for routines with one expr, `routine
I'd hazard that it's a limit in the internals of numpy for how they handle
broadcasting, but I can't be certain on that.
However, we can handle this, you just need to frame your problem in a
better way. You're trying to do optimization, so generally you'd frame your
optimal condition as a
value. I suspect there will be issues with how you define Routine, however,
since it looks like in the first example that the expression is returned
directly while in the second case it is not(?).
On Thursday, September 18, 2014 10:32:01 PM UTC-5, James Crist wrote:
I have a new SymPy type
.com/track/click?u=e2e180baf855ac797ef407fc7id=bf9baee8d0e=c9acab854b
In this talk, James Powell shows us his experiences with the worst and
most evil things he's encountered when coding in Python and gives demos
on how to get around those errors to turn them into good ideas.
Network
.
We could print out code for the expanded form of (x + a1)*(x + a2)*...*(x +
an), or we could factor it, and reduce the number of operations
significantly.
On Thu, Sep 18, 2014 at 6:03 PM, James Crist crist...@umn.edu wrote:
Oh hey, that's kinda neat.
On one of the slides you mention
And I suppose for matrix expressions, if a user tried to generate code for
x = A^-1*B, we could recognize that and replace it with a matrix solve
routine, rather than explicitly calculating the inverse.
On Thu, Sep 18, 2014 at 6:23 PM, James Crist crist...@umn.edu wrote:
Hmmm, it appears I
do this now, but I'd love to see some research in this area).
Aaron Meurer
On Thu, Sep 18, 2014 at 6:27 PM, James Crist cris...@umn.edu
javascript: wrote:
And I suppose for matrix expressions, if a user tried to generate code
for x = A^-1*B, we could recognize that and replace
I have a new SymPy type that serves to represent a unit of computation
(contains an expression/expressions that is/are being computed). I'd like
to be able to alias queries on the assumptions of the element to
assumptions of the underlying expression it represents. Example:
a, b, c =
.
- Jim
On Wednesday, September 3, 2014 10:49:03 AM UTC-5, Ondřej Čertík wrote:
Hi James,
On Mon, Sep 1, 2014 at 2:37 PM, James Crist cris...@umn.edu javascript:
wrote:
I prefer either to pass -ffast-math flag (setting compiler flags is
already an issue since we need
to indicate
no additional kwargs. If pow
is overridden by the user, it will not be included in the code printout.
Easy, simple, and requires no user facing changes.
Thoughts?
On Sat, Aug 30, 2014 at 3:04 PM, Björn Dahlgren bjo...@gmail.com wrote:
Cool! Great work!
On Friday, 29 August 2014 22:48:23 UTC+2, James
, James Crist cris...@umn.edu
javascript: wrote:
I was wondering about that. I wasn't sure if the overhead from looping
through the inputs multiple times would outweigh improvements from fast C
loops. Glad that in your case it does.
I've thrown a WIP PR up: https://github.com/sympy/sympy
wrote:
Jason
moorepants.info
+01 530-601-9791
On Fri, Aug 29, 2014 at 2:38 AM, James Crist crist...@umn.edu wrote:
I was planning on going to bed, but ended up working on this instead.
I
have no self control...
Anyway, I've uncovered some things:
1. Addition
For handling Pow? horner(x**11) results in x**11. Or were you recommending
applying horner to an entire expression tree?
On Fri, Aug 29, 2014 at 3:55 PM, Tim Lahey tim.la...@gmail.com wrote:
I recommend that you use the horner function in polys.
On 29 Aug 2014, at 16:48, James Crist wrote
I still need to do some cleanups and add tests, but I finally have this
working and thought I'd share. I'm really happy with this:
In [1]: from sympy import *
In [2]: a, b, c = symbols('a, b, c')
In [3]: expr = (sin(a) + sqrt(b)*c**2)/2
In [4]: from sympy.utilities.autowrap import ufuncify
on this today but it looks like you've by passed
what I had working. Do you have a PR with this?
Jason
moorepants.info
+01 530-601-9791
On Wed, Aug 27, 2014 at 11:11 PM, Matthew Rocklin mrock...@gmail.com
wrote:
Cool
On Wed, Aug 27, 2014 at 8:07 PM, James Crist crist...@umn.edu wrote:
I
tomorrow.
Jason
moorepants.info
+01 530-601-9791
On Wed, Aug 27, 2014 at 11:26 PM, James Crist crist...@umn.edu wrote:
Not yet. I wrote it this morning during an extremely boring meeting, and
haven't had a chance to clean it up. This doesn't solve your problem about
broadcasting a matrix
This is part of my continuing my work on making sympy capable of generating
code for evaluating matrix functions. So far the codeprinters and codegen
is done(ish); all that's left is autowrap.
Questions:
1. Should functions created by autowrap create matrices that would be
in-out parameters
routine.
On Tue, Aug 5, 2014 at 11:46 AM, Aaron Meurer asmeu...@gmail.com wrote:
On Mon, Aug 4, 2014 at 5:44 PM, Tim Lahey tim.la...@gmail.com wrote:
I've answered your questions below.
On 4 Aug 2014, at 18:27, James Crist wrote:
*1. Sympy Matrices are always 2 dimensional, should
I'm working on adding support ofr codegeneration with Matrix objects.
Currently an `indexed` type is supported that results in low-level
contiguous arrays. These are always converted into loops though (and I
don't really understand what they're for). In contrast, the intent here is
to provide
on this problem seems
sensible.
Thanks again,
-Matt
On Mon, Aug 4, 2014 at 4:13 PM, Jason Moore moorepa...@gmail.com wrote:
Jason
moorepants.info
+01 530-601-9791
On Mon, Aug 4, 2014 at 3:27 PM, James Crist crist...@umn.edu wrote:
I'm working on adding support ofr codegeneration with Matrix
for the outer Matrix as is already
shown in the docs?
On Mon, Aug 4, 2014 at 11:34 PM, Tim Lahey tim.la...@gmail.com wrote:
On 5 Aug 2014, at 0:14, James Crist wrote:
@Tim:
They're for representing tensors. Of course, the can be used for a number
of things, including calculating finite
On Friday, July 11, 2014 11:13:28 PM UTC-5, Richard Fateman wrote:
The obvious brute force method would be to use software floats in which
case you could increase the precision and the range of the
numbers involved. I'm assuming the NaNs come from division by zero where
the denominator is
* simplifying that
expression.
This way you are targeting only trouble expressions and not working with
the
whole expression tree.
On Friday, June 27, 2014 9:35:09 PM UTC-5, James Crist wrote:
signsimp didn't seem to help. I'll try to get a gist of the expression
up
tomorrow
On Wed, Jul 2, 2014 at 9:10 PM, James Crist cris...@umn.edu
javascript: wrote:
I wrote this up today. In physics.mechanics we often have to sub
symbols
for
values (or a smaller subset of symbols, i.e. the operating point).
For
the
huge expressions generated, `subs
code.
There's also the strategies submodule that Matthew added which looks
useful, but I couldn't figure out how to use it.
On Thursday, July 3, 2014 11:12:55 AM UTC-5, James Crist wrote:
Checking for these conditions in general would prove extremely tricky.
I've updated the gist with my
I wrote this up today. In physics.mechanics we often have to sub symbols
for values (or a smaller subset of symbols, i.e. the operating point).
For the huge expressions generated, `subs` is extremely slow. Also, it subs
inside derivatives, which is not ideal (we are currently using a hacky
signsimp didn't seem to help. I'll try to get a gist of the expression up
tomorrow, it's so huge I haven't succesfully ran simplify on it. However, I
had the same thing happen with a readable expression. The issue was
something like:
expr = sin(a)/tan(a)
expr.subs(a, 0)
nan
expr =
The code I'm working on in sympy.physics.mechanics often results in large
expressions (100,000 + operations, biggest I've seen was 31,000,000
operations). Once these expressions are obtained, we often need to
substitute in values/symbols (this is the operating point). However, due to
the
more natural to me too.
Should we make a convenience function that does this? I think this use
of jacobian would be lost on most people.
Aaron Meurer
On Sat, Jun 14, 2014 at 6:29 PM, James Crist cris...@umn.edu
javascript: wrote:
It's just the convention I'm most used to. Systems
-5, Aaron Meurer wrote:
That's a clever trick. I should have thought of that.
Is there any reason you let system = A*x - B instead of A*x + B? The
latter seems more natural.
Aaron Meurer
On Sat, Jun 7, 2014 at 12:28 AM, James Crist cris...@umn.edu
javascript: wrote:
I just answered
1 - 100 of 124 matches
Mail list logo