{"id":"203f9dc6-8abc-40da-8638-0cc33c602ae4","arxiv_id":"2608.08505","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A rank-r coarse space built from the top singular modes of the local scattering map makes the GOSM skeleton operator identity-plus-compact, giving superlinear GMRES convergence with explicit bounds.","lead":"For a specialized wave-equation solver called GOSM, this paper proves that the hard operator differs from a simple involution by a compact perturbation, and builds a coarse space preconditioner from that fact with proven superlinear convergence. The practical payoff awaits a promised follow-up paper with the implemented discrete version.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The compactness and Schatten claims rest on the unverified Section 4(a) hypothesis that T is a compact/Schatten perturbation of T+; for the standard local impedance T=ηId the hypothesis is false, and no admissible T other than the nonlocal T=T+ is exhibited.","rationale":"The reader's verdict is already CONDITIONAL, and my stress-test does not move it. I agree that Section 4(a) is load-bearing: Lemma 5.2 uses it directly, and every quantitative result (Lemma 5.6, Corollaries 5.7, 6.2, 6.5, Proposition 6.4) inherits it. I add a concrete sharpening: the standard local impedance T=ηId is not a compact perturbation of T+, because H^{1/2} compactly embeds into L² and L² continuously into H^{-1/2}, so ηId is compact while T+ is not; hence the assumption excludes the most common optimized Schwarz impedance. The paper notes in Section 5.2 that if T=T+ then A−A+ is Schatten via Weyl laws, so a theoretical admissible choice exists, but this is an implicit restriction to a nonlocal DtN exchange operator with potentially high application cost, and it is not stated as the scope of the theorem. Separately, Lemma 5.2's displayed factorization is not algebraically exact: using M^{-1}-N^{-1}=M^{-1}(N-M)N^{-1} gives an extra term 2i(T−T+)B N^{-1}B*; the compactness proof still works because this term is compact, but Lemma 5.6's Schatten norm estimate undercounts the Cp norm by at least the Cp norm of this term divided by β. This is patchable and does not affect the qualitative superlinearity claim. The remaining gap is the absence of a rank-scaling analysis with frequency and subdomain count: estimate (47) still carries (1/γ⋆)^n and the tail singular values σ_{r+j}(S−iId), so the advertised robustness against high frequency and trapping is not established by the displayed bounds. This, too, is a limitation rather than an internal inconsistency, and it is already reflected in the CONDITIONAL verdict. The analytical test proposed above would settle whether an admissible T exists and whether the Schatten constants need revision.","tokens_in":15070,"tokens_out":20298,"duration_ms":219474,"concrete_test":"Take the standard local impedance T=ηId on one subdomain and compute the essential spectrum, or the singular value decay, of T−T+; since T+ is an order-one noncompact pseudodifferential operator and ηId is compact as a map H^{1/2}→H^{-1/2}, the difference is not compact. Then, for the same geometry, repeat the check with T=T+ and verify from (22) that S−iId=2iT+B(A−iB*T+B)^{-1}(A+−A)(A+−iB*T+B)^{-1}B* belongs to the promised Schatten class; if the only admissible T is T+, the theorem must be restated with T=T+ as an explicit hypothesis rather than a generic compact perturbation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that S−iId is compact (Lemma 5.2) and that the coarse space built from its leading singular modes gives superlinear GMRES (Corollaries 5.3, 5.7, 6.2, 6.5, Proposition 6.4) has a single load-bearing premise: the standing assumption in Section 4(a) that the impedance operator T differs from the reference positive DtN operator T+ by a compact (later Schatten-p) perturbation. The paper never constructs, for the Helmholtz problem (1), an operator T satisfying this hypothesis. This is not a cosmetic gap: the natural local impedance used in optimized Schwarz methods, T=ηId (with the L² pairing), maps H^{1/2}(∂Ω) compactly into H^{-1/2}(∂Ω) via H^{1/2}⊂⊂L²⊂H^{-1/2}, whereas T+ is an order-one noncompact isomorphism; hence T−T+ is not compact and the entire construction fails for that choice. If the only admissible choice is the nonlocal operator T=T+ (mentioned in passing in Section 5.2), the paper should state this as an explicit scope restriction and address its computational cost. A secondary, patchable defect is that the factorization displayed in the proof of Lemma 5.2 drops the term 2i(T−T+)B(A+−iB*T+B)^{-1}B*; the qualitative compactness result survives, but Lemma 5.6 and Corollaries 5.7, 6.2, 6.5 miss an additive ∥T−T+∥_{C_p} contribution in their Schatten constants.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes a continuous-level skeleton formulation of the Generalized Optimized Schwarz Method (GOSM) for the Helmholtz problem (1) in a non-overlapping domain decomposition, written as (Id + ΠS)p = f. Under the standing assumption that the impedance operator T differs from the reference DtN operator T+ by a compact (eventually Schatten-p) perturbation, the paper proves that S − iId is compact, derives superlinear GMRES convergence bounds via Moret's theorem, and proposes two preconditioners: a first-order one based on (Id − iΠ)/2 and a refined coarse-space preconditioner built from a rank-r truncated SVD Sr = iId + Kr of S, inverted by the Woodbury formula. The main results are Proposition 6.4 and Corollary 6.5, which bound the GMRES residual of the preconditioned equation (46) by products of tail singular values σ_{r+j}(S − iId), yielding superlinear convergence controlled by the singular-value decay.","tokens_in":15199,"tokens_out":11456,"duration_ms":112130,"significance":"The paper is a credible and largely correct continuous-level analysis of a parametrix for GOSM. The algebraic chain from Lemma 5.1 to Corollary 5.3 is clean, the use of Moret's theorem is appropriate, and the Woodbury-based coarse-space inversion in Section 6.2 is a concrete algorithmic proposal. If the standing compactness assumption on T is satisfied, the proposed coarse space removes the subdomain-count and inf-sup dependence from the asymptotic convergence rate, and the explicit Schatten-class estimates give quantitative rates. The main weakness is that the key hypothesis is not instantiated: apart from the passing mention of T = T+, the paper gives no concrete nonlocal T satisfying the compactness/Schatten assumption, and it does not analyze the cost of applying such a T. This limits the applicability of the claimed construction as it stands, but the issue is local and fixable within the manuscript's scope.","major_comments":[{"comment":"The entire low-rank compressibility statement of the paper rests on the standing assumption in Section 4(a) that T − T+ is compact (or Schatten-p). The manuscript does not exhibit any concrete T for problem (1) other than T = T+, mentioned in passing in Section 5.2. For the standard local impedance T = ηId used in optimized Schwarz methods, T − T+ is not compact: T+ is an order-one isomorphism H^{1/2}(∂Ωj) → H^{-1/2}(∂Ωj), whereas ηId is compact in that scale. The construction therefore fails for the usual Robin-type transmission condition. Please either (i) prove that a class of admissible nonlocal impedance operators satisfies the compactness/Schatten assumption and discuss the cost of applying them, or (ii) state explicitly and prominently that the method is restricted to such T, with T = T+ as the only known example. As written, the reader cannot tell whether the proposed coarse space is applicable to the GOSM instances that motivated the paper.","section":"Section 4(a), Lemma 5.2, Corollaries 5.7, 6.2, 6.5"},{"comment":"The factorization displayed in the proof of Lemma 5.2 is incomplete. Setting R = (A − iB*TB)^{-1} and R+ = (A+ − iB*T+B)^{-1}, a correct derivation gives S − iId = 2i[(T − T+)B R+ B* + T B R(A+ − A − iB*(T+ − T)B)R+ B*], whereas the proof drops the term 2i(T − T+)B R+ B*. The compactness conclusion survives because this extra term is compact, but the Schatten bound in Lemma 5.6 and the subsequent quantitative estimates in Corollaries 5.7, 6.2, and 6.5 need an additional contribution involving ∥T − T+∥_{C_p} with a constant factor, and they are not supported by the proof as stated. Please correct the factorization and re-derive the constants.","section":"Lemma 5.2 proof and Lemma 5.6"}],"minor_comments":[{"comment":"The abbreviation 'GMRes' is used consistently but is nonstandard; the conventional spelling is 'GMRES' (see [30,31]). Please make the notation consistent.","section":"Throughout"},{"comment":"The sentence referring to 'the explicit expression (2)' should refer to formula (22) for the scattering operator S, not to the variational form (2).","section":"Section 5.2"},{"comment":"Typos: 'the later' should be 'the latter', 'developpement' should be 'development', 'reasonnable' should be 'reasonable', and 'equiped' should be 'equipped'. Also, 'Stecklov-Poincaré' should be 'Steklov-Poincaré'.","section":"Introduction and Section 6.2"},{"comment":"The inequality '⌊n/4⌋ ≤ n/4 ≤ ⌊n/4⌋ + 1' should have a strict inequality '⌊n/4⌋ ≤ n/4 < ⌊n/4⌋ + 1'; the subsequent argument is unaffected.","section":"Section 5 (after Eq. (31))"}],"recommendation":"major_revision","confidential_remarks":"The manuscript relies heavily on the author's prior work [3,4,6] for the existence and properties of the impedance operator T. The editor may wish to ask the authors for a precise pointer to where the compactness of T − T+ is established, or to require them to narrow the stated scope. The paper is explicitly Part I of a theoretical study, and the promised discrete counterpart is not included; novelty relative to the author's earlier papers should be clarified."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, this paper is the first to propose a coarse space construction specific to GOSM, and the key observation is genuinely new: Lemma 5.2 shows S - iId is compact, which leads to a parametrix and superlinear GMRES estimates. I checked the main algebraic chain - Lemma 5.1, Corollary 5.3, the unitary factor algebra, and the Woodbury argument in Proposition 6.4 - and those steps hold. Second, the whole construction hangs on an assumption about the impedance operator T that the paper never verifies for problem (1). Section 4(a) assumes T - T+ is compact, with a Schatten-p strengthening later. For the standard local impedance T = eta Id, T is compact as a map from H^{1/2} to H^{-1/2} while T+ is an order-one noncompact isomorphism, so T - T+ is not compact. The only admissible example mentioned is T = T+, which is nonlocal and computationally costly. That is a real scope restriction, not a cosmetic gap.\n\nThere is also a patchable defect: the factorization in the proof of Lemma 5.2 drops the term involving (T - T+) B (A+ - iB*T+B)^{-1} B*. That term is compact under the standing assumption, so the qualitative compactness survives, but Lemma 5.6 and Corollaries 5.7, 6.2, 6.5 miss an additive norm(T-T+)_{C_p} contribution. The coercivity bound (25) is imported from the author's prior work; it is load-bearing and not derived here, so the constants should be read with caution. Also, the claimed robustness against frequency and trapping outruns what (41) and (47) actually deliver, since those bounds still carry powers of 1/gamma* and the required rank growth is not analyzed.\n\nWho is this for? Specialists in non-overlapping DDM for Helmholtz. The paper deserves a serious referee. It is honest about being continuous-only and promises discrete implementation in Part II. My recommendation: accept for peer review, but require the corrected factorization, an explicit statement of admissible T with the local-impedance restriction addressed, and either a rank-scaling analysis or an honest deferral.","headline":"First coarse space for GOSM, built on a genuinely new compactness observation, but the central hypothesis on T is unverified and fails for the standard local impedance, so the scope is narrower than advertised.","tokens_in":696,"tokens_out":1517,"would_cite":true,"duration_ms":68724,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65N55","65F10","35J05","47B06"],"pacs":[],"model":"deepseek-v4-flash","headline":"The Generalized Optimized Schwarz Method's scattering operator is a compact perturbation of i times the identity, so its coarse-space-preconditioned GMRES converges superlinearly.","keywords":["generalized optimized Schwarz method","coarse space","domain decomposition","Helmholtz equation","GMRES","compact perturbation","Schatten class","parametrix"],"falsifier":"Take a two-subdomain Helmholtz problem in a square, choose T so that T−T+ is a rank-one smoothing perturbation, compute the singular values σ_j(S−iId) and the GMRES residuals for the preconditioned equation (46), and compare with the product bound (47): the central claim fails if residuals do not outpace the Elman geometric bound (27) while the singular-value product tends to zero. Conversely, taking T as a local Robin impedance (multiplication by a constant) should make σ_j(S−iId) plateau away from zero, showing that the compact-perturbation premise is what the construction needs.","tokens_in":14613,"feed_emoji":"🌊","tokens_out":11125,"duration_ms":111263,"temperature":0.7,"pith_summary":"The paper aims to show that the Generalized Optimized Schwarz Method (GOSM), a non-overlapping domain decomposition method for harmonic wave propagation, is low-rank compressible in a precise sense: its local scattering operator S differs from i times the identity by a compact operator. Working entirely at the continuous level, the author proves that the preconditioned skeleton equation is identity plus compact, so GMRES converges superlinearly with a rate controlled by the singular values of S−iId. The coarse space is the span of the leading singular modes of S−iId, and the Woodbury formula makes the coarse solve practical. If the argument holds, convergence of the preconditioned iteration depends on how fast those singular values decay, not directly on the number of subdomains. This first part is continuous; a discrete counterpart is announced in the paper.","feed_headline":"Compact tails make GOSM converge superlinearly","feed_subtitle":"A coarse space built from the leading singular modes of the scattering map controls the GMRES rate.","key_machinery":"The object doing the work is S−iId, the deviation of the subdomain scattering operator from i times the identity. Lemma 5.1 shows that replacing A and T by the reference operators A+ and T+ makes S exactly iId; Lemma 5.2 then shows that the same formula with the true operators is a compact perturbation. Corollary 5.3 lifts this to compactness of (Id+ΠS)^4+4Id using Π²=Id, which is what fits Moret's GMRES theorem. The refined preconditioner is built from the truncated singular value decomposition of S−iId, and the Woodbury formula (45) inverts the resulting rank-r update. The single identity that carries the argument is (Id+ΠS)^4 = −4(Id+X)^4 with X = (Π−iId)(S−iId)/2, which turns the compact tail of S−iId into a compact tail of the preconditioned operator.","core_discovery":"The central claim is Lemma 5.2: under the hypothesis that the impedance operator T used in the transmission condition differs from the reference positive-definite Dirichlet-to-Neumann operator T+ by a compact perturbation, S−iId is compact as a map from H(Σ)′ to H(Σ)′. The proof factors S−iId as 2iTB(A−iB*TB)^{-1}(A+−A−iB*(T+−T)B)(A+−iB*T+B)^{-1}B*, using Rellich compactness for A−A+ and the assumed compactness of T−T+. Because Π²=Id, this gives compactness of (Id+ΠS)^4+4Id, and via Moret's singular-value theorem it yields superlinear GMRES convergence for the unpreconditioned equation. The refined preconditioner replaces S by Sr=iId+Kr, where Kr is a rank-r truncation of S−iId; the preconditioned operator (Id+ΠSr)^{-1}(Id+ΠS) is then identity plus a residual whose GMRES factors are σ_{r+j}(S−iId)/γ⋆. The coarse space of dimension r is therefore the span of the leading singular modes of S−iId, and Proposition 6.4 converts this tail into a convergence bound.","pith_inferences":["This is effectively a parametrix construction: the coarse space is a finite-dimensional inverse of the compact tail, so the same template could apply to other wave domain-decomposition formulations whose scattering maps are compact perturbations of constants.","The continuous analysis suggests a discrete-hindsight prediction: on fine meshes, the convergence rate of the preconditioned iteration should track the decay of the discrete singular values of S−iId, so adaptively choosing r to capture all singular values above γ⋆ should make the method insensitive to frequency and to the number of subdomains.","The paper leaves the choice of T open; if a concrete non-local impedance can be constructed for which the Schatten p-norm ∥T−T+∥_{C_p} is small, the convergence bounds would become quantitative rather than existential.","Because the analysis is at the continuous level, a cheap numerical check of the singular-value tail on a family of subdomain partitions would test whether the predicted superlinear phase appears in practice."],"forward_implications":["GMRES on the unpreconditioned GOSM skeleton equation converges superlinearly, at a rate governed by the singular values of (Id+ΠS)^4+4Id divided by γ⋆^4 (Estimate (31)).","The simple preconditioner (Id−iΠ)/2 gives the cleaner bound ∏_{j=1}^n σ_j(S−iId)/γ⋆, removing the fourth power on γ⋆ (Lemma 6.1).","With a rank-r coarse space from the truncated SVD of S−iId, the GMRES residual bound shifts to the tail singular values σ_{r+j}(S−iId), so enlarging r can compensate a small γ⋆ (Proposition 6.4).","When A−A+ and T−T+ lie in a Schatten p-class, the product bounds become explicit rates involving (∥S−S_r∥_{C_p}/(γ⋆ n^{1/p}))^n (Corollaries 5.7, 6.2, 6.5).","The coarse-space solve is practical because the Woodbury formula reduces it to inverting an r×r matrix (Equation (45))."],"supporting_citations":[{"why":"Defines GOSM, the exchange operator Π, the local scattering operator S, and the skeleton equation (24); also supplies the equivalence Lemma 4.1.","marker":"[4]"},{"why":"Supplies the skeleton formulation and the coercivity estimate (25) that fixes γ⋆.","marker":"[3]"},{"why":"Its arguments are adapted to prove the coercivity bound for the GOSM operator Id+ΠS.","marker":"[6]"},{"why":"Provides Theorem 5.4, the singular-value product bound that turns compactness into superlinear GMRES convergence.","marker":"[21]"},{"why":"Gives the product-of-singular-values inequality used in Lemma 5.5 for Schatten-class rates.","marker":"[35]"},{"why":"Defines Schatten classes and the characterization of compact operators by decaying singular values.","marker":"[1]"},{"why":"Supplies the Weyl law used to place A−A+ in a Schatten class in the T=T+ case.","marker":"[23]"},{"why":"Gives the Poincaré-Steklov eigenvalue asymptotics used to place T−T+ in a Schatten class.","marker":"[27]"},{"why":"Provides the Woodbury formula that inverts the rank-r coarse-space correction.","marker":"[36]"}],"fun_headline_variants":["Compact tails make GOSM converge superlinearly","Singular modes build coarse space for GOSM preconditioner","Prefiltered GOSM converges superlinearly via compact tails","Coarse space from scattering map singular modes speeds GMRES"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the impedance operator T can be chosen so that T−T+ is compact (and Schatten-p for the quantitative rates), and that the coercivity constant γ⋆ from (25) is strictly positive; the paper assumes, rather than constructs, such a T for the Helmholtz problem (1). If T−T+ is not compact, S−iId need not have a decaying singular-value tail, and all the superlinear estimates lose their engine.","fun_headline_variants_meta":{"raw":{"variants":["Compact tails make GOSM converge superlinearly","Singular modes build coarse space for GOSM preconditioner","Prefiltered GOSM converges superlinearly via compact tails","Coarse space from scattering map singular modes speeds GMRES"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000156,"raw_usage":{"total_tokens":1198,"prompt_tokens":904,"completion_tokens":294,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":520,"completion_tokens_details":{"reasoning_tokens":224}},"tokens_in":520,"tokens_out":294,"duration_ms":3613,"temperature":1.0,"reasoning_tokens":224,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:38:34.548427+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a two-subdomain Helmholtz problem in a square, choose T so that T−T+ is a rank-one smoothing perturbation, compute the singular values σ_j(S−iId) and the GMRES residuals for the preconditioned equation (46), and compare with the product bound (47): the central claim fails if residuals do not outpace the Elman geometric bound (27) while the singular-value product tends to zero. Conversely, taking T as a local Robin impedance (multiplication by a constant) should make σ_j(S−iId) plateau away from zero, showing that the compact-perturbation premise is what the construction needs.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines GOSM, the exchange operator Π, the local scattering operator S, and the skeleton equation (24); also supplies the equivalence Lemma 4.1."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the skeleton formulation and the coercivity estimate (25) that fixes γ⋆."},{"cited_title":"Robusttreatmentofcross-pointsinoptimizedSchwarzmethods","cited_arxiv_id":null,"evidence_quote":"Its arguments are adapted to prove the coercivity bound for the GOSM operator Id+ΠS."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides Theorem 5.4, the singular-value product bound that turns compactness into superlinear GMRES convergence."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the product-of-singular-values inequality used in Lemma 5.5 for Schatten-class rates."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines Schatten classes and the characterization of compact operators by decaying singular values."},{"cited_title":"Netrusov and Y","cited_arxiv_id":null,"evidence_quote":"Supplies the Weyl law used to place A−A+ in a Schatten class in the T=T+ case."},{"cited_title":"Rozenblum","cited_arxiv_id":null,"evidence_quote":"Gives the Poincaré-Steklov eigenvalue asymptotics used to place T−T+ in a Schatten class."},{"cited_title":"Woodbury.Inverting modified matrices","cited_arxiv_id":null,"evidence_quote":"Provides the Woodbury formula that inverts the rank-r coarse-space correction."}],"review_version":1}