Example AlgEtQ_UniformizersCompletion (H42E8)

> _<x> := PolynomialRing(Integers());
> f := (x^8+16)*(x^8+81);
> A := EtaleAlgebra(f);
> O := MaximalOrder(A);
> // Consider a bunch of prime of O and their uniformizers in O.
> pp := PrimesAbove(2*3*5*7*O);
> unifs := Uniformizers(pp);
> // We now verify that each element is a uniformizer at the correct prime and a unit everywhere else
> assert Matrix([[ Valuation(mP(t)) : t in unifs ] where AP, mP := Completion(P) : P in pp]) eq 1;