{"id":"afae4634-4fcc-42d0-a218-ded95df1a61d","arxiv_id":"1908.09326","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"A new Log-Cholesky metric on SPD matrices gives closed-form Frechet means and parallel transport and avoids the swelling effect.","lead":"This paper defines a new Riemannian metric, the Log-Cholesky metric, on the space of symmetric positive definite matrices by taking Cholesky factors and applying a logarithmic scale to their diagonal entries. The metric has simple closed-form formulas for averages, geodesics and parallel transport, and it avoids the swelling effect.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: the Log-Cholesky construction is internally consistent; the only real gap is an implicit completeness step in Proposition 9, trivially filled by the global Euclidean chart.","rationale":"I read the paper as a self-contained differential-geometric construction; its acceptance hinges on the claims that (L+,*) with ~g is a bi-invariant flat abelian Lie group and that Cholesky pushforward transfers this to S+_m. These claims are internally consistent: the group operation is associative and commutative, left- and right-invariance of ~g checks out, the coordinate computation gives zero Christoffel symbols, and the isometry S transfers all objects. Proposition 9 is the only place where an external theorem is used, and its hypotheses are satisfied once completeness is made explicit; this is a one-line fix rather than a flaw. Lemma 6's use of the Lie group exponential similarly needs surjectivity, which is automatic in this vector-group chart. I therefore find no load-bearing concern about the central claim. The reader's CONDITIONAL verdict is reasonable because of unsupported numerical performance claims and minor omitted justifications, but those do not threaten the mathematical result. My concrete test would settle the completeness point definitively.","tokens_in":17164,"tokens_out":20435,"duration_ms":220835,"concrete_test":"Verify analytically that ψ(L)=⌊L⌋+log D(L) is a global isometry from (L+,~g) onto (L,⟨·,·⟩_F): compute Dψ_L(X)=⌊X⌋+D(L)^{-1}D(X) and confirm ⟨DψX,DψY⟩_F equals ~g_L(X,Y). If confirmed, completeness and Proposition 9's appeal to [8] are justified; also confirm in this chart that the Lie group exponential of (L+,*) is ψ^{-1}, so Lemma 6's Y exists for every pair.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central mathematical claim is sound. The weakest step is Proposition 9, which invokes Theorem 2.1 of [8] for existence and uniqueness of the Fréchet mean. That theorem needs a complete, simply connected, nonpositively curved manifold; the paper explicitly proves flatness and simple connectedness but leaves completeness implicit. Completeness is immediate: ψ(L)=⌊L⌋+log D(L) is a global isometry from (L+,~g) to the Euclidean space L with the Frobenius metric (the same coordinate system used in the proof of Proposition 3 has constant identity metric), so (L+,~g) is complete; since S is an isometry, (S+_m,g) is complete as well. The other implicit point, surjectivity of the Lie group exponential in Lemma 6, also holds because the group is an abelian vector group in ψ-coordinates. Hence no load-bearing failure of the central argument exists. The unsupported numerical timing and stability claims are peripheral and would not change the mathematical core, but they do justify a request for code if those claims are to be reused.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a new Riemannian metric, the Log-Cholesky metric, on the manifold of symmetric positive definite matrices. The construction first defines a commutative Lie group operation X⋄Y = ⌊X⌋+⌊Y⌋+D(X)D(Y) and a bi-invariant metric on the Cholesky space L+ of lower-triangular matrices with positive diagonal, then pushes both forward to SPD matrices via the Cholesky parametrization S(L)=LL^T. The paper derives closed-form expressions for geodesics, Riemannian exponential and logarithmic maps, geodesic distance, and parallel transport on both L+ and S+_m, and shows that (S+_m,⋄) is an abelian Lie group with bi-invariant Log-Cholesky metric. It then studies Fréchet means, proving existence and uniqueness under a finite second-moment condition and giving closed forms for the Log-Cholesky mean and average. The determinant identity log det(mean) = E log det(Q) is established, implying that the average's determinant lies between the minimum and maximum determinants of the input matrices, so the average avoids the swelling effect. The paper also makes computational-efficiency and numerical-stability claims relative to the affine-invariant and Log-Euclidean metrics.","tokens_in":17354,"tokens_out":15270,"duration_ms":155893,"significance":"The mathematical core is sound and the paper delivers a genuinely useful, parameter-free Riemannian structure on SPD matrices with closed-form geodesics, exponential and logarithmic maps, parallel transport, and Fréchet mean, together with a bi-invariant metric on an abelian Lie group. The determinant identity and the resulting no-swelling property are proved directly from the definitions rather than assumed. This gives the SPD community an additional computational tool, particularly for large-scale or computation-heavy applications where the simplicity of the formulas matters. The novelty—combining Cholesky coordinates with a logarithmic treatment of the diagonal—is modest but real, and the paper is careful to note the lack of congruence invariance relative to the affine-invariant metric. The main mathematical claims are checkable and, apart from the completeness gap noted below, correct; the numerical claims are less well supported but are peripheral to the mathematical contribution.","major_comments":[{"comment":"The proof invokes Theorem 2.1 of [8], which requires a complete, simply connected, nonpositively curved manifold. The proof establishes simple connectedness via the diffeomorphism ψ and zero sectional curvature via Proposition 8, but it never states or proves completeness of (L+,~g). This is load-bearing for the existence and uniqueness conclusion. The gap is easily closed: ψ(L)=⌊L⌋+log D(L) is a global isometry from (L+,~g) to Euclidean (L,‖·‖_F), since dψ_L(X)=⌊X⌋+D(L)^{-1}D(X), so (L+,~g) is complete; completeness of (S+_m,g) then follows from the isometry S. Please add this argument before invoking [8].","section":"§4.1, proof of Proposition 9"},{"comment":"The lemma as stated applies to an arbitrary abelian Lie group with a bi-invariant metric, but the proof assumes the existence of Y with exp(Y)=p^{-1}q. This requires the group exponential to be surjective, which holds for connected abelian Lie groups but should be stated explicitly; for a disconnected abelian Lie group, p^{-1}q need not lie in the identity component. The later application to L+ is safe because L+ is isomorphic to the vector group (L,+) via ψ, but the lemma should be restated with a connectedness hypothesis, or with the exponential-image condition made explicit, so that the proof is valid as written.","section":"§3.4, Lemma 6"}],"minor_comments":[{"comment":"In the row labeled 'Riemannian metric', the second argument of ~g is written as (D_P L)(W), but it should be (D_P L)(V); the displayed formula should read g_P(W,V) = ~g_{L(P)}((D_P L)(W), (D_P L)(V)).","section":"Table 2"},{"comment":"In the expression for τ_{P,Q}(W), the second bracket uses D(K)D(L^{-1})D(X) while the first uses D(K)D(L)^{-1}D(X). Since D(L^{-1})=D(L)^{-1} for L∈L+, the two are equal, but the notation should be made consistent.","section":"§3.4, Proposition 7"},{"comment":"The map x is defined as 'x : L+ → R', but the intended codomain is R^{m(m+1)/2}; please correct the typo.","section":"§3.1, proof of Proposition 3"},{"comment":"The phrase 'minimizes F2j(ex)' should read 'minimizes F2j(e^x)' or 'minimizes F2j(exp(x))' for clarity.","section":"§4.1, proof of Proposition 10"},{"comment":"The text says 'log(P) is a symmetric metric for an SPD matrix P'; this should be 'symmetric matrix' rather than 'symmetric metric'.","section":"§2.1"},{"comment":"The numerical timing and stability claims (9.3ms, 0.85ms, 0.2ms; expected relative difference 3.3×10^-2; largest-eigenvalue ratios 10^10 and 10^15) are not accompanied by a full experimental protocol or code. Since the abstract advertises computational efficiency and numerical stability, please supply a reproducible procedure or clearly mark these statements as indicative and soften the corresponding claims.","section":"§3.4 and §4.2"}],"recommendation":"minor_revision","confidential_remarks":"The numerical efficiency and stability comparisons in §3.4 and §4.2 are not reproducible from the manuscript alone, and the abstract and conclusion treat them as part of the contribution. The mathematical core, however, is correct and the completeness gap in Proposition 9 is a one-sentence fix. I would encourage the editor to request that the author either provide code and a detailed protocol for the numerical experiments or remove the specific numbers and keep only qualitative statements. This does not affect my positive assessment of the mathematical content."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The Log-Cholesky metric is a genuine new construction, not a repackaging. The Cholesky-space Lie group with the componentwise operation and the diagonal log scaling is natural, and pushing it forward to SPD matrices gives a bi-invariant metric with closed-form geodesics, parallel transport, and Fréchet mean. I checked the main derivations: the geodesic formula, the pushforward metric, the Lie group structure, and the mean formulas are all correct. The no-swelling property follows from the determinant identity, not from circular reasoning.\n\nThe novelty claim holds against the cited work. Dryden et al. use Euclidean distance on Cholesky factors, which swells; Grubisic and Pietersz are on correlation matrices with a different structure. The log-Cholesky operation on the diagonal is not in those papers.\n\nThe soft spots are minor. Proposition 9 invokes Bhattacharya-Patrangenaru for the Fréchet mean, which requires completeness. The paper proves flatness and simple connectedness but not completeness. As the stress test notes, completeness is immediate from the global isometry ψ(L)=⌊L⌋+log D(L), so the theorem applies. It should have been stated, but it is not a load-bearing flaw. Similarly, surjectivity of the group exponential in Lemma 6 is implicit but holds because the group is a vector group in those coordinates.\n\nThe bigger weakness is the numerical section. The timing numbers (9.3ms vs 0.85ms vs 0.2ms) are presented without code, without specifying how many terms of the Log-Euclidean series were used, and without a protocol. The stability claims about 10^10 and 10^15 condition numbers are also anecdotal. These claims are peripheral to the geometry, but they need to be reproducible if the paper is going to make computational promises.\n\nThere is a small typo in Table 2: the metric row repeats W instead of V in the second argument. Easily fixed.\n\nThe paper is for anyone doing statistical analysis or imaging with SPD matrices who wants a faster alternative to Log-Euclidean or affine-invariant metrics. The central contribution is solid, and the presentation is clear. It deserves a serious referee; the main referee asks should be to fill the completeness step, clean up the numerical experiments, and add code or a detailed protocol.","headline":"Log-Cholesky is a real, useful addition to SPD geometry; the central derivations hold up, and the gaps are small.","tokens_in":17868,"tokens_out":1540,"would_cite":true,"duration_ms":14592,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["47A64","26E60","53C35","22E99","32F45","53C22","15A22"],"pacs":[],"model":"deepseek-v4-flash","headline":"The Log-Cholesky metric gives SPD matrices a flat, bi-invariant geometry with closed-form averages and parallel transport.","keywords":["symmetric positive definite matrices","Cholesky decomposition","Log-Cholesky metric","bi-invariant metric","Frechet mean","swelling effect","parallel transport","Lie group"],"falsifier":"Take $P_1=\\begin{pmatrix}\\epsilon^2&0\\\\0&1\\end{pmatrix}$ and $P_2=\\begin{pmatrix}1&0\\\\0&\\epsilon^2\\end{pmatrix}$ and compute the Log-Cholesky geodesic midpoint $\\gamma(1/2)$. The claim predicts $\\det\\gamma(1/2)=\\sqrt{\\det P_1\\det P_2}=\\epsilon^2$; any computed determinant larger than $\\max\\{\\epsilon^2,\\epsilon^2\\}=\\epsilon^2$ would refute the swelling-free claim.","tokens_in":1409,"feed_emoji":"📐","tokens_out":3563,"duration_ms":78357,"temperature":0.7,"pith_summary":"Symmetric positive definite matrices encode covariance or diffusion information, but standard averaging can inflate determinants. This paper argues that pushing a carefully chosen metric from Cholesky space to SPD space yields a Riemannian metric with unusually concrete properties. If the construction is right, practitioners get a geometry where averages are easy to compute, geodesics and parallel transport have closed forms, and the determinant of any average stays between the determinants of the inputs. The paper also shows the SPD manifold becomes an abelian Lie group with a bi-invariant metric under this operation.","feed_headline":"Log-Cholesky metric stops covariance averages from swelling","feed_subtitle":"A new Riemannian metric gives closed-form means and parallel transport, with determinants staying in bounds.","key_machinery":"The central object is the Cholesky decomposition $S(L)=LL^\\top$, used as an isometry between Cholesky space $L^+$ and the SPD manifold. On $L^+$ the paper defines the metric $\\tilde g_L(X,Y)=\\sum_{i>j}X_{ij}Y_{ij}+\\sum_j X_{jj}Y_{jj}L_{jj}^{-2}$. In the coordinate chart that logs the diagonal, this metric has constant coefficients, so all Christoffel symbols vanish; geodesics, exponentials, and the group operation are then computed by simple formulas and pushed forward to SPD matrices.","core_discovery":"On the space of lower-triangular matrices with positive diagonal, the paper puts a metric on each tangent space by combining the Euclidean inner product on the strictly lower-triangular part with the squared relative change on the diagonal. Through the Cholesky map $P=LL^\\top$, this metric becomes the Log-Cholesky metric on SPD matrices. The paper proves that this metric is bi-invariant for a commutative group operation on SPD matrices, has identically zero sectional curvature, and admits closed-form geodesics, exponential and logarithmic maps, parallel transport, and Fréchet averages. It also proves that $\\log\\det$ of the average equals the average of $\\log\\det$'s, so the determinant of the Log-Cholesky average lies between the minimum and maximum determinant of the averaged matrices.","pith_inferences":["Beyond the paper: because the Log-Cholesky coordinates are globally Euclidean, standard Euclidean statistical tools could be applied to the log-Cholesky representation and mapped back, with an exact geodesic interpretation rather than an approximation.","Beyond the paper: the abelian Lie group structure suggests that Fourier-style analysis or convolution-type operations on SPD-valued data may be meaningful, though the paper does not develop this direction.","Beyond the paper: the closed-form Log-Cholesky average could serve as a fast initialization or surrogate for affine-invariant or Log-Euclidean means in large-scale pipelines, especially when computational cost is the bottleneck."],"forward_implications":["The Log-Cholesky average of $P_1,\\dots,P_n$ has a closed form involving only Cholesky factors and matrix logarithms of diagonal blocks, so no numerical optimization is needed.","The determinant identity $\\det E_n = \\left(\\prod_i \\det P_i\\right)^{1/n}$ means the average never inflates dispersion, for any number of matrices.","Parallel transport along geodesics has a closed formula that uses only Cholesky factors and inverses of triangular matrices, making it substantially cheaper than the Log-Euclidean approach.","The SPD manifold becomes flat and abelian, so geodesic interpolation between two SPD matrices keeps determinants between the endpoint determinants.","Because the Log-Cholesky average shares its determinant with the Log-Euclidean and affine-invariant averages, it offers a fast route to averages with the same determinant behavior."],"supporting_citations":[{"why":"Supplies the Log-Euclidean baseline and the swelling-effect criterion that the Log-Cholesky metric is designed to satisfy.","marker":"[4]"},{"why":"Supplies the theorem used to prove existence and uniqueness of the Fréchet mean on flat simply connected manifolds.","marker":"[8]"},{"why":"Supplies the Cholesky distance whose interpolation the paper shows still swells, motivating the new metric.","marker":"[11]"},{"why":"Supplies the uniqueness theorem for Cholesky factorization, making the map between SPD and Cholesky spaces bijective.","marker":"[16]"},{"why":"Supplies the Lie-group lemma connecting the bi-invariant connection to parallel transport, used to derive the closed form.","marker":"[29]"},{"why":"Supplies the affine-invariant metric against which the paper compares simplicity, speed, and congruence invariance.","marker":"[34]"}],"fun_headline_variants":["Log-Cholesky: a simple SPD metric with bounded determinants","New metric for SPD matrices: fast, stable, no swelling","Cholesky-based metric gives closed-form SPD averages","Log-Cholesky metric: closed-form parallel transport, no swelling","A better Riemannian metric for SPD matrices via Cholesky"],"cache_read_input_tokens":20096,"weakest_assumption_plain":"Existence and uniqueness of the Log-Cholesky average rests on a theorem that requires the manifold to be complete, connected, and nonpositively curved; the paper proves flatness and simple connectedness but leaves the completeness check implicit.","fun_headline_variants_meta":{"raw":{"variants":["Log-Cholesky: a simple SPD metric with bounded determinants","New metric for SPD matrices: fast, stable, no swelling","Cholesky-based metric gives closed-form SPD averages","Log-Cholesky metric: closed-form parallel transport, no swelling","A better Riemannian metric for SPD matrices via Cholesky"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000905,"raw_usage":{"total_tokens":3869,"prompt_tokens":896,"completion_tokens":2973,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":512,"completion_tokens_details":{"reasoning_tokens":2888}},"tokens_in":512,"tokens_out":2973,"duration_ms":21877,"temperature":1.0,"reasoning_tokens":2888,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:16:04.592754+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take $P_1=\\begin{pmatrix}\\epsilon^2&0\\\\0&1\\end{pmatrix}$ and $P_2=\\begin{pmatrix}1&0\\\\0&\\epsilon^2\\end{pmatrix}$ and compute the Log-Cholesky geodesic midpoint $\\gamma(1/2)$. The claim predicts $\\det\\gamma(1/2)=\\sqrt{\\det P_1\\det P_2}=\\epsilon^2$; any computed determinant larger than $\\max\\{\\epsilon^2,\\epsilon^2\\}=\\epsilon^2$ would refute the swelling-free claim.","supporting_citations":[{"cited_title":"Arsigny, P","cited_arxiv_id":null,"evidence_quote":"Supplies the Log-Euclidean baseline and the swelling-effect criterion that the Log-Cholesky metric is designed to satisfy."},{"cited_title":"Bhattacharya and V","cited_arxiv_id":null,"evidence_quote":"Supplies the theorem used to prove existence and uniqueness of the Fréchet mean on flat simply connected manifolds."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Cholesky distance whose interpolation the paper shows still swells, motivating the new metric."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the uniqueness theorem for Cholesky factorization, making the map between SPD and Cholesky spaces bijective."},{"cited_title":"Milnor, Morse Theory, Princeton University Press, 1963","cited_arxiv_id":null,"evidence_quote":"Supplies the Lie-group lemma connecting the bi-invariant connection to parallel transport, used to derive the closed form."},{"cited_title":"Pennec, P","cited_arxiv_id":null,"evidence_quote":"Supplies the affine-invariant metric against which the paper compares simplicity, speed, and congruence invariance."}],"review_version":1}