Re: [sage-support] publishing a program written with Sage in complement of a scientific article

2009-11-21 Thread Robert Bradshaw
On Nov 21, 2009, at 12:23 PM, Carlo Hamalainen wrote: On Sat, Nov 21, 2009 at 8:43 PM, Laurent moky.m...@gmail.com wrote: Thus I have two files .py (using Sage) to publish and add to my bibliography. Does one knows where and how I have to publish them ? Some website for sharing the code ?

Re: [sage-support] Can't run Cython's hello world example

2009-11-20 Thread Robert Bradshaw
On Nov 20, 2009, at 8:59 AM, Jorge E. ´Sanchez Sanchez wrote: Thank you very much Robert, I recently upgrade my Ubuntu 9.04 to 9.10 and because I have no problems with my binary Sage Version 4.1, Release Date: 2009-07-09 installation I thought I could stay with it, (last time I

Re: [sage-support] Can't run Cython's hello world example

2009-11-19 Thread Robert Bradshaw
On Nov 19, 2009, at 1:57 PM, Jorge E. ´Sanchez Sanchez wrote: Hi dear sage-support group: I was working a Differential equations problem with few elements in some regular 3D-mesh which needed not too much calculations but I am upgrading it to a bigger system. Then I shall need to do

Re: [sage-support] smoothness in plot3d

2009-11-18 Thread Robert Bradshaw
On Nov 18, 2009, at 10:09 AM, William Stein wrote: On Wed, Nov 18, 2009 at 9:54 AM, Laurent moky.m...@gmail.com wrote: Use the plot_points option. Type parametric_plot3d? for more details: - ``plot_points`` - (default: automatic, which is 75 for curves and [40,40]

Re: [sage-support] Complex Numbers in vectors and matrices

2009-11-18 Thread Robert Bradshaw
On Nov 18, 2009, at 6:25 PM, Ichnich wrote: Hi, I want to find the solution of a system of coupled differential equations. Therfore I need the eigenvalues and vectors of a matrix D and their complex conjugates. short example: sage: D=matrix([[0,1],[-1,0]]) sage:

Re: [sage-support] binary for ubuntu 8.04

2009-11-17 Thread Robert Bradshaw
On Nov 17, 2009, at 9:55 PM, William Stein wrote: On Tue, Nov 17, 2009 at 8:50 PM, Minh Nguyen nguyenmi...@gmail.com wrote: Hi Ajay, On Wed, Nov 18, 2009 at 3:28 PM, Ajay Rawat ajay.rawa...@gmail.com wrote: Hello, Is the binary for sage 4.2.1 is available for ubuntu 8.04 (64 bit).

Re: [sage-support] binary for ubuntu 8.04

2009-11-17 Thread Robert Bradshaw
On Nov 17, 2009, at 11:30 PM, Ajay Rawat wrote: Thanks for the prompt response I'm planning to compile from source. Any problem arises i'll inform you. I all ready had python installed does sage installation will clash with my python environment variable. thanking you No, Sage installs

[sage-support] Re: plot3d can't handle log(0)

2009-11-16 Thread Robert Bradshaw
On Nov 10, 2009, at 12:34 PM, kcrisman wrote: In 4.2.1.alpha0: sage: f(x,y)=ln(x) sage: P=plot3d(f,(x,0,1),(y,0,1)) sage: P ERROR: An unexpected error occurred while tokenizing input The following traceback may be corrupted or invalid The error message is: ('EOF in multi-line statement',

[sage-support] Re: plot3d can't handle log(0)

2009-11-16 Thread Robert Bradshaw
On Nov 16, 2009, at 7:23 PM, William Cauchois wrote: On Nov 16, 10:59 am, Robert Bradshaw rober...@math.washington.edu wrote: As a first pass, just surround the evaluation with a try/except, as you would in Python. Perhaps you could set the value in the except clause to nan. - Robert

[sage-support] Re: Out of memory allocating triangulation (implicit_plot3d)

2009-11-14 Thread Robert Bradshaw
On Nov 14, 2009, at 7:00 PM, Mike Witt wrote: Can anyone tell me if this is really an out of memory error? MemoryError: Out of memory allocating triangulation for type 'sage.plot.plot3d.implicit_surface.ImplicitSurface' Yes, see

[sage-support] Re: plotting a matrix with pixels corresponding to entries

2009-11-10 Thread Robert Bradshaw
On Nov 10, 2009, at 11:25 PM, Jonathan Bober wrote: Does anyone know how to plot a matrix so that each pixel in the output picture corresponds precisely to an entry in the matrix? To see the problem that I have, you might see the odd results that occur when you do: M =

[sage-support] Re: Gap Undeposited Implementations

2009-11-09 Thread Robert Bradshaw
On Nov 9, 2009, at 8:58 AM, Adam Sorkin wrote: I'd like to do it this way, running everything through sage, but it doesn't work. Here is what I'm getting: sage: g = SymmetricGroup(3) sage: t = [(2,1),(2,1)] sage: gap.Braid(g,t)

[sage-support] Re: subclass of hyperelliptic curve class

2009-11-05 Thread Robert Bradshaw
On Nov 5, 2009, at 1:09 PM, David Holmes wrote: Hi, I am trying to write a subclass of the class of hypereliptic curves over QQ, for example: class test (sage .schemes .hyperelliptic_curves .hyperelliptic_rational_field.HyperellipticCurve_rational_field): def __init__(self,C):

[sage-support] Re: Matrix inconsistency with decimal/fraction representations

2009-11-04 Thread Robert Bradshaw
On Nov 4, 2009, at 11:08 AM, Michael Orlitzky wrote: Jason Grout wrote: I don't think it's an issue of irrational versus rational. It's numerical precision and inexact floating point numbers. This matrix is terribly ill-conditioned. It is right on the border line between being

[sage-support] Re: Matrix inconsistency with decimal/fraction representations

2009-11-04 Thread Robert Bradshaw
On Nov 4, 2009, at 11:25 AM, Michael Orlitzky wrote: Robert Bradshaw wrote: I'd also like to point out that we don't just want to fall back and do everything over the rationals (even though any finite decimal expansion is rational) as things get much slower due to coefficient explosion

[sage-support] Re: Gap Undeposited Implementations

2009-11-04 Thread Robert Bradshaw
On Nov 4, 2009, at 7:27 PM, Adam Sorkin wrote: I am trying to load a Gap package for braid orbit computations. I would like to run everything through Sage if possible. This is an undeposited Implementation, so running sage -i gap_packages-4.4.10_4 doesn't retrieve this package. I'm

[sage-support] Re: Matrix inconsistency with decimal/fraction representations

2009-11-04 Thread Robert Bradshaw
On Nov 4, 2009, at 3:01 PM, Michael Orlitzky wrote: Simon King wrote: Hi Michael! On 4 Nov., 20:55, Michael Orlitzky mich...@orlitzky.com wrote: [...] it starts using floating point numbers internally. I didn't tell it to do that. You did. 0.5 is a floating point number. I guess it

[sage-support] Re: help about the function diff()

2009-11-04 Thread Robert Bradshaw
On Nov 4, 2009, at 7:43 PM, wxu...@sohu.com wrote: Hi, all I just saw that if I defined a function: f=f(e^t), How did you define f? Perhaps the f.diff(f,t) will give e^t*D[0](f)(e^t). and the question is what is the meaning of D[0](f)(e^t)? I can find that in the help of diff().

[sage-support] Re: Matrix inconsistency with decimal/fraction representations

2009-11-03 Thread Robert Bradshaw
On Nov 3, 2009, at 8:13 PM, Jason Grout wrote: Michael Orlitzky wrote: Jason Grout wrote: The problem here is that no one has really implemented a numerically stable version of echelon_form for RR. I thought we called scipy for rank calculations over RDF, but apparently we don't even

[sage-support] Re: Ubuntu 9.10

2009-11-02 Thread Robert Bradshaw
On Nov 2, 2009, at 3:22 PM, davedo2 wrote: OK, I downloaded Sage 4.2 and -testall seemed OK. However, I have a question: if I evaluate pari(7).isprime() in the notebook it returns True, but if I try it from the Sage command line it spawns a longish list of errors - what's up with that?

[sage-support] Re: Possible to view cython annotation from command line

2009-11-02 Thread Robert Bradshaw
On Nov 2, 2009, at 3:10 PM, Justin Domke wrote: Hello, When developing cython code, I find the annotation .html files (showing yellow lines where python calls are, etc.) useful. I can do this from the notebook interface, or (apparently) from command-line usage of cython, but I'm not

[sage-support] Re: Possible to view cython annotation from command line

2009-11-02 Thread Robert Bradshaw
On Nov 2, 2009, at 3:56 PM, Justin Domke wrote: Thanks, that works perfectly! One other question-- is there any way to get .html output when attaching a .spyx file? (Perhaps this is asking a lot...) No, I don't imagine it would be super hard to do though. - Robert

[sage-support] Re: [david.r.guich...@gmail.com: Re: sage server]

2009-10-29 Thread Robert Bradshaw
On Oct 29, 2009, at 3:39 PM, William Stein wrote: On Thu, Oct 29, 2009 at 3:27 PM, Jason Grout jason-s...@creativetrax.com wrote: William Stein wrote: Nothing is forbidden yet, but we should try not to load whatever we can get away with not loading in order to improve the import time.

[sage-support] Re: Python lists and removals

2009-10-27 Thread Robert Bradshaw
On Oct 27, 2009, at 3:15 AM, Francois Maltey wrote: Hello, I'm an old lisp-list user and python is my first use of dynamic array-list. Complexity for lisp-list is constant and fast o(1) when we add a new element at the head of the list. (cons e L) in Lisp or e::L in caml. Complexity

[sage-support] Re: powers of polynomials in characteristic p0

2009-10-20 Thread Robert Bradshaw
Interesting point, I've made http://trac.sagemath.org/sage_trac/ticket/7253 On Oct 19, 2009, at 1:49 PM, finotti wrote: Dear all, I need to compute some larger powers of polynomials in characteristic p0. I've noticed that Sage does not do it very efficiently, as even f^(p^n) takes a long

[sage-support] Re: Membership in the SAGE support group

2009-10-15 Thread Robert Bradshaw
On Oct 15, 2009, at 2:02 AM, domingo.sala...@syngenta.com wrote: Dear Colleague, I requested to become a member of the SAGE support group but I have not been contacted about it yet. I would like to post a question because I have had a problem installing SAGE in SUSE Linux (I believe

[sage-support] Re: Request for Ploting

2009-10-13 Thread Robert Bradshaw
On Oct 13, 10:09 am, Santanu Sarkar sarkar.santanu@gmail.com wrote: From some data set I plot the attached eps file using Matlab. How Can I do same thing using Sage only? Sage uses the extension to determine the file type. sage: g = plot(x^3-x, (x,2, -2)) sage: g.save(fig.eps) sage:

[sage-support] Re: fipy

2009-10-13 Thread Robert Bradshaw
On Oct 13, 2009, at 7:14 AM, Ajay Rawat wrote: My question is can fipy http://www.ctcms.nist.gov/fipy/ run on sage? -- I don't see any reason why not (though you would have to install it to see). - Robert --~--~-~--~~~---~--~~ To post to this group,

[sage-support] Re: Cython and SciPy special functions

2009-10-05 Thread Robert Bradshaw
On Oct 4, 2009, at 10:07 AM, mark mcclure wrote: I'm investigating the complex dynamics of Airy functions with Python and I wonder if I can speed up the process significantly with Cython. I've successfully sped up the code by about 50% but I might expect much greater speed improvement. I'm

[sage-support] Re: Cython and SciPy special functions

2009-10-05 Thread Robert Bradshaw
On Oct 5, 2009, at 1:26 PM, mark mcclure wrote: On Oct 5, 1:19 pm, Robert Bradshaw rober...@math.washington.edu wrote: you might want to call these functions in SciPy directly as C functions. (You'd have to look up the SciPy headers to see what to call them.) That was my first thought

[sage-support] Re: setuptools vs. Distribute

2009-09-29 Thread Robert Bradshaw
On Sep 29, 2009, at 3:25 AM, Flavio Coelho wrote: Hi, what the opinion of Sage developers regarding the fork from setuptools == Distribute? http://pypi.python.org/pypi/distribute/0.6.3 From what I have heard so far, setuptools will no longer be developed and Distribute already supports

[sage-support] Re: Symbolic derivative question

2009-09-29 Thread Robert Bradshaw
On Sep 28, 2009, at 11:30 AM, Tim Lahey wrote: On Sep 28, 2009, at 2:21 PM, kcrisman wrote: Dear support (and/or Burcin), How does Sage/Pynac support derivatives evaluated at a point (or does it)? E.g., sage: f = function('f',t) sage: h = f.diff(t,1) sage: h.subs(t=0) D[0](f)(0)

[sage-support] Re: twisted server

2009-09-29 Thread Robert Bradshaw
On Sep 29, 2009, at 4:00 PM, calcp...@aol.com wrote: Isn't dSage based on twisted? Yes, it is. It hasn't been maintained much these last couple of years though. - Robert --~--~-~--~~~---~--~~ To post to this group, send email to

[sage-support] Re: statistics in sage

2009-09-26 Thread Robert Bradshaw
Yes, there's some stuff in the finance.time_series for doing simple stats over the reals. (This really needs to be put somewhere more obvious, along with histogram, etc.) sage: sage.finance.time_series.TimeSeries([1,2,5]) [1., 2., 5.] sage: t =

[sage-support] Re: sage vmware image visible from Internet

2009-09-24 Thread Robert Bradshaw
On Sep 24, 2009, at 12:11 AM, ma...@mendelu.cz wrote: Hello all, I am trying to provide Sage server for my students. If I run Sage in vmware, I can access the vmware virtual machine from my PC only. The students can connect to my PC, but not to the virtual machine. I guess, some port

[sage-support] Re: sage vmware image visible from Internet

2009-09-24 Thread Robert Bradshaw
On Sep 24, 2009, at 1:17 AM, ma...@mendelu.cz wrote: On 24 zář, 09:17, Robert Bradshaw rober...@math.washington.edu wrote: This depends on what kind of PC you're running the server on. The way we do sagenb.org is have apache running on the host machine which forwards all incoming

[sage-support] Re: sagenb.org

2009-09-22 Thread Robert Bradshaw
On Sep 22, 2009, at 2:10 PM, William Stein wrote: On Tue, Sep 22, 2009 at 2:04 PM, Eric Jackson eric...@comcast.net wrote: Earlier today I signed up through sagenb.org so that I can publish my worksheet. The signing up process was simple, but whenever I try to login, I receive errors.

[sage-support] Re: Digamma Function in 4.1.1

2009-09-21 Thread Robert Bradshaw
On Sep 21, 2009, at 10:29 PM, The_Fool wrote: On Sep 21, 9:37 pm, Jason Grout jason-s...@creativetrax.com wrote: After modifying a file, do: sage -br which copies the modified files to the build directory and rebuilds Sage and then runs Sage. To just rebuild, just do sage -b I am

[sage-support] Re: sxrange / xsrange possible bug

2009-09-16 Thread Robert Bradshaw
On Sep 16, 2009, at 7:22 AM, Mariah wrote: Minh, Ok, how does the following look? Like what you want? The date looks funny. # HG changeset patch # User mariah.le...@gmail.com # Date 1253110290 14400 # Node ID

[sage-support] Re: Newbie question

2009-09-16 Thread Robert Bradshaw
On Sep 16, 2009, at 2:23 PM, Dan wrote: I've been using sagenb for a few weeks now and have the hang of most of the easy things. Is there a way to generate tables of a user defined function? I'm still switching over to Excel or the TI calculator to do that. TI table function lets you punch

[sage-support] Re: eval(string)

2009-09-14 Thread Robert Bradshaw
about [1, 2, 3, [4]*(130)]? How much memory do you have on your machine? - Robert On Sep 10, 1:56 pm, Robert Bradshaw rober...@math.washington.edu wrote: On Sep 10, 2009, at 12:24 PM, Robert Bradshaw wrote: On Sep 9, 2009, at 9:01 AM, Tim Dumol wrote: `eval(the_string, globals

[sage-support] Re: graphics broken in 10.6?

2009-09-12 Thread Robert Bradshaw
On Sep 11, 2009, at 8:31 PM, j9mosely wrote: I'm new to sage and I can't get any graphics, 2d or 3d, to draw. I upgraded to Snow Leopard shortly before downloading sage. When I run: circle((0,0), 1, rgbcolor=(1,1,0)) from the tutorial, I get the error:

[sage-support] Re: graphics broken in 10.6?

2009-09-12 Thread Robert Bradshaw
On Sep 11, 2009, at 11:25 PM, William Stein wrote: By the way, Robert Bradshaw just wrote: Almost certainly an incompatibility with 10.6--I would either wait until new binaries are available, or compile from scratch. (The latter is not as hard as it sounds--make sure you have xcode

[sage-support] Re: possible Polynomial_integer_dense_flint bug

2009-09-10 Thread Robert Bradshaw
On Sep 10, 2009, at 8:54 AM, Mariah wrote: Sorry, mistyped. Should be A = 2^(2^17+2^15) On Sep 10, 11:51 am, Mariah Lenox mariah.le...@gmail.com wrote: In sage-4.1.1 R.x = PolynomialRing(ZZ) A = 2^(2^17+s^15) a = A * x^31 b = (A * x) * x^30 a == b # prints False ??? I believe A

[sage-support] Re: old sagenb.org notebooks

2009-09-10 Thread Robert Bradshaw
it would be great to have a push/pull, or even just a sync feature of notebooks. That way, I could run my local notebook, sync to the public sagenb.org server, and therefore easily back up everything on sagenb.org in my account. Robert Bradshaw already implemented that, but it needs

[sage-support] Re: eval(string)

2009-09-10 Thread Robert Bradshaw
=[3,10,15,23,25,30,3,[5]*3] Need the repeated values for the 5. If I don't have repeated values your code works. I have done some error traping for eval. Users can not put something like rm. On Sep 8, 2:38 pm, Robert Bradshaw rober...@math.washington.edu wrote: On Sep 8, 2009, at 11:28 AM, Mikie wrote

[sage-support] Re: eval(string)

2009-09-10 Thread Robert Bradshaw
On Sep 10, 2009, at 12:24 PM, Robert Bradshaw wrote: On Sep 9, 2009, at 9:01 AM, Tim Dumol wrote: `eval(the_string, globals = {__builtins__:None}, locals = {})` should do it. This removes access from all functions. Add any functions that are needed by adding them to the locals dictionary

[sage-support] Re: eval(string)

2009-09-08 Thread Robert Bradshaw
On Sep 8, 2009, at 11:28 AM, Mikie wrote: Here is the function --- def BasicStats1a(exp1): v = exp1 v1 = eval(v);Count_=len(v1) sort_v1=sorted(v1) M1 =stats.mode(v1); v3=eval(str(M1[0])); v4=eval(str(M1[1])) R1 = stats.mean(v1);R2 =

[sage-support] Re: generating and exporting .dat file

2009-09-05 Thread Robert Bradshaw
I still have no idea what a .dat file is. If it's plain ascii, or a simple binary format, I'm sure Sage could do it with little difficulty. - Robert On Sep 5, 2009, at 1:26 PM, Afonso Henriques Silva Leite wrote: Does the mathematica file will help you? It is attached. Else, I also

[sage-support] Re: generating and exporting .dat file

2009-09-05 Thread Robert Bradshaw
Easy. sage: M = random_matrix(RDF, 10, 2) # get your data sage: s = M.str().replace('[','').replace(']','') sage: print s # is this what you want? 0.103609743105 -0.00797973769955 -0.96222133551 -0.341208831103 -0.904012926167-0.42972921542 -0.6744275756920.338249515207

[sage-support] Re: plotting and float

2009-09-02 Thread Robert Bradshaw
On Sep 1, 2009, at 8:58 PM, William Cauchois wrote: It seems to me that the error comes from feeding a function which uses i into fast_float (called by the plotting functions to compile the function to be plotted into an optimized form). I tried a simpler function using i and got the same

[sage-support] Re: Maxima Solve sys

2009-09-02 Thread Robert Bradshaw
GUI. - Robert On Sep 1, 10:02 pm, Robert Bradshaw rober...@math.washington.edu wrote: On Sep 1, 2009, at 12:16 PM, Mikie wrote: I took out the eval and for some reason it is working. Robert, this is function in my API (AlgCalc) http://pirsqrt.com:1843/ If I would give it to you would you

[sage-support] Re: Completing the square

2009-09-02 Thread Robert Bradshaw
On Wed, 2 Sep 2009, calcp...@aol.com wrote: That's a good question! I've written functions in MATLAB (well Octave actually) no problem. But I get confused where Python leaves off and Sage kicks in when witting functions here. The way to think about this is that Sage is just a huge Python

[sage-support] Re: Evaluation of time spent on each instruction...

2009-09-01 Thread Robert Bradshaw
On Sep 1, 2009, at 1:34 AM, Nathann Cohen wrote: Hello everybody !!! I would like to test the functions I write to learn which instructions take most of the time, and if possible find a way around.. There may not be some tool for this in Sage, but perhaps it already exists for

[sage-support] Re: A weird exception?

2009-09-01 Thread Robert Bradshaw
On Aug 29, 2009, at 11:05 PM, Rolandb wrote: Hi, Using math.log has a disadvantage; it is less accurate. sage: print n(math.log(2),100) sage: print n(log(2),100) 0.69314718055994528622676398300 0.69314718055994530941723212146 Yes, it is. Accuracy vs. speed is a common tradeoff one needs

[sage-support] Re: Maxima Solve sys

2009-09-01 Thread Robert Bradshaw
On Sep 1, 2009, at 9:17 AM, Mikie wrote: Here is function I am using to solve systems of linear equations. def MSolveSys(syss): eqns=eval(syss) solns=maxima.solve(syss) return solns Works great in the notebook, but when I put it in a Python script it rounds the coeficients of

[sage-support] Re: A weird exception?

2009-09-01 Thread Robert Bradshaw
On Sep 1, 2009, at 9:23 AM, Rolandb wrote: On 1 sep, 11:08, Robert Bradshaw rober...@math.washington.edu wrote: On Aug 29, 2009, at 11:05 PM, Rolandb wrote: Hi, Using math.log has a disadvantage; it is less accurate. sage: print n(math.log(2),100) sage: print n(log(2),100

[sage-support] Re: Maxima Solve sys

2009-09-01 Thread Robert Bradshaw
: Robert Bradshaw wrote: On Sep 1, 2009, at 9:17 AM, Mikie wrote: Here is function I am using to solve systems of linear equations. def MSolveSys(syss):    eqns=eval(syss)    solns=maxima.solve(syss)    return solns Works great in the notebook, but when I put it in a Python script it rounds

[sage-support] Re: Maxima Solve sys

2009-09-01 Thread Robert Bradshaw
busy. I have loaded load.js, but when it goes for another .js file it cannot find. You're probably thinking about another Robert--I don't know anything about jsmath. What I would do is look for instructions on their site. - Robert (Bradshaw) On Sep 1, 12:37 pm, Mikie thephantom6

[sage-support] Re: A weird exception?

2009-08-29 Thread Robert Bradshaw
On Aug 29, 2009, at 12:44 PM, Rolandb wrote: On 29 aug, 18:43, John H Palmieri jhpalmier...@gmail.com wrote: On Aug 29, 8:51 am, Simon King simon.k...@nuigalway.ie wrote: Note that expon uses Maxima, because you use the logarithm. So, I reckon that your problem is related

[sage-support] Re: English grammar of numbers

2009-08-29 Thread Robert Bradshaw
On Aug 29, 2009, at 5:25 PM, kcrisman wrote: On Aug 29, 6:45 pm, Robert Dodier robert.dod...@gmail.com wrote: Minh Nguyen wrote: Usually an comes before a word that starts with a vowel, i.e a, e, i, o, u. So one would say an eight o'clock meeting or an 8 o'clock meeting. More examples:

[sage-support] Re: How can I incorporate ready made graphics?

2009-08-27 Thread Robert Bradshaw
On Aug 27, 2009, at 9:56 PM, Ling Kah Jai wrote: I have found out a method of adding jpg or png graphics by: 1) First deposit a jpg or png file to a website server, example here is www.irnova.com.my/Sage/Proving%20Ring.PNG 2) Then incorporate this jpg or png file by using html: a)

[sage-support] Re: Weave broken in Sage 4.1 Mac intel 64 bit

2009-08-25 Thread Robert Bradshaw
On Aug 20, 2009, at 9:02 AM, felix wrote: Hi, this should have happened to other people, but I can't find some other post on this one. I'm not sure which update exactly caused the bug, since I didn't use weave since Sage 3.something. All I can say is, that weave doesn't work at all in Sage

[sage-support] Re: Loading pickles dumped from within a notebook

2009-08-24 Thread Robert Bradshaw
with a worksheet. - Robert On 24 ago, 18:59, Robert Bradshaw rober...@math.washington.edu wrote: On Mon, 24 Aug 2009, Flavio Coelho wrote: Hi, I am porting a python script to a notebook and in my script I save some data as pickles which  are loaded later for further processing. I noticed

[sage-support] Re: How does one define a real variable

2009-08-22 Thread Robert Bradshaw
On Aug 22, 2009, at 3:13 AM, Mani chandra wrote: Mani chandra wrote: Hi, How does one define a variable in sage and make sure the it's conjugate is the same as the variable itself? Thanks, Mani chandra Hi, Perhaps I wasn't specific. sage: var('a') sage: a.conjugate()

[sage-support] Re: How does one define a real variable

2009-08-22 Thread Robert Bradshaw
On Aug 22, 2009, at 1:23 PM, Golam Mortuza Hossain wrote: Hi, On Sat, Aug 22, 2009 at 2:28 PM, Robert Bradshawrober...@math.washington.edu wrote: You have to invoke simplify for it to look at assumptions. sage: var('a') sage: assume(a, 'real') sage: a.conjugate().simplify() a I

[sage-support] Re: How does one define a real variable

2009-08-22 Thread Robert Bradshaw
On Aug 22, 2009, at 2:41 PM, William Stein wrote: On Sat, Aug 22, 2009 at 2:39 PM, Robert Bradshaw rober...@math.washington.edu wrote: On Aug 22, 2009, at 1:23 PM, Golam Mortuza Hossain wrote: Hi, On Sat, Aug 22, 2009 at 2:28 PM, Robert Bradshawrober...@math.washington.edu wrote

[sage-support] Re: Optional package in Sage : Graphviz

2009-08-21 Thread Robert Bradshaw
On Fri, 21 Aug 2009, Jason Grout wrote: Nathann Cohen wrote: Hello As we are dealing in Sage with a patch concerning the plotting of Trees, I wondered why we were not using Graphviz to plot them.. The answer is the usual one : Graphviz is not GPL-Uncompatible... Even though, there

[sage-support] Re: force sage to print a matrix

2009-08-21 Thread Robert Bradshaw
On Fri, 21 Aug 2009, Simon King wrote: On 22 Aug., 00:57, William Stein wst...@gmail.com wrote: On Fri, Aug 21, 2009 at 3:18 PM, Simon King simon.k...@nuigalway.ie wrote: [...]  sage: print M.str() but there might be nicer (and more intuitive!) ways. I think that is the only way.  I can't

[sage-support] Re: change a string to sage expression by sage_eval

2009-08-19 Thread Robert Bradshaw
. Perhaps this come from the function the eval of phi[x,y], v[x,y], Ah, so this is a Mathematica string where a[x,y] is actually a function application. You might be able to get away with s.replace([, ().replace(], )) On 18 août, 11:02, Robert Bradshaw rober...@math.washington.edu wrote: On Aug

[sage-support] Re: How to properly define a 'Sage-function'?

2009-08-19 Thread Robert Bradshaw
On Wed, 19 Aug 2009, KvS wrote: Dear all, just started exploring Sage (via sagenb.org), I'm very enthousiastic about the concept and am very eager to leave 'black box' Mathematica asap. One issue however I can't seem to get my head around, namely what exactly is the 'right' way to think

[sage-support] Re: How to properly define a 'Sage-function'?

2009-08-19 Thread Robert Bradshaw
On Wed, 19 Aug 2009, KvS wrote: On Aug 19, 11:06 pm, Robert Bradshaw rober...@math.washington.edu wrote: On Wed, 19 Aug 2009, KvS wrote: Dear all, just started exploring Sage (via sagenb.org), I'm very enthousiastic about the concept and am very eager to leave 'black box' Mathematica

[sage-support] Re: Memory issues with maxima in SAGE

2009-08-18 Thread Robert Bradshaw
On Aug 17, 2009, at 2:21 PM, Robert Dodier wrote: On Aug 17, 12:06 pm, William Stein wst...@gmail.com wrote: My understanding is that one has to modify the maxima source code itself and recompile maxima. Not so -- EXT:SET-LIMIT can be called anytime after the Maxima session is launched.

[sage-support] Re: change a string to sage expression by sage_eval

2009-08-18 Thread Robert Bradshaw
On Aug 17, 2009, at 3:45 AM, Viny wrote: Hello! Is it possible to have the equivalent sage expression of the following string, by using sage_eval() or another function without getting this error in while the transcription : There was talk a while back about x[1] giving a sub-scripted x,

[sage-support] Re: Debugging; how to use sage -gdb?

2009-08-07 Thread Robert Bradshaw
On Aug 6, 2009, at 3:59 PM, Simon King wrote: Dear Minh, On 7 Aug., 00:46, Minh Nguyen nguyenmi...@gmail.com wrote: This is unrelated but: the version of GCC that is distributed with openSUSE 11.0 is $ gcc --version gcc (SUSE Linux) 4.3.1 20080507 (prerelease) [gcc-4_3-branch

[sage-support] Re: Debugging; how to use sage -gdb?

2009-08-07 Thread Robert Bradshaw
On Aug 7, 2009, at 1:16 AM, Simon King wrote: Hi Robert, On 7 Aug., 09:06, Robert Bradshaw rober...@math.washington.edu wrote: ... I came to asking about sage -gdb since sage's last words at crashing suggest to use sage -gdb. But then I would expect that it is explained in detail

[sage-support] Re: Y axis non at y=0.

2009-07-30 Thread Robert Bradshaw
On Jul 30, 2009, at 2:01 PM, Laurent wrote: and that sometimes the graph might be busy in any particular location one chose. Some users might also find the labeling of axes distracting. Anyway, please let us know what your ideas would be. * For me, it would be satisfactory to draw the

[sage-support] Re: The speed of random_element in finite fields versus Magma

2009-07-30 Thread Robert Bradshaw
On Jul 30, 2009, at 2:50 PM, VictorMiller wrote: I just did a test of SAGE versus Magma on the same computer. I had a finite field GF(19991^2), and timed generating a random element in SAGE and in Magma. I found, much to my surprise, that Magma was a factor of 7 times faster. Does anyone

[sage-support] Re: The speed of random_element in finite fields versus Magma

2009-07-30 Thread Robert Bradshaw
On Jul 30, 6:07 pm, Robert Bradshaw rober...@math.washington.edu wrote: On Jul 30, 2009, at 2:50 PM, VictorMiller wrote: I just did a test of SAGE versus Magma on the same computer. I had a finite field GF(19991^2), and timed generating a random element in SAGE and in Magma. I found, much

Re: [Cython] [sage-support] Re: changeing sage -b (cython)

2009-07-29 Thread Robert Bradshaw
On Jul 29, 2009, at 9:46 AM, Dag Sverre Seljebotn wrote: Robert Bradshaw wrote: I'm not sure if or when either of these will be available, but neither are trivial. Both questions are probably better asked on the cython lists. - Robert On Jul 22, 2009, at 12:17 PM, Ethan Van Andel wrote

[sage-support] Re: changeing sage -b (cython)

2009-07-29 Thread Robert Bradshaw
On Jul 29, 2009, at 1:14 PM, dagss wrote: On Jul 29, 7:22 pm, Robert Bradshaw rober...@math.washington.edu wrote: On Jul 29, 2009, at 9:46 AM, Dag Sverre Seljebotn wrote: Robert Bradshaw wrote: I'm not sure if or when either of these will be available, but neither are trivial. Both

[sage-support] Re: doctests and the long flag

2009-07-25 Thread Robert Bradshaw
On Jul 24, 2009, at 11:22 AM, Simon King wrote: I think it would be nice to have the possibility to do the following in doc tests: EXAMPLES:: sage: R.x,y,z = QQ[] #if testflag contains long sage: I = ... # some nasty ideal sage: G = I.groebner_basis() #if

[sage-support] Re: doctests and the long flag

2009-07-25 Thread Robert Bradshaw
On Jul 25, 2009, at 9:38 AM, William Stein wrote: On Fri, Jul 24, 2009 at 2:19 PM, Robert Bradshawrober...@math.washington.edu wrote: On Jul 24, 2009, at 11:22 AM, Simon King wrote: I think it would be nice to have the possibility to do the following in doc tests: EXAMPLES::

[sage-support] Re: changeing sage -b (cython)

2009-07-23 Thread Robert Bradshaw
I'm not sure if or when either of these will be available, but neither are trivial. Both questions are probably better asked on the cython lists. - Robert On Jul 22, 2009, at 12:17 PM, Ethan Van Andel wrote: Robert, Is there any prediction for when numpy complex types will work?

[sage-support] Re: changeing sage -b (cython)

2009-07-22 Thread Robert Bradshaw
After applying both http://trac.sagemath.org/sage_trac/ticket/6438 and http://trac.sagemath.org/sage_trac/ticket/4571 , the notebook seems to work just as well as the command line, or just raw Cython. Hopefully this makes it into 4.1.1. As for complex types, there's still some issues of

[sage-support] Re: changeing sage -b (cython)

2009-07-22 Thread Robert Bradshaw
On Jul 22, 2009, at 10:55 AM, Ethan Van Andel wrote: Robert, Should those double complex declarations give any speedup? Yes, lots (compared to using Python complex objects). The are complex analogues of int, double, float, etc. Are they declaring the numpy complex type or something else?

[sage-support] Re: changeing sage -b (cython)

2009-07-22 Thread Robert Bradshaw
On Jul 22, 2009, at 11:27 AM, William Stein wrote: On Wed, Jul 22, 2009 at 11:21 AM, Michael Madisonmadison.mich...@gmail.com wrote: Robert, Thanks for the feedback. I don't want to wait for 4.1.1 so I will attempt to follow the Sage Days 16 How to get started with developing Sage to apply

[sage-support] Re: compatibility mode in sage 4.x?

2009-07-21 Thread Robert Bradshaw
On Jul 20, 2009, at 7:49 AM, Jason Grout wrote: Stan Schymanski wrote: Thanks, Laurent. I figured that out, but it is quite tedious to have to go through each worksheet and change the soln[0][1] to soln[0].rhs(). There are other things that work differently now and of course there are

[sage-support] Re: How can calculate the time of execution of a command in sage

2009-07-16 Thread Robert Bradshaw
If you want to time running something just once, see cputime() and walltime(). - Robert On Jul 15, 2009, at 12:12 PM, Gustavo Rama wrote: Thanks, I'll try it. Cheers Gustavo On Jul 14, 10:42 pm, Simon King simon.k...@uni-jena.de wrote: Hi Gustavo! On 15 Jul., 03:17, Gustavo Rama

[sage-support] Re: Cython Modules

2009-07-16 Thread Robert Bradshaw
You need include_dirs = numpy_include_dirs in your Extension(...) declaration in module_list.py. - Robert On Jul 15, 2009, at 11:23 AM, Ethan Van Andel wrote: I want to include some of my cython code as a sage module. I followed the directions for adding the .pyx file to the sage library

[sage-support] Re: Is it possible to replace the list built-in type in Sage?

2009-07-16 Thread Robert Bradshaw
On Jul 14, 2009, at 3:35 PM, William Stein wrote: 2009/7/14 Carlos Córdoba ccordob...@gmail.com: Thanks John, I'd seen Python comprehensions before, but since I was trying to do all in a one-liner, I think I overlooked your elegant and simple solution. One comprehension at a time is

[sage-support] Re: assume()

2009-07-16 Thread Robert Bradshaw
On Jul 14, 2009, at 9:05 AM, Doug wrote: Hmm. I've also had trouble interpreting what assume() affects, and I'm glad to hear that I'm not the only one. What Robert says here helps a lot, but is there anything written anywhere else that goes into a bit more detail? I'm sure there's more to

[sage-support] Re: Cython typedef not working

2009-07-14 Thread Robert Bradshaw
There is a ticket in progress to fix this. I've started reviewing it, I'm hopefully it'll go into 4.1.1. On Jul 14, 2009, at 7:11 AM, Ethan Van Andel wrote: When I run this code: %cython %time import numpy as np cimport numpy as np COMPLEX = np.complex128 ctypedef np.complex128_t

[sage-support] Re: Newb basic algebra question

2009-07-14 Thread Robert Bradshaw
On Jul 14, 2009, at 9:18 AM, Doug wrote: Hi all; I'm brand new to sage and am finding it really fun and useful. Thanks. Of course, I don't know anyone else who uses it and thus need to rely on the general community when reading documentation and scratching my head leave me with open

[sage-support] Re: Nearest Integer Function in sage

2009-07-13 Thread Robert Bradshaw
On Jul 13, 2009, at 6:18 AM, mac8090 wrote: Hi A recent calculation I made in Sage (version 3.4) gave me a matrix where each value was out by a decimal less than 10 x e-16 or e-17. I have two questions: 1. Is this a usual rounding error in sage or is there an error in my calculations? I

[sage-support] Re: Is it possible to replace the list built-in type in Sage?

2009-07-13 Thread Robert Bradshaw
On Jul 12, 2009, at 2:49 PM, Kevin Horton wrote: Carlos Córdoba wrote: Thanks for your quick answers. Coming from Mathematica, I was expecting to add lists as vectors, multiply real numbers by lists, etc, without sub-classing or using another types (such as vectors in sage) Sage lists

[sage-support] Re: assume()

2009-07-11 Thread Robert Bradshaw
On Jul 11, 2009, at 4:39 PM, Minh Nguyen wrote: On Sun, Jul 12, 2009 at 9:12 AM, Nealneal.har...@gmail.com wrote: Hi everyone, I thought I'd share the following: sage: assume(x0) sage: solve([x^2-1],x) [x == -1, x == 1] Shouldn't it not give me the negative solution? Also: sage:

[sage-support] Re: scipy binomial pmf

2009-07-10 Thread Robert Bradshaw
. I am very impressed. Well, I figured it's about time someone did it :). I thought I was going to have to patch the numpy c sources myself, but shortly after diving into them I saw the much simpler solution. Kudos and thanks to Robert Bradshaw. If I ever meet you, I'd be happy to buy you

[sage-support] Re: scipy binomial pmf

2009-07-10 Thread Robert Bradshaw
itself? I would recommend creating a patch and trying to get it into Sage proper (there has been discussion about the lack of easy-to-use stats support). A lot has changed since 3.2.2. On Jul 10, 3:10 am, Robert Bradshaw rober...@math.washington.edu wrote: On Jul 9, 2009, at 1:57 PM, Kevin

[sage-support] Re: scipy binomial pmf

2009-07-07 Thread Robert Bradshaw
On Jul 6, 2009, at 2:21 PM, Mikie wrote: Yes, I am calling Scipy functions from a standard Python, thus no preparser. It does not recognize 1r, 10r, etc. Is there anyway to get it to work? I have also tried R, but r.binom does not work. An example, including the error message, would be

<    1   2   3   4   5   6   7   8   >