[sage-devel] Re: New Axiom package
On 8/14/07, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote: > > It build fine but the I/O is seriously bugged for me, at least for the > displaying of rational functions (but probably wrapping too)... > ... Oops, for the full axiom interface to work reliably you need sage-2.8.1 or the following patch: http://sage.math.washington.edu/home/page/axiom.py-0.3.patch (Not required just to run 'sage -axiom'.) On 8/13/07, William Stein <[EMAIL PROTECTED]> wrote: > Thanks. I've applied this for sage-2.8.1. > > On 8/13/07, Bill Page <[EMAIL PROTECTED]> wrote: > > On 8/13/07, William Stein <[EMAIL PROTECTED]> wrote: > > > Please resend your patch as an email attachment or post it somewhere > > > and send me a link. Your emailing + my email client mangled it. > > > > > > > Sorry. Here's a link: > > > > > > http://sage.math.washington.edu/home/page/axiom.py-0.3.patch > > > > > > > -- > William Stein > Associate Professor of Mathematics > University of Washington > http://www.williamstein.org > Regards, Bill Page. --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: ticket #59 (optimize elliptic curve arithmetic)
I think repeated squaring is not more efficient in all cases. For example over Z it is only more efficient if your multiplication algorithm is faster than the naive one (which is O(n^2) in the number of bits). So in the case of elliptic curves it really depends on the efficiency of the addition steps and on how fast the numerators/denominators grow. I am sure the elliptic curve experts on this forum can say more about this. Michel On Aug 14, 10:41 pm, David Harvey <[EMAIL PROTECTED]> wrote: > On Aug 14, 2007, at 4:19 PM, Alex Ghitza wrote: > > > > > Hi, > > > I've looked at ticket #59, in which David Kohel says: > > > William, my student noticed some slow performance with elliptic > > curves > > group law. I think there was a huge overhead in duplication: > > > sage: E = EllipticCurve([GF(101)(1),3]) > > sage: P = E([-1,1,1]) > > sage: timeit 2*P > > 100 loops, best of 3: 3.81 ms per loop > > sage: timeit P+P > > 1000 loops, best of 3: 1.81 ms per loop > > > Basically n*P was passing through all sorts of high-level layers for > > group schemes, abelian groups, and the like. > > > It turns out that n*P is inherited from multiplication in > > AdditiveGroupElement, which itself goes all the way to ModuleElement. > > But at that point, it is implemented efficiently (by repeated > > doubling). > > It seems to me that the loss of efficiency in comparison with P+P > > +...+P > > is all in the overhead. > > > I experimented with this and noticed that P+...+P is faster than > > n*P as > > long as n is smaller than 10, then the repeated doubling starts to > > bear > > fruit. So inside ell_point.py, I overloaded the multiplication > > operator, simply making it return P+...+P if n is less than 10, and > > return n*P if n is larger than 10 (see below). > > > This is not a very elegant solution. For one thing, the magic > > number 10 > > that works on my architecture (Intel Core Duo running Gentoo) might > > not > > be right for a completely different machine. Is there a way to > > automatically generate this number when the user compiles sage from > > source (by starting with a value like 100, say, and then doing a > > "binary > > search" until it hits the optimal number)? > > Even on a single architecture, the correct crossover point will vary > between different base fields. > > For example if the base field is Q, the best crossover will depend on > the size of the coordinates. If the coordinates are enormous, it's > probably best to do repeated doubling as much as possible. Similarly > if you are working over GF(p) for some very large p, you'll find a > different crossover from GF(101). > > This is going to be a difficult problem in general. I think a better > approach is to try to reduce the overhead, i.e. to get the n*P > pathway working faster, and/or to reduce the overhead in the addition > law itself. > > david --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: Getting Geometry in SAGE
On Tue, 14 Aug 2007, William Stein wrote: > > On 8/14/07, John Cremona <[EMAIL PROTECTED]> wrote: >> >> Whatever the fact of the case, I find the tone of this email >> offensive. I'm sure I am not the only one. I think an apology is >> called for. > > I'm very sorry for not carefully reading the email before moderating > it "send" -- I usually just check that postings aren't blatant spam. > In the future I'll be more careful when moderating messages. > > -- William i interprted john's statement to mean that arash owes you, william, an apology, and if that is what he meant, i heartily agree. -kyle > >> On 8/14/07, Arash Partow <[EMAIL PROTECTED]> wrote: >>> >>> Hi William, >>> >>> I believe you are the one that is confused about the GPL. >>> The GPL is a form of copyright that is provided with a >>> piece of work. In this case I permit people to use Wykobi >>> under the terms of the GPL open source license as opposed >>> to say the CPL, LGPL or the Mozilla license. >>> >>> As for commercial use, if you bothered to understand the >>> FSF's definition of derivative works you would know it is >>> a necessity for anyone wanting to use the code in a commercial >>> environment not to have to deal with the viral effects of the >>> GPL. >>> >>> I suggest in the future you actually "know" or at least become >>> familiar with what you are talking about before you make such >>> comments, as it does not bode well for someone that supposedly >>> has tenure. >>> >>> >>> Arash Partow >>> __ >>> Be one who knows what they don't know, >>> Instead of being one who knows not what they don't know, >>> Thinking they know everything about all things. >>> http://www.partow.net >>> >>> >>> >>> >>> On Aug 11, 9:47 am, "William Stein" <[EMAIL PROTECTED]> wrote: > Another GPL one is Wykobi > (www.wykobi.com). Wykobi has no documentation, is very new (version 0.0.03), and I'm pretty sure the author(s) are totally confused about what the GPL is. The claim to make Wykobi available under the GPL but only for non-commercial use -- which doesn't make sense because "GPL + non-commercial" is non-GPL compatible license. The GPL doesn't allow restrictions against commercial use. Anyway, Wykobi unfortunately isn't really GPL'd. -- William >>> >>> >>> >> >> >> -- >> John Cremona >> >>> >> > > > -- > William Stein > Associate Professor of Mathematics > University of Washington > http://www.williamstein.org > > > --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: New Axiom package
On 8/14/07, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote: > > It build fine but the I/O is seriously bugged for me, at least for the > displaying of rational functions (but probably wrapping too)... Thanks for the report. I will look into it. Do other things work ok? E.g. this script that posted earlier: > For example: > > sage: r=[0,3,32,375,5184,84035] > sage: R=axiom(r).guessExpRat(); R > >n > [[function= n (n + 2) ,order= 0]] > > sage: f=axiom('%s.1.function'%N.name()); f > >n > n (n + 2) > > sage: f.type() > Expression Integer > > sage: f.eval(axiom('n=5')) > 84035 > > sage: f.eval(axiom('n=6')) > 1572864 > > sage: Regards, Bill Page. --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: Getting Geometry in SAGE
On 8/14/07, David Joyner <[EMAIL PROTECTED]> wrote: > New members are moderated. I didn't post FlameGuy's > email but didn't delete it either. (Both William and I are moderators. > It was addressed to William and didn't think it appropriate > for me to delete emails adressed to him, but I did tell him > I would not post it for obvious reasons.) > I do think it is important to maintain decorum and Iapologize. > I'll try to improve communication with William so this doesn't > happen again. > Just to add to this, since sage-devel is a high traffic list, the spammers have found it and now probably about 1/3 of the posts by new members are spam emails. So by moderating new posts we are saving you from the trouble of getting spam via sage-devel. William --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: quaddouble in sage, and other floating point issues
On 8/14/07, cwitty <[EMAIL PROTECTED]> wrote: > On Aug 14, 12:59 am, Jonathan Bober <[EMAIL PROTECTED]> wrote: > > This is exactly what NTL does in its quad float class. Just about every > > function starts and ends with a macro to adjust the fpu, resulting in > > around 7 extra assembly instructions. In the following code, the > > overhead is quite significant - it takes around 21 seconds to execute on > > my machine, but only about 4 seconds without the START_FIX and END_FIX. > > Of course, this is not necessarily any sort of accurate test, but it > > does indicate that this can be an expensive operation. > > Yes, changing the floating-point modes is very slow on many (all?) x86 > processors. I believe it flushes the floating-point pipeline, which > takes many clock cycles. OK, how about this plan: (1) On systems with sse2, we do the option 3a (which is "If a processor supports sse2, then passing gcc -march=whatever -msse2 -mfpmath=sse (maybe the -march isn't needed) will cause gcc to use sse registers and instructions for doubles, and these have the proper precision.") (2) On systems without sse2 (old slow pentium 3's) we do the START_FIX and END_FIX. These computers are very slow anyways, so let them suffer (and the suffering is *only* for code that uses quaddouble, which is very little code anyways). William --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: Getting Geometry in SAGE
New members are moderated. I didn't post FlameGuy's email but didn't delete it either. (Both William and I are moderators. It was addressed to William and didn't think it appropriate for me to delete emails adressed to him, but I did tell him I would not post it for obvious reasons.) I do think it is important to maintain decorum and Iapologize. I'll try to improve communication with William so this doesn't happen again. On 8/14/07, Chris Chiasson <[EMAIL PROTECTED]> wrote: > > On Aug 14, 11:37 am, "William Stein" <[EMAIL PROTECTED]> wrote: > > On 8/14/07, John Cremona <[EMAIL PROTECTED]> wrote: > > > > > > > > > Whatever the fact of the case, I find the tone of this email > > > offensive. I'm sure I am not the only one. I think an apology is > > > called for. > > > > I'm very sorry for not carefully reading the email before moderating > > it "send" -- I usually just check that postings aren't blatant spam. > > In the future I'll be more careful when moderating messages. > > > > -- William > > > > This place is (pre-)moderated? > > > > > --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: New Axiom package
It build fine but the I/O is seriously bugged for me, at least for the displaying of rational functions (but probably wrapping too)... I paste this: -- k = var('k') tmp = [rising_factorial(x,k)/falling_factorial(x,k) for k in range(5)] pol1 = axiom(tmp[1]) pol1 pol1 pol1 pol4 = axiom(tmp[4]) pol4 pol4 pol4 pol4 pol4 pol4 pol4 pol4 pol1 = axiom(tmp[1]) pol1 pol1 pol1 pol1 pol1 pol1 pol1 pol1 -- and get -- -- | SAGE Version 2.8, Release Date: 2007-08-12 | | Type notebook() for the GUI, and license() for information.| -- Loading SAGE library. Current Mercurial branch is: paul sage: k = var('k') sage: tmp = [rising_factorial(x,k)/falling_factorial(x,k) for k in range(5)] sage: sage: pol1 = axiom(tmp[1]) sage: pol1 1 sage: pol1 1 sage: pol1 1 sage: pol4 = axiom(tmp[4]) sage: pol4 sage: pol4 sage: pol4 - (0)) * ((x) - (1))) * ((x) - (2))) * ((x) - (3))) 3 2 x + 6x + 11x + 6 -- 3 2 x - 6x + 11x - 6 sage: pol4 3 2 x + 6x + 11x + 6 -- 3 2 x - 6x + 11x - 6 sage: pol4 3 2 x + 6x + 11x + 6 -- 3 2 x - 6x + 11x - 6 sage: pol4 3 2 x + 6x + 11x + 6 -- 3 2 x - 6x + 11x - 6 sage: pol4 3 2 x + 6x + 11x + 6 -- 3 2 x - 6x + 11x - 6 sage: pol4 3 2 x + 6x + 11x + 6 -- 3 2 x - 6x + 11x - 6 sage: pol1 = axiom(tmp[1]) sage: pol1 3 2 x + 6x + 11x + 6 -- 3 2 x - 6x + 11x - 6 sage: pol1 3 2 x + 6x + 11x + 6 -- 3 2 x - 6x + 11x - 6 sage: pol1 3 2 x + 6x + 11x + 6 -- 3 2 x - 6x + 11x - 6 sage: pol1 1 sage: pol1 1 sage: pol1 1 sage: pol1 1 sage: pol1 1 sage: -- Notice how when I ask for "sage: pol4" I first get blanks a few times, then bits of the (presumably) internal representation, then what should be displayed. And that pol4 keeps on coming when I ask for "sage: pol1". That seems like a first serious bug. It also hangs when I try to do comvert SAGE rational functions to axiom rational functions without displaying them. I am running SAGE 2.8 on a MacBook under MacOSX Paul On Aug 14, 9:53 pm, "Bill Page" <[EMAIL PROTECTED]> wrote: > On 8/14/07, [EMAIL PROTECTED] wrote: > > > Worked for me too, MacOSX MacBook. > > (and I tried installing the publicly-available optional Axiom package > > yesterday and that didn't work) > > Paul > > The previous optional package 'axiom4sage-0.1.1' is compatible only > with versions of Sage prior to 2.8.1 and builds on a more restricted > range of machines using GCL rather than Clisp. > > The new 'axiom4sage-0.3.1' includes several important bug fixes as > well as the new GUESS package: > > http://arxiv.org/abs/math/0702086 > > Extended Rate, more GFUN > Authors: Martin Rubey > > Abstract: We present a software package that guesses formulae for > sequences of, for example, rational numbers or rational functions, > given the first few terms. Thereby we extend and complement Christian > Krattenthaler's program Rate and the relevant parts of Bruno Salvy and > Paul Zimmermann's GFUN. > > See also: > > http://wiki.axiom-developer.org/GuessingFormulasForSequences > > For example: > > sage: r=[0,3,32,375,5184,84035] > sage: R=axiom(r).guessExpRat(); R > >n > [[function= n (n + 2) ,order= 0]] > > sage: f=axiom('%s.1.function'%N.name()); f > >n > n (n + 2) > > sage: f.type() > Expression Integer > > sage: f.eval(axiom('n=5')) > 84035 > > sage: f.eval(axiom('n=6')) > 1572864 > > sage: > > Regards, > Bill Page. --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: New Axiom package
On 8/14/07, [EMAIL PROTECTED] wrote: > > Worked for me too, MacOSX MacBook. > (and I tried installing the publicly-available optional Axiom package > yesterday and that didn't work) > Paul > The previous optional package 'axiom4sage-0.1.1' is compatible only with versions of Sage prior to 2.8.1 and builds on a more restricted range of machines using GCL rather than Clisp. The new 'axiom4sage-0.3.1' includes several important bug fixes as well as the new GUESS package: http://arxiv.org/abs/math/0702086 Extended Rate, more GFUN Authors: Martin Rubey Abstract: We present a software package that guesses formulae for sequences of, for example, rational numbers or rational functions, given the first few terms. Thereby we extend and complement Christian Krattenthaler's program Rate and the relevant parts of Bruno Salvy and Paul Zimmermann's GFUN. See also: http://wiki.axiom-developer.org/GuessingFormulasForSequences For example: sage: r=[0,3,32,375,5184,84035] sage: R=axiom(r).guessExpRat(); R n [[function= n (n + 2) ,order= 0]] sage: f=axiom('%s.1.function'%N.name()); f n n (n + 2) sage: f.type() Expression Integer sage: f.eval(axiom('n=5')) 84035 sage: f.eval(axiom('n=6')) 1572864 sage: Regards, Bill Page. --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: ticket #59 (optimize elliptic curve arithmetic)
On Aug 14, 2007, at 4:19 PM, Alex Ghitza wrote: > Hi, > > I've looked at ticket #59, in which David Kohel says: > > > William, my student noticed some slow performance with elliptic > curves > group law. I think there was a huge overhead in duplication: > > sage: E = EllipticCurve([GF(101)(1),3]) > sage: P = E([-1,1,1]) > sage: timeit 2*P > 100 loops, best of 3: 3.81 ms per loop > sage: timeit P+P > 1000 loops, best of 3: 1.81 ms per loop > > Basically n*P was passing through all sorts of high-level layers for > group schemes, abelian groups, and the like. > > > It turns out that n*P is inherited from multiplication in > AdditiveGroupElement, which itself goes all the way to ModuleElement. > But at that point, it is implemented efficiently (by repeated > doubling). > It seems to me that the loss of efficiency in comparison with P+P > +...+P > is all in the overhead. > > I experimented with this and noticed that P+...+P is faster than > n*P as > long as n is smaller than 10, then the repeated doubling starts to > bear > fruit. So inside ell_point.py, I overloaded the multiplication > operator, simply making it return P+...+P if n is less than 10, and > return n*P if n is larger than 10 (see below). > > This is not a very elegant solution. For one thing, the magic > number 10 > that works on my architecture (Intel Core Duo running Gentoo) might > not > be right for a completely different machine. Is there a way to > automatically generate this number when the user compiles sage from > source (by starting with a value like 100, say, and then doing a > "binary > search" until it hits the optimal number)? Even on a single architecture, the correct crossover point will vary between different base fields. For example if the base field is Q, the best crossover will depend on the size of the coordinates. If the coordinates are enormous, it's probably best to do repeated doubling as much as possible. Similarly if you are working over GF(p) for some very large p, you'll find a different crossover from GF(101). This is going to be a difficult problem in general. I think a better approach is to try to reduce the overhead, i.e. to get the n*P pathway working faster, and/or to reduce the overhead in the addition law itself. david --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] ticket #59 (optimize elliptic curve arithmetic)
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi, I've looked at ticket #59, in which David Kohel says: William, my student noticed some slow performance with elliptic curves group law. I think there was a huge overhead in duplication: sage: E = EllipticCurve([GF(101)(1),3]) sage: P = E([-1,1,1]) sage: timeit 2*P 100 loops, best of 3: 3.81 ms per loop sage: timeit P+P 1000 loops, best of 3: 1.81 ms per loop Basically n*P was passing through all sorts of high-level layers for group schemes, abelian groups, and the like. It turns out that n*P is inherited from multiplication in AdditiveGroupElement, which itself goes all the way to ModuleElement. But at that point, it is implemented efficiently (by repeated doubling). It seems to me that the loss of efficiency in comparison with P+P+...+P is all in the overhead. I experimented with this and noticed that P+...+P is faster than n*P as long as n is smaller than 10, then the repeated doubling starts to bear fruit. So inside ell_point.py, I overloaded the multiplication operator, simply making it return P+...+P if n is less than 10, and return n*P if n is larger than 10 (see below). This is not a very elegant solution. For one thing, the magic number 10 that works on my architecture (Intel Core Duo running Gentoo) might not be right for a completely different machine. Is there a way to automatically generate this number when the user compiles sage from source (by starting with a value like 100, say, and then doing a "binary search" until it hits the optimal number)? Best, Alex this is the function I added to ell_point.py: +def _rmul_(self, n): +""" +Return n times the point, with respect to the group law on the +elliptic curve. + +INPUT: +self -- a point on an elliptic curve +n-- an integer + +OUTPUT: +the n-th multiple of the point + +EXAMPLES: +sage: E = EllipticCurve('389a') +sage: P = E([-1,1]) +sage: 3*P +(26/361 : -5720/6859 : 1) +sage: -2*P +(10/9 : 8/27 : 1) +""" +# for small values of n, return self+...+self +# for larger values, the generic multiplication is more efficient +if n<0: +result = (-n)*(-self) +elif n==0: +E = self.curve() +result = E(0) +elif n==1: +result = self +elif n==2: +result = self+self +elif n==3: +result = self+self+self +elif n==4: +result = self+self+self+self +elif n==5: +result = self+self+self+self+self +elif n==6: +result = self+self+self+self+self+self +elif n==7: +result = self+self+self+self+self+self+self +elif n==8: +result = self+self+self+self+self+self+self+self +elif n==9: +result = self+self+self+self+self+self+self+self+self +else: +result = AdditiveGroupElement._rmul_(self,n) +return result -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.7 (GNU/Linux) Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org iD8DBQFGwg5YdZTaNFFPILgRApFZAJ4voxanbGm9+DYnqRtVd9nHWLEUswCeNRJM AFdHE0B/Te/YWJXxtK0xI7I= =+Pn5 -END PGP SIGNATURE- --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: sage 2.8 build fails on G5 powerpc OS 10.4.10
On Aug 14, 2007, at 4:09 PM, Justin C. Walker wrote: >> Something goes wrong in LAPACK build: >> >> sage-spkg installed/lapack-20070723 2>&1 > [snip] >> make[2]: *** [lapack_install] Error 126 >> Error compiling lapack. > > A shot in the dark: can you check your install log for lines > containing "fork:"? I see plenty of weird failures on my G5 tower > that are caused by some kind of resource limit (I suspect its > 'maxprocperuid', but I don't know for sure). I've currently got it > bumped to 220 (which is the default value on the Core * Duo systems). Nope there is no "fork:". I have copied my whole install.log to http://sage.math.washington.edu/home/dmharvey/failed/install.log From memory I couldn't build 2.7.3 either, but not 100% sure. david --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: sage 2.8 build fails on G5 powerpc OS 10.4.10
On Aug 14, 2007, at 12:46 , David Harvey wrote: > > Something goes wrong in LAPACK build: > > sage-spkg installed/lapack-20070723 2>&1 [snip] > make[2]: *** [lapack_install] Error 126 > Error compiling lapack. A shot in the dark: can you check your install log for lines containing "fork:"? I see plenty of weird failures on my G5 tower that are caused by some kind of resource limit (I suspect its 'maxprocperuid', but I don't know for sure). I've currently got it bumped to 220 (which is the default value on the Core * Duo systems). Justin -- Justin C. Walker, Curmudgeon at Large Institute for the Absorption of Federal Funds --- If it weren't for carbon-14, I wouldn't date at all. --- --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: Getting Geometry in SAGE
On Aug 14, 11:37 am, "William Stein" <[EMAIL PROTECTED]> wrote: > On 8/14/07, John Cremona <[EMAIL PROTECTED]> wrote: > > > > > Whatever the fact of the case, I find the tone of this email > > offensive. I'm sure I am not the only one. I think an apology is > > called for. > > I'm very sorry for not carefully reading the email before moderating > it "send" -- I usually just check that postings aren't blatant spam. > In the future I'll be more careful when moderating messages. > > -- William > This place is (pre-)moderated? --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] sage 2.8 build fails on G5 powerpc OS 10.4.10
Something goes wrong in LAPACK build: sage-spkg installed/lapack-20070723 2>&1 lapack-20070723 Machine: Darwin George.local 8.10.0 Darwin Kernel Version 8.10.0: Wed May 23 16:50:59 PDT 2007; root:xnu-792.21.3~1/RELEASE_PPC Power Macintosh powerpc Deleting directories from past builds of previous/current versions of lapack-20070723 Extracting package /Users/david/sage-2.8/spkg/standard/ lapack-20070723.spkg ... -rw-r--r-- 1 david david 3442036 Jul 23 16:27 /Users/david/ sage-2.8/spkg/standard/lapack-20070723.spkg lapack-20070723/ lapack-20070723/.hg/ lapack-20070723/.hg/00changelog.i lapack-20070723/.hg/dirstate lapack-20070723/.hg/requires lapack-20070723/.hg/store/ . lapack-20070723/src/TESTING/zec.in lapack-20070723/src/TESTING/zed.in lapack-20070723/src/TESTING/zgbak.in lapack-20070723/src/TESTING/zgbal.in lapack-20070723/src/TESTING/zgd.in lapack-20070723/src/TESTING/zgg.in lapack-20070723/src/TESTING/zsb.in lapack-20070723/src/TESTING/zsg.in lapack-20070723/src/TESTING/ztest.in Finished extraction Host system uname -a: Darwin George.local 8.10.0 Darwin Kernel Version 8.10.0: Wed May 23 16:50:59 PDT 2007; root:xnu-792.21.3~1/RELEASE_PPC Power Macintosh powerpc GCC Version gcc -v Using built-in specs. Target: powerpc-apple-darwin8 Configured with: /private/var/tmp/gcc/gcc-5247.obj~4/src/configure -- disable-checking -enable-werror --prefix=/usr --mandir=/share/man -- enable-languages=c,objc,c++,obj-c++ --program-transform-name=/^[cg] [^.-]*$/s/$/-4.0/ --with-gxx-include-dir=/include/c++/4.0.0 -- build=powerpc-apple-darwin8 --host=powerpc-apple-darwin8 -- target=powerpc-apple-darwin8 Thread model: posix gcc version 4.0.1 (Apple Computer, Inc. build 5247) ( cd INSTALL; make; ./testlsame; ./testslamch; \ ./testdlamch; ./testsecond; ./testdsecnd; ./testversion ) sage_fortran -fPIC -c lsame.f -o lsame.o sage_fortran -fPIC -c lsametst.f -o lsametst.o sage_fortran -o testlsame lsame.o lsametst.o ld: table of contents for archive: /Users/david/sage-2.8/local/bin/../ lib/gcc-lib/powerpc-apple-darwin6.8/4.0.3//libf95.a is out of date; rerun ranlib(1) (can't load from it) ld: table of contents for archive: /Users/david/sage-2.8/local/bin/../ lib/gcc-lib/powerpc-apple-darwin6.8/4.0.3//libgcc.a is out of date; rerun ranlib(1) (can't load from it) ld: table of contents for archive: /Users/david/sage-2.8/local/bin/../ lib/gcc-lib/powerpc-apple-darwin6.8/4.0.3//libgcc_eh.a is out of date; rerun ranlib(1) (can't load from it) make[3]: *** [testlsame] Error 1 /bin/sh: line 1: ./testlsame: cannot execute binary file /bin/sh: line 1: ./testslamch: cannot execute binary file /bin/sh: line 1: ./testdlamch: cannot execute binary file /bin/sh: line 1: ./testsecond: No such file or directory /bin/sh: line 1: ./testdsecnd: cannot execute binary file /bin/sh: line 1: ./testversion: cannot execute binary file make[2]: *** [lapack_install] Error 126 Error compiling lapack. real0m0.419s user0m0.149s sys 0m0.163s sage: An error occured while installing lapack-20070723 Please email sage-devel http://groups.google.com/group/sage-devel explaining the problem and send the relevant part of of /Users/david/sage-2.8/install.log. Describe your computer, operating system, etc. If you want to try to fix the problem, yourself *don't* just cd to /Users/david/sage-2.8/spkg/build/lapack-20070723 and type 'make'. Instead (using bash) type "source local/bin/sage-env" from the directory /Users/david/sage-2.8 in order to set all environment variables correctly, then cd to /Users/david/sage-2.8/spkg/build/lapack-20070723 make[1]: *** [installed/lapack-20070723] Error 1 real60m26.672s user26m56.349s sys 17m18.560s david --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: [fricas-devel] New Axiom package
installs fine on suse 10.2 amd64 machine running sage 2.8. On 8/14/07, William Stein <[EMAIL PROTECTED]> wrote: > > Hi, > > Recently Waldek created a (forked) version of Axiom called "FriCAS", > which means "free CAS". He set things up so it builds and installs > in the SAGE environment in about 10-15 minutes (instead of hours > for the previous version of Axiom). I've created a SAGE package and > would like to request SAGE users to test out building this package on > OS X and Linux. To try the package just do this: > > sage -i axiom4sage-0.3.1.spkg > > After about 10-15 minutes you should then be able to do this: > > sage -axiom > > and see something like what is listed below. Thanks! > > fermat:~/sage-2.7.3 was$ ./sage -axiom > i i i i i i i ooooo o o > I I I I I I I 8 8 8 8 8 o 88 > I \ `+' / I 8 8 8 888 >\ `-+-' / 8 8 8 o 8 > `-__|__-'8 8 8 8 8 > |8 o 8 8 o 8 8 > --+-- o8oo ooo8ooo o 8 > > Copyright (c) Bruno Haible, Michael Stoll 1992, 1993 > Copyright (c) Bruno Haible, Marcus Daniels 1994-1997 > Copyright (c) Bruno Haible, Pierpaolo Bernardi, Sam Steingold 1998 > Copyright (c) Bruno Haible, Sam Steingold 1999-2000 > Copyright (c) Sam Steingold, Bruno Haible 2001-2006 > > FriCAS (AXIOM fork) Computer Algebra System > Version: FriCAS 2007-07-19 >Timestamp: Tuesday August 14, 2007 at 12:34:52 > - >Issue )copyright to view copyright notices. >Issue )summary for a summary of useful system commands. >Issue )quit to leave AXIOM and return to shell. > - > > (1) -> 2 + 3 > >(1) 5 > Type: PositiveInteger > (2) -> integrate(sin(x)*cos(x), x) > > > > > On 8/14/07, Gabriel Dos Reis <[EMAIL PROTECTED]> wrote: > > On Tue, 14 Aug 2007, Bill Page wrote: > > > > | > > | Waldek > > | > > | We haven't heard much from you on this list lately. :-( I presume you > > | have been busy with other things or perhaps on vacation? > > | > > | On 7/27/07, you wrote: > > | > ... > > | > I have put a tarball at: > > | > > > | > http://www.math.uni.wroc.pl/~hebisch/prog/fricas-25.nn4.tar.bz2 > > | > > > | > Configuration and dependencies are like in regular FriCAS build, > > | > the difference is that this tarball skips testsuite run (you > > | > can still run tests by hand typing "make all-input") and uses > > | > precompiled Lisp (and databases) for algebra build. > > | > > > | > Disclaimer: this is relatively quick hack, intended mostly as a > > | > proof of concept. > > | > > > | > > | Your "proof of concept" is about to be released to Sage users as a > > | package called 'axiom4sage-0.3.1'. If there are other things that you > > | think should be done on this build, please let's discuss it. (E.g. I > > | know you also experimented with a faster build that did not re-start > > | clisp for every algebra file. Does that work? Does it need more > > | testing?) > > | > > | Anyway, here is a very small patch needed to build this version of > > | FriCAS on MAC OSX: > > | > > | [EMAIL PROTECTED]:~/packages$ diff -au > > | axiom4sage-0.3/src/configure.ac.pamphlet > > | axiom4sage-0.3.1/src/configure.ac.pamphlet > > | --- axiom4sage-0.3/src/configure.ac.pamphlet2007-07-27 > > | 03:36:44.0 -0700 > > | +++ axiom4sage-0.3.1/src/configure.ac.pamphlet 2007-08-14 > > | 06:40:43.0 -0700 > > | @@ -530,7 +530,7 @@ > > | *solaris*) > > | AC_DEFINE([SUNplatform], [], [SunOS flavour]) > > | ;; > > | -powerpc*darwin*) > > | +*darwin*) > > | AC_DEFINE([MACOSXplatform], [], [MACOX flavour]) > > | CCF="-O2 -fno-strength-reduce -Wall -D_GNU_SOURCE \ > > | -I/usr/include -I/usr/include/sys" > > | [EMAIL PROTECTED]:~/packages$ > > | > > | --- > > | > > | Of course you also need to run 'build-setup.sh'. > > | > > | I think the change is self-explanatory: Not all MAC OSX systems > > | identify themselves as "powerpc...". I think '*darwin*' will be > > | sufficient, right? > > > > I believe that is OK. When I put that pattern in, I was being overly > > conservative at the time. > > > > -- Gaby > > > > > -- > William Stein > Associate Professor of Mathematics > University of Washington > http://www.williamstein.org > > > > --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs
[sage-devel] Re: quaddouble in sage, and other floating point issues
On Aug 14, 12:59 am, Jonathan Bober <[EMAIL PROTECTED]> wrote: > This is exactly what NTL does in its quad float class. Just about every > function starts and ends with a macro to adjust the fpu, resulting in > around 7 extra assembly instructions. In the following code, the > overhead is quite significant - it takes around 21 seconds to execute on > my machine, but only about 4 seconds without the START_FIX and END_FIX. > Of course, this is not necessarily any sort of accurate test, but it > does indicate that this can be an expensive operation. Yes, changing the floating-point modes is very slow on many (all?) x86 processors. I believe it flushes the floating-point pipeline, which takes many clock cycles. Carl --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: New Axiom package
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Installed successfully on Dell Latitude D620 (Core Duo T2400 1.83GHz) running Gentoo and sage-2.8. It's slow because the laptop was running on battery, and hence throttled down. real29m12.096s user26m1.034s sys 1m22.641s Successfully installed axiom4sage-0.3.1 Alex William Stein wrote: > Hi, > > Recently Waldek created a (forked) version of Axiom called "FriCAS", > which means "free CAS". He set things up so it builds and installs > in the SAGE environment in about 10-15 minutes (instead of hours > for the previous version of Axiom). I've created a SAGE package and > would like to request SAGE users to test out building this package on > OS X and Linux. To try the package just do this: > > sage -i axiom4sage-0.3.1.spkg > > After about 10-15 minutes you should then be able to do this: > > sage -axiom > > and see something like what is listed below. Thanks! > > fermat:~/sage-2.7.3 was$ ./sage -axiom > i i i i i i i ooooo o o > I I I I I I I 8 8 8 8 8 o 88 > I \ `+' / I 8 8 8 888 >\ `-+-' / 8 8 8 o 8 > `-__|__-'8 8 8 8 8 > |8 o 8 8 o 8 8 > --+-- o8oo ooo8ooo o 8 > > Copyright (c) Bruno Haible, Michael Stoll 1992, 1993 > Copyright (c) Bruno Haible, Marcus Daniels 1994-1997 > Copyright (c) Bruno Haible, Pierpaolo Bernardi, Sam Steingold 1998 > Copyright (c) Bruno Haible, Sam Steingold 1999-2000 > Copyright (c) Sam Steingold, Bruno Haible 2001-2006 > > FriCAS (AXIOM fork) Computer Algebra System > Version: FriCAS 2007-07-19 >Timestamp: Tuesday August 14, 2007 at 12:34:52 > - >Issue )copyright to view copyright notices. >Issue )summary for a summary of useful system commands. >Issue )quit to leave AXIOM and return to shell. > - > > (1) -> 2 + 3 > >(1) 5 > Type: PositiveInteger > (2) -> integrate(sin(x)*cos(x), x) > > > > -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.7 (GNU/Linux) Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org iD8DBQFGwehedZTaNFFPILgRAvbnAKCom1ymjtLZFKlFzWZHcSFbIY7uiQCfT+d+ 5nFtGHleE7G4N9fDcb3BuM4= =1dLp -END PGP SIGNATURE- --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: New Axiom package
Worked for me too, MacOSX MacBook. (and I tried installing the publicly-available optional Axiom package yesterday and that didn't work) Paul Axiom installation finished. real10m17.685s user8m5.767s sys 1m27.710s Successfully installed axiom4sage-0.3.1 Now cleaning up tmp files. Making SAGE/Python scripts relocatable... Making script relocatable paul-olivier-dehayes-computer:~ pdehaye$ On Aug 14, 6:15 pm, Jaap Spies <[EMAIL PROTECTED]> wrote: > William Stein wrote: > > for the previous version of Axiom). I've created a SAGE package and > > would like to request SAGE users to test out building this package on > > OS X and Linux. To try the package just do this: > > > sage -i axiom4sage-0.3.1.spkg > > > After about 10-15 minutes you should then be able to do this: > > 79 Axiom installation finished. > > real12m24.290s > user10m40.994s > sys 0m59.883s > Successfully installed axiom4sage-0.3.1 > Now cleaning up tmp files. > Making SAGE/Python scripts relocatable... > Making script relocatable > [EMAIL PROTECTED] sage-2.8]$ > > The installation could be more silent or if you prefer less verbose! > > Jaap --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: New Axiom package
William Stein wrote: > for the previous version of Axiom). I've created a SAGE package and > would like to request SAGE users to test out building this package on > OS X and Linux. To try the package just do this: > > sage -i axiom4sage-0.3.1.spkg > > After about 10-15 minutes you should then be able to do this: 79 Axiom installation finished. real12m24.290s user10m40.994s sys 0m59.883s Successfully installed axiom4sage-0.3.1 Now cleaning up tmp files. Making SAGE/Python scripts relocatable... Making script relocatable [EMAIL PROTECTED] sage-2.8]$ The installation could be more silent or if you prefer less verbose! Jaap --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: New Axiom package
Hello, The build worked perfectly with SAGE 2.7.3 on Ubuntu 7.04 on a Core 2 Duo. --Mike On 8/14/07, William Stein <[EMAIL PROTECTED]> wrote: > > Hi, > > Recently Waldek created a (forked) version of Axiom called "FriCAS", > which means "free CAS". He set things up so it builds and installs > in the SAGE environment in about 10-15 minutes (instead of hours > for the previous version of Axiom). I've created a SAGE package and > would like to request SAGE users to test out building this package on > OS X and Linux. To try the package just do this: > > sage -i axiom4sage-0.3.1.spkg > > After about 10-15 minutes you should then be able to do this: > > sage -axiom > > and see something like what is listed below. Thanks! > > fermat:~/sage-2.7.3 was$ ./sage -axiom > i i i i i i i ooooo o o > I I I I I I I 8 8 8 8 8 o 88 > I \ `+' / I 8 8 8 888 >\ `-+-' / 8 8 8 o 8 > `-__|__-'8 8 8 8 8 > |8 o 8 8 o 8 8 > --+-- o8oo ooo8ooo o 8 > > Copyright (c) Bruno Haible, Michael Stoll 1992, 1993 > Copyright (c) Bruno Haible, Marcus Daniels 1994-1997 > Copyright (c) Bruno Haible, Pierpaolo Bernardi, Sam Steingold 1998 > Copyright (c) Bruno Haible, Sam Steingold 1999-2000 > Copyright (c) Sam Steingold, Bruno Haible 2001-2006 > > FriCAS (AXIOM fork) Computer Algebra System > Version: FriCAS 2007-07-19 >Timestamp: Tuesday August 14, 2007 at 12:34:52 > - >Issue )copyright to view copyright notices. >Issue )summary for a summary of useful system commands. >Issue )quit to leave AXIOM and return to shell. > - > > (1) -> 2 + 3 > >(1) 5 > Type: PositiveInteger > (2) -> integrate(sin(x)*cos(x), x) > > > > > On 8/14/07, Gabriel Dos Reis <[EMAIL PROTECTED]> wrote: > > On Tue, 14 Aug 2007, Bill Page wrote: > > > > | > > | Waldek > > | > > | We haven't heard much from you on this list lately. :-( I presume you > > | have been busy with other things or perhaps on vacation? > > | > > | On 7/27/07, you wrote: > > | > ... > > | > I have put a tarball at: > > | > > > | > http://www.math.uni.wroc.pl/~hebisch/prog/fricas-25.nn4.tar.bz2 > > | > > > | > Configuration and dependencies are like in regular FriCAS build, > > | > the difference is that this tarball skips testsuite run (you > > | > can still run tests by hand typing "make all-input") and uses > > | > precompiled Lisp (and databases) for algebra build. > > | > > > | > Disclaimer: this is relatively quick hack, intended mostly as a > > | > proof of concept. > > | > > > | > > | Your "proof of concept" is about to be released to Sage users as a > > | package called 'axiom4sage-0.3.1'. If there are other things that you > > | think should be done on this build, please let's discuss it. (E.g. I > > | know you also experimented with a faster build that did not re-start > > | clisp for every algebra file. Does that work? Does it need more > > | testing?) > > | > > | Anyway, here is a very small patch needed to build this version of > > | FriCAS on MAC OSX: > > | > > | [EMAIL PROTECTED]:~/packages$ diff -au > > | axiom4sage-0.3/src/configure.ac.pamphlet > > | axiom4sage-0.3.1/src/configure.ac.pamphlet > > | --- axiom4sage-0.3/src/configure.ac.pamphlet2007-07-27 > > | 03:36:44.0 -0700 > > | +++ axiom4sage-0.3.1/src/configure.ac.pamphlet 2007-08-14 > > | 06:40:43.0 -0700 > > | @@ -530,7 +530,7 @@ > > | *solaris*) > > | AC_DEFINE([SUNplatform], [], [SunOS flavour]) > > | ;; > > | -powerpc*darwin*) > > | +*darwin*) > > | AC_DEFINE([MACOSXplatform], [], [MACOX flavour]) > > | CCF="-O2 -fno-strength-reduce -Wall -D_GNU_SOURCE \ > > | -I/usr/include -I/usr/include/sys" > > | [EMAIL PROTECTED]:~/packages$ > > | > > | --- > > | > > | Of course you also need to run 'build-setup.sh'. > > | > > | I think the change is self-explanatory: Not all MAC OSX systems > > | identify themselves as "powerpc...". I think '*darwin*' will be > > | sufficient, right? > > > > I believe that is OK. When I put that pattern in, I was being overly > > conservative at the time. > > > > -- Gaby > > > > > -- > William Stein > Associate Professor of Mathematics > University of Washington > http://www.williamstein.org > > > > --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this group at http://groups.go
[sage-devel] New Axiom package
Hi, Recently Waldek created a (forked) version of Axiom called "FriCAS", which means "free CAS". He set things up so it builds and installs in the SAGE environment in about 10-15 minutes (instead of hours for the previous version of Axiom). I've created a SAGE package and would like to request SAGE users to test out building this package on OS X and Linux. To try the package just do this: sage -i axiom4sage-0.3.1.spkg After about 10-15 minutes you should then be able to do this: sage -axiom and see something like what is listed below. Thanks! fermat:~/sage-2.7.3 was$ ./sage -axiom i i i i i i i ooooo o o I I I I I I I 8 8 8 8 8 o 88 I \ `+' / I 8 8 8 888 \ `-+-' / 8 8 8 o 8 `-__|__-'8 8 8 8 8 |8 o 8 8 o 8 8 --+-- o8oo ooo8ooo o 8 Copyright (c) Bruno Haible, Michael Stoll 1992, 1993 Copyright (c) Bruno Haible, Marcus Daniels 1994-1997 Copyright (c) Bruno Haible, Pierpaolo Bernardi, Sam Steingold 1998 Copyright (c) Bruno Haible, Sam Steingold 1999-2000 Copyright (c) Sam Steingold, Bruno Haible 2001-2006 FriCAS (AXIOM fork) Computer Algebra System Version: FriCAS 2007-07-19 Timestamp: Tuesday August 14, 2007 at 12:34:52 - Issue )copyright to view copyright notices. Issue )summary for a summary of useful system commands. Issue )quit to leave AXIOM and return to shell. - (1) -> 2 + 3 (1) 5 Type: PositiveInteger (2) -> integrate(sin(x)*cos(x), x) On 8/14/07, Gabriel Dos Reis <[EMAIL PROTECTED]> wrote: > On Tue, 14 Aug 2007, Bill Page wrote: > > | > | Waldek > | > | We haven't heard much from you on this list lately. :-( I presume you > | have been busy with other things or perhaps on vacation? > | > | On 7/27/07, you wrote: > | > ... > | > I have put a tarball at: > | > > | > http://www.math.uni.wroc.pl/~hebisch/prog/fricas-25.nn4.tar.bz2 > | > > | > Configuration and dependencies are like in regular FriCAS build, > | > the difference is that this tarball skips testsuite run (you > | > can still run tests by hand typing "make all-input") and uses > | > precompiled Lisp (and databases) for algebra build. > | > > | > Disclaimer: this is relatively quick hack, intended mostly as a > | > proof of concept. > | > > | > | Your "proof of concept" is about to be released to Sage users as a > | package called 'axiom4sage-0.3.1'. If there are other things that you > | think should be done on this build, please let's discuss it. (E.g. I > | know you also experimented with a faster build that did not re-start > | clisp for every algebra file. Does that work? Does it need more > | testing?) > | > | Anyway, here is a very small patch needed to build this version of > | FriCAS on MAC OSX: > | > | [EMAIL PROTECTED]:~/packages$ diff -au > | axiom4sage-0.3/src/configure.ac.pamphlet > | axiom4sage-0.3.1/src/configure.ac.pamphlet > | --- axiom4sage-0.3/src/configure.ac.pamphlet2007-07-27 > | 03:36:44.0 -0700 > | +++ axiom4sage-0.3.1/src/configure.ac.pamphlet 2007-08-14 > | 06:40:43.0 -0700 > | @@ -530,7 +530,7 @@ > | *solaris*) > | AC_DEFINE([SUNplatform], [], [SunOS flavour]) > | ;; > | -powerpc*darwin*) > | +*darwin*) > | AC_DEFINE([MACOSXplatform], [], [MACOX flavour]) > | CCF="-O2 -fno-strength-reduce -Wall -D_GNU_SOURCE \ > | -I/usr/include -I/usr/include/sys" > | [EMAIL PROTECTED]:~/packages$ > | > | --- > | > | Of course you also need to run 'build-setup.sh'. > | > | I think the change is self-explanatory: Not all MAC OSX systems > | identify themselves as "powerpc...". I think '*darwin*' will be > | sufficient, right? > > I believe that is OK. When I put that pattern in, I was being overly > conservative at the time. > > -- Gaby > -- William Stein Associate Professor of Mathematics University of Washington http://www.williamstein.org --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: Getting Geometry in SAGE
On 8/14/07, John Cremona <[EMAIL PROTECTED]> wrote: > > Whatever the fact of the case, I find the tone of this email > offensive. I'm sure I am not the only one. I think an apology is > called for. I'm very sorry for not carefully reading the email before moderating it "send" -- I usually just check that postings aren't blatant spam. In the future I'll be more careful when moderating messages. -- William > On 8/14/07, Arash Partow <[EMAIL PROTECTED]> wrote: > > > > Hi William, > > > > I believe you are the one that is confused about the GPL. > > The GPL is a form of copyright that is provided with a > > piece of work. In this case I permit people to use Wykobi > > under the terms of the GPL open source license as opposed > > to say the CPL, LGPL or the Mozilla license. > > > > As for commercial use, if you bothered to understand the > > FSF's definition of derivative works you would know it is > > a necessity for anyone wanting to use the code in a commercial > > environment not to have to deal with the viral effects of the > > GPL. > > > > I suggest in the future you actually "know" or at least become > > familiar with what you are talking about before you make such > > comments, as it does not bode well for someone that supposedly > > has tenure. > > > > > > Arash Partow > > __ > > Be one who knows what they don't know, > > Instead of being one who knows not what they don't know, > > Thinking they know everything about all things. > > http://www.partow.net > > > > > > > > > > On Aug 11, 9:47 am, "William Stein" <[EMAIL PROTECTED]> wrote: > > > > > > > Another GPL one is Wykobi > > > > (www.wykobi.com). > > > > > > Wykobi has no documentation, is very new (version 0.0.03), and > > > I'm pretty sure the author(s) are totally confused about what the > > > GPL is. The claim to make Wykobi available under the GPL but > > > only for non-commercial use -- which doesn't make sense because > > > "GPL + non-commercial" is non-GPL compatible license. The GPL > > > doesn't allow restrictions against commercial use. Anyway, Wykobi > > > unfortunately isn't really GPL'd. > > > > > > -- William > > > > > > > > > > > > > > > -- > John Cremona > > > > -- William Stein Associate Professor of Mathematics University of Washington http://www.williamstein.org --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: vpython package (boost unicode problem)
William Stein wrote: > I would definitely consider making this change, at least as long as nobody > sees any major problems with it in the next few days. One problem is > that doing "sage -upgrade" will require rebuilding every Python-related thing > i SAGE, which will be a bit painful. > > Anyway, comments? > Maybe not every Python related package. Definitely numpy and maybe some more using unicode, as is sage-2.8.spkg Jaap --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: Getting Geometry in SAGE
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi, About 30 seconds of google-ing "gpl commercial" yielded the following email written by Richard Stallman (see below). Of course, it is *theoretically* possible that he is also "confused about the GPL", but I think that's unlikely. Of course, I agree with the replies about the impropriety of the post's tone and wording, even if the statements made had been completely correct. Alex - - -Original Message- From: Richard Stallman [mailto:[EMAIL PROTECTED] Sent: Friday, January 14, 2000 8:41 PM To: [EMAIL PROTECTED] You will often hear people say that the GPL does not allow use of the code in commercial software. This is a subtle confusion. The GPL does not make any distinction between commercial and noncommercial software development. It allows businesses to do all the same things that individuals and schools are allowed to do. However, the GPL does forbid use of the code in *proprietary* (non-free) software. The GPL does not allow adding any additional restrictions to any program that includes the code. So any program that incorporates the GPL-covered code, whether it be commercial, academic, or avocational, must be released as free software. As a whole, it must be released under the GPL (though you can put a more lax GPL-compatible license such as X11 on the parts you write, if you wish). I think that people who say the GPL disallows commercial software are probably saying "commercial" but actually thinking "proprietary". But while they may understand the GPL rightly, lumping together commercial software and proprietary software is a grave mistake. There is plenty of noncommercial non-free software. There is also commercial free software--and we want to encourage more. If we want to encourage people to write commercial free software, first of all we should make sure they know it is a conceivable thing. So let's not use terminology which implies it is impossible! Note, the GNU Project and the Free Software Foundation are not associated with "open source". We are part of the Free Software Movement, not the Open Source Movement. For more explanation, see http://www.gnu.org/philosophy/free-software-for-freedom.html. - - Arash Partow wrote: > Hi William, > > I believe you are the one that is confused about the GPL. > The GPL is a form of copyright that is provided with a > piece of work. In this case I permit people to use Wykobi > under the terms of the GPL open source license as opposed > to say the CPL, LGPL or the Mozilla license. > > As for commercial use, if you bothered to understand the > FSF's definition of derivative works you would know it is > a necessity for anyone wanting to use the code in a commercial > environment not to have to deal with the viral effects of the > GPL. > > I suggest in the future you actually "know" or at least become > familiar with what you are talking about before you make such > comments, as it does not bode well for someone that supposedly > has tenure. > > > Arash Partow > __ > Be one who knows what they don't know, > Instead of being one who knows not what they don't know, > Thinking they know everything about all things. > http://www.partow.net > > > > > On Aug 11, 9:47 am, "William Stein" <[EMAIL PROTECTED]> wrote: >>> Another GPL one is Wykobi >>> (www.wykobi.com). >> Wykobi has no documentation, is very new (version 0.0.03), and >> I'm pretty sure the author(s) are totally confused about what the >> GPL is. The claim to make Wykobi available under the GPL but >> only for non-commercial use -- which doesn't make sense because >> "GPL + non-commercial" is non-GPL compatible license. The GPL >> doesn't allow restrictions against commercial use. Anyway, Wykobi >> unfortunately isn't really GPL'd. >> >> -- William >> > > > > > -BEGIN PGP SIGNATURE- Version: GnuPG v1.4.7 (GNU/Linux) Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org iD8DBQFGwdnBdZTaNFFPILgRAi0YAKCFMVmU27tJzTk8+fWRaIadWzr03wCfe5gl sWQUmUjL5F0c2BItKxk2Imw= =hKy7 -END PGP SIGNATURE- --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: file system slowdown on sage.math
On 8/14/07, Michael Abshoff <[EMAIL PROTECTED]> wrote: > The first thing I checked upon investigation was if somebody had pushed > the box into swap and that wasn't the case. > > An 'ls -1 > foo.txt' in the now infamous tmp-directory took in excess of > 10 seconds which is pretty bad and that was after an undetermined number > of files had been deleted. My script is now 1/10 of the way through deleted all the files in that directory :-). [EMAIL PROTECTED]:/home/rlmill/tmp$ tail -f nohup.out 1206640615.txt 0.107169889731 1206640758.txt 0.107170100043 1206640908.txt 0.107170310355 1206641547.txt 0.107170520667 --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] unicode UCS4
Hi, I changed the install script in python-2.5.1.p3, replaced: ./configure --prefix="$SAGE_LOCAL" with ./configure --prefix="$SAGE_LOCAL" --enable-unicode=ucs4 build sage-2.8 with no problems. -- All tests passed! Total time for all tests: 1670.5 seconds Hopefully I now can use precompiled libraries using unicode=ucs4 as e.g. libboost_python. Cheers, Jaap --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: vpython package (boost unicode problem)
On 8/14/07, Jaap Spies <[EMAIL PROTECTED]> wrote: > See also > > > http://mail.python.org/pipermail/distutils-sig/2006-July/006579.html > > for a discussion of Python and UCS2 vs. UCS4 > > Ubuntu and Fedora distribute Python compiled with UCS4. > Maybe the Python version of SAGE should comply with the distribution, > so for example libboost_python can be used. I would definitely consider making this change, at least as long as nobody sees any major problems with it in the next few days. One problem is that doing "sage -upgrade" will require rebuilding every Python-related thing i SAGE, which will be a bit painful. Anyway, comments? -- William --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: Getting Geometry in SAGE
On 8/14/07, Arash Partow <[EMAIL PROTECTED]> wrote: > I believe you are the one that is confused about the GPL. > The GPL is a form of copyright that is provided with a > piece of work. In this case I permit people to use Wykobi > under the terms of the GPL open source license as opposed > to say the CPL, LGPL or the Mozilla license. OK. > As for commercial use, if you bothered to understand the > FSF's definition of derivative works you would know it is > a necessity for anyone wanting to use the code in a commercial > environment not to have to deal with the viral effects of the > GPL. Redhat is a famous "commercial environment" that uses GPL'd code. The GPL license and "commercial" are not always incompatible. As another example, Star-P (http://www.interactivesupercomputing.com/products/) uses Octave as the compute backend for their product, and Octave is GPL'd. > I suggest in the future you actually "know" or at least become > familiar with what you are talking about before you make such > comments, as it does not bode well for someone that supposedly > has tenure. Thanks for your comments. However, this is not the correct level of civility for discussions on sage-devel. In fact, your lack of civility is another great argument against including Wykobi in SAGE. -- William --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: Getting Geometry in SAGE
Whatever the fact of the case, I find the tone of this email offensive. I'm sure I am not the only one. I think an apology is called for. John Cremona On 8/14/07, Arash Partow <[EMAIL PROTECTED]> wrote: > > Hi William, > > I believe you are the one that is confused about the GPL. > The GPL is a form of copyright that is provided with a > piece of work. In this case I permit people to use Wykobi > under the terms of the GPL open source license as opposed > to say the CPL, LGPL or the Mozilla license. > > As for commercial use, if you bothered to understand the > FSF's definition of derivative works you would know it is > a necessity for anyone wanting to use the code in a commercial > environment not to have to deal with the viral effects of the > GPL. > > I suggest in the future you actually "know" or at least become > familiar with what you are talking about before you make such > comments, as it does not bode well for someone that supposedly > has tenure. > > > Arash Partow > __ > Be one who knows what they don't know, > Instead of being one who knows not what they don't know, > Thinking they know everything about all things. > http://www.partow.net > > > > > On Aug 11, 9:47 am, "William Stein" <[EMAIL PROTECTED]> wrote: > > > > > Another GPL one is Wykobi > > > (www.wykobi.com). > > > > Wykobi has no documentation, is very new (version 0.0.03), and > > I'm pretty sure the author(s) are totally confused about what the > > GPL is. The claim to make Wykobi available under the GPL but > > only for non-commercial use -- which doesn't make sense because > > "GPL + non-commercial" is non-GPL compatible license. The GPL > > doesn't allow restrictions against commercial use. Anyway, Wykobi > > unfortunately isn't really GPL'd. > > > > -- William > > > > > > > -- John Cremona --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: Getting Geometry in SAGE
Hi William, I believe you are the one that is confused about the GPL. The GPL is a form of copyright that is provided with a piece of work. In this case I permit people to use Wykobi under the terms of the GPL open source license as opposed to say the CPL, LGPL, BSD or the Mozilla Public license. As for commercial use, if you bothered to understand the FSF's definition of derivative works you would know it is a necessity for anyone wanting to use the code in a commercial environment not to have to deal with the viral effects of the GPL. I suggest in the future you actually "know" or at least become familiar with what you are talking about before you make such comments, as it does not bode well for someone that supposedly has tenure. Arash Partow __ Be one who knows what they don't know, Instead of being one who knows not what they don't know, Thinking they know everything about all things. http://www.partow.net On Aug 11, 9:47 am, "William Stein" <[EMAIL PROTECTED]> wrote: > > > Another GPL one is Wykobi > > (www.wykobi.com). > > Wykobi has no documentation, is very new (version 0.0.03), and > I'm pretty sure the author(s) are totally confused about what the > GPL is. The claim to make Wykobi available under the GPL but > only for non-commercial use -- which doesn't make sense because > "GPL + non-commercial" is non-GPL compatible license. The GPL > doesn't allow restrictions against commercial use. Anyway, Wykobi > unfortunately isn't really GPL'd. > > -- William > --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: Getting Geometry in SAGE
Hi William, I believe you are the one that is confused about the GPL. The GPL is a form of copyright that is provided with a piece of work. In this case I permit people to use Wykobi under the terms of the GPL open source license as opposed to say the CPL, LGPL or the Mozilla license. As for commercial use, if you bothered to understand the FSF's definition of derivative works you would know it is a necessity for anyone wanting to use the code in a commercial environment not to have to deal with the viral effects of the GPL. I suggest in the future you actually "know" or at least become familiar with what you are talking about before you make such comments, as it does not bode well for someone that supposedly has tenure. Arash Partow __ Be one who knows what they don't know, Instead of being one who knows not what they don't know, Thinking they know everything about all things. http://www.partow.net On Aug 11, 9:47 am, "William Stein" <[EMAIL PROTECTED]> wrote: > > > Another GPL one is Wykobi > > (www.wykobi.com). > > Wykobi has no documentation, is very new (version 0.0.03), and > I'm pretty sure the author(s) are totally confused about what the > GPL is. The claim to make Wykobi available under the GPL but > only for non-commercial use -- which doesn't make sense because > "GPL + non-commercial" is non-GPL compatible license. The GPL > doesn't allow restrictions against commercial use. Anyway, Wykobi > unfortunately isn't really GPL'd. > > -- William > --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: build dependencies
On Tuesday 14 August 2007, Michael Abshoff wrote: > William Stein wrote: > > > > > Hello, > > > I'll wait for the answers to the above questions. I would be fine > > with including > > scons in SAGE. I think it's already in one of the optional packages > > (I can't remember > > for sure though). It's just a little python program, after all. > > > > It would take an absolutely massive amount of convincing to even begin > > to convince me to include the entire 17MB C++ boost library in SAGE. > > Well, PolyBori just about destroys any code out there for > GBasis-computations over boolean rings including Magma's F4 as well as > Faugere's F4 & F5. And it isn't only massively faster, but also needs > substancially less RAM. Just for the presentation they did in March or so, > I am not sure the MEGA paper had any benchmarks in them. Or ask Martin > about the CTC examples. It is definitely worth it, espcially if one could > cut down the amount of boost code needed to compile. I think we can benefit in two ways from PolyBori 1.) just use it for GB calculations. This way we can just use the existing Python-Boost wrapper, convert, calculate and convert back. In this scenario we would just have to investigate which pieces of boost we need. 2.) use PolyBori as the backend for the native SAGE datatype. So if a user constructs F_2[x_1,...,x_n]/FieldIdeal he/she gets QuotientRingElement_polybori just like he/she sometimes gets MPolynomial_libsingular now. If we are going to do that, we will need to wrap PolyBori again anyway to match our expected interface. However, in that case, we still want to allow/reimplement the native boost-power interface because some highlevel code of PolyBoRi is written in Python. If we can _replicate_ what is done using Boost with Cython -- I know we can, the question is how long it takes -- the highlevel-code would still work. So in this scenario we'll have to wrap twice. I strongly favor the second option. Martin -- name: Martin Albrecht _pgp: http://pgp.mit.edu:11371/pks/lookup?op=get&search=0x8EF0DC99 _www: http://www.informatik.uni-bremen.de/~malb _jab: [EMAIL PROTECTED] --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: file system slowdown on sage.math
Bill Hart wrote: > Hello, > To possibly add to the sage.math woes yesterday, I was testing a new > memory manager just before the slowdown. It had the unusual property > of allocating very many gigabytes of swap space in small chunks and > very little actual memory. My presumption was that the slowness had > something to do with sage.math trying to recover from the swap space > mess it created for itself after I killed the process that was causing > the problem. > Nope, when I checked there was less than 500MBs of swap used while there were more than 8GB of physical RAM available (even accounting for buffers). You also cannot allocate swap space directly because the allocs are buffered in system ram. The same thing holds for files mmaped into system RAM. The first thing I checked upon investigation was if somebody had pushed the box into swap and that wasn't the case. An 'ls -1 > foo.txt' in the now infamous tmp-directory took in excess of 10 seconds which is pretty bad and that was after an undetermined number of files had been deleted. > Needless to say, the new new memory manager doesn't have this > "feature". > One would hope so :) > Sorry if my tinkering contributed to the problem. > > Bill. > Cheers, Michael --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: vpython package (boost unicode problem)
Jaap Spies wrote: > I'm trying to make a vpython package for SAGE: > see http://www.vpython.org/index.html [...] Got an answer from the VPython mailing list: > Jaap Spies wrote: >> Bruce Sherwood wrote: >>> Dunno whether this might explain the problem that was seen with the Sage > [] >>> It seems that the option,"unicode=ucs4" ,is not used in compiling >>> "Python-2.5.1". > [] >> Is this option essential for vpython??? > > This whole problem has nothing to do with vpython, or only insofar as it uses > libboost_python for building. It is a libboost_python problem. > Libboost_python uses (by default) unicode, and therefore it needs to be > linked together with the same flavor of python, UCS2 or UCS4, that it was > built with. > > From the original post it seems clear that there were two pythons present on > the machine, the one from SAGE which believes in UCS2, and the one from the > linux distribution which believes in UCS4. The libboost_python was built > against the UCS4 python, so it cannot be used with UCS2 python. In other > words, if you want to build vpython with the python from SAGE, you first have > to build libboost_python with the same python. > > After reading up on this story, I still find it incredible that the python > developers consciously decided to allow python to adhere to one of two > incompatible persuasions, UCS2 and UCS4, so that compiled modules > (extensions) must adhere to the same persuasion as the python they were built > with, or they will crash. What is worse is that you don't see from the > outside which persuasion your python was built with, unless you run something > like > > python -c 'import sys; print sys.maxunicode' > > or you search for UCS2 or UCS4 in symbol names in the binary executable. > > AFAICT, UCS2 is still the default, i.e. you get this if you don't use any > --unicode configure flag. > > -- > Martin > > > > See also > http://mail.python.org/pipermail/distutils-sig/2006-July/006579.html for a discussion of Python and UCS2 vs. UCS4 Ubuntu and Fedora distribute Python compiled with UCS4. Maybe the Python version of SAGE should comply with the distribution, so for example libboost_python can be used. Jaap --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: file system slowdown on sage.math
To possibly add to the sage.math woes yesterday, I was testing a new memory manager just before the slowdown. It had the unusual property of allocating very many gigabytes of swap space in small chunks and very little actual memory. My presumption was that the slowness had something to do with sage.math trying to recover from the swap space mess it created for itself after I killed the process that was causing the problem. Needless to say, the new new memory manager doesn't have this "feature". Sorry if my tinkering contributed to the problem. Bill. On 14 Aug, 06:38, "Justin C. Walker" <[EMAIL PROTECTED]> wrote: > On Aug 13, 2007, at 18:51 , William Stein wrote: > > > > > > > On 8/13/07, mabshoff <[EMAIL PROTECTED] > > dortmund.de> wrote: > [snip] > > Robert Miller was cleaning up his home directories (we ran out > > of disk space yesterday) and ... > >8598 pts/16 DN+7:47 rm -rf tmp > > The problem is that the directory tmp contains almost 5 million files: > > The file system ain't made that can deal gracefully with that many > files in one directory. > > > sage:/home/rlmill/tmp# ls -1 > a > > sage:/home/rlmill/tmp# wc -l a > > 4754843 a > > How long did that 'ls' take? 'ls' has to read the whole directory, > sort it, and then print it :-}. Next time, use "ls -f1", which won't > sort the directory; it just takes the entries as they come in the read. > > > (Robert -- Why do you have so many tmp files? Maybe you should > > be using a database? Filesystems like ext3 aren't good at dealing > > with 5 million files in a directory...) > > Inquiring minds want to know... > > Justin > > -- > Justin C. Walker, Curmudgeon-at-Large > () The ASCII Ribbon Campaign > /\ Help Cure HTML Email --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] is_zero()
Hi, I have some problems with the way is_zero() is implemented (or maybe it is the documentation...) Is this behavior really desirable? sage: k = var('k') sage: pol = 1/(k-1) - 1/k -1/k/(k-1) sage: pol -1/((k - 1)*k) - (1/k) + 1/(k - 1) sage: pol.partial_fraction() 0 sage: pol.is_zero() False sage: pol.is_zero?? Type: builtin_function_or_method Base Class: String Form: Namespace: Interactive Source: def is_zero(self): """ Return True if self equals self.parent()(0). The default implementation is to fall back to 'not self.__nonzero__'. NOTE: Do not re-implement this method in your subclass but implement __nonzero__ instead. """ return not self sage: pol == pol.parent()(0) -1/((k - 1)*k) - (1/k) + 1/(k - 1) == 0 Paul --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: quaddouble in sage, and other floating point issues
On Tue, 2007-08-14 at 01:07 -0400, didier deshommes wrote: > 2007/8/14, William Stein <[EMAIL PROTECTED]>: > > Correct me if I'm wrong, but I don't think that makes sense. > > If I do this: > > > > sage: a = RQDF(5) > > sage: number_of_partitions(1000) > > 24061467864032622473692149727991 > > sage: del a > > > > then during the number_of_partitions call the CPU is set to > > the wrong mode. > > You are right, I did not think about this. > > 2007/8/13, Jonathan Bober <[EMAIL PROTECTED]>: > > So -ffloat-store doesn't work. So on cpus without sse2, the quaddouble > > wrapper needs to use the fpu fix, but probably should be rewritten so > > that it doesn't affect other things. > > > One possibility is that one wrap _add_c_impl, _sub_c_impl, > > etc., by setting the CPU mode, doing the op, then unsetting > > the cpu mode. This is probably the robust solution (??), > > but might have potential efficiency issues (which should be > > investigated with careful benchmarking). > > It's looking like the *only* solution. Looks fail-safe, and of course, > very ugly :). > This is exactly what NTL does in its quad float class. Just about every function starts and ends with a macro to adjust the fpu, resulting in around 7 extra assembly instructions. In the following code, the overhead is quite significant - it takes around 21 seconds to execute on my machine, but only about 4 seconds without the START_FIX and END_FIX. Of course, this is not necessarily any sort of accurate test, but it does indicate that this can be an expensive operation. In any case, the fpu fix should be probably be implemented as a macro/inline function in a header file (which would do nothing if quad double was compiled with sse2 support) rather than as a full blown function call, like it is now. #define START_FIX \ volatile unsigned short __old_cw, __new_cw; \ asm volatile ("fnstcw %0":"=m" (__old_cw)); \ __new_cw = (__old_cw & ~0x300) | 0x200; \ asm volatile ("fldcw %0": :"m" (__new_cw)); #define END_FIX asm volatile ("fldcw %0": :"m" (__old_cw)); int main() { int i,j; double x; double y = 1.2; double z = 1.3; for(i = 0; i < 1000; i++) { for(j = 0; j < 100; j++) { START_FIX x = y * z; END_FIX } } return 0; } > On a related note, it looks like GSL has a "-gsl" compilation flag > that "Use[s] GSL to control floating point rounding and precision". > > didier > > > > > --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---
[sage-devel] Re: sage-2.8
Hi On Mon, Aug 13, 2007 at 05:46:36AM -0700, William Stein wrote: > (3) Maybe (I haven't tested this) do >sage: notebook(reset=True) This works fine, thanks. It prompts for a new password and then continues. regards, Jan -- .~. /V\ Jan Groenewald /( )\www.aims.ac.za ^^-^^ --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@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-devel URLs: http://sage.scipy.org/sage/ and http://modular.math.washington.edu/sage/ -~--~~~~--~~--~--~---