FAQFAQ   SearchSearch   MemberlistMemberlist   UsergroupsUsergroups   RegisterRegister 
 ProfileProfile   PreferencesPreferences   Log in to check your private messagesLog in to check your private messages   Log inLog in 
Forum index » Science and Technology » Math » Symbolic
Magma Computation needed !
Post new topic   Reply to topic Page 1 of 1 [6 Posts] View previous topic :: View next topic
Author Message
Jürgen Böhm
science forum beginner


Joined: 10 Jun 2005
Posts: 15

PostPosted: Fri Jun 10, 2005 12:59 pm    Post subject: Magma Computation needed ! Reply with quote

Hello,

maybe someone with access to MAGMA can help me and try to execute the
following, probably quite expensive computation on his computer:

Z:=IntegerRing();
L<x, ac, pa, pb, pc, aa, bb, qa, qb, pd, qc, qd, ab, ad, ba, bc, bd, ra
, rb, rc, sb, sc, sd, ub, uc, ud, wa, wb, wc, A,
B>:=PolynomialRing(Z,31,"elim",29);
I:=ideal<L| pd*pb^2-pa*pc^2, qd*pb^2+2*pd*qb*pb-qa*pc^2-2*pa*qc*pc,
aa*pa+ba*qa-1
, ab*pb+bb*qb-1, ac*pc+bc*qc-1, ad*pd+bd*qd-1, ra*pa+sb*pb-1,
rb*pb+sc*pc-x^3-
A*x-B, rc*pc+sd*pd-1, wa*pa+ra*qa+ub*pb+sb*qb,
wb*pb+rb*qb+uc*pc+sc*qc-3*x^2-A
, wc*pc+rc*qc+ud*pd+sd*qd>;

Groebner(I);

I tried it with the MAGMA online calculator
http://modular.fas.harvard.edu/cgi-bin/calc/calc.py but it consumed too
much time and memory (limited there to 20 seconds). I would be very much
interested in the result, actually the question is whether n1 *
(4*A^3+27*B^2)^n2 with integers n1 and n2 appears in the Groebner-base.

Thanks in advance

Jürgen Böhm


-------------------------------------------------------------------------
Dipl.-Math. Jürgen Böhm e-mail: reverse: net dot gmx at jboehm
"At a time when so many scholars in the world are calculating, is it not
desirable that some, who can, dream ?" R. Thom
Back to top
Roman Pearce
science forum beginner


Joined: 03 May 2005
Posts: 37

PostPosted: Fri Jun 10, 2005 1:37 pm    Post subject: Re: Magma Computation needed ! Reply with quote

I'm getting down to a 120480 by 333717 matrix in the F4 algorithm.
Now, since this is a matrix over Z I'm going to guess that you're
doomed. However, I'll let it run for a few days (2.5 GHz PPC G5).

Roman
Back to top
Roman Pearce
science forum beginner


Joined: 03 May 2005
Posts: 37

PostPosted: Fri Jun 10, 2005 2:41 pm    Post subject: Re: Magma Computation needed ! Reply with quote

Sure enough I ran out of memory (4GB). You may be able to do this with
an 8GB machine. Does it have to be over the ring of integers, or did
you want rational number coefficients ?
Back to top
Jürgen Böhm
science forum beginner


Joined: 10 Jun 2005
Posts: 15

PostPosted: Fri Jun 10, 2005 6:45 pm    Post subject: Re: Magma Computation needed ! Reply with quote

Roman Pearce wrote:
Quote:

Sure enough I ran out of memory (4GB). You may be able to do this with
an 8GB machine. Does it have to be over the ring of integers, or did
you want rational number coefficients ?

Let me first thank you for your fast answer and your generous offer
of so much computer time.

Regarding the question of coefficients: This is a good question, that
I could not fully answer to myself. The goal of the computation is to
prove the squarefreeness of the so called division-polynomial of an
elliptic curve y^2=x^3+A*x+B over a field K with char K != 2, 3

If char K = 0 it would be enough to compute over the rational numbers,
but with characteristic p it would be useful to know that the hoped-for
element n1*(4*A^3+27*B^2)^n2 has only 2 and 3 as prime factors in n1.

At first I was quite sure that this had to be true, so it would be
worthwhile to find a groebner base over Z, but now it came to my mind
that [p] :E->E (multiplication by p on an elliptic curve E) has not
kernel Z/pZ x Z/pZ for p=char K. So the reasons to assume squarefreeness
of the division-polynomial psi[p] seem to break down in this case.

