[sage-devel] sage-2.8.7.alpha0

2007-10-13 Thread William Stein

Hi,

I've posted the first alpha release of sage-2.8.7 here:

   http://sage.math.washington.edu/tmp/

It is probably somewhat broken.  Building it, contributing
fixes, etc., would be greatly appreciated.

William

-- 
William Stein
Associate Professor of Mathematics
University of Washington
http://wstein.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: multivariate polynomial coefficients

2007-10-13 Thread Joel B. Mohler

Thanks for the discussion about this topic.  I send this mail to re-iterate 
and summarize.  It seems there are two things that you might want:
1)  Get the coefficient of a specific monomial in the multivariate polynomial 
ring.
2)  Get the coefficient of the polynomial in a tower of (two) polynomial 
extensions

I suggest that MPolynomial.coefficient method perform function #1 above.  It 
currently does some strange things over the base rings I've tested.  It does 
vary on base ring due to implementation differences.

I suggest that there be a new method which performs function #2 and I suggest 
it be called "polynomial_coefficient", but I really don't like that name.  I 
almost prefer "coefficient_polynomial" since it will sort beside coefficient 
in documentation and help people realize that if they are not happy 
with "coefficient" that they have an alternative.  Of course, each method's 
documentation should probably mention the other.

I'm wondering if we could have a vote on preferred syntax.  I'm not going to 
describe the parameters because if they are not clear enough from context, it 
probably isn't a good parameter choice :) :

sage: P.=ZZ[]
sage: f=(1-v)*(1-2*w)*(1-3*x); f
-6*v*w*x + 2*v*w + 3*v*x + 6*w*x - v - 2*w - 3*x + 1
##
#  Alternative Number 1
##
sage: f.polynomial_coefficient({w:0,v:1})
3*x - 1
##
#  Alternative Number 2
##
sage: f.polynomial_coefficient([x,y,z],v)
3*x - 1
##
#  Alternative Number 3
#  Overloading of coefficient to handle both functionalities
##
sage: f.coefficient(v,base=ZZ['x,y,z'])
3*x - 1
sage: f.coefficient(v)  # here's an example of functionality #1 above
-1

Those are listed in the order of my preference for the interface.  It feels to 
me that "Alternative Number 3" would not be much fun to implement (and 
involve coercions that some view as suspect).  I basically 
implemented "Alternative Number 1" as a function in my own program and it 
feels like a straightforward implementation just iterating through the 
polydict.  I'd love to have some other voters.

NOTE:  It is NOT sufficient to simply provide a monomial expression to 
the "polynomial_coefficient" method because it doesn't provide for the case 
when you explicitly want an exponent to be 0 -- this point seems to me to 
have caused a great deal of confusion in the current implementation.

--
Joel

--~--~-~--~~~---~--~~
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] Braid Goups / Coxeter Groups / Root Systems / ...

2007-10-13 Thread Mike Hansen

I just saw Robert's ticket about implementing / wrapping braid groups
in SAGE ( http://www.sagetrac.org/sage_trac/ticket/876 ) which is
pretty coincidental since I was thinking about the exact same thing
before I went to bed last night.  I wrote an email to Nikos
Apostolakis who had expressed interest in this on sage-support in June
to see if he was still interested.

GAP3 had support for braid groups provided by the "share package"
CHEVIE.  Unfortunately, when the change to GAP4 was made, this was
never ported over.  As David Joyner mentioned,  one possibility is to
have a (optional) GAP3 package from which the functionality could be
wrapped.  Another possibility is to port the GAP3 chevie code over to
GAP4.  I do not know how much of an undertaking this would be because
my familiarity with GAP programming is almost non-existent.  The third
possibility would be to port chevie directly over to SAGE; I looked at
the code for chevie, and this is definitely doable.  It would however
take a bit of time -- it's about 6000 lines of GAP code and about 2000
lines of comments.

It would be good to have a framework for root systems and the like
that all of these things can fit into.  This may be helpful in that
aspect: http://mupad-combinat.sourceforge.net/Wiki/RootSystems

--Mike

--~--~-~--~~~---~--~~
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.7.alpha0

2007-10-13 Thread mabshoff



On Oct 13, 10:10 am, "William Stein" <[EMAIL PROTECTED]> wrote:
> Hi,
>
> I've posted the first alpha release of sage-2.8.7 here:
>
>http://sage.math.washington.edu/tmp/
>
> It is probably somewhat broken.  Building it, contributing
> fixes, etc., would be greatly appreciated.

Hello,

building it went without a hitch, but there is trouble with the
doctests:


sage -t  const.tex
sage -t  devel/sage-main/sage/matrix/matrix_integer_dense.pyx
sage -t  devel/sage-main/sage/schemes/elliptic_curves/
lseries_ell.py
sage -t  devel/sage-main/sage/schemes/elliptic_curves/
ec_database.py
sage -t  devel/sage-main/sage/graphs/graph.py
sage -t  devel/sage-main/sage/sets/set.py
sage -t  devel/sage-main/sage/dsage/tests/testdoc.py
sage -t  devel/sage-main/sage/combinat/sloane_functions.py
sage -t  devel/sage-main/sage/rings/morphism.pyx
sage -t  devel/sage-main/sage/rings/polynomial/real_roots.pyx
sage -t  devel/sage-main/sage/rings/integer_mod.pyx
sage -t  devel/sage-main/sage/rings/residue_field.pyx


There are issues with coercion and also some fallout from the l-series
split up as well as some small ntl issues with the occasional
precision issue thrown in. If you fix any of the following please open
a trac ticket and attach a patch. Target the ticket for 2.8.7.

Cheers,

Michael

In Detail:

sage -t  const.tex
**
File "const.py", line 749:
sage: vals = E.Lseries_values_along_line(1-I, 1+10*I, 100) #
critical line
Exception raised:
Traceback (most recent call last):
  File "/tmp/Work-mabshoff/sage-2.8.7-alpha0/local/lib/python2.5/
doctest.py", line 1212, in __run
compileflags, 1) in test.globs
  File "", line 1, in 
