18
$\begingroup$

Hot to plot the field of algebraic numbers in the complex plane?

enter image description here

In this picture, the color of a point indicates the degree of the polynomial of which it’s a root:

red = rational numbers

green = roots of quadratic polynomials,

blue = roots of cubic polynomials

yellow = roots of quartic polynomials, and so on

I tried first few steps:

data = Table[(-b + Sqrt[b^2 - 4 a c])/(2 a), {a, 1, 100}, {b, -100, 100}, {c, -100, 100}];

ListPlot[{Re[#], Im[#]} & /@ data, AxesOrigin -> {0, 0}, 
PlotStyle -> [Green, PointSize[.02]]]

but it doesn't work..

Thank you in advance to any one who may be able to give me some ideas

$\endgroup$
4
  • 1
    $\begingroup$ The techniques in this question might be useful. $\endgroup$ May 23, 2016 at 17:57
  • $\begingroup$ Your link has the relevant subroutines in c++ code, but sadly the author left out a few required file: lset.c and rnd/frnd.c $\endgroup$
    – M.R.
    May 23, 2016 at 18:38
  • 1
    $\begingroup$ I've done similar algebraic number renderings before, see flickr.com/photos/104348204@N05/11637086656/sizes/o (very large image) and flickr.com/photos/104348204@N05/16570026460/sizes/o (moderately large image). The first image is cubics (with a special restriction on the linear coefficient) and the second is an image of the quartic roots. The typical recipe is to generate a large image array of zeros, lay down numbers at various points, blur the image, and colorize it. $\endgroup$ May 23, 2016 at 19:20
  • $\begingroup$ Do you know what determines the size of these points in that picture? $\endgroup$
    – Greg Hurst
    May 24, 2016 at 15:01

1 Answer 1

25
$\begingroup$

Your code works fine, but it's missing half the roots, and a Flattening of the list of numbers prior to applying Re and Im helps. Adding those in:

data = Flatten[
   Table[{(-b + Sqrt[b^2 - 4 a c])/(2 a), (-b - 
        Sqrt[b^2 - 4 a c])/(2 a)}, {a, 1, 20}, {b, -20, 20}, {c, -20, 
     20}]];
ListPlot[{Re[#], Im[#]} & /@ data, PlotRange -> {{-3, 3}, {-3, 3}}, 
 AspectRatio -> 1]

gets you this:

enter image description here

Which is pretty nice!

If you'd like to generate images more like the one in the Wiki article or my examples, you can do the following steps:

  1. Generate a list of complex roots of polynomials of some order N, where the vector of polynomial coefficients $\mathbf{v}=(a_0,\ldots,a_N)$ satisfies $|\mathbf{v}|\leq R$ for some cutoff radius $R$.

  2. Convert each root into a 2d coordinate by taking the real and imaginary parts.

  3. Create a sparse array, and increment the element at each root's coordinate by $\frac{1}{|\mathbf{v}|^2}$. That way, roots of simple polynomials (ie, polynomials whose coefficient vectors lie near the origin in $\mathbb{R}^{N+1}$) are "brightest", and roots of complicated polynomials (ie, polynomials whose coefficient vectors lie near the surface of the hypersphere) are dimmer. You'll need to translate, scale, and round the coordinates to be positive integers, as matrix indices must be positive integers.

  4. Blur the sparse array. The raw sparse array isn't pretty, and blurring will reveal make it obvious which points are bright, and which are dim. Using the built-in GaussianFilter is the quickest way, but it's not the best-looking. My preferred way is to FFT-convolve with a Lorentz-like distribution $\left(\frac{\gamma^2}{\gamma^2+x^2+y^2}\right)^\alpha$, where $\alpha$ determines the tail heaviness; $\alpha=1.15$ works pretty good.

  5. Colorize the resulting array $M$ and render it using Image. One easy way is to feed Image the tensor-product $M\otimes\mathbf{c}$ where $\mathbf{c}$ is some 3-element RGB vector. In Mathematica 10, TensorProduct preserves PackedArrays, whereas in older versions it doesn't, so beware of RAM consumption if you're using an older version and want large images. Incidentally, that link has some examples of algebraic number plots.

Here's an example you can get from the above:

enter image description here

IIRC the green points are $N=3$ (cubics), and the yellow ones are $N=2$ (quadratics).

Similar with just quartics (I added additional red lighting for points near the origin):

enter image description here

And a large image of monic cubics where the linear coefficient is restricted to be an 8th root of unity:

https://www.flickr.com/photos/104348204@N05/11637086656/sizes/o/

Example Code

Here's some code to do the above:

https://www.dropbox.com/s/ushalt5feb19wal/Root%20Generator.nb?dl=0 https://www.dropbox.com/s/wo4405fb73yrxgp/Root%20Visualizer.nb?dl=0

The "Root Generator" generates roots, converts them into a sparse array, and exports it to a file (steps 1-3). The "Root Visualizer" imports the files, and renders them (steps 4-5). It's got decent comments/instructions inside, but you may need to tweak the parameters to get what you want.

$\endgroup$
5
  • $\begingroup$ Hello DD, can you maybe include an example implementation of your idea in this answer? I know that you've probably shown some of it in the other answers, but a user might benefit from seeing your description and sample code side-by-side to get the idea. $\endgroup$ May 24, 2016 at 4:38
  • $\begingroup$ @J.M. Added some links to sample code (rather than pasting the code into the answer, as it's a little long). $\endgroup$ May 24, 2016 at 5:02
  • $\begingroup$ Very nice. I'll test it out as soon as I'm in front of a proper machine. :) $\endgroup$ May 24, 2016 at 5:13
  • $\begingroup$ @DumpsterDoofus +1 very nice. thanks :) $\endgroup$
    – vito
    May 25, 2016 at 9:45
  • $\begingroup$ I guess I am probably confused, but who is responsible for creating the great graphics above? I would like to ask permission to use them in papers on number theory and metamathematics. I would also love to give due credit and acknowledgment to the maths artist/s. Uh...what I mean is, what is the proper legal name I should use in published docs. $\endgroup$ Jul 12, 2017 at 0:33

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.