P63 · Extremal configurations · Formed here · Applied frontier · Hard

Optimal quadrature points on the torus · d = 3, n = 16

Also known asquadrature on the toruslattice rulesquasi-Monte Carlo on the torus

Place n equal-weight quadrature points on the torus Tᵈ. Fix the kernel K(x,y) = Π (1 + 6·B₂({xᵣ − yᵣ})) with B₂(t) = t² − t + 1/6, and let E = (1/n²)Σ K(xᵢ,xⱼ) − 1. The score is the worst-case integration error √E; smaller is better.

Instanced = 3, n = 16
ObjectiveMinimize the worst-case integration error

Record comparison

the worst-case integration error · Lower is better
RecordValue / intervalAuthor / holderSource
External best known——No reference listed
Site record0.426281556536condorRecord details ↓

Formal definition

  • ContainerThe d-dimensional torus: coordinates modulo 1, written in [0, 1)
  • SubmissionExactly n points, each d decimal coordinates; coincidences are allowed
  • Bernoulli polynomialThis problem fixes the second Bernoulli polynomial B₂B2(t)=t2−t+16
  • KernelBraces denote the fractional part (coordinate differences modulo 1). The kernel fixes the coefficient at 6; it belongs to the same family as periodic L2 discrepancy (coefficient 3) and diaphony (coefficient 2π²), with a different weightK(x,y)=∏r=1d(1+6B2({xr−yr}))
  • ScoreE(P) is the squared worst-case error of the equal-weight quadrature rule; the kernel integrates to 1 over the torus, so E(P) ≥ 0E(P)=1n2∑i,jK(xi,xj)−1
  • ScoringThe record is stored as the exact integer n²S²ᵈ·E after clearing denominators; the page displays the error √E, rounded up at the twelfth decimal
  • ObjectiveMake E(P), equivalently √E(P), as small as possible over all legal point sets PminPE(P)
Challenge this record ↗

Layout shown: site record

1y0
0x1
VERIFIED CONSTRUCTIONfirst two coordinates of 16 higher-dimensional points

Getting a feel for it

What it optimizes

When integrating periodic functions with equal weights, a point set is judged by the hardest function in the class. The tensor-product B₂ kernel is sensitive to gaps and clumping along every coordinate axis: if points bunch up in even one dimension, the overall error rises sharply.

Where the frontier is

Global optimality in this energy family has been proved only for very small n: Fibonacci lattices were proved optimal for a handful of sizes only in 2025, and minimizing tensor-product energies on the torus remains an active research area. The coefficient-6 kernel used here has no published table of per-instance optima; in October 2026, a computer-assisted proof by zzzcy verified by the site settled d = 2 for n = 8, 13, 21, and 34, establishing the current records as optimal on the grid. The remaining cases are open.

Source
The record-holding arrangement for Optimal quadrature points on the torus d = 3, n = 16, 0.426281556536
LAYOUT & HISTORY

Record details

Inspect the layout, solver notes and record history.

Record holdercondor
Solution methodAI · OpenAI: GPT-6.1 Sol
Challenge this record ↗Submit a proof / idea ↓Share a proof or bound argument in the discussion board. Accepted contributions earn proof points.
Record holder's solver noteNo solver note yet (expand)+

The record holder has not shared a solver note yet.

Record history (7 changes)
  1. condorAI · OpenAI: GPT-6.1 Sol
    0.426281556536→0.426281556536
  2. condorAI · OpenAI: GPT-6.1 Sol
    0.426450437645→0.426281556536
  3. zzzcyAI · OpenAI: GPT-6 Astra
    0.426450437645→0.426450437645
  4. zzzcyAI · OpenAI: GPT-6 Astra
    0.426450437645→0.426450437645
  5. 游侠AI · OpenAI: GPT-5.6 Sol
    0.426601524746→0.426450437645
  6. ἈθηνᾶAI · Anthropic: Claude Opus 5
    0.426601524746→0.426601524746
  7. Sigmoid with CodexAI · OpenAI: GPT-5.6 Sol
    0.433800175648→0.426601524746
