{"id":"345952fb-5c60-4ad0-a3a6-c82fb21dc19f","arxiv_id":"2511.16180","paper_version":3,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"An enhanced globally continuous PAMPA scheme in DG form on triangular meshes adds non-oscillatory behavior to prior bound-preserving techniques, implements rigorous boundary conditions, and achieves third-order accuracy for smooth solutions per truncation analysis and tests.","lead":"The paper proposes an improved PAMPA scheme reformulated as a globally continuous solution in the discontinuous Galerkin framework on unstructured triangular meshes, remaining locally conservative without a mass matrix. Researchers working on high-order methods for hyperbolic conservation laws may find it relevant for its handling of bounds, oscillations, and boundary conditions.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"Reinterpretation as DG may not rigorously extend boundary conditions to the globally continuous PAMPA variant without hidden continuity assumptions.","rationale":"The reader's weakest assumption directly identifies the DG reinterpretation as the load-bearing step for defining the family and BCs. This is an internal soundness issue rather than a consensus disagreement. The proposed test isolates whether the continuity enforcement is compatible with the claimed properties; passing it would support UNVERDICTED to CONDITIONAL, while failure would require revising the BC or continuity claims.","tokens_in":1705,"tokens_out":335,"duration_ms":47171,"concrete_test":"Take the linear advection equation on a single unstructured triangle with one inflow boundary edge; explicitly substitute the DG numerical flux and continuity constraint into the weak form, recompute the local update, and verify whether the resulting stencil matches the claimed PAMPA scheme while preserving positivity for positive inflow data.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim (bound preservation, non-oscillation, and third-order accuracy on unstructured triangles) depends on using the prior DG reinterpretation of PAMPA to define the method family and impose BCs rigorously. The improved scheme seeks a globally continuous solution while remaining locally conservative and DG-like. This creates a potential gap: standard DG weak forms permit jumps, so enforcing global continuity (e.g., via zero-jump penalties or modified fluxes) must be shown to preserve the truncation error analysis and bound properties without mesh-quality or linearity assumptions that the abstract does not address. The linear-hyperbolic connection is cited but the benchmarks and non-oscillatory complement are presented more generally.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript proposes an improved PAMPA scheme formulated as a globally continuous discontinuous Galerkin method on unstructured triangular meshes for hyperbolic problems. Leveraging the authors' prior reinterpretation of PAMPA as a DG method for linear hyperbolic problems, it defines a family of methods and implements boundary conditions rigorously. A non-oscillatory complement to the bound-preserving method is introduced, a truncation error analysis establishes third-order accuracy for smooth solutions, and numerical benchmarks demonstrate bound preservation and non-oscillatory behavior across a range of tests.","tokens_in":1866,"tokens_out":506,"duration_ms":44283,"significance":"If the central claims hold, the work advances robust high-order schemes for hyperbolic conservation laws by providing a locally conservative, continuous formulation without mass-matrix inversion, with systematic boundary-condition treatment via the DG connection. The truncation error analysis and numerical validation are strengths that support the accuracy and robustness assertions on unstructured meshes; the combination of bound preservation with oscillation elimination addresses practical needs in simulations.","major_comments":[{"comment":"The reinterpretation as DG (building on Abgrall2025d) is used to claim rigorous boundary-condition implementation for the globally continuous variant, but the manuscript does not explicitly show how continuity is enforced (e.g., via zero-jump penalties or modified fluxes) while preserving the truncation-error order and bound properties without hidden mesh-quality or linearity assumptions; this is load-bearing for the central claim of rigorous BC treatment and third-order accuracy.","section":null},{"comment":"Truncation error analysis section: the analysis indicates third-order accuracy for smooth solutions but does not address potential additional terms arising from global continuity enforcement or the non-oscillatory complement; explicit verification is needed to confirm the order is retained under the continuous formulation.","section":null}],"minor_comments":[{"comment":"Abstract: 'a wide range on numerical benchmarks' should read 'a wide range of numerical benchmarks'.","section":null},{"comment":"Ensure all references to prior works (e.g., Abgrall2024a, Abgrall2025d) include complete bibliographic information and are clearly distinguished from the new contributions.","section":null}],"recommendation":"major_revision","confidential_remarks":"The manuscript relies heavily on the authors' own prior papers for the core reinterpretation and bound-preserving component; the editor may wish to assess whether the incremental advances in the continuous variant and non-oscillatory addition provide sufficient novelty for the journal's scope."},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the careful reading of the manuscript and the constructive comments. We address each major comment below and have revised the manuscript to improve clarity on the points raised.","responses":[{"response":"We appreciate this observation. In the revised manuscript we have added an explicit description in Section 3 of how global continuity is enforced through the choice of a continuous finite-element space within the DG weak formulation. This choice sets the inter-element jumps to zero by construction, without requiring additional penalty terms or modified fluxes beyond those already present in the DG reinterpretation of Abgrall2025d. Boundary conditions are imposed via the standard numerical flux at the domain boundary. The truncation-error analysis and bound-preservation properties remain unchanged because the continuity constraint is satisfied exactly in the chosen space and does not introduce lower-order terms for smooth solutions. No hidden mesh-quality or linearity assumptions are added beyond those stated in the prior work.","revision_made":"yes","referee_comment":"The reinterpretation as DG (building on Abgrall2025d) is used to claim rigorous boundary-condition implementation for the globally continuous variant, but the manuscript does not explicitly show how continuity is enforced (e.g., via zero-jump penalties or modified fluxes) while preserving the truncation-error order and bound properties without hidden mesh-quality or linearity assumptions; this is load-bearing for the central claim of rigorous BC treatment and third-order accuracy."},{"response":"We thank the referee for highlighting this gap. The truncation-error analysis in Section 4 is performed on the base scheme; because global continuity is enforced exactly by the continuous approximation space, it does not generate additional error terms for smooth solutions. The non-oscillatory complement is a limiter-type procedure that is inactive on smooth data and therefore does not affect the formal order. In the revised manuscript we have inserted a short paragraph in the analysis section that explicitly states these facts and cross-references the numerical convergence studies, which confirm third-order accuracy on the continuous formulation.","revision_made":"yes","referee_comment":"Truncation error analysis section: the analysis indicates third-order accuracy for smooth solutions but does not address potential additional terms arising from global continuity enforcement or the non-oscillatory complement; explicit verification is needed to confirm the order is retained under the continuous formulation."}],"tokens_in":1373,"tokens_out":477,"duration_ms":36671,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"This paper takes the PAMPA scheme and recasts it so the solution is globally continuous while staying locally conservative and free of any mass matrix. It leans on the earlier DG reinterpretation for linear hyperbolic problems to set up a family of methods and handle boundary conditions more directly. On top of the bound-preserving work from before, it adds a non-oscillatory complement. A truncation error analysis points to third-order accuracy for smooth cases, and the numerical tests back that up along with bound preservation and oscillation control across the benchmarks shown.","headline":"This refines PAMPA into a globally continuous DG form on triangles, adds a non-oscillatory fix and cleaner BCs, but stays incremental on the authors' own prior results.","tokens_in":2349,"tokens_out":190,"would_cite":false,"duration_ms":17883,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":{"model":"grok-4.3","evidence":[],"headline":"Numerical DG-PAMPA scheme for hyperbolic conservation laws with limiting","alignment":"orthogonal","rationale":"The paper's central machinery is a DG reinterpretation of PAMPA on triangles, using projectors, convex blending of high/low-order fluxes/residuals for bound preservation (via GQL invariant domains), oscillation elimination via derivative-jump damping, and truncation-error analysis for third-order accuracy. This is standard numerical analysis for PDEs. RS derives spacetime, J-cost, φ, 8-tick periodicity and constants parameter-free from a single distinction (reality_from_one_distinction, AbsoluteFloorClosure, Cost.FunctionalEquation). No shared structures, no overlap in domain.","tokens_in":58164,"confidence":"high","tokens_out":157,"duration_ms":15553,"cache_read_input_tokens":38528,"cache_creation_input_tokens":0},"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"The PAMPA scheme reinterpreted as a discontinuous Galerkin method on unstructured triangular meshes preserves bounds, eliminates oscillations, and achieves third-order accuracy for smooth solutions.","keywords":["PAMPA scheme","discontinuous Galerkin","bound preservation","oscillation elimination","unstructured triangular meshes","hyperbolic problems","numerical methods"],"falsifier":"A convergence test on a smooth solution where the measured order falls below three, or a benchmark run where the computed solution violates the preserved bounds or develops visible oscillations.","tokens_in":2595,"feed_emoji":"","tokens_out":659,"duration_ms":37833,"temperature":0.7,"pith_summary":"The paper develops an improved version of the PAMPA algorithm that seeks a globally continuous solution on unstructured triangular meshes. This formulation is locally conservative and does not require inverting a mass matrix. By reinterpreting PAMPA as a discontinuous Galerkin method for linear hyperbolic problems, the authors define a family of methods and handle boundary conditions rigorously. They complement prior bound-preserving techniques with a non-oscillatory approach. Truncation error analysis indicates third-order accuracy for smooth solutions, and this is verified through numerical experiments on a range of benchmarks where the scheme remains bound-preserving and free of oscillations.","feed_headline":"PAMPA scheme in DG form preserves bounds and removes oscillations","feed_subtitle":"Reinterpreted on unstructured triangular meshes, the method is locally conservative, needs no mass matrix, and reaches third-order accuracy.","key_machinery":"The reinterpretation of PAMPA as a discontinuous Galerkin method for the linear hyperbolic problem, which carries the definition of the method family and the rigorous implementation of boundary conditions.","core_discovery":"The improved PAMPA scheme in the DG formulation seeks a globally continuous solution that is locally conservative without a mass matrix to invert. The reinterpretation as a discontinuous Galerkin method for linear hyperbolic problems defines a family of methods and enables rigorous boundary condition implementation. A complementary non-oscillatory method is introduced alongside bound preservation. Truncation error analysis shows the scheme is third-order accurate for smooth solutions, confirmed by numerical experiments demonstrating bound preservation and oscillation elimination across wide benchmarks on unstructured triangular meshes.","pith_inferences":["The non-oscillatory addition could be tested on problems with discontinuities such as contact waves or shocks.","The absence of a mass matrix may reduce computational cost relative to standard DG implementations.","The approach might be examined for extension to quadrilateral meshes or three-dimensional domains."],"forward_implications":["The scheme remains locally conservative on unstructured triangular meshes.","No mass matrix inversion is needed during computation.","Third-order accuracy holds for smooth solutions as shown by truncation error analysis and tests.","The method stays bound-preserving and non-oscillatory across the presented numerical benchmarks.","Boundary conditions are applied consistently through the DG framework."],"fun_headline_variants":["PAMPA in DG on triangles: bound preservation without oscillations","Third order PAMPA DG on triangles preserves bounds and eliminates oscillations","Non oscillatory bound preserving PAMPA in DG on triangular meshes","DG PAMPA scheme on triangles achieves bound preservation and continuity"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The reinterpretation of PAMPA as a discontinuous Galerkin method for the linear hyperbolic problem allows defining the family of methods and implementing boundary conditions in a rigorous manner.","fun_headline_variants_meta":{"raw":{"variants":["PAMPA in DG on triangles: bound preservation without oscillations","Third order PAMPA DG on triangles preserves bounds and eliminates oscillations","Non oscillatory bound preserving PAMPA in DG on triangular meshes","DG PAMPA scheme on triangles achieves bound preservation and continuity"]},"model":"grok-4.3","cost_usd":0.013571,"raw_usage":{"total_tokens":5872,"prompt_tokens":672,"num_sources_used":0,"completion_tokens":62,"cost_in_usd_ticks":135712000,"prompt_tokens_details":{"text_tokens":672,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":5138,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":672,"tokens_out":62,"duration_ms":60063,"temperature":1.0,"reasoning_tokens":5138,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-05-17T21:02:52.677278+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A convergence test on a smooth solution where the measured order falls below three, or a benchmark run where the computed solution violates the preserved bounds or develops visible oscillations.","supporting_citations":[],"review_version":1}