[sage-support] ask.sagemath.org

2014-06-28 Thread Mike Hansen
Hello all, We've recently updated http://ask.sagemath.org to the newest version of askbot and installed some additional tools to help with spam. New users should be able to register and post questions, but they will have to complete a reCAPTCHA until they reach karma 10. Let me know if you have

Re: [sage-support] Re: Completely confused about simple Sage matrix operation error

2013-07-01 Thread Mike Hansen
I think problem is actually due to the inverse using a non-numerically stable echelon form algorithm for inexact fields. For example, if you using matrices over RDF: M=matrix(RDF,[[7,3,10,13],[1,1,2,2],[1,2,3,4],[1,3,5,7]]);det(M);invM=M^(-1);invM*M;det(invM) you don't get this problem. My gues

Re: [sage-support] format for output to file to be used by fortran

2012-08-20 Thread Mike Hansen
On Mon, Aug 20, 2012 at 1:21 PM, bsmile wrote: > BTW, I would like to ask a basic question. SymmetricGroupRepresentation has > three options, "specht", "orthogonal", "seminormal". Could you please > explain briefly their differences, and when which one should normally be > used. It seems "orthogon

Re: [sage-support] Re: Group action

2012-05-14 Thread Mike Hansen
On Mon, May 14, 2012 at 11:20 PM, Nathann Cohen wrote: >> One thing to watch out for is that the generators returned by >> automorphism_group contain symbols that may not be the actual vertices. I >> realised this once after several frustrating hours of bizarre results from >> my program. I'm not

Re: [sage-support] Re: Use lazy_attribute in cython?

2012-05-02 Thread Mike Hansen
On Wed, May 2, 2012 at 12:08 AM, P Purkayastha wrote: > Currently, in my Sage (python) code, I used lazy_attribute to achieve this. > I tested the computation using cython and without lazy attribute and it > speeds up by a factor of about 3 compared to the python code. So, I am > looking for somet

Re: [sage-support] Sage container class?