vals = E.Lseries_values_along_line(Integer(1)-I,
Integer(1)+Integer(10)*I, Integer(100)) # critical line###line 749:
sage: vals = E.Lseries_values_along_line(1-I, 1+10*I, 100) #
critical line
AttributeError: 'EllipticCurve_rational_field' object has no
attribute 'Lseries_values_along_line'
**
File "const.py", line 750:
sage: L = [(z[1].real(), z[1].imag()) for z in vals]
Exception raised:
Traceback (most recent call last):
  File "/tmp/Work-mabshoff/sage-2.8.7-alpha0/local/lib/python2.5/
doctest.py", line 1212, in __run
compileflags, 1) in test.globs
  File "", line 1, in 
L = [(z[Integer(1)].real(), z[Integer(1)].imag()) for z in
vals]###line 750:
sage: L = [(z[1].real(), z[1].imag()) for z in vals]
NameError: name 'vals' is not defined
**
File "const.py", line 751:
sage: p = line(L, rgbcolor=(3/4,1/2,5/8))
Exception raised:
Traceback (most recent call last):
  File "/tmp/Work-mabshoff/sage-2.8.7-alpha0/local/lib/python2.5/
doctest.py", line 1212, in __run
compileflags, 1) in test.globs
  File "", line 1, in 
p = line(L, rgbcolor=(Integer(3)/Integer(4),Integer(1)/
Integer(2),Integer(5)/Integer(8)))###line 751:
sage: p = line(L, rgbcolor=(3/4,1/2,5/8))
NameError: name 'L' is not defined
**

sage -t  devel/sage-main/sage/matrix/matrix_integer_dense.pyx
**
File "matrix_integer_dense.pyx", line 1520:
sage: A.lll()
Exception raised:
Traceback (most recent call last):
  File "/tmp/Work-mabshoff/sage-2.8.7-alpha0/local/lib/python2.5/
doctest.py", line 1212, in __run
compileflags, 1) in test.globs
  File "", line 1, in 
A.lll()###line 1520:
sage: A.lll()
  File "matrix_integer_dense.pyx", line 1550, in
matrix_integer_dense.Matrix_integer_dense.lll
det = ZZ(det2.sqrt_approx())
  File "integer_ring.pyx", line 242, in
integer_ring.IntegerRing_class.__call__
  File "integer.pyx", line 319, in integer.Integer.__init__
  File "real_mpfr.pyx", line 1001, in
real_mpfr.RealNumber._integer_
ValueError: Attempt to coerce non-integral RealNumber to Integer
**

sage -t  devel/sage-main/sage/schemes/elliptic_curves/lseries_ell.py
**
File "lseries_ell.py", line 59:
sage: L.taylor_series(series_prec=3)
Expected:
-1.28158145691931e-23 + (7.26268290635587e-24)*z +
0.759316500288427*z^2 + O(z^3)
Got:
-2.69129566562797e-23 + (1.52514901968783e-23)*z +
0.759316500288427*z^2 + O(z^3)
**


sage -t  devel/sage-main/sage/schemes/elliptic_curves/ec_database.py
**
File "ec_database.py", line 21:
sage: elliptic_curves.rank(n=5, rank=3, tors=2, 

[sage-devel] Re: sage-2.8.7.alpha0

2007-10-13 Thread William Stein

Hi,

I just landed in Albuquerque, NM (free wifi, excellent power
outlets, etc., at their airport!).

Anyway, I get similar doctest failures on most architecture.
As Michael says, if anybody wants to help on this, report the
failure as a trac ticket, assign it to sage-2.8.7, and fix it.
Then hopefully I can just apply a bunch of patches and be ready
to go :-).  Please don't add new unrelated patches to 2.8.7.

By the way, the failures on PPC os x are an exception --
though the build works, there are way *more* failures there, e.g.,

**
File "set.py", line 1196:
sage: sqrt(RR(2)) in X # since sqrt(RR(2)) is currently approximated
Exception raised:
Traceback (most recent call last):
  File "/Users/was/sage-2.8.7-alpha0/local/lib/python2.5/doctest.py",
line 1212, in __run
compileflags, 1) in test.globs
  File "", line 1, in 
sqrt(RR(Integer(2))) in X # since sqrt(RR(2)) is currently
approximated###line 1196:
sage: sqrt(RR(2)) in X # since sqrt(RR(2)) is currently approximated
  File 
