{
 "artifacts": null,
 "category": "training",
 "description": "Train the PDE network against a learned family of normalized weak test functions instead of minimizing a pointwise strong-form residual. For the p-Laplace equation, this directly approximates the W^{-1,p'} residual norm and should remain effective when the solution or its second derivatives are not smooth enough for reliable collocation-point differentiation.",
 "formulas_latex": [
  "$$\\|f\\|_{W^{-1,p^{\\prime}}(\\Omega)}:=\\sup_{\\varphi\\in W^{1,p}_{0}(\\Omega)\\setminus\\{0\\}}\\frac{\\langle f,\\varphi\\rangle_{W^{-1,p^{\\prime}}(\\Omega)\\times W^{1,p}_{0}(\\Omega)}}{\\|\\varphi\\|_{W^{1,p}(\\Omega)}}.$$",
  "$$B_{\\theta}(\\varphi)=\\int_{\\Omega}|\\nabla u_{\\theta}|^{p-2}\\nabla u_{\\theta}\\cdot\\nabla\\varphi\\,dx-\\int_{\\Omega}f\\varphi\\,dx,$$",
  "$$L_{\\mathrm{dual}}(\\theta)=\\max_{\\eta}\\frac{|B_{\\theta}(\\psi_{\\eta})|}{\\left(\\int_{\\Omega}(|\\psi_{\\eta}|^{p}+|\\nabla\\psi_{\\eta}|^{p})\\,dx\\right)^{1/p}+\\epsilon}.$$",
  "$$\\|u_{\\theta}-g\\|_{W^{s,p}(\\partial\\Omega)}^{p}=\\int_{\\partial\\Omega}|u_{\\theta}(x)-g(x)|^{p}\\,d\\sigma(x)+\\int_{\\partial\\Omega}\\int_{\\partial\\Omega}\\frac{|(u_{\\theta}-g)(x)-(u_{\\theta}-g)(y)|^{p}}{|x-y|^{m+sp}}\\,d\\sigma(x)d\\sigma(y),\\quad s=1-\\frac{1}{p},\\quad m=d-1.$$"
 ],
 "id": 128,
 "implementation": "(1) Replace the interior strong-form PINN loss, which requires evaluating div(|grad u_theta|^{p-2} grad u_theta), with the weak residual above. Use an MLP u_theta(x) for the solution and a second MLP psi_eta(x) for test functions. Enforce psi_eta = 0 on the boundary by setting psi_eta(x) = b(x) times tilde_psi_eta(x), where b(x) is a known boundary-distance factor or a product of coordinate factors on a box. (2) At each outer update, draw interior points x_i and boundary points z_j. Compute a_i = |grad u_theta(x_i)|^{p-2} grad u_theta(x_i), then estimate B with B_hat = |Omega|/N times the sum over i of [a_i dot grad psi_eta(x_i) - f(x_i) psi_eta(x_i)]. Estimate the denominator with D_hat = (|Omega|/N times the sum over i of [|psi_eta(x_i)|^p + |grad psi_eta(x_i)|^p])^{1/p} + epsilon. Ascend eta on |B_hat|/D_hat for 3 to 10 critic steps, then descend theta on L_dual + lambda_B L_B. Stop gradients through the critic during the solution-network update if joint optimization becomes unstable. (3) Estimate L_B from the fractional boundary formula by sampling random boundary pairs (z_j, z_k), using the pairwise term |g_theta(z_j)-g_theta(z_k)|^p divided by |z_j-z_k|^{m+sp}. The dual-norm supremum is estimated empirically by the critic; the exponent relation p' = p/(p-1) and the weak-form construction come directly from the paper. Clip critic gradients and optionally normalize its output after each critic update so that D_hat remains near one. (4) First test on a two-dimensional unit-square manufactured problem with p = 3 and both smooth and low-regularity target solutions. Compare against a four-layer tanh strong-form PINN and a weak Galerkin-style baseline at equal collocation counts. Track relative W^{1,p} error, fractional boundary error, dense-grid PDE residual, seed variance, and wall-clock convergence. Success means lower variance and lower W^{1,p} error at equal training steps, particularly for targets with sharp gradients or insufficient second-derivative regularity.",
 "math_summary": "The paper defines the negative Sobolev norm by equation (2.1): ||f||_{W^{-1,p'}(Omega)} = sup_{varphi in W^{1,p}_0(Omega)\\{0}} \u003cf,varphi\u003e / ||varphi||_{W^{1,p}(Omega)}, where Omega is a subset of R^d, p is at least 1, p' = p/(p-1) is the Holder-conjugate exponent, W^{1,p}_0(Omega) contains test functions with zero boundary trace, and the pairing is between W^{-1,p'}(Omega) and W^{1,p}_0(Omega). For a network u_theta solving -div(|grad u|^{p-2} grad u) = f, define A_p(q) = |q|^{p-2}q and the weak residual pairing B_theta(varphi) = integral over Omega of A_p(grad u_theta) dot grad varphi minus integral over Omega of f varphi. The transferable loss is L_dual(theta) = max_eta |B_theta(psi_eta)| / (||psi_eta||_{W^{1,p}(Omega)} + epsilon), where psi_eta is a test network constrained to have approximately zero boundary trace and epsilon is positive. A boundary mismatch g_theta = u_theta - g can additionally be measured in W^{s,p}(partial Omega), with s = 1 - 1/p, using the Slobodeckij expression shown in the formulas.",
 "math_tags": [
  "pde",
  "functional-analysis",
  "convex-analysis"
 ],
 "ml_areas": [
  "loss",
  "training",
  "regularization"
 ],
 "paper": {
  "arxiv_id": "2608.24205",
  "arxiv_url": "https://arxiv.org/abs/2608.24205",
  "summary_what_math_gives_to_ml": "The paper's transferable contribution is to replace pointwise PDE residual penalties with residuals measured in the dual space W^{-1,p'}, which is better aligned with weak solutions and limited regularity. The dual norm turns residual training into an adversarial weak-form test: the network is penalized according to the worst normalized test function rather than the largest pointwise derivative error. Its boundary treatment also suggests matching traces in the fractional space W^{1-1/p,p}(partial Omega), avoiding an unnecessarily strong pointwise boundary penalty. A practical first transfer is a PINN loss with alternating neural test-function maximization and a Monte Carlo approximation of the fractional boundary seminorm.",
  "title": "Robust training and rigorous error analysis of physics-informed neural networks for the $p$-Laplace equation",
  "year": "2026"
 },
 "ratings": {
  "difficulty": 5,
  "novelty": 7,
  "usefulness": 7
 },
 "solves": [
  "stability",
  "accuracy",
  "sample-efficiency"
 ],
 "title": "Adversarial negative-Sobolev PINN residual",
 "url": "https://synthcore.org/idea/128/adversarial-negative-sobolev-pinn-residual",
 "verification": {
  "peer_reviewed": false,
  "status": "unverified",
  "status_label": "Unverified",
  "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": null,
    "tested": false
   }
  }
 }
}