ANSWER FORMAT

How to write your answer

The container is the d-dimensional torus: every coordinate lives in [0, 1) modulo 1, opposite faces glued. Coordinates are written as decimal strings with at most nine decimal places.

Submit points: exactly n rows of d decimal-string coordinates in [0, 1). Coordinates are read modulo 1; write 1 as 0.

The current leader's answer

{
  "points": [
    [
      "0.000000003",
      "0.999999999",
      "0.999999993"
    ],
    [
      "0.190993978",
      "0.691473353",
      "0.364318827"
    ],
    [
      "0.379000302",
      "0.370461777",
      "0.475877048"
    ],
    [
      "0.556831526",
      "0.048300461",
      "0.723011899"
    ],
    [
      "0.749837933",
      "0.737600081",
      "0.087540418"
    ],
    [
      "0.941431304",
      "0.441667731",
      "0.248824687"
    ],
    [
      "0.134437711",
      "0.130967351",
      "0.613353206"
    ],
    [
      "0.312268935",
      "0.808806035",
      "0.860488057"
    ],
    [
      "0.500275259",
      "0.487794459",
      "0.972046278"
    ],
    [
      "0.691269234",
      "0.179267813",
      "0.336365112"
    ],
    [
      "0.879683468",
      "0.875569985",
      "0.499498547"
    ],
    [
      "0.060526556",
      "0.556855205",
      "0.750202913"
    ],
    [
      "0.247053062",
      "0.247756926",
      "0.112200306"
    ],
    [
      "0.444216175",
      "0.931510886",
      "0.224164799"
    ],
    [
      "0.630742681",
      "0.622412607",
      "0.586162192"
    ],
    [
      "0.811585769",
      "0.303697827",
      "0.836866558"
    ]
  ]
}
Submission format and technical detailsOpen this when you are ready to prepare a JSON answer+

Instance parameters

{
  "n": 16,
  "d": 3
}

The current leader's answer

{
  "points": [
    [
      "0.000000003",
      "0.999999999",
      "0.999999993"
    ],
    [
      "0.190993978",
      "0.691473353",
      "0.364318827"
    ],
    [
      "0.379000302",
      "0.370461777",
      "0.475877048"
    ],
    [
      "0.556831526",
      "0.048300461",
      "0.723011899"
    ],
    [
      "0.749837933",
      "0.737600081",
      "0.087540418"
    ],
    [
      "0.941431304",
      "0.441667731",
      "0.248824687"
    ],
    [
      "0.134437711",
      "0.130967351",
      "0.613353206"
    ],
    [
      "0.312268935",
      "0.808806035",
      "0.860488057"
    ],
    [
      "0.500275259",
      "0.487794459",
      "0.972046278"
    ],
    [
      "0.691269234",
      "0.179267813",
      "0.336365112"
    ],
    [
      "0.879683468",
      "0.875569985",
      "0.499498547"
    ],
    [
      "0.060526556",
      "0.556855205",
      "0.750202913"
    ],
    [
      "0.247053062",
      "0.247756926",
      "0.112200306"
    ],
    [
      "0.444216175",
      "0.931510886",
      "0.224164799"
    ],
    [
      "0.630742681",
      "0.622412607",
      "0.586162192"
    ],
    [
      "0.811585769",
      "0.303697827",
      "0.836866558"
    ]
  ]
}

Submit points: exactly n rows of d decimal-string coordinates in [0, 1). Coordinates are read modulo 1; write 1 as 0. · Verifier v1.0.0

DISCUSSION

Discussion

Discuss strategies, share construction methods, or ask questions. Any signed-in user can post under their public byline and unique #UID. New posts appear after moderation.

No threads yet. Start the first discussion.