"/Users/was/sage-2.8.7-alpha0/local/lib/python2.5/site-packages/sage/sets/set.py",
line 1205, in __contains__
return (x in self.__X and x not in self.__Y) \
  File 
"/Users/was/sage-2.8.7-alpha0/local/lib/python2.5/site-packages/sage/sets/primes.py",
line 56, in __contains__
if not x in ZZ:
  File "parent.pyx", line 121, in parent.Parent.__contains__
  File "integer_ring.pyx", line 242, in
integer_ring.IntegerRing_class.__call__
  File "integer.pyx", line 319, in integer.Integer.__init__
  File "real_mpfr.pyx", line 1001, in real_mpfr.RealNumber._integer_
ValueError: Attempt to coerce non-integral RealNumber to Integer
**


On 10/13/07, mabshoff <[EMAIL PROTECTED]> wrote:
>
>
>
> On Oct 13, 10:10 am, "William Stein" <[EMAIL PROTECTED]> wrote:
> > Hi,
> >
> > I've posted the first alpha release of sage-2.8.7 here:
> >
> >http://sage.math.washington.edu/tmp/
> >
> > It is probably somewhat broken.  Building it, contributing
> > fixes, etc., would be greatly appreciated.
>
> Hello,
>
> building it went without a hitch, but there is trouble with the
> doctests:
>
>
> sage -t  const.tex
> sage -t  devel/sage-main/sage/matrix/matrix_integer_dense.pyx
> sage -t  devel/sage-main/sage/schemes/elliptic_curves/
> lseries_ell.py
> sage -t  devel/sage-main/sage/schemes/elliptic_curves/
> ec_database.py
> sage -t  devel/sage-main/sage/graphs/graph.py
> sage -t  devel/sage-main/sage/sets/set.py
> sage -t  devel/sage-main/sage/dsage/tests/testdoc.py
> sage -t  devel/sage-main/sage/combinat/sloane_functions.py
> sage -t  devel/sage-main/sage/rings/morphism.pyx
> sage -t  devel/sage-main/sage/rings/polynomial/real_roots.pyx
> sage -t  devel/sage-main/sage/rings/integer_mod.pyx
> sage -t  devel/sage-main/sage/rings/residue_field.pyx
>
>
> There are issues with coercion and also some fallout from the l-series
> split up as well as some small ntl issues with the occasional
> precision issue thrown in. If you fix any of the following please open
> a trac ticket and attach a patch. Target the ticket for 2.8.7.
>
> Cheers,
>
> Michael
>
> In Detail:
>
> sage -t  const.tex
> **
> File "const.py", line 749:
> sage: vals = E.Lseries_values_along_line(1-I, 1+10*I, 100) #
> critical line
> Exception raised:
> Traceback (most recent call last):
>   File "/tmp/Work-mabshoff/sage-2.8.7-alpha0/local/lib/python2.5/
> doctest.py", line 1212, in __run
> compileflags, 1) in test.globs
>   File "", line 1, in 
> vals = E.Lseries_values_along_line(Integer(1)-I,
> Integer(1)+Integer(10)*I, Integer(100)) # critical line###line 749:
> sage: vals = E.Lseries_values_along_line(1-I, 1+10*I, 100) #
> critical line
> AttributeError: 'EllipticCurve_rational_field' object has no
> attribute 'Lseries_values_along_line'
> **
> File "const.py", line 750:
> sage: L = [(z[1].real(), z[1].imag()) for z in vals]
> Exception raised:
> Traceback (most recent call last):
>   File "/tmp/Work-mabshoff/sage-2.8.7-alpha0/local/lib/python2.5/
> doctest.py", line 1212, in __run
> compileflags, 1) in test.globs
>   File "", line 1, in 
> L = [(z[Integer(1)].real(), z[Integer(1)].imag()) for z in
> vals]###line 750:
> sage: L = [(z[1].real(), z[1].imag()) for z in vals]
> NameError: name 'vals' is not defined
> **
> File "const.py", line 751:
> sage: p = line(L, rgbcolor=(3/4,1/2,5/8))
> Exception raised:
> Traceback (most recent call last):
>   File "/tmp/Work-mabshoff/sage-2.8.7-alpha0/local/

[sage-devel] my plenary ams talk

2007-10-13 Thread William Stein

Hi,

Here is a link to my slides and demo for my ams plenary talk:

  http://sagemath.org/why/alb/

-- 
William Stein
Associate Professor of Mathematics
University of Washington
http://wstein.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: sage-2.8.7.alpha0

2007-10-13 Thread Jaap Spies

William Stein wrote:

> 
> It is probably somewhat broken.  Building it, contributing
> fixes, etc., would be greatly appreciated.

A lot of test failures on Fedora 7:

The following tests failed:


 sage -t  devel/sage-main/sage/calculus/wester.py
 sage -t  devel/sage-main/sage/libs/pari/gen.pyx
 sage -t  devel/sage-main/sage/matrix/matrix2.pyx
 sage -t  devel/sage-main/sage/matrix/matrix_integer_dense.pyx
 sage -t  
devel/sage-main/sage/schemes/hyperelliptic_curves/jacobian_generic.py
 sage -t  
devel/sage-main/sage/schemes/elliptic_curves/ell_number_field.py
 sage -t  devel/sage-main/sage/schemes/elliptic_curves/sha.py
 sage -t  
