{"id":"de1e748b-7fc5-4c60-ac63-3ef78263fc6c","arxiv_id":"2506.02429","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A coordinate-free, arbitrary-dimension Rust implementation of Finite Element Exterior Calculus with first-order Whitney forms, validated on Hodge-Laplace eigenvalue and source problems.","lead":"A BSc thesis presents formoniq, a Rust library that solves Hodge-Laplace problems on arbitrary-dimensional meshes using differential forms, with geometry stored intrinsically as edge lengths instead of global coordinates. It is verified on two problems: the 1-form spectrum on a torus and a manufactured Poisson-like source problem on cubes in 2D and 3D.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: the manufactured solution satisfies the actual natural boundary conditions in star-trace form, and the observed first-order L2 convergence is consistent with lowest-order Whitney 1-forms.","rationale":"I read the manuscript in good faith as a BSc thesis whose central claim is that formoniq solves Hodge-Laplace eigenvalue and source problems for first-order Whitney forms on coordinate-free simplicial complexes. The strongest evidence is the recovery of two zero eigenvalues on the torus (b1 = 2) and the clean O(h) convergence for both the solution and its exterior derivative in the manufactured source problem. I specifically stress-tested the reader's selected weakest assumption: that the manufactured solution's zero-trace assertion is unverified and that the RHS may omit boundary contributions. Direct analytic inspection shows this concern is unfounded once the boundary conditions are read in their correct star-trace form from Eq. 82. On each boundary face, the normal component of u vanishes and the normal-tangential components of du vanish, so the mixed variational formulation (Eq. 83) is genuinely boundary-free. The paper's Eq. 93 is sloppy in omitting the Hodge star, but the actual conditions needed by the weak form are satisfied. I also examined the O(h) L2 rate. The paper worries that FEEC theory predicts O(h^2) via Aubin-Nitsche, but for lowest-order Whitney 1-forms the local space contains all constant forms yet not all linear 1-forms; a simple reference-triangle interpolant of a linear field such as u = x dx yields a constant, so L2 approximation is first-order. Thus the measured O(h) is exactly what this element should give, and the implementation's convergence results are consistent with correctness. The remaining caveats (no codifferential error, no k > 1 or n > 3 runs, flat-torus baseline for a curved torus, no commit-hash-pinned reproducibility) are real limitations on the strength of the validation, but they do not constitute a load-bearing flaw in the central mathematical or computational claim. I therefore find no significant objection that would change the reader's conditional verdict.","tokens_in":49309,"tokens_out":56655,"duration_ms":613963,"concrete_test":"Evaluate tr ⋆ u and tr ⋆ du numerically at sample points on each boundary facet of the [0,π]^n meshes for n = 2 and n = 3, using the analytic expressions Eq. 88 and Eq. 90; confirm both quantities vanish to round-off. As an additional check, rerun the 2D source convergence study with a higher-order quadrature rule for the RHS assembly; if the L2 rate changes materially, the reported convergence was quadrature-limited, whereas if it stays at O(h), the rate is the intrinsic Whitney approximation rate.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reader's weakest assumption does not land. The strong Hodge-Laplace source problem (Eq. 82) imposes tr ⋆ u = 0 and tr ⋆ du = 0, not the plain traces Tr u = 0 and Tr du = 0 written in Eq. 93. For the manufactured 1-form u_i = sin^2(x_i) ∏_{j≠i} cos(x_j), on the boundary face x_l = 0 or π the normal component u_l vanishes, so tr ⋆ u = 0. For du, the coefficient of any normal-tangential term dx_l ∧ dx_i contains the factor sin(x_l), which vanishes on that face; the remaining tangent-tangent components do not contribute to the pullback of ⋆ du. Hence tr ⋆ du = 0 as well. The mixed formulation (Eq. 83) therefore has no missing boundary contributions and the assembled RHS is complete. The O(h) L2 rate flagged by the authors is also not a correctness red flag: lowest-order Whitney 1-forms reproduce constant forms but not general linear 1-forms, so first-order L2 approximation is the expected intrinsic rate for this element, not O(h^2). The remaining limitations (no codifferential error measurement, no validation for k ≠ 1 or n > 3, and the idealized flat-torus comparison) reduce the breadth of the validation but do not undermine the central claim that the library solves the stated Hodge-Laplace problems.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a Rust library, formoniq, implementing Finite Element Exterior Calculus (FEEC) on abstract, coordinate-free simplicial complexes in arbitrary dimensions. The geometry is encoded through edge lengths via a Regge-calculus-derived Riemannian metric. The library provides first-order Whitney basis functions, closed-form element mass matrices, discrete exterior derivatives via incidence matrices, and solvers for Hodge-Laplace eigenvalue and source problems using mixed formulations. Numerical validation includes a 1-form eigenvalue problem on a torus (where the two zero eigenvalues correctly reflect b1 = 2) and a manufactured-solution convergence study for the 1-form source problem in 2D and 3D, showing O(h) convergence for both the solution and its exterior derivative.","tokens_in":49495,"tokens_out":32736,"duration_ms":317247,"significance":"If the claims hold, this is a useful contribution as a working reference implementation of coordinate-free FEEC in arbitrary dimensions, with a clean separation between topology and metric geometry. Concrete strengths include a machine-checked unit test for the defining Whitney integral property (Section 4.2.1), closed-form integration of Whitney mass matrices without quadrature (Eq. 70), and external validation targets: an analytic manufactured solution and the topological requirement of two harmonic 1-forms on the torus. The observed first-order convergence is consistent with lowest-order Whitney elements. The main weaknesses are the narrow validation scope (only 1-forms, n = 2,3, trivial topology for the source problem, and only a heuristic flat-torus baseline for the eigenvalue problem) and several presentation inaccuracies, but these do not, in this referee's reading, invalidate the central derivation.","major_comments":[],"minor_comments":[{"comment":"The boundary conditions are written as Tr u = 0 and Tr du = 0, but the strong form in Eq. (82) imposes tr ⋆ u = 0 and tr ⋆ du = 0. Although the manufactured solution in Eq. (88) does satisfy the correct star-trace conditions (the normal component of u and the normal-tangential components of du vanish on the boundary), the notation should be corrected and a brief verification or explicit statement added, since the current text asserts the wrong trace conditions.","section":"Section 7.2, Eq. (93)"},{"comment":"The discussion of the L2 convergence rate is misleading. For lowest-order Whitney 1-forms the L2 approximation order is O(h) because the space contains constant 1-forms but not all linear 1-forms; an Aubin-Nitsche duality argument does not generally raise the L2 error to O(h^2) for these elements. The observed O(h) rate for the solution value is therefore expected and not a sign of a flaw; the paragraph should be revised to state this.","section":"Section 7.2, convergence discussion"},{"comment":"The block matrix in Eq. (81) is described as 'symmetric indefinite', but with D = C^T the off-diagonal blocks are -C and D, which are not transposes of each other. The generalized eigenvalue problem is equivalent to a symmetric problem via the Schur complement (C^T M^{-1} C + L), but the matrix itself is non-symmetric. This statement should be corrected.","section":"Chapter 6.1, after Eq. (81)"},{"comment":"The comparison to the idealized flat torus spectrum is only heuristic, since the mesh is a curved torus. The first non-zero group (4.116 vs 4) is close, but the second group (about 14.45 vs 16) deviates by roughly 10%. The language 'excellent agreement' and 'closely matches' overstates the agreement; the authors should either soften these claims or compute reference eigenvalues for the specific piecewise-flat discretized torus.","section":"Section 7.1"},{"comment":"The source problem is validated only for the trivial cohomology case (homology_dim = 0 on the cube). The implementation of the harmonic-form constraints in Eq. (83), which is a central part of the mixed formulation for general topology, is not exercised in any numerical test. As the paper explicitly states that verification is partial, this is a limitation rather than an error, but a test on a domain with non-trivial topology would considerably strengthen the central claim about the source problem.","section":"Section 7.2"},{"comment":"Minor typographical and language issues: 'combing' should be 'combining' (Section 4.2.1), 'irregardless' should be 'regardless' (Section 3.1), and the bibliography lists 'mathemtical' (reference [10]). These do not affect the technical content.","section":"Throughout"}],"recommendation":"minor_revision","confidential_remarks":"The manuscript is clearly a BSc thesis and reads as a software report rather than a focused journal article. The central numerical approach is sound and the code-level details are consistent with the FEEC theory, but the validation is narrow. I recommend minor revision for the technical content; the editorial fit and the amount of implementation detail (full code listings, package installation instructions) may need to be considered separately by the editor."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nShort version: this BSc thesis delivers a working, arbitrary-dimension FEEC library in Rust with coordinate-free simplicial complexes (metric from edge lengths via Regge calculus), Whitney 1-forms, and mixed Hodge-Laplace eigenvalue and source solvers. The mathematical content is standard FEEC, but the software combination is genuinely new and useful. The reader's weakest assumption about missing boundary terms in the manufactured-solution validation does not land: the strong problem imposes tr ⋆ u = 0 and tr ⋆ du = 0, and the chosen u satisfies those in star-trace form, so no boundary contributions are missing. The O(h) L2 rate is also not a red flag; lowest-order Whitney 1-forms only reproduce constant forms, so first-order L2 convergence is the intrinsic expected rate.\n\nWhat the paper does well: the assembly follows the standard mixed FEEC derivation, the Whitney basis property is unit-tested for dimensions 0 through 4, and the torus eigenvalue problem recovers the two zero modes (b1 = 2) cleanly. The manufactured source problem shows clean O(h) convergence in the H(d) seminorm in 2D and 3D. The author is candid about scope: no codifferential error, only k = 1, only n = 2, 3, and an idealized flat-torus baseline. Calling a 10% deviation in the second eigenvalue group \"excellent agreement\" is generous, but the topological and spectral validation still holds.\n\nSoft spots, in proportion: the validation is deliberately partial, and the thesis says so. For a software contribution, that is the main limitation. I would want a commit hash and archived test meshes before reproducing the numbers; as posted, the build depends on specific PETSc/SLEPc versions. The \"arbitrary dimensions\" claim is demonstrated only up to 3, but the data structures, quadrature, and assembly are dimension-generic, so I credit it.\n\nWho this is for: anyone building FEEC or DEC software, or working with intrinsic geometry discretizations (numerical relativity is an obvious audience). The comparison to other libraries is brief but useful.\n\nVerdict: this deserves a serious referee. The central claims hold; the open items are breadth and reproducibility, not correctness. If I were refereeing, I would accept with requests for versioned code and one extra validation, say a 2-form problem or a curved simplicial surface.","headline":"A solid, honest BSc thesis that ships a working arbitrary-dimension, coordinate-free FEEC library in Rust; the math is known, the software is the contribution, and the reader's boundary-trace concern does not land on the actual equations.","tokens_in":50113,"tokens_out":3609,"would_cite":true,"duration_ms":38290,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65N30","58A10","65N25"],"pacs":[],"model":"deepseek-v4-flash","headline":"A Rust library called formoniq solves Hodge-Laplace eigenvalue and source problems on coordinate-free n-dimensional simplicial complexes whose only geometry is edge lengths; the torus's two zero eigenvalues and O(h) convergence validate it.","keywords":["finite element exterior calculus","Hodge-Laplace operator","Whitney forms","coordinate-free geometry","Regge calculus","simplicial complex","discrete differential forms","Rust"],"falsifier":"Evaluate, by high-order quadrature, the boundary integrals that the variational formulation silently drops: on each face of the cube $[0,\\pi]^n$, pull back the exact solution $u$ and the trace of $\\star du$ to the face and integrate them against Whitney test forms. Any nonzero result means the assembled right-hand side $\\mathbf{b}$ is missing boundary contributions, and the reported $O(h)$ convergence describes a source problem whose data differ from those stated in Eq. (82).","tokens_in":48963,"feed_emoji":"🦀","tokens_out":16148,"duration_ms":147577,"temperature":0.7,"pith_summary":"The paper reports formoniq, a finite element library written in Rust that solves elliptic Hodge-Laplace eigenvalue and source problems on $n$-dimensional simplicial complexes whose geometry consists only of edge lengths, with no vertex coordinates anywhere. The claim is that a full Finite Element Exterior Calculus pipeline—Whitney $k$-form basis functions, Hodge mass matrices, the discrete exterior derivative, and the mixed weak formulation—works on abstract coordinate-free manifolds, with each cell's Riemannian metric reconstructed intrinsically from its edges via the law of cosines. The supporting evidence is the 1-form spectrum of a torus, where the two zero eigenvalues correctly reproduce the first Betti number $b_1 = 2$, and manufactured-source tests in 2D and 3D showing errors in both the solution and its exterior derivative decreasing at rates close to $O(h)$. A reader should care because the approach cuts the usual dependence of FEEC implementations on an ambient embedding, opening the same solver to curved and high-dimensional domains.","feed_headline":"Two zero eigenvalues reveal the torus in a coordinate-free solver","feed_subtitle":"A Rust library solves Hodge-Laplace problems on n-dimensional meshes whose only geometry is edge lengths.","key_machinery":"The load-bearing construction is the pairing of the Whitney complex with an intrinsic metric. A Whitney $k$-form is the piecewise-linear differential form attached to a $k$-simplex, with degrees of freedom given by integrals over $k$-simplices; its local shape functions are built from barycentric coordinates, $\\lambda_{i_0\\dots i_k} = k!\\sum_{l=0}^{k}(-1)^l \\lambda_{i_l}\\, d\\lambda_{i_0}\\wedge\\cdots\\wedge\\widehat{d\\lambda_{i_l}}\\wedge\\cdots\\wedge d\\lambda_{i_k}$. The geometry enters only through each cell's constant metric tensor, recovered from edge lengths by the law of cosines in the spirit of Regge calculus, so the Hodge inner product on forms is computed from the inverse Gramian $G^{-1}$ with the standard determinant rule $\\langle e_I, e_J\\rangle = \\det[\\langle dx^{I_i}, dx^{J_j}\\rangle]$. All four Galerkin operators needed for the mixed Hodge-Laplace formulation reduce to Hodge mass matrices combined with the purely combinatorial coboundary matrices $d_k = \\partial_{k+1}^{\\mathsf{T}}$: $D_k = M_k d_{k-1}$, $C_k = d_{k-1}^{\\mathsf{T}} M_k$, and $L_k = d_k^{\\mathsf{T}} M_{k+1} d_k$. The eigenvalue problem is then the sparse symmetric-indefinite generalized eigenproblem $\\left[\\begin{smallmatrix} M_{k-1} & -C \\\\ D & L\\end{smallmatrix}\\right]\\left[\\begin{smallmatrix}\\sigma\\\\u\\end{smallmatrix}\\right] = \\lambda \\left[\\begin{smallmatrix}0&0\\\\0&M_k\\end{smallmatrix}\\right]\\left[\\begin{smallmatrix}\\sigma\\\\u\\end{smallmatrix}\\right]$, and the source problem augments this system with the harmonic-form block that enforces $u \\perp \\mathfrak{H}^k$.","core_discovery":"The central claim is that the discrete de Rham complex built from first-order Whitney forms can be fully assembled from just two inputs: the combinatorics of a simplicial complex and the length of every edge. Given those lengths, Regge calculus supplies each cell's constant metric tensor $\\mathbf{G}$ through the law of cosines, $G_{ij} = \\frac{1}{2}(d_{0i}^2 + d_{0j}^2 - d_{ij}^2)$, and the Hodge mass matrices, source vectors, and the mixed Hodge-Laplace system are then computed by quadrature and assembly without ever consulting an embedding. The paper's validation is two-fold: on a torus with major radius $0.5$ and minor radius $0.2$, the computed spectrum of $\\Delta_1$ contains two zero eigenvalues, matching the two harmonic 1-forms predicted by Hodge theory and the topology ($b_1 = 2$), followed by eigenvalue groups near the flat-torus values $4$, $16$, and $25$; and for the manufactured 1-form solution $u = \\sum_i \\sin^2(x_i)\\prod_{j\\neq i}\\cos(x_j)$ on $[0,\\pi]^n$, the $\\mathrm{L}^2$ errors of $u$ and $du$ both converge at roughly $O(h)$ in $n = 2$ and $n = 3$. The author presents these results as confirmation that intrinsic, coordinate-free discretization preserves both the topological and the geometric content of the Hodge-Laplacian.","pith_inferences":["The paper's own convergence tables show the $\\mathrm{L}^2$ rate stuck at $O(h)$ instead of the $O(h^2)$ an Aubin-Nitsche duality argument suggests; a natural next experiment is refining a curved, non-flat mesh and watching whether the rate stays at $O(h)$, which would implicate the piecewise-flat metric approximation rather than the discrete forms themselves.","Because the formulation separates topology from geometry, the same code path could compute harmonic representatives on any closed manifold given only edge lengths, effectively turning the eigensolver into a largely combinatorial cohomology tool once the metric is treated as arbitrary input.","The paper leaves the codifferential error $\\|\\delta(u-u_h)\\|_{\\mathrm{L}^2}$ unmeasured, so the claimed convergence in the full $H\\Lambda$ norm is incomplete; measuring it on the same manufactured solutions is the cheapest direct check of the remaining half of the Hodge-Laplacian.","The torus eigenvalue groups sit slightly above the flat-torus predictions (about 4.116 versus 4, and 14.45 versus 16); refining the torus mesh and tracking whether these groups converge downward to the predicted values would separate pure discretization error from the geometric approximation of the curved torus."],"forward_implications":["The solver doubles as a homology detector: on a closed manifold, the number of zero eigenvalues of the discrete Hodge-Laplacian on $k$-forms equals the $k$-th Betti number, so the same code that produced $b_1 = 2$ on the torus gives the harmonic-form dimension on any mesh whose edge lengths are supplied.","Measured convergence of $\\|du - du_h\\|_{\\mathrm{L}^2}$ at rate near $O(h)$ matches the theoretical expectation for first-order Whitney elements in the $H\\Lambda$ norm, which the author takes as evidence that the exterior-derivative discretization and its assembly are correct.","Because only edge lengths define the geometry, any simplicial mesh with an admissible (Cayley-Menger-realizable) metric can be solved on, regardless of dimension or curvature, as long as the input comes as abstract simplices plus lengths.","The same mixed formulation covers every grade of the de Rham complex—grade 0 reduces to the scalar Laplacian with empty auxiliary blocks—so the implementation is not restricted to the 1-form problems used in the validation."],"supporting_citations":[{"why":"Supplies the mixed weak formulation of the Hodge-Laplace problem and the homological FEEC framework the thesis builds on.","marker":"[3]"},{"why":"Provides the Hodge-Laplace theory, the harmonic-form/cohomology link, the Galerkin matrix structure, and the expected first-order convergence rates in the H-Lambda norm.","marker":"[4]"},{"why":"Defines the Whitney basis forms and the integration and interpolation maps from which the discrete spaces are constructed.","marker":"[10]"},{"why":"Origin of the edge-length-based metric reconstruction (Regge calculus) used as the intrinsic geometry of each cell.","marker":"[9]"},{"why":"Justifies the first-order piecewise-flat geometry as an admissible geometric variational crime.","marker":"[11]"},{"why":"Supplies the Cayley-Menger realizability test ensuring that given edge lengths define a genuine Euclidean simplex.","marker":"[24]"},{"why":"The SLEPc eigensolver that computes the torus spectrum and the harmonic 1-forms in the eigenvalue validation.","marker":"[17]"},{"why":"The mesh generator used to produce the torus mesh for the non-trivial-topology eigenvalue experiment.","marker":"[21]"}],"fun_headline_variants":["Rust FEEC solver zeroes in on torus topology from edge lengths alone","Coordinate-free Rust solver recovers torus's two zero modes","Edge lengths alone: Rust FEEC finds torus harmonic 1-forms","Intrinsic FEEC in Rust: edge metrics capture torus topology","Rust library computes Hodge-Laplace from edge lengths only"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The convergence study rests on the assertion that the manufactured solution $u$ and its exterior derivative $du$ have vanishing traces on the boundary of $[0,\\pi]^n$, which is what removes all boundary terms from the mixed weak form; the paper states this assertion but never verifies it by computation, so if any of those boundary integrals is actually nonzero, the reported $O(h)$ rates describe a different boundary-value problem than the one stated.","fun_headline_variants_meta":{"raw":{"variants":["Rust FEEC solver zeroes in on torus topology from edge lengths alone","Coordinate-free Rust solver recovers torus's two zero modes","Edge lengths alone: Rust FEEC finds torus harmonic 1-forms","Intrinsic FEEC in Rust: edge metrics capture torus topology","Rust library computes Hodge-Laplace from edge lengths only"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000583,"raw_usage":{"total_tokens":2744,"prompt_tokens":949,"completion_tokens":1795,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":565,"completion_tokens_details":{"reasoning_tokens":1711}},"tokens_in":565,"tokens_out":1795,"duration_ms":10758,"temperature":1.0,"reasoning_tokens":1711,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T11:25:01.071243+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate, by high-order quadrature, the boundary integrals that the variational formulation silently drops: on each face of the cube $[0,\\pi]^n$, pull back the exact solution $u$ and the trace of $\\star du$ to the face and integrate them against Whitney test forms. Any nonzero result means the assembled right-hand side $\\mathbf{b}$ is missing boundary contributions, and the reported $O(h)$ convergence describes a source problem whose data differ from those stated in Eq. (82).","supporting_citations":[{"cited_title":"Whitney, Geometric integration theory","cited_arxiv_id":null,"evidence_quote":"Defines the Whitney basis forms and the integration and interpolation maps from which the discrete spaces are constructed."},{"cited_title":"Geometric Variational Crimes: Hilbert Complexes, Finite Element Exterior Calculus, and Problems on Hypersurfaces,","cited_arxiv_id":null,"evidence_quote":"Justifies the first-order piecewise-flat geometry as an admissible geometric variational crime."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Cayley-Menger realizability test ensuring that given edge lengths define a genuine Euclidean simplex."}],"review_version":1}