{"id":"db5eccf6-b579-4dde-a09e-07fecf1b481b","arxiv_id":"1908.01796","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"An equilibration-based cut-point discretization allows ordinary planar Cartesian grid finite differences to solve PDEs on closed level-set surfaces with observed second-order accuracy.","lead":"A new numerical method solves partial differential equations on curved closed surfaces by mapping the surface onto flat Cartesian grids and using standard finite differences. Tests on spheres and other shapes show second-order accuracy for diffusion, advection, and shallow water flow.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unproven equilibration consistency: O(h^3) secondary-point interpolation could degrade Δ_h E_h to O(h), and no truncation test isolates this risk.","rationale":"The reader correctly identified the equilibration operator E_h as the load-bearing assumption, but the verdict should be conditional rather than an unconditional accept: the numerical tests show second-order behavior in selected examples, yet no direct consistency check of Δ_h E_h rules out an O(h) interpolation error. The paper itself acknowledges the difficulty of proving convergence (Section 5), and the only positivity proof is for curves in R^2 (Appendix A), so the general R^3 claim lacks analytic support. The proposed truncation test would settle whether the interpolation error cancels. Until that check is performed or a convergence proof is supplied, the general second-order claim should be treated as plausible but not fully established; hence CONDITIONAL rather than UNCHANGED.","tokens_in":19140,"tokens_out":14803,"duration_ms":171637,"concrete_test":"Compute the exact truncation error of the reduced operator: for a smooth non-polynomial function u (e.g., u = Y_3^2 on the unit sphere, or u = cos(x−y+z) on the Cassini oval), form u^p as the restriction of u to primary points, compute E_h u^p via (2.4), evaluate Δ_h(E_h u^p) at every primary point using the divergence-form stencil (3.11)/(3.15), and compare with the exact Δu at those points. Measure max and L2 norms for N = 80, 160, 320, 640. If the rate is O(h), the equilibration interpolation is a first-order bottleneck and the second-order claim is not general; if the rate is O(h^2), the concern does not land. Repeating with a higher-frequency smooth u, such as u = sin(8x), would expose any missing cancellation more sharply.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central claim of second-order accuracy for general surface PDEs rests on the composition Δ_h E_h (eqs. 3.23–3.25) being a consistent O(h^2) approximation of the Laplace–Beltrami operator. The equilibration map E_h fixes secondary values by quadratic interpolation (2.1)–(2.4), so for a smooth function the pointwise interpolation error is O(h^3). A primary-point finite-difference stencil divides by h^2; if a neighboring cut point is secondary, that O(h^3) error can enter the discrete operator as O(h) unless cancellations occur across the stencil or the global equilibration system enforces them. The author explicitly states that proving convergence 'seems difficult' (Section 5), and the only resolvent-positivity proof is for curves in R^2 (Appendix A), giving no analytical protection in R^3. The convergence tables (e.g., Table 3.1) are consistent with second order, but they do not separate the equilibration contribution from the primary stencil truncation; the shallow-water max-norm data (Tables 4.2–4.3) are also less clean. The most load-bearing assumption is therefore that this interpolation error cancels or is otherwise benign, and the paper provides no direct consistency measurement of Δ_h E_h.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a finite-difference method for solving PDEs on smooth closed surfaces, built on the Ying--Wang cut-point representation of a level-set surface. Unknowns are carried at primary cut points, values at secondary cut points are obtained by solving a quadratic-interpolation equilibration system, and differential operators are written in each of the three coordinate planes after projection. The method is applied to surface diffusion and Poisson problems with the Laplace--Beltrami operator, to a linear advection equation with an exact solution on the sphere, and to the shallow water equations on a sphere using a Lax--Wendroff-type scheme with artificial viscosity. The paper reports second-order convergence in extensive numerical tables, computes eigenvalues of the discrete spherical Laplacian, proves a resolvent-positivity result for the analogous discretization on curves in two dimensions, and explicitly states that a full convergence proof for surfaces is not yet available.","tokens_in":19448,"tokens_out":11716,"duration_ms":128024,"significance":"If the observed second-order accuracy holds in the advertised generality, the method is a simple and attractive alternative to triangulated, closest-point, or RBF-based surface PDE solvers for smooth closed surfaces given by level sets: it reduces surface computations to standard finite differences on planar Cartesian grids without boundary conditions. The central discretization is parameter-free in the sense that the finite-difference coefficients are obtained from Taylor expansions rather than from data fitting, and the paper provides unusually complete numerical evidence, including convergence tables for diffusion, advection, and shallow water, eigenvalue tests for the discrete Laplace--Beltrami operator on the sphere, and a high-order quadrature rule. The author's explicit acknowledgment of the missing convergence proof and the candid reporting of less regular shallow-water maximum-norm results are also strengths of the presentation. The main reservation is that the abstract's phrasing of the accuracy claim is broader than what the shallow-water tables support in the maximum norm.","major_comments":[],"minor_comments":[{"comment":"The paper states in Section 5 that proving convergence 'seems difficult,' and the consistency analysis in Section 3 is performed for the primary stencil assuming exact values at all stencil points. Since the O(h^3) quadratic-interpolation error at secondary points could in principle enter a second-difference stencil as O(h), it would be informative to include a direct numerical consistency test of the composed operator Delta_h E_h (for example, apply Delta_h E_h to a known smooth function on a nontrivial surface and tabulate the discrete Laplace--Beltrami error as a function of h). This would isolate the equilibration contribution from the primary stencil truncation and would further substantiate the observed second-order convergence.","section":"Sections 3 and 5, Eqs. (2.1)-(2.4), (3.23)-(3.25)"},{"comment":"The abstract's claim of 'second order accuracy in examples with a variety of equations, including ... shallow water equations' is broader than the reported data. In Table 4.2 with nu=1, the maximum-norm errors for Phi v do not decrease between N=320 and N=640 (e.g., day 1 values 2.26e-3 and 2.17e-3), and the nu=0.5 results in Table 4.3 are non-monotone, with the day-1 maximum error for Phi v increasing from 4.48e-3 at N=320 to 1.37e-2 at N=640. The L2 errors do show roughly second-order behavior, so I suggest qualifying the accuracy statement for the shallow water example as L2 accuracy and restating the feasibility-only goal for nonlinear hyperbolic problems.","section":"Abstract, Section 4, Tables 4.2 and 4.3"},{"comment":"In the BDF2 L2 column for N=320, the entry printed as '2.33-5' should presumably be '2.33e-5'; please correct the typographical formatting.","section":"Table 3.2, Cassini oval row"},{"comment":"The artificial-viscosity parameter nu takes the values 0.5 and 1 in the shallow water experiments, and the convergence behavior depends noticeable on this choice. A brief sentence on how these values were selected, and on the sensitivity of the reported rates to nu, would improve reproducibility.","section":"Section 4, Eqs. (4.33)-(4.34)"}],"recommendation":"minor_revision","confidential_remarks":"The paper is well written and the numerical evidence for the linear problems is strong. I do not regard the absence of a convergence proof as a blocker, since the paper explicitly identifies this as an open challenge and the convergence tables already exercise the full composition Delta_h E_h. The main revision request is to align the abstract and conclusions with the actual shallow-water convergence behavior and, if convenient, to add the direct consistency test mentioned in the minor comments."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version: the paper does something genuinely new — it takes Ying-Wang's cut-point surface discretization, built for boundary integral equations, and turns it into a general finite-difference toolbox for surface PDEs. The Laplace–Beltrami, gradient, divergence, and shallow-water formulations are derived cleanly in three coordinate patches, and the numerical evidence is consistently second order across diffusion, Poisson, advection, and the shallow-water equations. The computations look honest: no fitted parameters in the central claim, Taylor expansions are given for the discrete operators, and the convergence tables are detailed across multiple surfaces including a nonconvex Cassini oval. The author explicitly states that a convergence proof 'seems difficult' and gives a resolvent-positivity proof only for curves in R2. That is the soft spot, and it is real: the equilibration step uses quadratic interpolation at secondary cut points, and a naive error count suggests the O(h^3) interpolation error could enter the discrete Laplacian as O(h) after division by h^2. The paper does not isolate and measure this contribution directly, so the observed second-order accuracy relies on cancellations that are not analyzed. This is a genuine gap, but not a fatal one: the empirical convergence is consistent across many runs, and the method is simple enough that others can and probably will test it. The shallow-water max-norm data are less clean — for ν=.5, the max error in Φv stops decreasing between N=320 and N=640 — but the L2 errors continue to improve, and the authors themselves present this as a feasibility study, not a mature scheme. Who is this for? Anyone working on surface PDEs who wants an alternative to closest-point, level-set extension, or point-cloud methods, especially if they already use Cartesian grids and want something simpler. It is not a finished theory paper — the missing convergence proof matters — but it is a solid, clearly written methods paper that deserves a serious referee and, I think, publication after revision. My recommendation: send it to review. It will be read, used, and probably extended.","headline":"A solid, honest methods paper: the cut-point equilibration idea is adapted to surface PDEs with convincing second-order numerics, though the load-bearing interpolation step remains unproved.","tokens_in":19926,"tokens_out":1539,"would_cite":true,"duration_ms":16249,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M06","65M50","58J35","58J45","35Q86"],"pacs":[],"model":"deepseek-v4-flash","headline":"Representing a closed surface by cut points on a planar grid, with values at secondary points fixed by quadratic interpolation, lets ordinary finite differences solve Laplace-Beltrami diffusion, advection, and shallow-water equations on…","keywords":["Cartesian grids","surface PDEs","Laplace-Beltrami operator","equilibration","level set method","shallow water equations","finite difference methods"],"falsifier":"Run the diffusion test of Section 3 on a Cassini oval with a near-self-intersecting value such as b/a = 1.02, where the surface has large curvature; if halving the grid size h no longer reduces the relative L2 error by roughly a factor of four, the claimed second-order accuracy fails.","tokens_in":18965,"feed_emoji":"🧮","tokens_out":7057,"duration_ms":62344,"temperature":0.7,"pith_summary":"Most PDE solvers on curved surfaces need a mesh or an extension of the equation into a neighborhood of the surface. This paper proposes a simpler route: if the surface is given as a level set, it is represented by the points where the surface crosses the intervals of a regular three-dimensional Cartesian grid. Some of these cut points carry independent unknowns (primary points); the rest (secondary points) get their values from the primary ones through a quadratic-interpolation 'equilibration' step. The paper shows that with this representation, the Laplace-Beltrami operator, surface gradients, divergence, advection, and the shallow-water equations can be discretized by ordinary finite differences in the three coordinate planes, and the observed accuracy is second order on spheres, ellipsoids, and a nonconvex Cassini oval.","feed_headline":"Flat grids solve curved-surface PDEs at second order","feed_subtitle":"Cut points plus quadratic interpolation let ordinary finite differences handle spheres, ellipsoids, and ovals.","key_machinery":"The machinery is the equilibration operator E_h. Given values at primary cut points, E_h fills in values at secondary cut points by solving the system u_s = Π_sp u_p + Π_ss u_s, whose rows are the quadratic interpolation condition (2.1) along the surface; since I − Π_ss is strictly diagonally dominant, the system is invertible. Composing E_h with the finite-difference discretization of the surface differential operator — e.g., the nine-point Laplace-Beltrami stencil Δ_h in (3.11) or (3.15), or the coordinate-plane formulas for surface gradient and divergence — produces the reduced operator Δ_h E_h on primary points. The equilibration is what lets each primary point use the coordinate plane in which the surface normal is nearly aligned with a coordinate axis, and it is what transfers updated values back to secondary points between time steps.","core_discovery":"The central claim is that a single, fixed discretization of any smooth closed level-set surface reduces surface PDEs to standard planar finite-difference calculations. The surface is sampled at cut points where it meets grid intervals; primary points are those closest to grid points, and secondary points are tied to primary points by solving a small interpolation system. The discrete Laplace-Beltrami operator is built as a nine-point stencil in the appropriate coordinate plane, and the equilibration operator E_h passes values from primary to secondary points so that stencils can be applied seamlessly as the coordinate system changes across the surface. The paper demonstrates second-order convergence for diffusion and Poisson problems, for a linear advection equation with a known solution, and for the shallow-water equations on a sphere using a MacCormack Lax-Wendroff scheme with artificial viscosity; it also computes low eigenvalues of the discrete spherical Laplacian and shows they converge to the exact values at second order. The author is careful to state that convergence is observed but not proved, and a proof is left for the future.","pith_inferences":["A natural next step, not taken in the paper, is to use higher-order interpolation in the equilibration step; if the accuracy barrier is really the quadratic interpolation, this would open the way to fourth-order surface PDE solvers with the same cut-point structure.","The paper's proof of resolvent positivity covers only closed curves in the plane; if a similar argument can be devised for surfaces in R^3, it would supply the missing maximum-principle stability theory for implicit time stepping.","The surface-discretization machinery could plausibly be combined with a semi-Lagrangian time step by interpolating from neighboring primary points, a possibility the paper mentions but does not develop.","The restriction to smooth level-set surfaces is explicit; applying the method to point clouds or to surfaces with rapidly varying curvature would require additional machinery, and those settings are outside the paper's claim."],"forward_implications":["Reaction-diffusion systems on surfaces are a direct extension, since the paper's Laplace-Beltrami discretization can be combined with pointwise reaction terms (Section 5).","Poisson's equation on a closed surface can be solved by augmenting the reduced Laplace-Beltrami operator with one scalar constraint, and the paper shows second-order convergence for this formulation on the sphere.","The method's coordinate-plane formulas for surface gradient, divergence, and material derivative mean that planar finite-volume or finite-difference schemes can be ported to smooth closed surfaces with only the equilibration step added.","Surface integrals can be computed directly from the cut points with a high-order quadrature rule, allowing conserved quantities such as mass and energy to be monitored in long-time simulations without building a separate surface mesh."],"supporting_citations":[{"why":"Introduces the cut-point surface discretization with primary/secondary points and the equilibration interpolation that the paper builds on.","marker":"[32]"},{"why":"Provides the standard shallow-water test problem on the sphere used to validate the method's accuracy.","marker":"[29]"},{"why":"Supplies the cut-point quadrature rule on the sphere that the paper adapts for surface integrals.","marker":"[31]"},{"why":"Presents the simple quadrature method for singular and near-singular integrals on closed surfaces that underlies Appendix B.","marker":"[2]"},{"why":"Supplies the M-matrix theory used in Appendix A to prove positivity of the resolvent for the discrete Laplace-Beltrami operator on closed curves.","marker":"[15]"}],"fun_headline_variants":["Curved surfaces, flat grids: second-order PDE solver","Level-set surface PDEs solved via planar Cartesian stencils","Cut-point equilibration brings surface PDEs down to flat grid","No surface mesh needed: planar grids hit second-order accuracy","Projection-based discretization handles spheres, ellipsoids, ovals"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The scheme's accuracy depends on the quadratic interpolation that fills in secondary-point values being accurate enough not to spoil the second-order finite-difference stencils it is composed with, and this is verified numerically but not proven for surfaces in three dimensions.","fun_headline_variants_meta":{"raw":{"variants":["Curved surfaces, flat grids: second-order PDE solver","Level-set surface PDEs solved via planar Cartesian stencils","Cut-point equilibration brings surface PDEs down to flat grid","No surface mesh needed: planar grids hit second-order accuracy","Projection-based discretization handles spheres, ellipsoids, ovals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000214,"raw_usage":{"total_tokens":1425,"prompt_tokens":946,"completion_tokens":479,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":562,"completion_tokens_details":{"reasoning_tokens":393}},"tokens_in":562,"tokens_out":479,"duration_ms":5128,"temperature":1.0,"reasoning_tokens":393,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:03:12.017863+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the diffusion test of Section 3 on a Cassini oval with a near-self-intersecting value such as b/a = 1.02, where the surface has large curvature; if halving the grid size h no longer reduces the relative L2 error by roughly a factor of four, the claimed second-order accuracy fails.","supporting_citations":[{"cited_title":"Ying and W.-C","cited_arxiv_id":null,"evidence_quote":"Introduces the cut-point surface discretization with primary/secondary points and the equilibration interpolation that the paper builds on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the standard shallow-water test problem on the sphere used to validate the method's accuracy."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Presents the simple quadrature method for singular and near-singular integrals on closed surfaces that underlies Appendix B."},{"cited_title":"Hackbusch, Elliptic Differential Equations: Theory and Numerical Trea tment, sec- ond ed., Springer, Berlin, 2017","cited_arxiv_id":null,"evidence_quote":"Supplies the M-matrix theory used in Appendix A to prove positivity of the resolvent for the discrete Laplace-Beltrami operator on closed curves."}],"review_version":1}