devel/sage-main/sage/schemes/elliptic_curves/ell_rational_field.py
 sage -t  devel/sage-main/sage/schemes/elliptic_curves/ec_database.py
 sage -t  devel/sage-main/sage/schemes/generic/scheme.py
 sage -t  devel/sage-main/sage/modular/modsym/space.py
 sage -t  devel/sage-main/sage/modular/modsym/modsym.py
 sage -t  devel/sage-main/sage/modular/dirichlet.py
 sage -t  devel/sage-main/sage/rings/morphism.pyx
 sage -t  devel/sage-main/sage/rings/number_field/number_field_ideal.py
 sage -t  
devel/sage-main/sage/rings/number_field/number_field_element_quadratic.pyx
 sage -t  
devel/sage-main/sage/rings/number_field/number_field_element.pyx
 sage -t  devel/sage-main/sage/rings/number_field/maps.py
 sage -t  devel/sage-main/sage/rings/number_field/class_group.py
 sage -t  
devel/sage-main/sage/rings/number_field/number_field_ideal_rel.py
 sage -t  devel/sage-main/sage/rings/number_field/number_field_base.pyx
 sage -t  devel/sage-main/sage/rings/number_field/order.py
 sage -t  devel/sage-main/sage/rings/number_field/number_field.py
 sage -t  devel/sage-main/sage/rings/number_field/morphism.py
 sage -t  devel/sage-main/sage/rings/power_series_ring.py
 sage -t  devel/sage-main/sage/rings/power_series_ring_element.pyx
 sage -t  devel/sage-main/sage/rings/residue_field.pyx
 sage -t  devel/sage-main/sage/rings/rational_field.py
 sage -t  devel/sage-main/sage/rings/algebraic_real.py
 sage -t  devel/sage-main/sage/rings/integer_mod.pyx
 sage -t  
devel/sage-main/sage/rings/polynomial/polynomial_quotient_ring.py
 sage -t  devel/sage-main/sage/rings/polynomial/polynomial_element.pyx
 sage -t  
devel/sage-main/sage/rings/polynomial/polynomial_singular_interface.py
 sage -t  
devel/sage-main/sage/rings/polynomial/multi_polynomial_ring_generic.pyx
 sage -t  devel/sage-main/sage/rings/polynomial/real_roots.pyx
 sage -t  
devel/sage-main/sage/rings/polynomial/polynomial_element_generic.py
 sage -t  devel/sage-main/sage/sets/set.py
 sage -t  devel/sage-main/sage/modules/free_module.py
 sage -t  devel/sage-main/sage/combinat/sloane_functions.py
 sage -t  devel/sage-main/sage/categories/category_types.py
Total time for all tests: 1755.2 seconds
Please see /home/jaap/downloads/sage-2.8.7-alpha0/tmp/test.log for the complete 
log from this test.
[EMAIL PROTECTED] sage-2.8.7-alpha0]$

I can send details if you want.

Jaap

In sloane_functions.py sequence can now be indexed by reals!?
a(1.0) returns a value. Formerly an error was raised.


--~--~-~--~~~---~--~~
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: my plenary ams talk

2007-10-13 Thread Timothy Clemans

Why is one of the screen shots with something about number theory of
the old notebook?

On 10/13/07, William Stein <[EMAIL PROTECTED]> wrote:
>
> Hi,
>
> Here is a link to my slides and demo for my ams plenary talk:
>
>   http://sagemath.org/why/alb/
>
> --
> William Stein
> Associate Professor of Mathematics
> University of Washington
> http://wstein.org
>
> >
>

--~--~-~--~~~---~--~~
You received this message because you are subscribed to the Google Groups 
"sage-devel" group.
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?hl=en
-~--~~~~--~~--~--~---



[sage-devel] Re: multivariate polynomial coefficients

2007-10-13 Thread didier deshommes

2007/10/13, Joel B. Mohler <[EMAIL PROTECTED]>:
>
> Thanks for the discussion about this topic.  I send this mail to re-iterate
> and summarize.  It seems there are two things that you might want:
> 1)  Get the coefficient of a specific monomial in the multivariate polynomial
> ring.
> 2)  Get the coefficient of the polynomial in a tower of (two) polynomial
> extensions
>
> I suggest that MPolynomial.coefficient method perform function #1 above.  It
> currently does some strange things over the base rings I've tested.  It does
> vary on base ring due to implementation differences.
>
> I suggest that there be a new method which performs function #2 and I suggest
> it be called "polynomial_coefficient", but I really don't like that name.  I
> almost prefer "coefficient_polynomial" since it will sort beside coefficient
> in documentation and help people realize that if they are not happy
> with "coefficient" that they have an alternative.  Of course, each method's
> documentation should probably mention the other.

+1,
Those should be 2 separate methods, as they try to solve 2 (slightly)
different problems.

>
> I'm wondering if we could have a vote on preferred syntax.  I'm not going to
> describe the parameters because if they are not clear enough from context, it
> probably isn't a good parameter choice :) :
>
> sage: P.=ZZ[]
> sage: f=(1-v)*(1-2*w)*(1-3*x); f
> -6*v*w*x + 2*v*w + 3*v*x + 6*w*x - v - 2*w - 3*x + 1
> ##
> #  Alternative Number 1
> ##
> sage: f.polynomial_coefficient({w:0,v:1})
> 3*x - 1