2012-03-21 Thread Mike Hansen
On Wed, Mar 21, 2012 at 4:17 PM, Emil wrote: > Alternatively, if I just make a dictionary of SageObjects and pickle, > would that work? Thanks, A normal dictionary should work: sage: loads(dumps(dict(a=2, b=random_matrix(ZZ, 4 {'a': 2, 'b': [ 0 -1 -2 0] [-2 1 -1 -1] [ 3 1 -1 1] [ 4 6 -1

[sage-support] Re: [sage-combinat-devel] Re: Is there a way to get rid of the new google groups?

2012-03-19 Thread Mike Hansen
On Mon, Mar 19, 2012 at 2:11 PM, Nicolas M. Thiery wrote: > On Mon, Mar 19, 2012 at 05:26:16PM +0800, Dima Pasechnik wrote: >> Hi sage-combinat-devel folks, >> could you advice on how you manage your group? > > As far as I remember, I haven't done anything specific.

Re: [sage-support] __call__() much slower than a "real" sage function ?

2012-03-14 Thread Mike Hansen
On Wed, Mar 14, 2012 at 12:08 AM, Laurent wrote: > My question in the following : is it normal that the second one is *much* > slower than the first one ? The reason why the second one is much faster is Sage has something called "fast_callable" which can optimize the evaluation of the function ov

Re: [sage-support] Polynomials with Irrational coefficients

2012-02-19 Thread Mike Hansen
Hello, On Sun, Feb 19, 2012 at 12:29 PM, ObsessiveMathsFreak wrote: > Sage doesn't seem to like these > > f(x,y)=pi*x^10*y+3*x > B=f(x,y).polynomial(SR) > print B.coefficients() > > Traceback (click to the left of this block for traceback) > ... > TypeError > > What's going wrong here? Is there a

Re: [sage-support] Re: Coefficients of multivariate polynomials along with degrees

2012-02-19 Thread Mike Hansen
On Sun, Feb 19, 2012 at 4:38 AM, ObsessiveMathsFreak wrote: > Basically, I want to restrict the variables of a multivariate > polynomial to a certain set of symbols. > > B.variables() should return > [x,y] > instead of > [a,x,y] > > How can this be done. I have a patch at http://trac.sagemath.org

Re: [sage-support] Coefficients of multivariate polynomials along with degrees

2012-02-08 Thread Mike Hansen
On Wed, Feb 8, 2012 at 3:29 PM, ObsessiveMathsFreak wrote: > Is there no way of getting sage to give back the degree's of the > corresponding multivariate polynomials as well. That is to return > something like > > [[10,[100,1]],[3,[1,0]]] exponents() is what you need: sage: zip(B.coefficients()

Re: [sage-support] Re: transforming numbers into lists.

2011-12-08 Thread Mike Hansen
On Thu, Dec 8, 2011 at 5:24 PM, Eric Kangas wrote: > Now that I have code ready. There is another issue that I found. I > need to figure out how to find a set sequence of numbers in the list. > I know in mathematica there is a function called take which allows for > you to display a partial list f

Re: [sage-support] Interrupted Sphinx docstrings

2011-11-18 Thread Mike Hansen
Hello, On Fri, Nov 18, 2011 at 1:48 AM, Nathann Cohen wrote: > Would anyone know if it is possible in a module's docstring to > "interrupt" the documentation ? > > I would like to write something like that : > > "This module deals with []. Here is the first block of functions" > > def functio

Re: [sage-support] why doesn't [lambda x: bessel_J(n, x) for n in [0,1]] work?

2011-10-04 Thread Mike Hansen
On Tue, Oct 4, 2011 at 9:09 PM, Dan Drake wrote: > It seems unreasonably annoying to plot a bunch of Bessel functions > together. How can I work around this? This is due to the way binding works in list comprehensions. For a good explanation and workaround, see http://stackoverflow.com/questions

Re: [sage-support] How Maxima is called from Python ?

2011-07-19 Thread Mike Hansen
On Tue, Jul 19, 2011 at 8:00 AM, Christophe BAL wrote: > Hello, > I would like to know how Sage calls Maxima from Python ? Which files have I > to look for in the source ? Sage (as of 4.7) uses pexpect to communicate with Maxima. The main places you'll want to look are in $SAGE_ROOT/devel/sage/s

Re: [sage-support] infinite polynomial ring as a free commutative algebra

2011-07-06 Thread Mike Hansen
On Wed, Jul 6, 2011 at 9:51 AM, jeremy chabot wrote: > Hi, I am trying to implement  the infinite polynomial ring as a free > commutative algebra. > > I am unsure of exactly how to do this. So far I have: > > --- > X.=

Re: [sage-support] Change language for a cell.

2011-05-24 Thread Mike Hansen
On Mon, May 23, 2011 at 11:59 PM, mankoff wrote: > Is there a way to change the language for one cell? For example, if in > a worksheet I'd like to have just one paragraph of LaTeX or a chunk of > R. It is awkward to wrap it all in a latex('...') function call. If a > command exists that means "ev

Re: [sage-support] Decomposing polynomials from other polynomials using Gröbner bases

2011-04-05 Thread Mike Hansen
On Tue, Apr 5, 2011 at 1:24 PM, Johan S. R. Nielsen wrote: > Let's say that I have a multivariate polynomial ring R which contains > the polynomials p, f1, ..., fn. I also know that p is in the ideal J = > . Now I wish to write p as a polynomial in the f- > polynomials. How can I do that with Sage

Re: [sage-support] Absolute value of matrices

2011-04-04 Thread Mike Hansen
Hello, On Tue, Apr 5, 2011 at 12:00 AM, pong wrote: > By that I simply mean a function that on input a real matrix M returns > the matrix N such that n[i][j] = abs(m[i][j]). > > This can be achieve by something like: > > n = len(M.rows()); m =len(M.columns()); N = matrix(n,m,lambda i,j: > abs(M[i

Re: [sage-support] Re: Arithmetic Progession syntax

2011-04-01 Thread Mike Hansen
On Fri, Apr 1, 2011 at 8:08 PM, ObsessiveMathsFreak wrote: >> sage: [0,0.2,..,1] >> [0.000, 0.200, 0.400, >> 0.600, 0.800, 1.00] > > May sage installation appears to be having trouble with this syntax > > OK Here > sage: [0.0,

Re: [sage-support] Re: Arithmetic Progession syntax

2011-03-28 Thread Mike Hansen
On Mon, Mar 28, 2011 at 1:57 PM, ObsessiveMathsFreak wrote: > Thank you. That is more or less what I was looking for. > > However, is there any shorthand way of getting srange to include the > final endpoint? While > srange(0,1.0,0.1,include_endpoint=True) > works, it is somewhat more verbose than

Re: [sage-support] accented characters

2011-03-26 Thread Mike Hansen
On Sat, Mar 26, 2011 at 10:53 AM, clodemil wrote: > Hi All, > > Sage does not accept accented characters like: é è à give a syntax > error. In Python 2.x, you can only use the following characters for identifiers (things like function names, variable names, etc.): the uppercase and lowercase let

Re: [sage-support] Random Permutation

2011-03-04 Thread Mike Hansen
On Sat, Mar 5, 2011 at 6:38 AM, Santanu Sarkar wrote: > How one can generate random permutation over the set > A={1,2,...,50}? sage: Permutations(50).random_element() [23, 40, 6, 27, 33, 43, 4, 34, 8, 39, 10, 35, 14, 41, 37, 21, 50, 44, 13, 46, 7, 16, 1, 2, 31, 38, 49, 36, 12, 28, 9, 5, 3, 19, 32

Re: [sage-support] A problem with assume real?

2011-02-17 Thread Mike Hansen
On Thu, Feb 17, 2011 at 1:30 PM, Jeff wrote: > To be consistent with the assumption, though, I think that iimag(a) > should return zero. > > Can anyone comment? Unfortunately, the assumptions system (still) is not well integrated with the Pynac code. When you do assume(a, 'real'), it sends that

Re: [sage-support] abs(pi*i) Bug

2011-02-07 Thread Mike Hansen
On Mon, Feb 7, 2011 at 1:27 PM, Laurent wrote: > Le 07/02/2011 13:05, Clemens Heuberger a écrit : >> >> I encountered the following bug: >> >> sage: abs(pi*I) >> I*pi > > I do not know how we define absolute value in Sage. Even in math in general, > I'm not sure of what means the absolute value of

Re: [sage-support] edges in graphs

2011-02-04 Thread Mike Hansen
On Fri, Feb 4, 2011 at 11:56 PM, Johannes wrote: > Hi list > how can i add a new edge (a->b) to a given graph G (n.n. connected), > just in the case that there is no path (a -> ... -> b) before? You should use "shortest_path" which uses Dijkstra's algorithm under the hood. if g.shortest_path('a'

Re: [sage-support] plotting box function

2011-02-04 Thread Mike Hansen
On Fri, Feb 4, 2011 at 11:37 AM, Renato Budinich wrote: > hello, > why is the below code plotting a flat function rather than a box one? When you do, > plot(box(x,1),(x,-3,3)) it evaluates box(x,1) which returns 0 because the variable x is not always less than 1. You need to delay the evaluati

Re: [sage-support] Re: [sage-devel] simplify_full on matrices.

2011-01-15 Thread Mike Hansen
On Sat, Jan 15, 2011 at 12:57 PM, Ivan Andrus wrote: > Perhaps Matrix (and other classes) should have a method .map() that does this, They already have an .apply_map() method which does this. --Mike -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this

Re: [sage-support] Catching a DeprecationWarning

2011-01-10 Thread Mike Hansen
On Mon, Jan 10, 2011 at 10:34 AM, Simon King wrote: > Is there a way to catch a deprecation warning (as one would catch an > error)? Nope, there's no way to catch it because there is no exception that is actually thrown. > Even stranger: When I catch the BaseException then the print statement >

Re: [sage-support] How to reinitialize a generator.

2011-01-09 Thread Mike Hansen
2011/1/9 José Díaz de Greñu y de Pedro : > Anything in the doc about that? http://docs.python.org/library/itertools.html#itertools.tee --Mike -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroup

Re: [sage-support] How to reinitialize a generator.

2011-01-09 Thread Mike Hansen
2011/1/9 José Díaz de Greñu y de Pedro : > But that's in Phyton is it exactly the same for SAGE? Yep! --Mike -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit thi

Re: [sage-support] Root finding of a real number

2010-12-22 Thread Mike Hansen
On Wed, Dec 22, 2010 at 9:15 PM, Santanu Sarkar wrote: > How one can find 10-th root of a real number in Sage? You can use the nth_root() method on real numbers: sage: f = 3.1415; f 3.141500 sage: f.nth_root(10) 1.12127904623998 --Mike -- To post to this group, send email to sage-supp

Re: [sage-support] sage-4.3.4 cygwin binary : pni flag missing complain

2010-12-07 Thread Mike Hansen
2010/12/7 Jorge E. ´Sanchez Sanchez : > Then I give again the ./sage command, ocurring something different: >    "Setting permissions of DOT_SAGE directory so only you can read and write > it ... > ... /home/ssgrge/sage-4.3.4/local/bin/sage-sage: line 206:3304 Illegal > instruction (core dumped) >

Re: [sage-support] bool(arcsin(x) == 2*arctan(x/(1+sqrt(1-x^2)))) returns false !!!

2010-12-06 Thread Mike Hansen
> Here's the same sort of thing in Mathematica. > > In[3]:= 12 == 2 > > Out[3]= False > > In[4]:= 1 == 1 > > Out[4]= True > > In[5]:= AcrSin[x] == 2 ArcTan[x/(1+Sqrt[1+x^2])] > >                                     x > Out[5]= AcrSin[x] == 2 ArcTan[] >                              

Re: [sage-support] Re: Bug (?) in full_simplify()

2010-11-20 Thread Mike Hansen
Hello, On Sat, Nov 20, 2010 at 7:57 AM, pevzi wrote: > Thank you for your reply. > But so is there any way to check if the expressions are equivalent? You could do this: sage: a = (x/(2*sqrt(x+1)) + 1/(2*sqrt(x+1)*(sqrt(x+1)+1))) sage: b = ((x*(sqrt(x+1)+1)+1)/(2*sqrt(x+1)*(sqrt(x+1)+1))) sage:

Re: [sage-support] Heavy memory usage in a basic example: what happens? How to manage it?

2010-11-16 Thread Mike Hansen
Hello, On Tue, Nov 16, 2010 at 12:03 AM, Rolandb wrote: > 3) I could not find a Sage standard function to deliver this output. > There are so many functions around... > Is there such a function? An iterator is even better. N.B.: The above > example "splits" can easily transformed to an iterator,

Re: [sage-support] sylvester matrix

2010-11-05 Thread Mike Hansen
On Fri, Nov 5, 2010 at 12:14 PM, andrew ewart wrote: > looking at the code in ticket 7458 (link:http://trac.sagemath.org/ > sage_trac/attachment/ticket/7458/sylvester.patch) > > however sage doesnt seem to recognise the sylvester_matrix command > any help Did you apply the patch and run "sage -b"

Re: [sage-support] Re: sage won't simplify something that is clearly 0

2010-11-03 Thread Mike Hansen
On Wed, Nov 3, 2010 at 9:16 AM, Rob H. wrote: > Here's another simple example of basic simplifications that aren't > processed: > > Does anyone have any suggestions on how to fix/circumvent these > problems? Don't use polynomial rings over the "symbolic ring" -- the "symbolic ring" is not really

Re: [sage-support] problem with a limit containing a fraction of square roots

2010-10-29 Thread Mike Hansen
On Fri, Oct 29, 2010 at 11:06 PM, Laurent wrote: > sage: f(x)=sqrt(x)/( sqrt(4-x)-sqrt(4+x)  ) > sage: f.limit(x=0) > Infinity > sage: limit(f,x=0) > Infinity > sage: f(0.0001) > -199.99983591 > > A computation by hand shows that the limit is actually -Infinity. Why does > Sage say Infinity ?

Re: [sage-support] Sage and skew-schurs

2010-10-27 Thread Mike Hansen
On Wed, Oct 27, 2010 at 10:42 PM, vasu wrote: > Hi all > I have tried searching all over the combinatorics sections of the > sagemath wiki, but I could not find if there is an implementation of > the skew schur functions. More specifically, here is what I am looking > for: > given partitions 'lamb

Re: [sage-support] Digits of precision in modified source code not works

2010-10-18 Thread Mike Hansen
On Mon, Oct 18, 2010 at 3:54 AM, gerrob wrote: > What am I doing in a wrong way? There are two differences. 1) In the .py file, there is no preparsing so the "1" in the function is of type "int" rather than Integer. This only causes a problem with #2. 2) You are not importing the exp function

Re: [sage-support] Array

2010-10-10 Thread Mike Hansen
On Sun, Oct 10, 2010 at 7:26 AM, Santanu Sarkar wrote: > How one can define an array of length 1? > For example an array of length 4 is X = [1, 2, 1, 3]. In that example, you are actually defining a Python *list* rather than an array. If you want to create a list of length 1, one way to

Re: [sage-support] modulo operator for reals

2010-10-09 Thread Mike Hansen
On Mon, Oct 4, 2010 at 4:03 PM, jpc wrote: > Why does the plot >   plot(k%1,k,0,5) For a symbolic variable 'k', k%1 returns an error since the notion of a symbolic mod operation hasn't been implemented. > What can be used instead ? You can just delay the application of the mod operation: plot(

Re: [sage-support] Re: Segfault when printing a copy of symbolic object

2010-10-07 Thread Mike Hansen
On Thu, Oct 7, 2010 at 6:25 PM, kcrisman wrote: > LOL - I should have checked that!  Trac isn't working now - can you or > drkirkby put that on the ticket when it comes back up? I've already posted a patch. --Mike -- To post to this group, send email to sage-support@googlegroups.com To unsubsc

Re: [sage-support] Re: Segfault when printing a copy of symbolic object

2010-10-07 Thread Mike Hansen
It doesn't have anything to do with the compiler warning messages -- the problem is much simpler than that. The implementation of Element.__copy__ is wrong for Expression objects; namely, it doesn't set the _gobj attribute. To fix it, we just add a proper implementation of copy: def __copy__

Re: [sage-support] Obtaining data from a computation in process

2010-09-13 Thread Mike Hansen
On Mon, Sep 13, 2010 at 12:41 PM, Alastair Irving wrote: > It depends on precisely what form your code takes.  If you're running the > loop at the top level with S as a global variable then you should be able to > do ctrl-c to terminate the computation and then look at S or any other > global vari

Re: [sage-support] field of symbolic expressions

2010-09-06 Thread Mike Hansen
On Mon, Sep 6, 2010 at 11:58 AM, Michael Beeson wrote: > I want to create a vector space of dimension three,  over the field > whose elements are symbolic expressions.   (The reason is that then I > can do vector calculus on such objects, which represent surfaces in 3- > space if the expressions d

Re: [sage-support] Nested if statements inside sum

2010-09-03 Thread Mike Hansen
One option is: sage: def F(m): :return sum([(m*n+2 if m*n < 0 else m*n+8) for n in range(1, 10)]) sage: F(2) 162 sage: F(-2) -72 --Mike -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegro

Re: [sage-support] compiling filename.sage generated by LaTeX on my Mac

2010-08-27 Thread Mike Hansen
On Fri, Aug 27, 2010 at 8:05 AM, Matt Boelkins wrote: > Matthew-Boelkinss-MacBook-Pro:z_sagetex_files boelkinm$ sage sagetest.sage > -bash: sage: command not found Your terminal does not know where to find the "sage" command. Suppose Sage in installed in /Applications/sage-4.5.2, then you should

Re: [sage-support] Octave and SAGE

2010-08-26 Thread Mike Hansen
On Thu, Aug 26, 2010 at 12:53 PM, calcp...@aol.com wrote: > Just wondering, since Octave is bundled with SAGE, is MPITB also > included? Octave is not bundled with Sage. If you have Octave already installed, then Sage can interface with it via pexpect. --Mike -- To post to this group, send em

Re: [sage-support] subtraction with dates?

2010-08-24 Thread Mike Hansen
On Tue, Aug 24, 2010 at 3:44 PM, Dr. David Kirkby wrote: > I was thinking of suggesting they use Sage to find out the number of days > between two dates, and so would realise that there is more than a month > between 5th Jan and the 16th Feb! sage: from datetime import date sage: date(day=16, mon

Re: [sage-support] Fixed Point arithmetic

2010-08-24 Thread Mike Hansen
On Tue, Aug 24, 2010 at 12:14 PM, Jeroen Demeyer wrote: > Is there any support in Sage for fixed point arithmetic?  That is, > computing with real numbers with a fixed number of bits after the > decimal point? Fredrik Johansson has done some work in this area -- see http://code.google.com/p/fastf

Re: [sage-support] Re: How can I make implicit_multiplication default?

2010-08-21 Thread Mike Hansen
On Sat, Aug 21, 2010 at 5:21 PM, Oscar Lazo wrote: > Unfortunately this is not available in the command line, so editing > ~/.sage/init.sage will not work. Can this too be set as a default? You can add the following to your init.sage: from sagenb.misc.support import automatic_names automatic_nam

Re: [sage-support] How can I make implicit_multiplication default?

2010-08-21 Thread Mike Hansen
On Sat, Aug 21, 2010 at 9:31 AM, Oscar Lazo wrote: > implicit_multiplication is very time-saving, and a much more natural way to > write things. I'd love it to be default. In your "~/.sage/init.sage" file, add the line "implicit_multiplication(True)". Whenever you start up Sage, file will be run

Re: [sage-support] Error installing cadabra.spkg

2010-08-20 Thread Mike Hansen
On Fri, Aug 20, 2010 at 3:29 PM, Oscar Gerardo Lazo Arjona wrote: > But I do have pcre installed on my ubuntu 10.04 system. If I type man pcre I > get a man page for perl compatible regular expressions, but if I type pcre I > get a pcre: command not found It is looking for pcre.h so my guess is t

Re: [sage-support] Re: How to deal with wrapper_descriptor / slot wrapper

2010-08-01 Thread Mike Hansen
Hello, > You say that "f.__pow__.__doc__ is a special read-only attribute", > apparently even at the time when the code is generated. But for what > reason? Is there really no way around? There isn't really any way to make f.__pow__.__doc__ work for your object. The reason is that there is reall

Re: [sage-support] Re: How to deal with wrapper_descriptor / slot wrapper

2010-07-31 Thread Mike Hansen
Hello, On Tue, Jul 27, 2010 at 12:30 PM, Simon King wrote: > So, the fact that RingElement uses the category framework means that > both the code *and* the documentation are hidden behind layers of > abstraction, to the extent that they can actually not be accessed *at > all*, unless one untars t

Re: [sage-support] scipy v0.8

2010-07-28 Thread Mike Hansen
Hello, On Wed, Jul 28, 2010 at 5:58 AM, maldun wrote: > When will the 0.8 release of scipy come into sage, the new fft > features would come quite in handy for me :) My guess maybe in about a month. We need to upgrade to Numpy 1.4.1 first and resolve and potential problems with both numpy and s

Re: [sage-support] Functional Programming doubt

2010-07-23 Thread Mike Hansen
On Thu, Jul 22, 2010 at 10:04 PM, vasu wrote: > I am iterating through all partitions p  of a fixed length k ( let's > say using Partitions (size, length=k) ) and I have written a function > which return true if the partition p has atleast one part repeated. I > can use a for loop and get the thin

Re: [sage-support] combinatoric polynomials

2010-07-20 Thread Mike Hansen
On Tue, Jul 20, 2010 at 4:26 PM, tvn wrote: > Is there a built in function that given  n variables and d degree, > compute all the terms from degree 0 to d  of all those n variables > combined ? The # of terms =  bionomial(n+d,d) Not built in exactly, but here is one way: sage: vars = (1,) + va

Re: [sage-support] Re: Importing a Python module

2010-07-20 Thread Mike Hansen
On Tue, Jul 20, 2010 at 9:27 AM, KvS wrote: > Thanks for the quick reply. I tried putting "from sage.all import *" Could you post the code to http://sage.pastebin.com . It is better to avoid using "import *" and explicitly list the things that you want to import. That way you know what objects

Re: [sage-support] Importing a Python module

2010-07-20 Thread Mike Hansen
On Tue, Jul 20, 2010 at 8:58 AM, KvS wrote: > So I guess my problem is that I don't understand the namespaces > involved somehow, since sgn() lives somewhere in a Sage namespace and > importing doesn't place the code in the right name space or something? > Does anybody have a hint how I may do thi

Re: [sage-support] needs *FAST* polynomial multiplication

2010-07-14 Thread Mike Hansen
On Wed, Jul 14, 2010 at 11:51 PM, Nathann Cohen wrote: > As such a polynomial will grow very large, very quickly, I would like > to know which would be the best way to compute it. I already know > things like the degree of the final polynomial, if it can help. Have you tried using the default pol

Re: [sage-support] types for functions

2010-07-14 Thread Mike Hansen
On Wed, Jul 14, 2010 at 3:04 PM, Johannes wrote: > scalar(v,w) = v * m * w This is not what you want to do -- http://sagemath.org/doc/tutorial/tour_functions.html explains what is going on here. > otherwise, if i do it this way, directly on the console it works: > v1 = vector([1,1,1]) > v2 = vec

Re: [sage-support] Re: how to change show() defaults for 2D plots

2010-07-11 Thread Mike Hansen
On Sun, Jul 11, 2010 at 5:48 AM, Johannes wrote: > Sorry for posting twice. I just noticed that > my first message got through all right, though I did > not receive it myself! Is this a common behaviour > of google groups? Yes, this is a common behavior of Google Groups and Gmail. You can set de

Re: [sage-support] ssse3 chipset instruction error

2010-07-09 Thread Mike Hansen
On Fri, Jul 9, 2010 at 11:19 PM, emil wrote: > Any help on how to solve this? Is a rebuild on the mentioned eeepc701 > necessary? > How do I build sage with the ssse3 flag disabled on a ssse3 machine > (for the next live CD / Puppy Linux sage package) to avoid this > problem? Did you build Sage w

Re: [sage-support] Notebooks "forgets" assignments

2010-07-09 Thread Mike Hansen
2010/7/9 Jorge E. ´Sanchez Sanchez : >> Which could be the problem? Thanks in advance, It sounds like the issue could be that the notebook closes the worksheet session after that amount of inactivity. Then, when you execute the next command it starts up a new session in which "a" hasn't been defi

Re: [sage-support] Re: Notebook appearance: change font and syntax highlighting?

2010-07-08 Thread Mike Hansen
On Thu, Jul 8, 2010 at 7:24 PM, Oscar Lazo wrote: > There is currently no syntax highlighting in the notebook interface. > Although it has been discussed. If you are interested in that feature, > you are encouraged to develop it. You will find the Sage community > very welcomming, lot's of people

Re: [sage-support] Re: sagetex question

2010-07-08 Thread Mike Hansen
On Thu, Jul 8, 2010 at 4:10 PM, Ra wrote: > Hi again, > > I have looked into my problem further, and it seems that the latex > command is at fault. As you see, only when I multiply the fraction by > -2 (or anything other than -1), do I get a negative fraction back from > the latex command. > > Any

Re: [sage-support] Strange (non-)substitution of variables whose name is not equal to their representation

2010-07-07 Thread Mike Hansen
On Wed, Jul 7, 2010 at 1:49 PM, David Sanders wrote: > I now want to substitute eps=1, so I do > > a.subs(eps = 1) > > but the response is still 3*epsilon ! This is due to the way Python functions work. Basically, doing a.subs(eps=1) is the same as doing a.subs(**{'eps': 1}) When you use key

Re: [sage-support] Unexact computation with roots of unity.

2010-06-22 Thread Mike Hansen
On Tue, Jun 22, 2010 at 11:13 AM, cjung wrote: > 0.?e-18 + 0.?e-19*I # <- should be zero I think that this is actually zero -- it just doesn't know it for sure yet. For example, look at sage: q = QQbar(sqrt(2)) sage: p = q^2; p 2.000? sage: p == 2 True sage: p 2 and sage: q = QQba

Re: [sage-support] Re: Computing in the Eisenstein Integers

2010-06-17 Thread Mike Hansen
On Thu, Jun 17, 2010 at 6:28 AM, kcrisman wrote: >> Would it be very difficult to implement a GaussianIntegers class which >> used the symbolic I instead but otherwise used number fields behind >> the curtain, or would that lead to too much confusion elsewhere? >> The symbolic I uses number field

Re: [sage-support] strange behavior in matrix substitution

2010-06-12 Thread Mike Hansen
Hello, On Sat, Jun 12, 2010 at 5:27 PM, Byungchul Cha wrote: > > Please tell me if this is a bug, or, I'm missing something obvious... > > sage: a = 3 # Assign a value to a variable a > sage: b = a # Create a copy of a This does not create a copy of a. When you do "a = 3", this creates new obje

Re: [sage-support] Data list

2010-06-10 Thread Mike Hansen
On Thu, Jun 10, 2010 at 12:03 PM, Marco Boretto wrote: > I'm tryng to to simple thing like this: > > m=[0.6158, 0.5893, 0.5682, 0.51510, 0.4980, 0.4750, 0.5791, > 0.5570,0.5461, 0.4970, 0.4920, 0.4358, 0.422, 0.420] > m.count > > i want to know the number of the object but the answer is You use l

Re: [sage-support] intallation fédora 13 fail

2010-06-08 Thread Mike Hansen
On Tue, Jun 8, 2010 at 7:08 AM, Jaap Spies wrote: > > This is now: http://trac.sagemath.org/sage_trac/ticket/9189 > There's an spkg at http://sage.math.washington.edu/home/mhansen/libgcrypt-1.4.4.p3.spkg which needs to be tested. --Mike -- To post to this group, send email to sage-support@goog

Re: [sage-support] Using _ is a lot less helpful than it used to be

2010-06-02 Thread Mike Hansen
On Wed, Jun 2, 2010 at 4:49 PM, Dan Drake wrote: > And yes, I get this problem too. It seems a little unpredictable, and is > quite irritating. It happens after you view a docstring's help from the command line. I made a patch a few weeks ago which fixes this: http://trac.sagemath.org/sage_trac

Re: [sage-support] question about what is in vector spaces/free modules

2010-06-02 Thread Mike Hansen
On Wed, Jun 2, 2010 at 6:43 AM, kcrisman wrote: > I guess I'm a little surprised that even the 'vector' explicitly > designated to be over ZZ isn't in there.  I think I'm mathematically > ok, too - taking the echelon matrix, isn't my vector e_1-(1/3)*e_2- > (1/3)*e_3 ? Those coefficients are in Q

Re: [sage-support] Retrieving xy data from implicit plots

2010-06-02 Thread Mike Hansen
On Wed, Jun 2, 2010 at 12:22 AM, William Stein wrote: > Can you elaborate?   I read the above and thought. > > Claim: your proposal is worse. > Proof: argument that my proposal is better. We already do most of what you suggested. In contour_plot.py, we have the following: g, ranges = setup_for_

Re: [sage-support] Retrieving xy data from implicit plots

2010-06-02 Thread Mike Hansen
On Wed, Jun 2, 2010 at 12:05 AM, William Stein wrote: > Then at least the data can be manipulated.  It seems to me that this > should be at least as good as what we already have.    The student > suggested that one could even "connect the dots" using some sort of > optimization technique. > > But

Re: [sage-support] Sage to run automatically in CLI

2010-05-30 Thread Mike Hansen
On Sun, May 30, 2010 at 10:33 AM, clinton bowen wrote: > How do i get sage to run automatically in ubuntu?  I would like to be > able to just type 'sage' into command line and it run automatically. I > used to be able to do this but not anymore and I forgot how to do > this. I would just add the

Re: [sage-support] set graphic viewer Linux

2010-05-30 Thread Mike Hansen
On Sun, May 30, 2010 at 11:13 AM, bb wrote: > Sage Graphic will always open Gimp on my Linux System (debian lenny). I > would like to display Sage graphics for instance with Kview or some program > else. I am not shure if this is a Sage issue, may be it is a debian default. On your system, Sage i

Re: [sage-support] error occurred while installing libgcrypt-1.4.4.p2: undefined reference to symbol 'gpg_strerror'

2010-05-28 Thread Mike Hansen
On Fri, May 28, 2010 at 12:36 AM, domingo.domingogo...@gmail.com wrote: > I have all the necessary packages and additionally libgpg-error and > libgpg-error-devel, could someone point > what I am doing wrong? There seems to be a problem with libgcrypt on Fedora 13. It looks like someone found a

Re: [sage-support] Re: Sage Server - displaying Matlab plots

2010-05-27 Thread Mike Hansen
On Thu, May 27, 2010 at 5:42 PM, William Stein wrote: > Mike Hansen just told me that he will create such a patch tonight. > > http://trac.sagemath.org/sage_trac/ticket/9072 The patch is up, but it needs testing. --Mike -- To post to this group, send email to sage-support@googlegrou

Re: [sage-support] Re: Sage Server - displaying Matlab plots

2010-05-27 Thread Mike Hansen
On Thu, May 27, 2010 at 12:20 PM, Vinod wrote: > The simple matlab code I wrote is functional and the sage server > evaluates it properly and generates the .png file in the directory. > Only issue that need to be addressed is to get the sage server display > the plots in the web browser itself. We

Re: [sage-support] Real exponential of a matrix (rotations with sine/cosine)

2010-05-25 Thread Mike Hansen
Hello, On Tue, May 25, 2010 at 8:11 AM, Laurent wrote: > By the way, the following works : > > sage : a=R[0][0] > sage: (a.real_part()+a.imag_part()).simplify_full() > cos(x) > > How can automatize it ? Is it a way to perform this simplification to each > element of the matrix at once ? You can

Re: [sage-support] Re: Are there statisticians using Sage?

2010-05-20 Thread Mike Hansen
On Thu, May 20, 2010 at 11:33 AM, Nathann Cohen wrote: > Hello ! > >>    * As a language, Python is vastly superior to R.   Python has good >> support for object oriented programming, a very wide selection of >> existing programs and libraries, and supports threads for handling >> realtime data.  

Re: [sage-support] Racional coercion of very small decimal numbers.

2010-05-07 Thread Mike Hansen
On Fri, May 7, 2010 at 10:18 AM, Rodrigo wrote: > I would expect both results to be the same. Am I missing something? This is due to the way floating point numbers are stored. In particular, sage: a = 6.62606896 sage: b = 66.2606896 sage: a.exact_rational() * 10 == b.exact_rational() False See

Re: [sage-support] How to add grid lines to plot?

2010-05-06 Thread Mike Hansen
Hello, On Thu, May 6, 2010 at 2:14 PM, Eugene Goldberg wrote: > There is plot: > p=plot(sin,0,2*pi) > How to add grid lines to this plot? plot(sin,0,2*pi,gridlines=True) You can look at the documentation of plot for more examples. --Mike -- To post to this group, send email to sage-support@g

Re: [sage-support] enumerating vectors in Z^d

2010-05-06 Thread Mike Hansen
On Wed, May 5, 2010 at 8:38 PM, Eva wrote: > Problem background: I want to find the smallest (in l_2 norm) vectors > that satisfy a certain condition - creating a minimal set of coset > representatives for Z^d/A(Z^d) for a dilation matrix A.  My idea is to > start with 0 and enumerate vectors in Z

Re: [sage-support] coersion

2010-05-05 Thread Mike Hansen
On Wed, May 5, 2010 at 2:46 PM, Matt Bainbridge wrote: > Sage knows how to coerce from Frac(ZZ[x]) to Frac(QQ[x]).  There is no > coercion going the other way, though there should be one, since these > two rings are equivalent.  Is there a reasonable way for me to define > my own coercion? This d

Re: [sage-support] rational function arithmetic

2010-05-05 Thread Mike Hansen
Hello, On Wed, May 5, 2010 at 5:03 AM, Matt Bainbridge wrote: > I wrote a sage program which does a lot of arithmetic in the field of > rational functions Frac(Q[x,y,z]).  The problem is that sage doesn't > check for common divisors of the numerator and denominator, so after > doing a lot of arit

Re: [sage-support] Elliptic curves over non-fields

2010-05-03 Thread Mike Hansen
Hello, On Mon, May 3, 2010 at 1:34 AM, Jeroen Demeyer wrote: > I wonder if there is a way to make Sage compute with points on elliptic > curves over rings which are not fields.  I would like to use it for a > very simple implementation of the ECM factoring algorithm for > educational purposes.  T

Re: [sage-support] Eigenvectors in radical format

2010-05-01 Thread Mike Hansen
On Sat, May 1, 2010 at 7:33 AM, Laurent wrote: > How do I have to compute eigenspaces of a matrix whose eigenvalues could > have a radical form ? > Or, alternatively, how do I get the radical form of a in the following ? > sage: sageMatrix = matrix(QQ,[[1,4],[4,2]]) > sage: eig=sageMatrix.eigenval

Re: [sage-support] Re: plot

2010-04-22 Thread Mike Hansen
On Thu, Apr 22, 2010 at 4:02 AM, pallab wrote: > Thanks for your answer, it seems complicated :). I think it would have > been much easier if there is a way to specify the plot variables in > plots. Plotting has basically the exact same syntax as Mathematica: plot(4*x, (x,0,2)) Just lower case

Re: [sage-support] plot

2010-04-22 Thread Mike Hansen
On Thu, Apr 22, 2010 at 2:46 AM, pallab wrote: > How to specify variables in a plot. Say I want something  like, > > At=S*x > S=4 > plot(At,0,2) Note that when you assign "S=4" it doesn't change the value of S in At. > > Question is that how to instruct Sage that I want to plot against 'x'. > In

Re: [sage-support] non-commutative variable?

2010-04-21 Thread Mike Hansen
On Wed, Apr 21, 2010 at 8:06 PM, Elizabeth wrote: > I'd like to know if there is a way to define a variable in Sage as non- > commuting. > I am trying to multiply a few matrices together and take the trace, > *but* the matrix entries are non-commuting elements.  So, I need Sage > to leave a term l

Re: [sage-support] Finite fields and tab completion?

2010-04-20 Thread Mike Hansen
On Tue, Apr 20, 2010 at 9:59 PM, Alasdair wrote: > does nothing.  I'm sure this is new behaviour (I'm using version > 4.3.2) - but how do I create a field whose name allows tab completion > of its methods? I believe this is specific to 4.3.2. The patch which fixes this and was merged in 4.3.3 is

Re: [sage-support] power function with runtime error

2010-04-19 Thread Mike Hansen
On Mon, Apr 19, 2010 at 11:52 AM, bb wrote: > I get a runtime error, but just would expect infinity! Is there something > wrong or any explanation? This is because when you do 2^3^4^5 you are computing that number exactly as an integer, and that number is definitely not infinity. If you wanted t

Re: [sage-support] Re: solve polynomial limit?

2010-04-16 Thread Mike Hansen
On Fri, Apr 16, 2010 at 9:57 AM, Owen wrote: >> Is there another method approach I could take?  I'd like to reach 6 at >> least.  My homework depends on it!  :) There's also the to_poly_solve argument to solve which cause non-exact answers to be returned: sage: solve([x^5 - x^4 - x^3 - x^2 - x -

Re: [sage-support] ReST raw format

2010-04-15 Thread Mike Hansen
On Thu, Apr 15, 2010 at 5:08 AM, bb wrote: > Is it possible to get the raw format of the documentation texts? > Just now one may get > 1st a copy from the text presentet in the browser window > 2nd a copy of the HTML-source On any of the documentation pages, there is a "Show Source" link. For ex

  1   2   3   4   5   >