/* G is a positive-definite integral Gram matrix, e.g. A2. */
G = [2,-1; -1,2];
qfauto(G)[1] /* 12 */Every value is an exact integer. The cubic and root-lattice orders are computed by integer products, with the diagram-symmetry factors included; duals use the same orders. The laminated and Coxeter-Todd orders are computed from their integral catalogue Gram matrices by PARI qfauto [4], implementing the Plesken-Souvignier algorithm [2]. The attached generator checks that each returned matrix has determinant $\pm1$ and preserves the Gram matrix.
Independent Gram-matrix computations verify all 116 cubic, root and dual entries, with known controls of orders $8$ and $12$ for the square and hexagonal lattices. Contragredient matrices preserve the inverse Gram forms for all 44 listed duals. All 17 catalogue matrices agree entry by entry with the published catalogue and with the attachment of the classical-lattice packing table. The 26 applicable stored Coxeter Poincaré polynomials evaluated at $1$ agree with the root and cubic orders after the diagram-symmetry factors. Further checks compare the exceptional root orders and the Leech order with the catalogue, the Barnes-Wall order with OEIS A014116, the real Coxeter-Todd order with the index-two relation on page $429$ of [1], and the dimension-$23$ and dimension-$24$ orders with twice the stored orders of $\mathrm{Co}_2$ and $\mathrm{Co}_1$. The Coxeter-Todd source has a printed final factor $8$ where its index-two argument and independent Gram computation give $7$.