{
 "artifacts": [
  {
   "name": "bench_report.json",
   "url": "https://synthcore.org/code/1144/bench_report.json"
  },
  {
   "name": "custom_beam_track.py",
   "url": "https://synthcore.org/code/1144/custom_beam_track.py"
  },
  {
   "name": "experiment.py",
   "url": "https://synthcore.org/code/1144/experiment.py"
  },
  {
   "name": "report.md",
   "url": "https://synthcore.org/code/1144/report.md"
  },
  {
   "name": "results.json",
   "url": "https://synthcore.org/code/1144/results.json"
  },
  {
   "name": "results2.json",
   "url": "https://synthcore.org/code/1144/results2.json"
  },
  {
   "name": "run_bench.py",
   "url": "https://synthcore.org/code/1144/run_bench.py"
  }
 ],
 "category": "training",
 "description": "Use a spectral eigenvalue-counting function to bracket each target mode before neural optimization. The network then solves only within an interval containing exactly one eigenfrequency, preventing optimization from repeatedly collapsing to the lowest mode or jumping between modes.",
 "download_zip": "https://synthcore.org/download/1144",
 "formulas_latex": [
  "$$j(\\omega^*)=j_0+\\operatorname{sgn}\\!\\left(\\det\\mathbf W(\\omega^*)\\right)n(\\omega^*)$$",
  "$$j_0=\\sum_m j_m$$",
  "$$j_m=\\operatorname{highest\\ integer}\\!\\left[\\frac{\\alpha}{\\pi}\\right]+\\frac{1}{2}\\left(2-\\operatorname{sign}(\\beta L)-\\operatorname{sign}\\!\\left(\\beta L-\\frac{(\\gamma L)^2}{\\beta L}\\right)\\right)$$",
  "$$\\mathcal L_k(\\theta,\\omega)=\\|\\mathbf B\\mathbf W(\\omega)\\mathbf A q_\\theta\\|_2^2+\\lambda_n(q_\\theta^T\\mathbf Mq_\\theta-1)^2+\\lambda_{\\mathrm{box}}[\\max(0,\\omega_k^- -\\omega)^2+\\max(0,\\omega-\\omega_k^+)^2]$$"
 ],
 "id": 2950,
 "implementation": "1. Integration point: add a preprocessing and training wrapper around a neural eigenmode solver. The wrapper receives an assembled symmetric dynamic stiffness matrix \\(W(\\omega)\\), a boundary-condition-reduced matrix, and a requested mode index \\(k\\). 2. Pseudocode: scan an increasing frequency grid \\(\\omega_i\\); at each point perform an LDLᵀ factorization of the reduced matrix and count negative diagonal pivots \\(n(\\omega_i)\\); compute the element correction \\(j_0=\\sum_mj_m\\) using the supplied \\(\\alpha,\\beta,\\gamma\\) formula; evaluate \\(j(\\omega_i)\\). Find adjacent grid points \\([\\omega_k^-,\\omega_k^+]\\) with \\(j(\\omega_k^-)=k-1\\) and \\(j(\\omega_k^+)=k\\), then bisect until the interval width is below tolerance. Initialize a separate MLP head for each bracket, optimize its parameters and frequency using the DSM residual and box penalty, and optionally freeze frequency for the first 20% of training. 3. LDLᵀ inertia and bisection are deterministic; the network estimates only the eigenvector. Near poles or tiny pivots, use a threshold \\(|d_i|\u003c10^{-10}\\|W\\|_\\infty\\) and repeat with a small frequency perturbation. 4. First experiment: use clamped-free and clamped-clamped beams with 16 elements and target modes 1–8. Compare unconstrained DSM-PINN, DSM-PINN with orthogonality penalties, and WW-bracketed DSM-PINN using equal 10,000-step budgets. Evaluate frequency error, modal assurance criterion against FEM eigenvectors, and the number of runs converging to an already-found lower mode. 5. Pre-register mechanism tests: every accepted bracket must contain exactly one reference eigenfrequency; the count must be monotone nondecreasing on a frequency grid except for numerical tolerance; and the fraction of runs converging to a lower mode must fall below 5%, versus at least 25% for unconstrained optimization. Measure bracket correctness against a high-accuracy generalized eigensolver and mode identity using \\(\\mathrm{MAC}(u,v)=|u^Tv|^2/[(u^Tu)(v^Tv)]\\). Falsification is systematic count non-monotonicity, more than 2% incorrect brackets, or no reduction in lower-mode collapse.",
 "math_summary": "The paper uses the Wittrick–Williams relation \\(j(\\omega^*)=j_0+\\operatorname{sgn}(\\det\\mathbf W(\\omega^*))n(\\omega^*)\\), where \\(j(\\omega^*)\\) is the number of structure eigenfrequencies below \\(\\omega^*\\), \\(j_0=\\sum_mj_m\\) counts element frequencies when all nodal displacements are restrained, and \\(n(\\omega^*)\\) is the number of negative pivots in an LDLᵀ factorization of the assembled boundary-conditioned dynamic stiffness matrix. For the beam element, the extracted formula is \\(j_m=\\operatorname{highest\\ integer}[\\alpha/\\pi]+0.5(2-\\operatorname{sign}(\\beta L)-\\operatorname{sign}(\\beta L-(\\gamma L)^2/(\\beta L)))\\). The count is piecewise constant and increases by one at each eigenfrequency. The transferable mechanism is inertia-based eigenvalue counting of an indefinite frequency-dependent operator, which produces certified brackets for mode isolation.",
 "math_tags": [
  "spectral-theory",
  "linear-algebra",
  "numerical-analysis",
  "dynamical-systems"
 ],
 "ml_areas": [
  "initialization",
  "training-dynamics",
  "loss"
 ],
 "paper": {
  "arxiv_id": "2608.28683",
  "arxiv_url": "https://arxiv.org/abs/2608.28683",
  "summary_what_math_gives_to_ml": "The paper's transferable asset is not the beam-specific matrix itself, but the replacement of derivative-heavy strong-form PINN residuals by an assembled, frequency-dependent operator whose element basis already satisfies the governing differential equation. This creates a derivative-free physics loss and exposes the spectral structure of the problem through the nonlinear matrix pencil \\(\\mathbf W(\\omega)\\). The Wittrick–Williams count supplies a constructive way to bracket and enumerate higher eigenfrequencies, avoiding mode collapse and repeated convergence to the lowest eigenvalue. The strongest ML transfer is therefore a dynamic-stiffness operator layer combined with certified frequency bracketing for neural eigenmode models on discretized linear PDEs.",
  "title": "A Framework Integrating the Dynamic Stiffness Matrix with Physics-Informed Neural Networks for Solving Eigenvalue Problems and Analysing Dynamic Response",
  "year": "2026"
 },
 "ratings": {
  "difficulty": 5,
  "novelty": 8,
  "usefulness": 7
 },
 "solves": [
  "accuracy",
  "stability",
  "sample-efficiency"
 ],
 "title": "Wittrick–Williams Mode Enumerator",
 "url": "https://synthcore.org/idea/2950/wittrick-williams-mode-enumerator",
 "verification": {
  "peer_reviewed": false,
  "stage1_mechanism_check": {
   "worked": true,
   "confidence": 7,
   "verdict": "Built a compact Wittrick–Williams-style inertia enumerator for a clamped-free FEM beam, with frequency scanning, bisection bracketing, and bracketed neural optimization. The count was monotone and the target bracket contained exactly one reference eigenfrequency. In eight fixed-seed trials, lower-mode collapse decreased from 12.5% to 0%, but exact target-mode identification remained 12.5% and MAC stayed weak; therefore the isolation mechanism worked, while an end-to-end accuracy win was not demonstrated.",
   "metrics": {
    "baseline": "8 seeds: lower-mode rate 12.5%, exact target-mode identification 12.5%",
    "idea": "8 seeds: lower-mode rate 0%, exact target-mode identification 12.5%; bracket width 6.16e-4 and exactly 1 eigenfrequency contained"
   },
   "how_to_run": "python3 experiment.py \u003e results2.json",
   "files": [
    "experiment.py",
    "results2.json",
    "REPORT.txt"
   ],
   "limitations": "Only one 8-element clamped-free beam and target mode were tested; the FEM inertia count substitutes for the full beam-element j0 correction, CUDA/CPU parity and LDL pivot perturbation behavior were not independently benchmarked, and the neural optimizer was a simple direct-vector prototype rather than the paper's full DSM-PINN."
  },
  "status": "unverified",
  "status_label": "Unverified",
  "updated_at": "2026-09-02T12:51:41",
  "verdict_source": "deterministic test code (paired-seed permutation statistics)",
  "verification_axes": {
   "benchmark_mechanism": {
    "confirmed": null,
    "tested": false
   },
   "practical_benchmark": {
    "beats_baseline": null,
    "tested": false
   },
   "toy_mechanism_gate": {
    "confirmed": true,
    "tested": true
   }
  }
 }
}