[sage-support] Re: problem with plot3d on ubuntu 64bit
I never got Java to work period in Ubuntu for 64bit with the Firefox built for it. I later reinstalled Ubuntu but the 32bit edition, and Java worked fine in the notebook. On Jan 15, 6:20 pm, David Joyner [EMAIL PROTECTED] wrote: Hi: I have an old 64 bit machine with 64bit ubuntu fiesty fawn loaded on it. I just noticed a problem with plot3d on it (plot3d runs fine on my intel macbook): In sage 2.9.3: sage: x = var(x) sage: y = var(y) sage: p = plot3d(x^2-y^2,(-1,1),(-1,1)) --- type 'exceptions.NameError' Traceback (most recent call last) /mnt/drive_hda1/sagefiles/sage-2.9.alpha5/ipython console in module() type 'exceptions.NameError': name 'plot3d' is not defined In sage 2.10.alpha1: sage: x = var(x) sage: y = var(y) sage: plot3d(x^2-y^2,(-1,1),(-1,1)).show() ## long time but nothing happens sage: p = plot3d(x^2-y^2,(-1,1),(-1,1)) sage: p ## long time but nothing happens sage: type(p) type 'sage.plot.plot3d.parametric_surface.ParametricSurface' sage: show(p)## nothing happens sage: I'll try running sage -testall on both of these to see if something pops up. - David Joyner --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://sage.math.washington.edu/sage/ and http://sage.scipy.org/sage/ -~--~~~~--~~--~--~---
[sage-support] Re: associativity of addition on ell. curves
Dear David, one can construct the apropriate quotient P=QQ['a,b,x1,x2,x3'] K.a,b,x1,x2,x3=FractionField(P) R.y1,y2,y3=K['y1,y2,y3'] I = R.ideal(y1^2 - x1^3 - a*x1 - b, y2^2 - x2^3 - a*x2 - b, y3^2 - x3^3 - a*x3 - b) S=quotient(R,I) sage: S Quotient of Multivariate Polynomial Ring in y1, y2, y3 over Fraction Field of Multivariate Polynomial Ring in a, b, x1, x2, x3 over Rational Field by the ideal (y1^2 - x1^3 - a*x1 - b, y2^2 - x2^3 - a*x2 - b, y3^2 - x3^3 - a*x3 - b) but E = EllipticCurve(S, [a,b]) does not work (in the present implementation) Andrzej On 15 Sty, 08:35, achrzesz [EMAIL PROTECTED] wrote: Dear Carl, I like your code; it is elegant and realy quick but it seems that finishing your code with Paul Zimmermann's approach I2 = singular(I).groebner() #print I.reduce(n12); print singular.reduce((n12), I2) (althout less elegant) is a little bit faster (0.06 - 0.05 on my comp. :) I want also add a question to David Harvey questions In experimenting with Lenstra factorization method one needs multiplication on ell.curv over the Ring Z_(p*q). GP-Pari allows for that so I'm doing it using gp interface. Is it difficult to implement a similar functionality in Sage? Andrzej Chrzeszczyk On 15 Sty, 03:28, David Harvey [EMAIL PROTECTED] wrote: On Jan 14, 2008, at 10:09 PM, Carl Witty wrote: Here is a more idiomatic way to do this computation in Sage. We work in the fraction field of a multivariate polynomial ring; this means that our polynomial arithmetic is handled by libSingular instead of by maxima, and that we can get the numerator directly with numerator, since fraction field elements are always normalized. Also, we use Sage's wrapper of ideals and Groebner bases (which I believe is implemented with libSingular), rather than calling Singular. (Avoiding the call to factor(s1-s2) means that this version is much faster.) sage: R.x1,y1,x2,y2,x3,y3,a,b = QQ[] sage: eq1 = y1^2 -(x1^3+a*x1+b) sage: eq2 = y2^2 -(x2^3+a*x2+b) sage: eq3 = y3^2 -(x3^3+a*x3+b) sage: lambda12 = (y1 - y2)/(x1 - x2) sage: x4 = (lambda12*lambda12 - x1 - x2) sage: nu12 = (y1 - lambda12*x1) sage: y4 = (-lambda12*x4 - nu12) sage: lambda23 = ((y2 - y3)/(x2 - x3)) sage: x5 = (lambda23*lambda23 - x2 - x3) sage: nu23 = (y2 - lambda23*x2) sage: y5 = (-lambda23*x5 - nu23) sage: s1 =(x1 - x5)*(x1 - x5)*((y3 - y4)*(y3-y4) - (x3+x4)*(x3-x4)* (x3- x4)) sage: s2 =(x3 - x4)*(x3 - x4)*((y1 - y5)*(y1-y5) - (x1+x5)*(x1-x5)* (x1- x5)) sage: n12 = numerator(s1-s2) sage: I = ideal([eq1,eq2,eq3]) sage: I.reduce(n12) 0 What would be *really* nice is if we could work directly in the fraction field of the quotient of R.x1,y1,x2,y2,x3,y3,a,b by the appropriate ideal. (Does that even make sense? Is the ideal prime?) I tried to do this but Sage gave up pretty quickly on me. A nice encore would be to do this using Sage's elliptic curve class to do the actual arithmetic. After all EllipticCurves can be defined over any field Here's my dream session: sage: R.x1,y1,x2,y2,x3,y3,a,b = QQ[] sage: I = R.ideal(y1^2 - x1^3 - a*x1 - b, y2^2 - x2^3 - a*x2 - b, y3^2 - x3^3 - a*x3 - b) sage: S = FractionField(R.quotient(I)) # currently barfs sage: E = EllipticCurve(S, [a, b]) sage: P1 = E(x1, y1) sage: P2 = E(x2, y2) sage: P3 = E(x3, y3) sage: (P1 + P2) + P3 == P1 + (P2 + P3) True Here's the traceback in the FractionField line: /Users/david/sage-2.9/local/lib/python2.5/site-packages/sage/rings/ fraction_field.py in FractionField(R, names) 104 if not ring.is_Ring(R): 105 raise TypeError, R must be a ring -- 106 if not R.is_integral_domain(): 107 raise TypeError, R must be an integral domain. 108 return R.fraction_field() /Users/david/sage-2.9/local/lib/python2.5/site-packages/sage/rings/ quotient_ring.py in is_integral_domain(self) 226 227 -- 228 return self.defining_ideal().is_prime() 229 230 def cover_ring(self): /Users/david/sage-2.9/local/lib/python2.5/site-packages/sage/rings/ ideal.py in is_prime(self) 275 276 def is_prime(self): -- 277 raise NotImplementedError 278 279 def is_principal(self): This suggests maybe the only barrier here is checking primality of the ideal? After that, the fraction field magic should just work right? But surely there is code somewhere to check primality, isn't this in singular or something? I don't know anything about the implementation of multivariate polynomial rings, maybe someone else can help out here. david --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this group at
[sage-support] Re: Doc Day 1 announcement: January 17th, 2008, 9am-5pm PST
Someone at the joint meetings looked over the tutorial and came back the next day with a suggestion: he wanted a list of all the pre- defined words in sage, e.g. list, Random, RationalField, etc. I told him it would be a long list, and showed him the auto-completion. He said maybe a short list of the most important commands would still be nice; one of his points was that while reading the tutorial he became confused about what was already defined vs. what was being defined. I'm not sure what the best way to address that would be, but I promised him I would bring it up. -Marshall On Jan 16, 10:52 am, mabshoff [EMAIL PROTECTED] dortmund.de wrote: Hello folks, after doing 8 Bug Days the time has come to spend some more time on the documentation, which also needs a lot of work. So after some discussion in IRC we decided to get together in IRC at the above date and time and work on the documentation. Since we have never done a doc day it isn't 100% clear to me how the whole thing will go down, but I assume we will just go with the flow. Thoughts? Suggestions? Cheers, Michael --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://sage.math.washington.edu/sage/ and http://sage.scipy.org/sage/ -~--~~~~--~~--~--~---
[sage-support] jmol and Camino
Hi, I get this error on Camino (... don't ask): ReferenceError: _jmolInitCheck is not defined While trying to plot this: var('v') f1 = (x, x*sin(v), x*cos(v)) parametric_plot3d(f1, (x,0,5), (v,0,2*pi)) I must say that this works great in Firefox and Safari, I just thought you would like to know. Best, -- Hector --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://sage.math.washington.edu/sage/ and http://sage.scipy.org/sage/ -~--~~~~--~~--~--~---
[sage-support] Elementary symmetric function expansion (bug?)
Hi, First of all, I don't know whether to post about this here, or file a bug report, or something else. I'm new to community developed software. Let me know if this isn't the right place to report the following. I'm using sage 2.9.2 and 2.9.3 on different (both PPC) machines, both running OS X 10.4. Consider the following commands which should expand an elementary symmetric function as a polynomial in 3 variables: # Code sage: k=SFAElementary(QQ) sage: f=k([2]) sage: f.expand(3) --- type 'exceptions.TypeError' Traceback (most recent call last) /Users/bjones/ipython console in module() /Applications/sage/local/lib/python2.5/site-packages/sage/combinat/ sfa.py in expand(self, n, alphabet) 1041 #if len(part) n: 1042 #continue - 1043 res += self_mc[part] * resPR(e(part, n, alphabet)) 1044 return res 1045 /Users/bjones/multi_polynomial_libsingular.pyx in sage.rings.polynomial.multi_polynomial_libsingular.MPolynomialRing_libsingular.__call__() /Users/bjones/multi_polynomial.pyx in sage.rings.polynomial.multi_polynomial.MPolynomial._polynomial_() /Users/bjones/multi_polynomial_libsingular.pyx in sage.rings.polynomial.multi_polynomial_libsingular.MPolynomialRing_libsingular.__call__() /Applications/sage/local/lib/python2.5/site-packages/sage/rings/ rational_field.py in __call__(self, x, base) 180 if isinstance(x, sage.rings.rational.Rational): 181 return x -- 182 return sage.rings.rational.Rational(x, base) 183 184 def construction(self): /Users/bjones/rational.pyx in sage.rings.rational.Rational.__init__() /Users/bjones/rational.pyx in sage.rings.rational.Rational.__set_value() /Users/bjones/polynomial_element.pyx in sage.rings.polynomial.polynomial_element.Polynomial._rational_() type 'exceptions.TypeError': cannot coerce nonconstant polynomial # End Similar commands, but on objects from SFASchur, SFAPower, SFAMonomial, etc.., work perfectly. It seems to be only SFAElementary that has a problem. Thanks, BFJ --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://sage.math.washington.edu/sage/ and http://sage.scipy.org/sage/ -~--~~~~--~~--~--~---
[sage-support] Re: Elementary symmetric function expansion (bug?)
On Jan 16, 10:30 pm, mabshoff [EMAIL PROTECTED] dortmund.de wrote: On Jan 16, 10:23 pm, BFJ [EMAIL PROTECTED] wrote: Hi, First of all, I don't know whether to post about this here, or file a bug report, or something else. I'm new to community developed software. Let me know if this isn't the right place to report the following. I'm using sage 2.9.2 and 2.9.3 on different (both PPC) machines, both running OS X 10.4. Consider the following commands which should expand an elementary symmetric function as a polynomial in 3 variables: SNIP Hi BJF, I can reproduce the problem with 2.9.3, but it is fixed in 2.10.alpha3: -- | SAGE Version 2.10.alpha3, Release Date: 2008-01-14 | | Type notebook() for the GUI, and license() for information.| -- sage: sage: k=SFAElementary(QQ) sage: sage: f=k([2]) sage: sage: f.expand(3) x0*x1 + x0*x2 + x1*x2 sage: Mike Hansen did post a patch that was merged early on in 2.10.alpha2 or so. He might be able to tell you if the bug was fixed on purpose or by accident. 2.10 should be out by the weekend or maybe Monday. Alpha3 is available at http://sage.math.washington.edu/home/mabshoff/release-cycles-2.10/sag... Mike: could you add a doctest that tests this behavior? Oops, looking at the actual failure this is a libSingular bug. I can't find the exact ticket number at the moment, but I remember merging that fix. Sorry for the double post. Similar commands, but on objects from SFASchur, SFAPower, SFAMonomial, etc.., work perfectly. It seems to be only SFAElementary that has a problem. Thanks, BFJ Cheers, Michael Cheers, Michael --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://sage.math.washington.edu/sage/ and http://sage.scipy.org/sage/ -~--~~~~--~~--~--~---
[sage-support] Re: Doc Day 1 announcement: January 17th, 2008, 9am-5pm PST
On Wednesday 16 January 2008, Ted Kosan wrote: I have 2 documentation requests: 1) If someone could add documentation to each of the javascript functions in the notebook, that would be helpful. 2) Add a short .txt file in each of the directories in the Sage distribution which describes what the purpose of the directory is (perhaps named purpose.txt?). I think this would be helpful for people like me who are in the process of learning how to contribute code and packages to Sage. Thanks :-) Ted I'd like to second that. I don't think I'll be able to participate in this doc day due to the time difference (I think PST is 8 hours behind GMT). Also, I've never used IRC before so it will be a learning exercise if I do get to join in. Bill -- +---+ | Bill Purvis, Amateur Mathematician| | email: [EMAIL PROTECTED] | | http://bil.members.beeb.net | +---+ --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://sage.math.washington.edu/sage/ and http://sage.scipy.org/sage/ -~--~~~~--~~--~--~---
[sage-support] Re: creating users for a sage notebook
On Jan 16, 2008 2:26 PM, William F Hammond wrote: William -- You may recall that we met briefly in San Diego. After seeing Sage there, I've moved to set it up, and I'm rather pleased with what I've seen so far. I gather that it's supposed to be possible to associate users other than admin with a sage notebook (in a secure way), and I suppose that the user names have nothing to do with usernames in the operating system, but I don't see directions for creating such notebook users aside from the following bit in the sage command line output for notebook?: accounts -- (default: False) if True, any visitor to the website will be able to create a new account. If False, only the admin can create accounts (currently, this can only be done by running with accounts=True for a few minutes, or on the command line with, e.g., nb = load('sage_notebook/nb.sobj') nb.set_accounts(True) nb.add_user(username, password, [EMAIL PROTECTED], user) nb.save() I assume the command line sequence should be run in the sage command line interface before launching the affected notebook. But when I execute the first of those four lines I see: sage: nb = load('sage_notebook/nb.sobj') --- type 'exceptions.IOError' Traceback (most recent call last) /home/hammond/ipython console in module() /home/hammond/sage_object.pyx in sage.structure.sage_object.load() type 'exceptions.IOError': [Errno 2] No such file or directory:\ 'sage_notebook/nb.sobj' I cannot make sense out of these traceback lines. (/home/hammond/.sage/sage_notebook/nb.sobj does exist.) First type cd $HOME/.sage I also do not understand how to create accounts when running notebook() for a short while with accounts=True. Click on the link that says Sign up for a new SAGE Notebook account right below the login box. William Thanks for any clues you may be able to provide. -- Bill -- William Stein Associate Professor of Mathematics University of Washington http://wstein.org --~--~-~--~~~---~--~~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://sage.math.washington.edu/sage/ and http://sage.scipy.org/sage/ -~--~~~~--~~--~--~---