What is the conclusion of all this ? I think it would be worthwhile to
do the calculation over Q if this is easier, and I will have to do some
experimental calculations to see, if my conjecture is actually true in
char K = p.

By the way: What machine did you use ? Was it an Apple with PPC or do
other workstations with PPC also exist ? (The magma online calculator
says to use a dual Opteron Sun-workstation - I would be very happy if I
could afford such a system and one of these expensive Magma licenses
too..)

Jürgen

-------------------------------------------------------------------------
Dipl.-Math. Jürgen Böhm e-mail: reverse: net dot gmx at jboehm
"At a time when so many scholars in the world are calculating, is it not
desirable that some, who can, dream ?" R. Thom
Back to top
Roman Pearce
science forum beginner


Joined: 03 May 2005
Posts: 37

PostPosted: Sat Jun 11, 2005 12:18 am    Post subject: Re: Magma Computation needed ! Reply with quote

Well, I'm a graduate student so things like computer time and Magma
licenses are much cheaper Smile I used an Apple PPC machine. You can
buy PPC machines from IBM as well, but they are generally quite
expensive (with server-grade hardware, etc).

As for your problem, computing Groebner bases over a finite field is
much faster than either the rational numbers or over Z, especially in
Magma. Will it help you to do this problem over GF(p^n) ?

Roman
Back to top
Jürgen Böhm
science forum beginner


Joined: 10 Jun 2005
Posts: 15

PostPosted: Sat Jun 11, 2005 7:38 pm    Post subject: Re: Magma Computation needed ! Reply with quote

Roman Pearce wrote:
Quote:

.....
As for your problem, computing Groebner bases over a finite field is
much faster than either the rational numbers or over Z, especially in
Magma. Will it help you to do this problem over GF(p^n) ?


Unfortunately it seems to be necessary to compute at least over Q for
the characteristic 0 case.

In characteristic p the division-polynomial psi[p] has - at least for
the cases I checked with experimental calculations - multiple roots
(contrary to what I first conjectured). To get a full understanding of
the cases with char K <> 0 it would be necessary to compute over Z. As I
actually don't need a groebner base, but only polynomials fi with
sum(i=1..12) fi * hi = n1*(4*A^3+27*B^2)^n2, where fi are polynomials in
the variables x, ac, pa, pb, pc,... A, B and hi are the polynomials
from the ideal I:

h1= pd*pb^2-pa*pc^2,
h2= qd*pb^2+2*pd*qb*pb-qa*pc^2-2*pa*qc*pc
.....

one could try an ansatz for the fi and solve a system of linear
equations. I already tried to do this with maple and got a system that
looked quite easy on first sight, but Maple didn't seem to make any
progress on it.

Jürgen

-------------------------------------------------------------------------
Dipl.-Math. Jürgen Böhm e-mail: reverse: net dot gmx at jboehm
"At a time when so many scholars in the world are calculating, is it not
desirable that some, who can, dream ?" R. Thom
Back to top
Google

Back to top
Display posts from previous:   
Post new topic   Reply to topic Page 1 of 1 [6 Posts] View previous topic :: View next topic
The time now is Wed Jan 07, 2009 10:54 pm | All times are GMT
Forum index » Science and Technology » Math » Symbolic
Jump to:  

Similar Topics
Topic Author Forum Replies Last Post
No new posts TURMEL: Parker Crown ducks exam he sa... John Turmel Engineering 1 Thu Jul 13, 2006 2:07 am
No new posts help with symbolic matrix computation... csviks Symbolic 0 Fri Jul 07, 2006 8:58 pm
No new posts help needed with spatial dependence niraj doshi Math 0 Fri Jul 07, 2006 2:56 am
No new posts Quality Manager needed : Full Time : ... gitproinc@yahoo.com Chem 0 Wed Jul 05, 2006 11:55 pm
No new posts Quality Manager needed : Full Time : ... gitproinc@yahoo.com Chem 0 Wed Jul 05, 2006 11:48 pm

The eBay Song | Mortgages | Money News | Problem Mortgage | Magazine Subscriptions
Copyright © 2004-2005 DeniX Solutions SRL
Other DeniX Solutions sites: Electronics forum |  Medicine forum |  Unix/Linux blog |  Unix/Linux documentation |  Unix/Linux forums


Powered by phpBB © 2001, 2005 phpBB Group
[ Time: 0.2132s ][ Queries: 16 (0.1259s) ][ GZIP on - Debug on ]