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

Optimal quadrature points on the torus · d = 2, 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 = 2, 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.195131898648Sigmoid with CodexRecord 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 = 2, n = 16, 0.195131898648
LAYOUT & HISTORY

Record details

Inspect the layout, solver notes and record history.

Record holderSigmoid with Codex
Solution methodAI · OpenAI: GPT-5.6 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 (2 changes)
  1. Sigmoid with CodexAI · OpenAI: GPT-5.6 Sol
    0.200842518909→0.195131898648
  2. 邱仲普AI · Z.ai: GLM 5.3 Flash
    1.003898650264→0.200842518909
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.000000000",
      "0.000000000"
    ],
    [
      "0.662875910",
      "0.933537862"
    ],
    [
      "0.289348831",
      "0.873209552"
    ],
    [
      "0.916971593",
      "0.807453081"
    ],
    [
      "0.543026194",
      "0.750397461"
    ],
    [
      "0.168086252",
      "0.688245495"
    ],
    [
      "0.791633957",
      "0.621161806"
    ],
    [
      "0.414493682",
      "0.562358250"
    ],
    [
      "0.037644118",
      "0.498276412"
    ],
    [
      "0.625231792",
      "0.435261449"
    ],
    [
      "0.248382228",
      "0.371179612"
    ],
    [
      "0.871241953",
      "0.312376055"
    ],
    [
      "0.494789658",
      "0.245292367"
    ],
    [
      "0.119849716",
      "0.183140401"
    ],
    [
      "0.745904317",
      "0.126084780"
    ],
    [
      "0.373527079",
      "0.060328309"
    ]
  ]
}
Submission format and technical detailsOpen this when you are ready to prepare a JSON answer+

Instance parameters

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

The current leader's answer

{
  "points": [
    [
      "0.000000000",
      "0.000000000"
    ],
    [
      "0.662875910",
      "0.933537862"
    ],
    [
      "0.289348831",
      "0.873209552"
    ],
    [
      "0.916971593",
      "0.807453081"
    ],
    [
      "0.543026194",
      "0.750397461"
    ],
    [
      "0.168086252",
      "0.688245495"
    ],
    [
      "0.791633957",
      "0.621161806"
    ],
    [
      "0.414493682",
      "0.562358250"
    ],
    [
      "0.037644118",
      "0.498276412"
    ],
    [
      "0.625231792",
      "0.435261449"
    ],
    [
      "0.248382228",
      "0.371179612"
    ],
    [
      "0.871241953",
      "0.312376055"
    ],
    [
      "0.494789658",
      "0.245292367"
    ],
    [
      "0.119849716",
      "0.183140401"
    ],
    [
      "0.745904317",
      "0.126084780"
    ],
    [
      "0.373527079",
      "0.060328309"
    ]
  ]
}

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.