Re: [sage-support] Echelon Form

2011-07-19 Thread amaseam
2011/7/19 Santanu Sarkar :
> Size of the matrix is (30,16). Entries are at least 3000 bit integer.

Try using pastebin.com.

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


[sage-support] Re: Lattice Reduction of two matrix

2011-07-19 Thread Maarten Derickx
Some matrix procedures also have an option which gives you the 
transformation matrix.

sage: M=MatrixSpace(ZZ,4,4)
sage: A=M.random_element()
sage: A.echelon_form(transformation=True)
(
[   101 3032]  [  85   31  4847]
[   011 2176]  [  61   22  3475]
[   002 5636]  [ 158   57  899   13]
[   000 9095], [ 255   92 1452   21]
)

So an efficient way would be just to multiply with that matrix.

Currently there is sadly enough no such option for LLL reduction so you have 
to find an explicit base change yourself.

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


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

2011-07-19 Thread Christophe BAL
Thanks a lot for this informations.

Christophe.

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


Re: [sage-support] Echelon Form

2011-07-19 Thread Santanu Sarkar
Size of the matrix is (30,16). Entries are at least 3000 bit integer.

On 17 July 2011 19:21, William Stein  wrote:

> On Thu, Jul 14, 2011 at 11:26 PM, Santanu Sarkar
>  wrote:
> > I want to find Echelon Form using following:
> > E,U=M3.echelon_form(transformation=True)
> >
> > But it terminates with the following message:
> >
> > ***   the PARI stack overflows !
> >   current stack size: 1600 (15.259 Mbytes)
> >   [hint] you can increase GP stack with allocatemem()
> >
> > Traceback (click to the left of this block for traceback)
> > ...
> > RuntimeError.
> >
> >
> > Similar problem  for Hermite Normal Form
>
> Post the *exact* input that replicates this problem.
>
> William
>
> >
> >
> >
> > --
> > To post to this group, send email to sage-support@googlegroups.com
> > To unsubscribe from this group, send email to
> > sage-support+unsubscr...@googlegroups.com
> > For more options, visit this group at
> > http://groups.google.com/group/sage-support
> > URL: http://www.sagemath.org
> >
>
>
>
> --
> William Stein
> Professor of Mathematics
> University of Washington
> http://wstein.org
>
> --
> To post to this group, send email to sage-support@googlegroups.com
> To unsubscribe from this group, send email to
> sage-support+unsubscr...@googlegroups.com
> For more options, visit this group at
> http://groups.google.com/group/sage-support
> URL: http://www.sagemath.org
>

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


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

2011-07-19 Thread Mike Hansen
On Tue, Jul 19, 2011 at 8:00 AM, Christophe BAL  wrote:
> Hello,
> I would like to know how Sage calls Maxima from Python ? Which files have I
> to look for in the source ?

Sage (as of 4.7) uses pexpect to communicate with Maxima.  The main
places you'll want to look are in
$SAGE_ROOT/devel/sage/sage/interfaces/maxima.py and
$SAGE_ROOT/devel/sage/sage/interfaces/expect.py .

In an 4.7.1, Sage will use ECL to use Maxima as "library" for it's
internal use of Maxima. See
http://trac.sagemath.org/sage_trac/ticket/7377 for more information on
this.

--Mike

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


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

2011-07-19 Thread Christophe BAL
Hello,
I would like to know how Sage calls Maxima from Python ? Which files have I
to look for in the source ?

Best regards.
Christophe

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org


Re: [sage-support] Re: How to calculate the affine coordinats of a point

2011-07-19 Thread Johannes
that sounds comprehensibly. explains some more errormessages i got a few
times (like 'bool has no attribute len' when using solve)

greatz

Am 18.07.2011 16:54, schrieb luisfe:

> For me it looks like:
> 
> In solve, when writting a == b you assume that a and b are expressions
> involving several variables. If a and b are expressions, then a == b
> is also an expression.
> 
> However, p and x*p1+y*p2 are NOT expressions, but vectors. And
> equality of vectors is not the same as equality of expressions
> 
> sage: x == y
> x == y
> sage: type(_)
> 
> sage: vector([x]) == vector([y])
> False
> sage: type(_)
> 
> 
> So, in fact, you are passing the following command:
> 
> sage: [x * p1 + y * p2 == p]
> [False]
> 
> sage: solve([False],x,y)
> 
> Which has no solution.  It is subtle, but I would not consider it a
> bug. If you really want to use solve, you may try the following:
> 
> sage: solve(x * p1 + y * p2 - p,x,y)
> [[x == (1/4), y == (3/4)]]
> 
> In this case, the input is a vector, that is an iterable, so solve
> extracts its components and equals them to zero.
> 
> Being said that, I recommend you to use the linear algebra
> interpretation that I suggested, since it will probably be much more
> efficient.
> 

-- 
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org