[sage-support] Re: Different behaviour of "sage -gap" and "gap pexpect interface"

2014-03-03 Thread Simon King
Hi! A partial reply to my own question: On 2014-03-03, Simon King wrote: > Doing the same in the pexpect interface, I get > sage: gap.RequirePackage('"ModIsom"') > Traceback (most recent call last): > ... > RuntimeError: Gap produced error output > Error, Variable: 'WeightVector' is re

[sage-support] Re: bug on error reporting in sage notebook?

2014-03-03 Thread kcrisman
On Monday, March 3, 2014 4:20:22 PM UTC-5, Fred Gruber wrote: > > In the sage 6 notebook , if I type the following R code (in R mode) > myfunc2<-function(x){ > print(x) > if (x>1){ > print('x larger than 1') > } > > > where there is a missing bracket I do not get any errors from sage. It

[sage-support] Re: plotting matrix transformations of unit ball in the p-norm

2014-03-03 Thread P Purkayastha
We used something like this for a class once. You can follow which point is transformed where def vector_plotter(vec_list, points=False, **kwds): pic = Graphics() for c, v in zip(rainbow(len(vec_list), 'rgbtuple'), vec_list): if points: pic += point(v, rgbcolor=c, *

[sage-support] bug on error reporting in sage notebook?

2014-03-03 Thread Fred Gruber
In the sage 6 notebook , if I type the following R code (in R mode) myfunc2<-function(x){ print(x) if (x>1){ print('x larger than 1') } where there is a missing bracket I do not get any errors from sage. It just does not define the function myfunc2 This is very serious because if I have

[sage-support] Different behaviour of "sage -gap" and "gap pexpect interface"

2014-03-03 Thread Simon King
Hi! I thought that sage's pexpect interface to GAP just launches "sage -gap" and then sends strings that are evaluated in the subprocess. Hence, I am very surprised by the following. I have installed GAP's ModIsom package, and in "sage -gap" I see that it really is available: gap> RequirePackag

Re: [sage-support] Re: sage->nauty->directed graphs isomorphism vs. magma -- smth is very wrong

2014-03-03 Thread Volker Braun
AFAIK both implement the same basic algorithm. My guess is that nauty quickly realizes that it has a huge problem and starts constructing graph invariants in an attempt to answer it to the negative. If you were to compare two graphs that are actually isomorphic I'm pretty sure that you wouldn't

Re: [sage-support] Re: sage->nauty->directed graphs isomorphism vs. magma -- smth is very wrong

2014-03-03 Thread Tom Boothby
Once the graph is constructed, is_isomorphic throws away the vertex labels and just works with pointers to ints. Constructing the graph itself happens in a blink. Sadly, the bulk of the time is spent in is_isomorphic. The answer to Aleksandr's question is: Yes, this is a known issue -- nauty has

[sage-support] Re: plotting matrix transformations of unit ball in the p-norm

2014-03-03 Thread Jason Grout
On 3/2/14, 9:52, martyall wrote: I would like to plot the "effects" of an 2 by 2 matrix over the reals on the unit sphere in the p-norm. By effects I mean I would like to plot the unit sphere and it's image under a given matrix. In clearest case when p=2 I found the following code from trac #9728

[sage-support] Re: sage->nauty->directed graphs isomorphism vs. magma -- smth is very wrong

2014-03-03 Thread Chris Godsil
I would not be surprised it it was the finite field arithmetic that is causing the difference. On Friday, February 28, 2014 4:18:44 PM UTC-5, Aleksandr Kodess wrote: > > As far as I know both sage and magma utilize Brendan McKay's program nauty > in order to check whether two given graphs (direc

Re: [sage-support] sage.libs.pari.gen.PariError: ZZ_123...1[y]/(y^2 + ...) is not a field in FpX_ffintersect

2014-03-03 Thread martin . vgagern
Reported as http://trac.sagemath.org/ticket/15886 Thanks for listening. -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-support+unsubscr...@googlegroups.co

Re: [sage-support] sage.libs.pari.gen.PariError: ZZ_123...1[y]/(y^2 + ...) is not a field in FpX_ffintersect

2014-03-03 Thread Jeroen Demeyer
On 2014-03-03 18:03, martin.vgag...@gmx.net wrote: I guess now is a good time for the bug report Absolutely, go for it! -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an e

Re: [sage-support] sage.libs.pari.gen.PariError: ZZ_123...1[y]/(y^2 + ...) is not a field in FpX_ffintersect

2014-03-03 Thread martin . vgagern
On Monday, March 3, 2014 5:31:18 PM UTC+1, martin@gmx.net wrote: > Nevertheless, here is a minimal example: > > https://gist.githubusercontent.com/gagern/9320350/raw/d1896f7a5a05b24075098941c9c3ff156ca6c139/MinimalReprodcing6.sage > > It is minimal in the following sense: For the whole express

Re: [sage-support] sage.libs.pari.gen.PariError: ZZ_123...1[y]/(y^2 + ...) is not a field in FpX_ffintersect

2014-03-03 Thread martin . vgagern
On Monday, March 3, 2014 11:30:50 AM UTC+1, Jeroen Demeyer wrote: > > > I have no clue what that error message means, and more importantly, how > to avoid it. > The concrete problem (I guess) is that 123633238138850861 is not a prime > number. But now the question becomes: where does this number

[sage-support] Re: Most efficient way to enumerate group elements

2014-03-03 Thread Karl Möller
I have been able to get ~1.7x with this somewhat ugly hack: %cython from sage.all import * def permutationMatrices(n): import numpy cdef int deg = SymmetricGroup(n).degree() cdef int i_max = SymmetricGroup(n).order() cdef int chunk = 1 cdef int i = 0 step = range(0,i_ma

[sage-support] Polynomial rings and homogenization

2014-03-03 Thread Tristan
Hi all, I'm trying to use Sage to implement the number field sieve factorization algorithm, but I've reached a sticking point with defining polynomials and some of their in built functions. I'm running Sage-5.13 through VirtualBox 4.3.8 and homogenization does not work as expected. Following th

Re: [sage-support] sage.libs.pari.gen.PariError: ZZ_123...1[y]/(y^2 + ...) is not a field in FpX_ffintersect

2014-03-03 Thread Jeroen Demeyer
On 2014-03-03 09:11, martin.vgag...@gmx.net wrote: As you can see in that file, I've encountered quite a number of places where the current behavior of sage surprised or even annoyed me. I think as a general rule, you should try to avoid symbolic computations for any serious work. Try to conver

Re: [sage-support] sage.libs.pari.gen.PariError: ZZ_123...1[y]/(y^2 + ...) is not a field in FpX_ffintersect

2014-03-03 Thread Jeroen Demeyer
On 2014-03-03 09:11, martin.vgag...@gmx.net wrote: sage.libs.pari.gen.PariError: ZZ_123633238138850861[y]/(y^2 + 114948053005503479*y + 71680486823734693) is not a field in FpX_ffintersect I have no clue what that error message means, and more importantly, how to avoid it. The concrete problem

[sage-support] sage.libs.pari.gen.PariError: ZZ_123...1[y]/(y^2 + ...) is not a field in FpX_ffintersect

2014-03-03 Thread martin . vgagern
I'm trying to find a formula concerned with some computation about the relative position of two ellipses. But some result along the way looks strange when evaluated with RDF numbers, and trying to plug in algebaric numbers instead, I encounter the following error message: sage.libs.pari.gen.Par

[sage-support] Most efficient way to enumerate group elements

2014-03-03 Thread Karl Möller
Hi I am using SAGE to interface with GAP. (I don't really know GAP, just playing with the SAGE interface to it). My problem is that I am trying to efficiently enumerate the group elements (and get their matrix representations) of various permutation groups. Now, of course this problem quickly

[sage-support] plotting matrix transformations of unit ball in the p-norm

2014-03-03 Thread martyall
I would like to plot the "effects" of an 2 by 2 matrix over the reals on the unit sphere in the p-norm. By effects I mean I would like to plot the unit sphere and it's image under a given matrix. In clearest case when p=2 I found the following code from trac #9728 which does exactly what I want.

[sage-support] Group and subgroup have different domains

2014-03-03 Thread Fredrik
Hi, This might be more appropriate at the ask.sage site, but it is down at the moment. The problem is this: I have a permutation group G defined as the group of permutations on {0..14}. Then I ask for subgroups, but the subgroups are indexed with {1..15}. Are the indices just translated 1 to t