+1 for this notation.

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.7.alpha0

2007-10-13 Thread Jason Grout

mabshoff wrote:
> 
> 
> On Oct 13, 10:10 am, "William Stein" <[EMAIL PROTECTED]> wrote:
>> Hi,
>>
>> I've posted the first alpha release of sage-2.8.7 here:
>>
>>http://sage.math.washington.edu/tmp/
>>
>> It is probably somewhat broken.  Building it, contributing
>> fixes, etc., would be greatly appreciated.
> 
> Hello,
> 
> building it went without a hitch, but there is trouble with the
> doctests:

[snip]

> **
> 
> sage -t  devel/sage-main/sage/graphs/graph.py
> **
> File "graph.py", line 2116:
> sage: h.edges()
> Expected:
> [((2, 3, None), (3, 4, None), None),
> ((1, 2, None), (2, 4, None), None),
> ((1, 2, None), (2, 3, None), None),
> ((1, 3, None), (3, 4, None), None)]
> Got:
> [((1, 3, None), (3, 4, None), None), ((2, 3, None), (3, 4, None),
> None), ((1, 2, None), (2, 4, None), None), ((1, 2, None), (2, 3,
> None), None)]
> **

This is odd.  The test passes on my Ubuntu 32-bit P4 system.  It looks
like the edges are returned in a different order for you, so the answer
is correct, just output differently.

I can't build the alpha right now, but someone that has the alpha built,
can you see if the following code returns the indicated result?  If so,
we can change the doctest to the following.

