[{"data":1,"prerenderedAt":-1},["ShallowReactive",2],{"doc-detail-83252-en":3,"doc-seo-83252-105":30,"detail-sidebar-cat-0-en-105":91},{"code":4,"msg":5,"data":6},0,"success",{"doc_id":7,"user_id":8,"nickname":9,"user_avatar":10,"doc_module":4,"category_id":11,"category_name":12,"doc_title":13,"doc_description":14,"doc_content":15,"file_id":16,"file_url":17,"file_type":18,"file_size":19,"view_count":20,"is_deleted":4,"is_public":21,"is_downloadable":21,"audit_status":21,"page_count":22,"language":23,"language_code":24,"site_id":25,"html_lang":24,"table_of_contents":26,"faqs":27,"seo_title":13,"seo_description":14,"update_tm":28,"read_time":29},83252,962075114101,"Seraphina","https://ap-avatar.wpscdn.com/avatar/e000253a75eb197efd?x-image-process=image/resize,m_fixed,w_180,h_180&k=1780044092746381165",8,"Research & Report","A Black-Box, Multilevel Algebraic Preconditioning Framework for Conforming Finite Elements","A black-box multilevel algebraic preconditioning framework is presented for conforming finite-element discretizations based on the least-squares algebraic-multigrid domain-decomposition (LS-AMG-DD) method. The approach clarifies when a Gram structure A = G⊤G arises, proving that on a prescribed degree-of-freedom cover C, a C-local Gram representation exists exactly when A admits a C-local symmetric positive semidefinite (SPSD) splitting. Elementwise SPSD energies then yield element-local Gram representations by appropriate choice of local factors, enabling LS-AMG-DD without detailed geometric hierarchy. Numerical tests show robustness and faster convergence than classical AMG, with errors reduced by 2–5 orders of magnitude in comparisons.","arXiv :2607 .07485v1 [math .NA] 8 Jul 2026  \nA BLACK-BOX, MULTILEVEL ALGEBRAIC PRECONDITIONING FRAMEWORK FOR CONFORMING FINITE ELEMENTS ∗  \nOLIVER A. KRZYSIK†, BEN S. SOUTHWORTH‡, AND GOLO A. WIMMER§  \nAbstract. Recently we introduced the least-squares algebraic-multigrid domain-decomposition (LS-AMG-DD) method as a multilevel, algebraic preconditioner for sparse symmetric positive definite (SPD) matrices that admit a Gram representation A = G⊤ G [37] . The factor G induces a local symmetric positive semidefinite (SPSD) splitting of A used to define local spectral problems from which an interpolation P is built, and a coarse-level Gram operator induced under Galerkin coarsening, Ac = G⊤cGc, for Gc := GP. This paper clarifies when this Gram structure arises, showing that, on a prescribed degree-of-freedom cover C, a C-local Gram representation of A exists if and only if A admits a C-local SPSD splitting. We then connect this viewpoint to conforming finite-element discretizations, where bilinear forms are naturally assembled from elementwise SPSD energies and therefore admit element-local Gram representations after choosing local factors (e.g., via algebraic factorizations of element blocks) . Taken together, these observations provide an essentially black-box route for applying LS-AMG-DD to conforming finite-element problems. Numerical tests illustrate the robustness of the method on several problems for which classical AMG methods require more than 105 iterations to converge, including high-order discretizations of grad–div in H(div), anisotropic hyperdiffusion in H2 , and linear elasticity in vector H 1 . Moreover, in some comparisons with existing AMG methods, LS-AMG-DD produces errors that are 2–5 orders of magnitude smaller, even when all methods are stopped at the same relative residual tolerance.  \nKey words. least squares, Gram matrix, algebraic multigrid, domain decomposition, spectral coarse spaces, conforming finite elements  \nMSC codes. 65F08, 65N30, 65N55  \n1. Introduction. Robust algebraic solvers for large sparse symmetric positive definite (SPD) systems remain a central need in scientific computing. In real application workflows, often an algebraic solver is necessary, as a user or application code may be unable or unwilling to provide more detailed information, such as a full geometric mesh hierarchy or a specific physical problem formulation a priori. Moreover, robustness is critical to ensure that a full multiphysics simulation does not fail due to sensitivity or lack of robustness of a single inner linear solver. For scalar elliptic problems posed in H 1 , algebraic multigrid (AMG) is often the method of choice in high-performance codes, and classical AMG and smoothed aggregation can deliver near-optimal complexity on isotropic, mildly anisotropic, and heterogeneous problems [44, 41] . However, even restricted to SPD systems that arise in the context of numerical partial differential equations (PDEs), AMG can quickly fail when confronted with complexities such as coupled systems of PDEs, finite-element spaces aside from H 1 , strong anisotropy and directional behavior, or high-order discretizations.  \nThese challenges all relate to standard assumptions underlying many AMG methods around near-nullspace modes, i.e., modes e for which Ae ≈ 0. Under typical pointwise smoothing, these are the slowest modes to converge, and also the modes  \n∗  \nFunding: This work was supported by the Laboratory Directed Research and Development program of Los Alamos National Laboratory, under project number 20240261ER. Los Alamos National Laboratory report number LA-UR-26-25292 .  \n†Theoretical Division, Los Alamos National Laboratory ([okrzysik@lanl.gov](okrzysik@lanl.gov), [https://orcid.org/](https://orcid.org/)[ ](https://orcid.org/)0000-0001-7880-6512.  \n‡Theoretical Division, Los Alamos National Laboratory ([southworth@lanl.gov](southworth@lanl.gov), [https://orcid.org/](https://orcid.org/)[ ](https://orcid.org/)0000-0002-02","cbCaibIoYcZbvz6G","https://ap.wps.com/l/cbCaibIoYcZbvz6G","pdf",2979966,3,1,25,"English","en",105,"# Introduction\n## Robust algebraic solvers for large SPD systems\n## Limitations of classical AMG assumptions\n## Near-nullspace modes and failure cases\n# Gram structure and SPSD splitting\n## When C-local Gram representations exist\n## Constructing interpolation from local spectral problems\n# Connection to conforming finite elements\n## Elementwise SPSD energies and local Gram factors\n# Numerical results","[{\"question\":\"What does the LS-AMG-DD framework aim to accomplish for conforming finite elements?\",\"answer\":\"It provides a multilevel algebraic preconditioner that can be applied in a largely black-box manner to conforming finite-element problems, improving robustness for challenging SPD systems.\"},{\"question\":\"Under what condition does a C-local Gram representation of A exist?\",\"answer\":\"A C-local Gram representation exists on a prescribed degree-of-freedom cover C if and only if A admits a C-local SPSD splitting.\"},{\"question\":\"Why can classical AMG fail in certain PDE-discretization settings?\",\"answer\":\"Failures arise when key assumptions about near-nullspace modes—such as locally smooth behavior across strong connections or representability by a direct sum of local coarse basis functions—break down under issues like strong anisotropy or higher-order and non-H1 finite-element spaces.\"}]",1784186269,63,{"code":4,"msg":31,"data":32},"ok",{"site_id":25,"language":24,"slug":33,"title":13,"keywords":34,"description":14,"schema_data":35,"social_meta":86,"head_meta":88,"extra_data":90,"updated_unix":28},"a-black-box-multilevel-algebraic-preconditioning-framework-for-conforming-finite-elements","",{"@graph":36,"@context":85},[37,53,68],{"@type":38,"itemListElement":39},"BreadcrumbList",[40,44,48,50],{"item":41,"name":42,"@type":43,"position":21},"https://docshare.wps.com","Home","ListItem",{"item":45,"name":46,"@type":43,"position":47},"https://docshare.wps.com/document/","Document",2,{"item":49,"name":12,"@type":43,"position":20},"https://docshare.wps.com/document/research-report/",{"item":51,"name":13,"@type":43,"position":52},"https://docshare.wps.com/document/a-black-box-multilevel-algebraic-preconditioning-framework-for-conforming-finite-elements/83252/",4,{"url":51,"name":13,"@type":54,"author":55,"headline":13,"publisher":57,"fileFormat":60,"inLanguage":24,"description":14,"dateModified":61,"datePublished":62,"encodingFormat":60,"isAccessibleForFree":63,"interactionStatistic":64},"DigitalDocument",{"name":9,"@type":56},"Person",{"url":41,"name":58,"@type":59},"DocShare","Organization","application/pdf","2026-07-25","2026-07-16",true,{"@type":65,"interactionType":66,"userInteractionCount":20},"InteractionCounter",{"@type":67},"ViewAction",{"@type":69,"mainEntity":70},"FAQPage",[71,77,81],{"name":72,"@type":73,"acceptedAnswer":74},"What does the LS-AMG-DD framework aim to accomplish for conforming finite elements?","Question",{"text":75,"@type":76},"It provides a multilevel algebraic preconditioner that can be applied in a largely black-box manner to conforming finite-element problems, improving robustness for challenging SPD systems.","Answer",{"name":78,"@type":73,"acceptedAnswer":79},"Under what condition does a C-local Gram representation of A exist?",{"text":80,"@type":76},"A C-local Gram representation exists on a prescribed degree-of-freedom cover C if and only if A admits a C-local SPSD splitting.",{"name":82,"@type":73,"acceptedAnswer":83},"Why can classical AMG fail in certain PDE-discretization settings?",{"text":84,"@type":76},"Failures arise when key assumptions about near-nullspace modes—such as locally smooth behavior across strong connections or representability by a direct sum of local coarse basis functions—break down under issues like strong anisotropy or higher-order and non-H1 finite-element spaces.","https://schema.org",{"og:url":51,"og:type":87,"og:title":13,"og:site_name":58,"og:description":14},"article",{"robots":89,"canonical":51},"index,follow",{"doc_id":7,"site_id":25},{"code":4,"msg":5,"data":92},[93,97,101,105,110,115,120,123,128,131,135],{"id":21,"doc_module":4,"doc_module_name":46,"category_name":94,"show_sort_weight":95,"slug":96},"Story & Novel",90,"story-novel",{"id":47,"doc_module":4,"doc_module_name":46,"category_name":98,"show_sort_weight":99,"slug":100},"Literature",80,"literature",{"id":52,"doc_module":4,"doc_module_name":46,"category_name":102,"show_sort_weight":103,"slug":104},"Exam",70,"exam",{"id":106,"doc_module":4,"doc_module_name":46,"category_name":107,"show_sort_weight":108,"slug":109},5,"Comic",60,"comic",{"id":111,"doc_module":4,"doc_module_name":46,"category_name":112,"show_sort_weight":113,"slug":114},6,"Technology",50,"technology",{"id":116,"doc_module":4,"doc_module_name":46,"category_name":117,"show_sort_weight":118,"slug":119},7,"Healthcare",40,"healthcare",{"id":11,"doc_module":4,"doc_module_name":46,"category_name":12,"show_sort_weight":121,"slug":122},30,"research-report",{"id":124,"doc_module":4,"doc_module_name":46,"category_name":125,"show_sort_weight":126,"slug":127},9,"Religion & Spirituality",20,"religion-spirituality",{"id":126,"doc_module":4,"doc_module_name":46,"category_name":129,"show_sort_weight":126,"slug":130},"World Cup","world-cup",{"id":132,"doc_module":4,"doc_module_name":46,"category_name":133,"show_sort_weight":132,"slug":134},10,"Lifestyle","lifestyle",{"id":136,"doc_module":4,"doc_module_name":46,"category_name":137,"show_sort_weight":106,"slug":138},19,"General","general"]