← LibraryComputing the Maximal Order: the Round 2 AlgorithmEngineering · MathematicsLesson 1/4← PrevNext →
GuidePublished 6 Aug 20265 min readBy Kevin JoginComputational Number TheoryNumber Fields IIRound 2 AlgorithmMaximal Order
Skip to the main content

MathematicsNumber Fields II

Computing the Maximal Order: the Round 2 Algorithm

Enlarging an order one prime at a time until it is maximal — the standard method, and the factoring obstruction it cannot escape.

Executive summary

Maximal at every prime implies maximal

An order is the maximal order exactly when it is maximal at every prime dividing the index, and only primes whose square divides the polynomial discriminant can do so. Round 2 therefore factors that discriminant, tests each candidate prime with the Dedekind criterion, and for each failure computes the p-radical and its ring of multipliers — a strictly larger order. Iterating at each prime terminates because the index drops by a factor of p each time.

Learning objectives

  • State the Pohst–Zassenhaus theorem and its role.
  • Apply the Dedekind criterion as a cheap maximality test.
  • Compute the p-radical and the ring of multipliers.
  • Trace the Round 2 loop and prove its termination.
  • Report an order correctly when the discriminant is not fully factored.

Section 01The Pohst–Zassenhaus theorem

An order O is maximal if and only if it is p-maximal for every prime p with p2 dividing disc(O). This localises the problem: instead of one global search, a finite list of independent local computations.

  1. Stage 01Compute disc(T)By the subresultant algorithm.
  2. Stage 02Factor the square partOnly primes whose square divides the discriminant can divide the index — the squarefree part is discarded immediately.
  3. Stage 03Test each candidateApply the Dedekind criterion. Most primes pass and need no further work.
  4. Stage 04Enlarge where neededFor each failing prime, compute the p-maximal overorder; merge all results into a single integral basis by Hermite normal form.
The factoring obstruction

Step 2 is the bottleneck, and it is not an implementation weakness — no algorithm is known that determines the maximal order without it. Where the discriminant cannot be fully factored, the result must be reported as maximal only at the primes that were resolved.

Section 02The p-radical and the ring of multipliers

The p-radical of O is the set of elements some power of which lies in pO — equivalently, the nilradical of O/pO pulled back to O.

AlgorithmComputing the p-radicalin: order O, prime p  →  out: the p-radical Ip
  1. Let q be the smallest power of p with q ≥ n, where n is the degree.
  2. Compute the matrix of the map x ↦ xq on O/pO, an Fp-linear map. Linear because the Frobenius is additive in characteristic p.
  3. The p-radical Ip is the preimage in O of the kernel of that map.
  4. Return Ip as an HNF module containing pO.
One kernel computation over Fp. The exponent q ensures every nilpotent element is captured, since nilpotency degree is bounded by the dimension.

The ring of multipliers of the radical is then the candidate larger order:

O′ = { xK : x IpIp }
Why the ring of multipliers is larger

If O is not p-maximal, the radical fails to be invertible, and elements outside O multiply it into itself. Those elements are algebraic integers not in O, so O′ strictly contains O. When O is already p-maximal the construction returns O itself, which is the loop's termination signal.

Section 03The Round 2 algorithm

AlgorithmRound 2 (p-maximal overorder)in: order O, prime p  →  out: the p-maximal overorder of O
  1. Apply the Dedekind criterion. If O is p-maximal, return O. Cheap test first — it succeeds for most primes.
  2. Compute the p-radical Ip.
  3. Compute the ring of multipliers O′ of Ip, as the kernel of a linear map expressed via the regular representation.
  4. If O′ = O, then O is p-maximal; return O.
  5. Otherwise set O ← O′ and return to step 2. The index [O′ : O] is a positive power of p, so disc drops by at least p2.
Termination: each iteration divides the p-part of the index by at least p, and the index is finite. The number of iterations is bounded by the p-adic valuation of the index, typically very small.
Round 4 and the alternatives

Round 2 is the standard because it is simple and robust. The Round 4 algorithm of Ford and Letard uses Newton polygons and Montes-style techniques to obtain the same result faster for high-degree fields and for primes with large ramification, at the cost of considerably more intricate implementation.

Section 04Reporting the result

How to record a maximal order computation
SituationCorrect claim
disc(T) fully factoredK determined unconditionally
Unfactored part proved squarefreeK determined unconditionally
Unfactored part remainsOrder maximal at all primes below the factoring bound; field discriminant known only up to a square factor
Only trial division performedMaximal at small primes; treat the result as provisional
Downstream consequences of a non-maximal order

Every later computation — prime decomposition, class group, units — assumes maximality. Running them on a non-maximal order gives well-formed and incorrect answers. The conditionality must therefore propagate with the data, not be recorded once and forgotten.

ReferenceFrequently asked questions

How often does the Dedekind criterion succeed immediately?

For most primes. It fails only when the factorisation of the defining polynomial modulo p has repeated factors and an additional gcd condition holds, which is uncommon — so the expensive radical computation typically runs for only a handful of primes.

Can the maximal order be computed for large degree fields?

Yes, subject to the discriminant factoring. The linear algebra is polynomial in the degree; the practical limit is almost always the size of the discriminant rather than the degree itself.

What if the defining polynomial has huge coefficients?

Reduce it first. Polynomial reduction — using LLL on the lattice of the order to find a generator of small height — can shrink the discriminant enormously and may convert an infeasible factorisation into an easy one.

NavigateContinue in this stream

Curated next steps from this page. The site also surfaces algorithmically related reading below.

ProvenanceSources and further reading

This page is an original KEVOS explanatory article. It presents the underlying mathematics — definitions, algorithms, complexity results and selection criteria — in KEVOS editorial voice. No text is reproduced from any copyrighted source. Where numerical tables are relevant, KEVOS links to live authoritative databases rather than republishing static values.

Page ID
KV-MATH-0036
Taxonomy
ENG-MATH — Engineering / Mathematics
Collection
COL-CANT-001
Topic stream
CANT-ADVANCED-FIELDS
Version
1.1.0 / content 2026.08
Last reviewed
2026-08-06

Continue learning

NEXT LESSON →Prime Decomposition: the Buchmann–Lenstra MethodGuide · MathematicsComputing Galois Groups of Number FieldsGuide · MathematicsClass Group and Unit Computation in General Number FieldsGuide · MathematicsDiscriminants and Integral BasesGuide · Mathematics