MathGroup Archive 2008

[Date Index] [Thread Index] [Author Index]

Search the Archive

Glitch in NumberFieldFundamentalUnits

  • To: mathgroup at
  • Subject: [mg92373] Glitch in NumberFieldFundamentalUnits
  • From: "Doug Hensley" <dahensley at>
  • Date: Mon, 29 Sep 2008 07:08:46 -0400 (EDT)
  • References: <D0D555E836994F10AE26728DDF2AECB2@DougPC> <>

Try this:

Timing[NumberFieldFundamentalUnits[Root[x^3+4 x^2-10 x-1,1]]]

You'll get something like this:

Prepend::normal : Nonatomic expression expected at position 1 in 
Prepend[System`NumberFieldDump`pvec,AlgebraicNumber[Root[-1-10 #1+4 
Slot[â? 1â?¡]2+#13&,1,0],{104,-39,-10}]].

{8491.51,NumberFieldFundamentalUnits[Root[-1-10 #1+4 #12+#13&,1]]}

Now try it with another package that offers this functionality, just to get 
a peek at the answer:


You'll get something like this: 

in far less than one second. (Meaning, that if alpha is one of the roots of 
p=x^3+4x^2-10x-1, then the set {alpha,126alpha^2-184 alpha-81} comprises a 
fundamental system of units for Q(alpha). )

Mathematica's NumberFieldFundamentalUnits gets fundamental unit systems for 
hundreds of comparably complicated inputs, and the outputs it obtains are 
not ever so much simpler in the other cases. Therefore, there's a bug. It's 
the better kind of bug, because it doesn't return a wrong answer. It just 
hangs and then supplies an error message.

  • Prev by Date: Re: Launch notebook in FE and evaluate it immediately
  • Next by Date: Re: ClearAll error message
  • Previous by thread: Graph Theory / Cycle Matrix
  • Next by thread: mathematica newbie ?