{"id":"9d0df9d9-fd20-476e-ac98-764459c0e90d","arxiv_id":"2411.15565","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":1,"one_line_summary":"The paper derives a stability bound for a Galerkin/least-squares method with L2-optimal test functions on B-splines and reports numerical comparisons with SUPG, but benchmark and claim inconsistencies undermine the stated conclusions.","lead":"An isogeometric finite element paper derives and tests a Galerkin/least-squares stabilization built from L2-optimal test functions for advection-dominated diffusion. Its central practical claim is that this parameter-free method matches a standard tuned stabilization, but the paper's own error tables contradict parts of its abstract.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The first model problem's benchmark is invalid: the exact solution (2) does not satisfy the stated homogeneous PDE, and the weak forms omit the y-advection term for beta=(1,1), so the reported errors are meaningless.","rationale":"The paper's main contribution is a numerical comparison, and the first model problem is one of two benchmarks. Both the exact solution and the weak forms are wrong for this problem: the exact solution does not satisfy the homogeneous PDE, and the weak forms drop the y-advection term, so the implementation solves a different PDE. This is a correctness flaw at the foundation of the numerical evidence, not a mere exposition problem. The stability proof's condition epsilon <= h/(2C^2) also fails for the reported experiments (e.g., epsilon=0.003, h=0.1), but that is a secondary concern; even if the proof were extended, the numerical comparison would still rest on an invalid benchmark. The reader's verdict of REJECT is appropriate because the manuscript's own evidence does not support the central claim for the first model problem, and the issue is fixable only by redoing the formulation and the experiments.","tokens_in":12673,"tokens_out":5953,"duration_ms":51527,"concrete_test":"Re-run the first model problem with the correct strong form u_x + u_y - eps*Laplacian(u) = 0 and with the y-advection term included in the weak forms, using a manufactured nontrivial solution (e.g., u_mfg = sin(pi*x)*sin(pi*y) with the corresponding right-hand side). Compare GLS, SUPG, and optimal-test-function solutions on the same ten-element uniform mesh and recompute the relative L2 and H1 errors. If the ranking changes or if the new errors differ substantially from Table 11, the central claim is unsupported as written.","verdict_should_be":"REJECT","load_bearing_attack":"The central claim that the Galerkin/least-squares method is a parameter-free stabilizer competitive with SUPG rests on the numerical comparisons in Section 3. For the first model problem, beta=(1,1), the stated PDE is u_x + u_y - eps*Laplacian(u) = 0, but every weak form used (Eqs. (4), (5), (12), and the GLS derivation in Section 2.4) contains only d/dx and the diffusion terms, omitting d/dy. Thus the method actually solves a different boundary value problem (u_x - eps*Laplacian(u) = 0) on the unit square with zero Dirichlet data, which has only the trivial solution. The 'exact solution' in Eq. (2) is a product of two 1D factors each solving -eps*f'' + f' = 1; applying the stated 2D operator gives 2*f(x)*f(y), not zero, so it is not a solution to the stated homogeneous PDE. Consequently, the relative L2 and H1 errors in Tables 1, 3, 5, 7, 9, and 11 for the first problem do not measure error against a true solution of the problem being solved. This invalidates the direct evidence for the paper's headline comparison, because one of the two model problems is faulty at the level of both the continuous problem and the discrete implementation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies stabilization of isogeometric finite element methods for advection-diffusion problems. It derives a least-squares method from Petrov-Galerkin optimal testing in the L2 norm, combines it with the standard Galerkin term to obtain a Galerkin/least-squares method, proves a coercivity bound for B-spline spaces, and compares the resulting schemes with SUPG on two model problems. The theoretical stability analysis in Section 4 is a self-contained argument based on an inverse inequality for B-splines. However, the numerical part is seriously compromised: the first benchmark problem is mis-specified, with an exact solution that does not satisfy the stated PDE, weak forms that omit the y-advection term for beta=(1,1), and an abstract and conclusions that contradict the reported error tables.","tokens_in":12919,"tokens_out":8107,"duration_ms":67640,"significance":"If the stability analysis and the numerical comparisons were correct, the paper would offer a parameter-free stabilized isogeometric variant with a rigorous coercivity bound. The coercivity proof for the Galerkin/least-squares form in Section 4 is a useful and largely self-contained contribution. However, the paper's central comparative claim—that the proposed Galerkin/least-squares method is competitive with or better than SUPG without tuning—rests on Tables 11 and 12, and for the first model problem those tables are obtained from an incorrectly posed benchmark. The second model problem is correctly posed, but it alone does not establish the headline conclusions, which are also in direct contradiction with the table values.","major_comments":[{"comment":"The exact solution in Eq. (2) is not a solution of the stated homogeneous PDE (1) with beta=(1,1). For f(t)=t+(e^{t/epsilon}-1)/(1-e^{1/epsilon}) one has -epsilon f''+f'=1, so for u=f(x)f(y) the left-hand side of (1) equals f(x)+f(y), which is not zero. With zero Dirichlet data, the unique solution of the homogeneous equation (1) is u identically zero. Consequently, the relative errors reported in Tables 1, 3, 5, 7, 9, and 11 for the first model problem are computed against a function that is not the exact solution of the problem being solved, and the comparison for that benchmark is meaningless.","section":"Section 2.1.1, Eq. (2)"},{"comment":"The weak forms used throughout the paper contain only the x-advection term and the two diffusion terms; the y-advection term is absent. For the first model problem with beta=(1,1), this means the implemented method solves a different boundary value problem, namely u_x - epsilon(u_xx+u_yy)=0 with zero Dirichlet data, whose only solution is zero. The non-zero computed solutions reported in the tables indicate that the actual right-hand side or boundary data used in the implementation differ from the equations stated in the paper, so the experiments are not reproducible from the written formulation.","section":"Section 2.2, Eqs. (4), (5), (12), and (14)-(16)"},{"comment":"The abstract states that least-squares (Petrov-Galerkin with L2-optimal test functions) outperforms the other stabilization methods for small Peclet numbers, and the conclusion states that the optimal test functions behave better than SUPG for epsilon=0.1 and 0.01. Tables 11 and 12 show the opposite: for example, Table 11 at epsilon=0.01 gives L2 errors 86.70 (PG-optimal), 28.20 (GLS), and 33.71 (SUPG); Table 12 at epsilon=0.01 gives 52.87, 17.10, and 20.88, respectively. The headline claims are thus not supported by the paper's own data.","section":"Abstract and Section 5 vs. Tables 11-12"}],"minor_comments":[{"comment":"Equation (1) writes 'beta_y d u/d x' for the second advection term; presumably this should be 'beta_y d u/d y'.","section":"Section 2.1, Eq. (1)"},{"comment":"The caption reads 'epsilon in {0.1, 0.001, 0.0003}' while the text and Table 3 use {0.1, 0.01, 0.003}.","section":"Figure 7 caption"},{"comment":"Table 8 includes a row for epsilon=0.1 for the second problem, although Section 2.1.2 defines the second problem for epsilon in {0.01, 0.001, 0.0001}.","section":"Section 2.4, Table 8"},{"comment":"The continuity bound for the term h^{-1}(beta·grad u_h, v_h) is K|beta|/epsilon, which grows as epsilon tends to zero; the phrase 'provided sufficiently small epsilon > 0' is therefore misleading, since the inequality itself holds for any fixed epsilon with a constant that depends on epsilon.","section":"Section 4, Lemma 2"}],"recommendation":"reject","confidential_remarks":"The stability analysis in Section 4 is a reasonable and largely independent contribution, but the numerical comparison—the main evidence for the proposed method—is invalidated by the first model problem's mis-specification, and the abstract and conclusions directly contradict the reported tables. These are load-bearing issues that would require redoing the experiments and reinterpreting the results, so I cannot recommend revision within the current scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The first thing you should know: the first model problem doesn't hold together. The stated PDE for beta=(1,1) is u_x + u_y - eps(u_xx+u_yy)=0, but every weak form used in the paper, from (4) through (16), keeps only the x-advection term. The method actually solves u_x - eps Δu = 0 with zero Dirichlet data, whose only solution is zero. The \"exact solution\" in (2) is a product of two 1D factors, each solving -eps f'' + f' = 1; applying the stated 2D operator to that product does not give zero. So the errors reported in Tables 1,3,5,7,9,11 are relative to a function that is not a solution. That's a load-bearing flaw, and the stress-test note is right about it.\n\nNow the credit. The stability analysis in Section 4 is a real piece of work. Theorem 1 proves coercivity of a GLS-type bilinear form on B-spline spaces under the condition eps ≤ h/(2C^2), where C comes from a standard inverse inequality. The proof is short, and it holds up on inspection. The explicit formula (11) for the L2-optimal test functions is clean, and the paper correctly acknowledges that this is the least-squares method in disguise. The second model problem, Eriksson-Johnson, is standard, and its weak forms are consistent; the results for it in Tables 6,8,10,12 are probably meaningful.\n\nThe soft spots beyond the broken benchmark: the conclusions contradict the paper's own tables. The abstract and Section 5 say L2-optimal test functions outperform SUPG for eps = 0.1 and 0.01, but in every row of Tables 11 and 12 the PG-optimal method has larger L2 and H1 errors than SUPG. That is not a matter of interpretation. The claim that GLS is competitive with SUPG is actually supported by the tables, which is useful, but the headline comparison is not. Also minor: the \"exact solution\" (2) is numerically unstable for small eps, and the eps=0.003 experiments violate the condition in Theorem 1.\n\nWho gets value from this: someone working on stabilized isogeometric methods might find the coercivity bound and the formulation comparison worth a look, but only after the first problem is re-specified and the experiments rerun. As it stands, I would not cite the numerical results. A serious editor should send this to peer review, because the stability proof is formally grounded and the flaw is identifiable and fixable. But the current version should be rejected, with an invitation to resubmit after a major revision.","headline":"A broken first model problem invalidates the headline numerical comparison; the stability proof is a plausible separate contribution.","tokens_in":13491,"tokens_out":5331,"would_cite":false,"duration_ms":44124,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65N30","65N12","65D07"],"pacs":[],"model":"deepseek-v4-flash","headline":"By combining standard Galerkin testing with L2-optimal test functions, this paper derives a stabilization for isogeometric B-spline discretizations that needs no tuning parameters and matches SUPG on advection-dominated diffusion.","keywords":["isogeometric analysis","B-spline basis functions","advection-diffusion equation","Petrov-Galerkin optimal test functions","least-squares finite element method","Galerkin/least-squares stabilization","SUPG","coercivity"],"falsifier":"Substitute the exact solution in Eq. (2) into the PDE in Eq. (1) with $\\beta=(1,1)$, including the y-advection term, and evaluate the residual; if the residual is not zero, the benchmark and all error tables do not measure the model problem the paper intends to solve.","tokens_in":48,"feed_emoji":"🌊","tokens_out":11885,"duration_ms":219090,"temperature":0.7,"pith_summary":"Advection-dominated diffusion is a standard test case where plain Galerkin finite element methods produce spurious oscillations unless the mesh is carefully adapted. This paper argues that a stabilization scheme for isogeometric (B-spline) discretizations can be built directly from residual minimization: using the $L_2$ norm, the optimal test function for each trial basis function is just the advection-diffusion operator applied to that basis function, which turns the method into a least-squares formulation. Combining that with an ordinary Galerkin term yields a Galerkin/least-squares method whose coercivity is proved for B-spline spaces under a mesh-size-diffusion relation, with no tuning parameter. The paper's numerical comparison on two benchmark problems with coarse uniform meshes indicates this combined scheme is as accurate as SUPG and sometimes better, especially for small diffusion. If the claim is right, stabilized simulations on B-spline meshes can be set up without hand-tuning stabilization constants.","feed_headline":"No-tune stabilized isogeometric method rivals SUPG","feed_subtitle":"Combining Galerkin and least-squares test functions stabilizes B-spline solutions on coarse uniform meshes.","key_machinery":"The central object is the $L_2$-optimal test function. For each trial basis function $e_i$, minimizing the residual in $L_2$ yields the test function $\\hat{e}_i = \\partial e_i/\\partial x - \\varepsilon(\\partial^2 e_i/\\partial x^2 + \\partial^2 e_i/\\partial y^2)$, the image of the trial function under the advection-diffusion operator; this explicit formula turns the Petrov-Galerkin method into a least-squares method without solving auxiliary problems. The stabilized Galerkin/least-squares variant adds a $1/h$-weighted Galerkin term to this optimal testing, and the coercivity proof rests on an inverse inequality for B-splines, $\\|\\Delta v_h\\|_{L_2} \\le C h^{-1}\\|\\nabla v_h\\|_{L_2}$, whose constant $C$ enters the stability condition $\\varepsilon \\le h/(2C^2)$.","core_discovery":"The paper's central claim is that $L_2$-norm residual minimization provides an explicit, analytic family of optimal test functions for advection-diffusion problems, and that combining these with standard Galerkin test functions produces a stabilized isogeometric method with proven stability. In the $L_2$ setting, the optimal test function for a B-spline basis function $e_i$ is the image of $e_i$ under the differential operator, $\\partial e_i/\\partial x - \\varepsilon(\\partial^2 e_i/\\partial x^2 + \\partial^2 e_i/\\partial y^2)$; the resulting Petrov-Galerkin method is exactly a least-squares formulation. Adding a mesh-weighted Galerkin term with a $1/h$ factor gives the Galerkin/least-squares scheme, and the paper proves a coercivity estimate $b(v_h, v_h) \\ge \\frac{1}{2}((\\varepsilon/h)\\|\\nabla v_h\\|^2 + \\|\\beta\\cdot\\nabla v_h\\|^2)$ for $\\varepsilon \\le h/(2C^2)$, where $C$ comes from a B-spline inverse inequality. The numerical experiments on two advection-diffusion benchmarks with ten uniform elements show the combined scheme matches or beats SUPG in $L_2$ accuracy for the smallest diffusion coefficients tested, and it requires no stabilization parameter to tune.","pith_inferences":["If the coercivity proof transfers to other B-spline spaces satisfying the same inverse inequality, the no-tune stabilization could extend to locally refined, trimmed, or higher-dimensional isogeometric meshes without new parameter choices.","Because the $L_2$-optimal test function is known in closed form, the same operator-action identity could be reused as a cheap residual or error indicator for adaptive mesh refinement; the paper does not explore this.","The mesh-dependent $1/h$ coupling suggests the stabilization strength is encoded in the discretization, which may make the method easier to port to time-dependent or nonlinear advection-diffusion problems where SUPG constants are usually re-tuned."],"forward_implications":["On a uniform mesh of ten elements, the combined Galerkin/least-squares formulation gives the lowest $L_2$ errors among all tested methods for the smallest diffusion parameters in both benchmark problems.","The pure $L_2$-optimal-test (least-squares) method is unconditionally stable but overly diffusive on coarse meshes for very small $\\varepsilon$, so the combination, not the pure form, is the practically useful scheme.","The stability theorem guarantees coercivity whenever $\\varepsilon \\le h/(2C^2)$, so the mesh-to-diffusion ratio, not a user-tuned constant, determines stability.","The paper's comparisons indicate the combined method requires no tuning of stabilization parameters, unlike SUPG, which needs the parameter $\\tau$ to be chosen."],"supporting_citations":[{"why":"supplies the two benchmark advection-diffusion problems and the exact solutions used in every error table.","marker":"[2]"},{"why":"supplies the inverse inequality for B-splines from which the coercivity constant C and the condition ε ≤ h/(2C²) are derived.","marker":"[3]"},{"why":"defines the SUPG method that serves as the main baseline in the numerical comparisons.","marker":"[4]"},{"why":"provides the least-squares finite element framework and its stability properties that the L2-optimal-test formulation builds on.","marker":"[10]"},{"why":"introduces the Galerkin/least-squares formulation that the combined test-function scheme is designed to recover.","marker":"[15]"},{"why":"supplies the stability theory used to conclude well-posedness of the discrete formulation from continuity and coercivity.","marker":"[16]"}],"fun_headline_variants":["L2 residual minimization yields tune-free stabilized isogeometric FEM","No-parameter stabilization for isogeometric methods beats SUPG at low Pe","Optimal L2 test functions stabilize isogeometric FEM without tuning","Tune-free isogeometric stabilization from residual minimization","Provably stable isogeometric method from L2 residual minimization"],"cache_read_input_tokens":15616,"weakest_assumption_plain":"The load-bearing premise is that the benchmark exact solutions stated in the paper really do solve the model problems as written, so the reported relative errors compare each method against the correct target.","fun_headline_variants_meta":{"raw":{"variants":["L2 residual minimization yields tune-free stabilized isogeometric FEM","No-parameter stabilization for isogeometric methods beats SUPG at low Pe","Optimal L2 test functions stabilize isogeometric FEM without tuning","Tune-free isogeometric stabilization from residual minimization","Provably stable isogeometric method from L2 residual minimization"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001107,"raw_usage":{"total_tokens":4665,"prompt_tokens":1049,"completion_tokens":3616,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":665,"completion_tokens_details":{"reasoning_tokens":3528}},"tokens_in":665,"tokens_out":3616,"duration_ms":23329,"temperature":1.0,"reasoning_tokens":3528,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:10:10.184014+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Substitute the exact solution in Eq. (2) into the PDE in Eq. (1) with $\\beta=(1,1)$, including the y-advection term, and evaluate the residual; if the residual is not zero, the benchmark and all error tables do not measure the model problem the paper intends to solve.","supporting_citations":[{"cited_title":"Eriksson, C","cited_arxiv_id":null,"evidence_quote":"supplies the two benchmark advection-diffusion problems and the exact solutions used in every error table."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the inverse inequality for B-splines from which the coercivity constant C and the condition ε ≤ h/(2C²) are derived."},{"cited_title":"Hughes, Leopoldo P","cited_arxiv_id":null,"evidence_quote":"defines the SUPG method that serves as the main baseline in the numerical comparisons."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the least-squares finite element framework and its stability properties that the L2-optimal-test formulation builds on."},{"cited_title":"J., Franca, L","cited_arxiv_id":null,"evidence_quote":"introduces the Galerkin/least-squares formulation that the combined test-function scheme is designed to recover."},{"cited_title":"Error bounds for finite element method","cited_arxiv_id":null,"evidence_quote":"supplies the stability theory used to conclude well-posedness of the discrete formulation from continuity and coercivity."}],"review_version":1}