sage: g=DiGraph([[1..4], lambda i,j: ihttp://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.7.alpha0

2007-10-13 Thread Jaap Spies

Jason Grout wrote:

> This is odd.  The test passes on my Ubuntu 32-bit P4 system.  It looks
> like the edges are returned in a different order for you, so the answer
> is correct, just output differently.
> 
> I can't build the alpha right now, but someone that has the alpha built,
> can you see if the following code returns the indicated result?  If so,
> we can change the doctest to the following.
>

On Fedora 7:
[EMAIL PROTECTED] sage-2.8.7-alpha0]$ ./sage -t  
devel/sage-main/sage/graphs/graph.py
sage -t  devel/sage-main/sage/graphs/graph.py
  [25.8 s]

--
All tests passed!
Total time for all tests: 25.8 seconds


> sage: g=DiGraph([[1..4], lambda i,j: i sage: h=g.line_graph()
> sage: h.vertices()
> [(1, 2, None),
>  (1, 3, None),
>  (1, 4, None),
>  (2, 3, None),
>  (2, 4, None),
>  (3, 4, None)]
> sage: h.am()
> [0 0 0 1 1 0]
> [0 0 0 0 0 1]
> [0 0 0 0 0 0]
> [0 0 0 0 0 1]
> [0 0 0 0 0 0]
> [0 0 0 0 0 0]
> 

sage: g=DiGraph([[1..4], lambda i,j: ihttp://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.7.alpha0

2007-10-13 Thread Jason Grout

Jaap Spies wrote:
> Jason Grout wrote:
> 
>> This is odd.  The test passes on my Ubuntu 32-bit P4 system.  It looks
>> like the edges are returned in a different order for you, so the answer
>> is correct, just output differently.
>>
>> I can't build the alpha right now, but someone that has the alpha built,
>> can you see if the following code returns the indicated result?  If so,
>> we can change the doctest to the following.
>>
> 
> On Fedora 7:
> [EMAIL PROTECTED] sage-2.8.7-alpha0]$ ./sage -t  
> devel/sage-main/sage/graphs/graph.py
> sage -t  devel/sage-main/sage/graphs/graph.py
>   [25.8 s]
> 
> --
> All tests passed!
> Total time for all tests: 25.8 seconds


Great, thanks.  I filed a patch at #878.

Jason


--~--~-~--~~~---~--~~
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.7.alpha0

2007-10-13 Thread Justin C. Walker


On Oct 13, 2007, at 11:22 AM, Jason Grout wrote:
mabshoff wrote:
>> On Oct 13, 10:10 am, "William Stein" <[EMAIL PROTECTED]> wrote:
[snip]
> This is odd.  The test passes on my Ubuntu 32-bit P4 system.  It looks
> like the edges are returned in a different order for you, so the  
> answer
> is correct, just output differently.
>
> I can't build the alpha right now, but someone that has the alpha  
> built,
> can you see if the following code returns the indicated result?  If  
> so,
> we can change the doctest to the following.
>
> sage: g=DiGraph([[1..4], lambda i,j: i sage: h=g.line_graph()
> sage: h.vertices()
> [(1, 2, None),
>  (1, 3, None),
>  (1, 4, None),
>  (2, 3, None),
>  (2, 4, None),
>  (3, 4, None)]
> sage: h.am()
> [0 0 0 1 1 0]
> [0 0 0 0 0 1]
> [0 0 0 0 0 0]
> [0 0 0 0 0 1]
> [0 0 0 0 0 0]
> [0 0 0 0 0 0]

Built from scratch on a Mac OS X (10.4.10) Core Duo.  The tests  
aren't complete, but many are failing.  For this, I get:

sage: g=DiGraph([[1..4], lambda i,j: ihttp://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.7.alpha0

2007-10-13 Thread cwitty

On Oct 13, 10:30 am, Jaap Spies <[EMAIL PROTECTED]> wrote:
> In sloane_functions.py sequence can now be indexed by reals!?
> a(1.0) returns a value. Formerly an error was raised.

In SloaneSequence.__call__() (which is shared among all sequences) the
first thing it does is "m = ZZ(n)", to coerce its argument to an
integer.  As of 2.8.7, coercing reals->integer succeeds if the real is
exactly equal to an integer.  a(1.1) will still fail.

So: this is a feature, not a bug.  (And I will be submitting a patch
to sloane_functions.py accordingly.)

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: sage-2.8.7.alpha0

2007-10-13 Thread Jaap Spies

cwitty wrote:
> On Oct 13, 10:30 am, Jaap Spies <[EMAIL PROTECTED]> wrote:
>> In sloane_functions.py sequence can now be indexed by reals!?
>> a(1.0) returns a value. Formerly an error was raised.
> 
> In SloaneSequence.__call__() (which is shared among all sequences) the
> first thing it does is "m = ZZ(n)", to coerce its argument to an
> integer.  As of 2.8.7, coercing reals->integer succeeds if the real is
> exactly equal to an integer.  a(1.1) will still fail.
> 
> So: this is a feature, not a bug.  (And I will be submitting a patch
> to sloane_functions.py accordingly.)
> 

Sorry, but I don't know yet I consider that a feature! In the context
of the OEIS is seems silly to allow index 20.0, but maybe there are
arguments I do not understand. Deleting my examples seems the best way
to go.

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: sage-2.8.7.alpha0

2007-10-13 Thread mabshoff

Carl Witty has opened tickets for all the issues in trac and there are
already some patches. So if you want to help out joins us in IRC, so
that way we do not duplicate efforts.

The additional number of doctests failing on FC7 in 32 bit mode might
be related to coercion issues, see #833 for the blocker and some back
traces. Valgrind also seems to indicate the same/similar issue.

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: andlinux

2007-10-13 Thread Jaap Spies

Bill Page wrote:
> William,
> 
> I am building Sage from source under andLinux right now. Everything
> looks good so far. I will report when it completes. It is also
> possible that an existing linux binary could work, perhaps provided
> some prerequisite libraries are added.
> 
> I really like andLinux so far and I would very interested in other
> reactions and comments.
> 
> Regards,
> Bill Page.
> 
> On 10/12/07, William Stein <[EMAIL PROTECTED]> wrote:
>> Hi,
>>
>> Does anybody want to try to make sage run on Windows via this:
>>
>>http://www.andlinux.org/
>>
>> Bill Page just pointed out that they're using it for FriCAS (=Axiom fork).
>>
>> It's definitely not a panacea, but it might be a good idea to support it for
>> some users.
>>
>>  -- William
>>

Yes, I could install it, download sage-2.8.6.tar, compile it
and run it. But I've problems for instance to run
install_scripts('/usr/local/bin')
because there the "'" is not accepted.

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: sage-2.8.7.alpha0

2007-10-13 Thread Justin C. Walker


On Oct 13, 2007, at 1:10 AM, William Stein wrote:

>
> Hi,
>
> I've posted the first alpha release of sage-2.8.7 here:
>
>http://sage.math.washington.edu/tmp/
>
> It is probably somewhat broken.  Building it, contributing
> fixes, etc., would be greatly appreciated.

Here're the test results for a Mac OS X/Core Duo run with the alpha0:

The following tests failed:


 sage -t  devel/sage-main/sage/combinat/sloane_functions.py
 sage -t  devel/sage-main/sage/dsage/tests/testdoc.py
 sage -t  devel/sage-main/sage/graphs/graph.py
 sage -t  devel/sage-main/sage/matrix/matrix_integer_dense.pyx
 sage -t  devel/sage-main/sage/rings/integer_mod.pyx
 sage -t  devel/sage-main/sage/rings/morphism.pyx
 sage -t  devel/sage-main/sage/rings/polynomial/real_roots.pyx
 sage -t  devel/sage-main/sage/rings/residue_field.pyx
 sage -t  devel/sage-main/sage/schemes/elliptic_curves/ 
ec_database.py
 sage -t  devel/sage-main/sage/sets/set.py
Total time for all tests: 2613.1 seconds

Justin

--
Justin C. Walker, Curmudgeon-At-Large
Institute for the Enhancement of the Director's Income

When LuteFisk is outlawed,
Only outlaws will have LuteFisk





--~--~-~--~~~---~--~~
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.7.alpha0

2007-10-13 Thread Jaap Spies

Hi Michael,

mabshoff wrote:
> Carl Witty has opened tickets for all the issues in trac and there are
> already some patches. So if you want to help out joins us in IRC, so
> that way we do not duplicate efforts.
> 

No, sorry but I will not join the party. I need a good night rest!
Do you ever sleep?

The next days my sailing season will end.
See: http://www.jaapspies.nl/sailing.html

> The additional number of doctests failing on FC7 in 32 bit mode might
> be related to coercion issues, see #833 for the blocker and some back
> traces. Valgrind also seems to indicate the same/similar issue.
> 

Wednesday I'll see what has been done.

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: sage-2.8.7.alpha0

2007-10-13 Thread cwitty

On Oct 13, 1:56 pm, Jaap Spies <[EMAIL PROTECTED]> wrote:
> cwitty wrote:
> > On Oct 13, 10:30 am, Jaap Spies <[EMAIL PROTECTED]> wrote:
> >> In sloane_functions.py sequence can now be indexed by reals!?
> >> a(1.0) returns a value. Formerly an error was raised.
>
> > In SloaneSequence.__call__() (which is shared among all sequences) the
> > first thing it does is "m = ZZ(n)", to coerce its argument to an
> > integer.  As of 2.8.7, coercing reals->integer succeeds if the real is
> > exactly equal to an integer.  a(1.1) will still fail.
>
> > So: this is a feature, not a bug.  (And I will be submitting a patch
> > to sloane_functions.py accordingly.)
>
> Sorry, but I don't know yet I consider that a feature! In the context
> of the OEIS is seems silly to allow index 20.0, but maybe there are
> arguments I do not understand. Deleting my examples seems the best way
> to go.

Well, I could be convinced otherwise.  I think it's a good thing that
ZZ(3.0) == 3, but Sloane sequences could do more input validation than
just calling ZZ().  I don't actually care whether an index of 20.0
works or not.

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: sage-2.8.7.alpha0

2007-10-13 Thread Jaap Spies

cwitty wrote:

> Well, I could be convinced otherwise.  I think it's a good thing that
> ZZ(3.0) == 3, but Sloane sequences could do more input validation than
> just calling ZZ().  I don't actually care whether an index of 20.0
> works or not.
> 

Yes, you are right, there should be more input validation.

You don't care about the indexing, but in the context of integer
sequences the index must be an integer, Most entries of the OEIS
are supposed to be indexed by positive integers, some are by
non negative integers and quite a few allow a negative index,
but Neil Sloane is removing this option.

Regards,

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: Braid Goups / Coxeter Groups / Root Systems / ...

2007-10-13 Thread David Joyner

On 10/13/07, Mike Hansen <[EMAIL PROTECTED]> wrote:
>
> I just saw Robert's ticket about implementing / wrapping braid groups
> in SAGE ( http://www.sagetrac.org/sage_trac/ticket/876 ) which is
> pretty coincidental since I was thinking about the exact same thing
> before I went to bed last night.  I wrote an email to Nikos
> Apostolakis who had expressed interest in this on sage-support in June
> to see if he was still interested.

Please also check out the GAP package braid
http://sage.math.washington.edu/home/wdj/gap/braid.tar.gz
written by K. Magaard, S. Shpectorov and H. Voelklein,
and HAP, which I think has some briad-related stuff but I can't find
it right now.

>
> GAP3 had support for braid groups provided by the "share package"
> CHEVIE.  Unfortunately, when the change to GAP4 was made, this was
> never ported over.  As David Joyner mentioned,  one possibility is to
> have a (optional) GAP3 package from which the functionality could be
> wrapped.  Another possibility is to port the GAP3 chevie code over to
> GAP4.  I do not know how much of an undertaking this would be because
> my familiarity with GAP programming is almost non-existent.  The third
> possibility would be to port chevie directly over to SAGE; I looked at
> the code for chevie, and this is definitely doable.  It would however
> take a bit of time -- it's about 6000 lines of GAP code and about 2000
> lines of comments.


I spoke with Frank Luebeck while in Germany last month. While a nice
guy and a very talented mathematician and programmer, Frank is
very very particular about the code he releases. The Chevie port
to GAP4 is basically finished and in my opinion should be released.
It would be a complete waste of time to work on GAP3 stuff.
I told Frank that it should be released now and he simply disagrees.
(It is his work on this for the most part, though he does have co-authors.)
He said it is high on his priority list to go through the code and carefully
bug-check it as soon as possible.

However, he has not ported the Hecke algebra section of Chevie and has no
immediate plans to do so.


>
> It would be good to have a framework for root systems and the like
> that all of these things can fit into.  This may be helpful in that
> aspect: http://mupad-combinat.sourceforge.net/Wiki/RootSystems



>
> --Mike
>
> >
>

--~--~-~--~~~---~--~~
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: my plenary ams talk

2007-10-13 Thread David Joyner

On 10/13/07, William Stein <[EMAIL PROTECTED]> wrote:
>
> Hi,
>
> Here is a link to my slides and demo for my ams plenary talk:
>
>   http://sagemath.org/why/alb/

Looks like it will be a great talk. Wish I could make it!

>
> --
> William Stein
> Associate Professor of Mathematics
> University of Washington
> http://wstein.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: andlinux

2007-10-13 Thread Bill Page

On 10/13/07, Jaap Spies wrote:
> ...
> > On 10/12/07, William Stein <[EMAIL PROTECTED]> wrote:
> >> Hi,
> >>
> >> Does anybody want to try to make sage run on Windows via this:
> >>
> >>http://www.andlinux.org/
> >>
> >> Bill Page just pointed out that they're using it for FriCAS (=Axiom fork).
> >>
> >> It's definitely not a panacea, but it might be a good idea to support it 
> >> for
> >> some users.
> >>
> >>  -- William
> >>
>
> Yes, I could install it, download sage-2.8.6.tar, compile it
> and run it.

Great! What version of Windows are you using?

> But I've problems for instance to run
> install_scripts('/usr/local/bin')
> because there the "'" is not accepted.
>

When I run install_scripts I get the following:

andLinux:~/sage-2.8.6# ./sage
--
| SAGE Version 2.8.6, Release Date: 2007-10-06   |
| Type notebook() for the GUI, and license() for information.|
--

sage: install_scripts('/usr/local/bin')

Checking that SAGE has the command 'gap' installed
Created script '/usr/local/bin/gap'

Checking that SAGE has the command 'gp' installed
Created script '/usr/local/bin/gp'

Checking that SAGE has the command 'singular' installed
Created script '/usr/local/bin/singular'

Checking that SAGE has the command 'maxima' installed
Created script '/usr/local/bin/maxima'

Checking that SAGE has the command 'M2' installed
The command 'M2' is not available; not adding shortcut

Checking that SAGE has the command 'kash' installed
The command 'kash' is not available; not adding shortcut

Checking that SAGE has the command 'mwrank' installed
Created script '/usr/local/bin/mwrank'

Checking that SAGE has the command 'ipython' installed
Created script '/usr/local/bin/ipython'

Checking that SAGE has the command 'hg' installed
Created script '/usr/local/bin/hg'

Checking that SAGE has the command 'hgmerge' installed
Created script '/usr/local/bin/hgmerge'
Finished creating scripts.
You need not do this again even if you upgrade or move SAGE.
The only requirement is that the command 'sage' is in the PATH.
sage: exit
Exiting SAGE (CPU time 0m0.08s, Wall time 2m33.33s).
andLinux:~/sage-2.8.6# ls /usr/local/bin
axiom  gap  gp  hg  hgmerge  ipython  maxima  mwrank  singular
andLinux:~/sage-2.8.6#



What do you mean by: "the "'" is not accepted"?

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: andlinux

2007-10-13 Thread Jaap Spies

Bill Page wrote:
> On 10/13/07, Jaap Spies wrote:
>> ...
>>> On 10/12/07, William Stein <[EMAIL PROTECTED]> wrote:
 Hi,

 Does anybody want to try to make sage run on Windows via this:

http://www.andlinux.org/

 Bill Page just pointed out that they're using it for FriCAS (=Axiom fork).

 It's definitely not a panacea, but it might be a good idea to support it 
 for
 some users.

  -- William

>> Yes, I could install it, download sage-2.8.6.tar, compile it
>> and run it.
> 
> Great! What version of Windows are you using?

Windows XP, etcetera.

> 
>> But I've problems for instance to run
>> install_scripts('/usr/local/bin')
>> because there the "'" is not accepted.
>>

> 
> 
> What do you mean by: "the "'" is not accepted"?
> 

I could not type the command: install_scripts('/usr/local/bin')
the first ' was replaced by 
which makes no sense to me!

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: multivariate polynomial coefficients

2007-10-13 Thread cwitty

On Oct 13, 6:51 am, "Joel B. Mohler" <[EMAIL PROTECTED]> wrote:
> I'm wondering if we could have a vote on preferred syntax.  I'm not going to
> describe the parameters because if they are not clear enough from context, it
> probably isn't a good parameter choice :) :
>
> sage: P.=ZZ[]
> sage: f=(1-v)*(1-2*w)*(1-3*x); f
> -6*v*w*x + 2*v*w + 3*v*x + 6*w*x - v - 2*w - 3*x + 1
> ##
> #  Alternative Number 1
> ##
> sage: f.polynomial_coefficient({w:0,v:1})
> 3*x - 1

As you might guess, I vote for #1.

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: my plenary ams talk

2007-10-13 Thread Timothy Clemans

In number four on page 11: "SAGE has a Huge, active, and well rounded developer
community: sage-devel mailing list has over 200 subscribers,
working very hard on everything from highly optimized arithmetic,
to high school education, to computing modular forms. Usually
about 30 people get patches accepting into SAGE every month."

accepting should be accepted

On 10/13/07, David Joyner <[EMAIL PROTECTED]> wrote:
>
> On 10/13/07, William Stein <[EMAIL PROTECTED]> wrote:
> >
> > Hi,
> >
> > Here is a link to my slides and demo for my ams plenary talk:
> >
> >   http://sagemath.org/why/alb/
>
> Looks like it will be a great talk. Wish I could make it!
>
> >
> > --
> > William Stein
> > Associate Professor of Mathematics
> > University of Washington
> > http://wstein.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: andlinux

2007-10-13 Thread Bill Page

On 10/13/07, Jaap Spies wrote:
> ...
> I could not type the command: install_scripts('/usr/local/bin')
> the first ' was replaced by 
> which makes no sense to me!
>

It doesn't make any sense to me either. What happens when you type
something else with a quoted string? E.g.

  maxima('1+1')

Does that work or do you get the  in that case too?

What about if you type " instead of ' ? Like this:

install_scripts("/usr/local/bin")

Are you using any special (non-US) keyboard configuration?

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/
-~--~~~~--~~--~--~---