pith. sign in

arxiv: 2510.25252 · v2 · submitted 2025-10-29 · 🧮 math.NA · cs.NA

Spectral analysis of the stiffness matrix sequence in the approximated Stokes equation

Samuele Ferri , Chiara Giraudo , Valerio Loi , Miroslav Kuchta , Stefano Serra-Capizzano This is my paper

Pith reviewed 2026-05-18 03:35 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords spectral analysisStokes problemTaylor-Hood elementsvariable viscositymatrix sequencesfinite element discretizationlocalization
0
0 comments X

The pith

The stiffness matrices from Taylor-Hood approximations of the 2D Stokes problem with variable viscosity exhibit localized eigenvalues and a limiting spectral distribution under weak regularity assumptions on the viscosity.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper examines the eigenvalues of the matrix sequences that arise when the two-dimensional Stokes equations with variable viscosity are discretized by Taylor-Hood P2-P1 finite elements. It establishes that the spectra localize in intervals that do not depend on the mesh size and that the matrices possess a well-defined limiting eigenvalue distribution as the discretization is refined. These properties matter because they directly affect the conditioning of the resulting linear systems and therefore the design of efficient iterative solvers for incompressible flow problems. The analysis is carried out under only weak regularity requirements on the viscosity coefficient.

Core claim

Localization and distributional spectral results are provided for the matrix sequences arising from the Taylor-Hood P2-P1 approximation of variable viscosity for the 2d Stokes problem under weak assumptions on the regularity of the diffusion. These results are accompanied by numerical tests and visualizations. A preliminary study of the impact of the findings on the preconditioning problem is also presented.

What carries the argument

The sequence of stiffness matrices generated by the Taylor-Hood P2-P1 finite-element discretization of the variable-viscosity two-dimensional Stokes equations, whose eigenvalues are shown to localize and to possess a limiting distribution.

If this is right

  • Eigenvalues of the discrete Stokes matrices remain bounded independently of the mesh size under the stated weak regularity assumptions.
  • The limiting spectral distribution permits asymptotic prediction of the convergence behavior of Krylov solvers as the mesh is refined.
  • Matrix-free or symbol-based preconditioners can be constructed that exploit the known spectral information.
  • The same localization framework applies to other low-regularity variable-coefficient Stokes-type problems.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Analogous localization may hold for three-dimensional or higher-order mixed finite-element Stokes discretizations.
  • The distributional result could be combined with existing symbol analysis to produce explicit bounds on the number of iterations required by common iterative solvers.
  • The approach might extend to the linearized Navier-Stokes equations when the convective term is treated as a lower-order perturbation.

Load-bearing premise

The localization and distributional results depend on the diffusion coefficient satisfying only weak regularity assumptions.

What would settle it

If eigenvalues computed on successively refined meshes for a viscosity with a jump discontinuity fail to remain inside the predicted localization interval, the claimed spectral localization would be disproved.

Figures

Figures reproduced from arXiv: 2510.25252 by Chiara Giraudo, Miroslav Kuchta, Samuele Ferri, Stefano Serra-Capizzano, Valerio Loi.

Figure 1
Figure 1. Figure 1: Effect of variable viscosity (4) on conditioning of the Stokes model (1) using the preconditioner (3). Lagrangian (AL) preconditioners. The Schur complement-based precondition￾ers have been developed e.g. in [40, 41, 17]. In these works the inverse of the Schur complement is approximated using BFBT approach (derived from a Least-Squares Commutator [20]) as S −1 h ≈ (BhD −1 h B T h ) −1 (BhD −1 h AhD −1 h B… view at source ↗
Figure 2
Figure 2. Figure 2: Ordering for the cells/viscosity (left), [PITH_FULL_IMAGE:figures/full_fig_p017_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Adherence of the eigenvalues of Ax,n to its symbol for n = 16. The first line displays the plots corresponding to Group 1 (on the left) and Group 2 (on the right), while the second line shows Group 3; from left γ = 1, γ = 10, γ = 100. Absence of Outliers Throughout all tests, we do not observe eigenvalues of An outside the essential range of the symbol. The reason is both numerical and theoretical. Indeed,… view at source ↗
Figure 4
Figure 4. Figure 4: Adherence of the singular values to its symbol for [PITH_FULL_IMAGE:figures/full_fig_p032_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Adherence of the eigenvalues of Mn to its symbol for n = 16. The first line displays the plots corresponding to Group 1 (on the left) and Group 2 (on the right), while the second line shows Group 3; from left γ = 1, γ = 10, γ = 100. converge within 1000 iterations. Tables 1 and 2 show the number of iterations used by PGMRES to achieve the tolerance of 10−5 , while figure 6 displays the successful clusterin… view at source ↗
Figure 6
Figure 6. Figure 6: Cluster around 1 of the singular values of the preconditioned system. [PITH_FULL_IMAGE:figures/full_fig_p034_6.png] view at source ↗
read the original abstract

In the present paper, we analyze in detail the spectral features of the matrix sequences arising from the Taylor-Hood $\mathbb{P}_2$-$\mathbb{P}_1$ approximation of variable viscosity for $2d$ Stokes problem under weak assumptions on the regularity of the diffusion. Localization and distributional spectral results are provided, accompanied by numerical tests and visualizations. A preliminary study of the impact of our findings on the preconditioning problem is also presented. A final section with concluding remarks and open problems ends the current work.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript analyzes the spectral features of matrix sequences from the Taylor-Hood P2-P1 finite-element discretization of the 2D Stokes problem with variable viscosity. It establishes localization and distributional spectral results under weak regularity assumptions on the diffusion coefficient, supported by numerical tests, visualizations, and a preliminary discussion of implications for preconditioning.

Significance. If the central claims hold, the work would usefully extend spectral localization techniques to variable-coefficient indefinite problems under minimal regularity, aiding analysis of iterative solvers. The combination of analytical arguments with numerical validation and the preconditioning study is a positive feature.

major comments (1)
  1. The distributional spectral result under merely measurable (L^∞) diffusion coefficients is load-bearing for the main claim, yet the argument that the symbol remains well-defined and that the essential spectrum localization carries over to the indefinite Stokes block structure is not fully detailed; standard GLT multiplication by a discontinuous coefficient requires additional approximation or density arguments that are not automatic here.
minor comments (2)
  1. Notation for the block structure of the discrete Stokes operator should be introduced earlier and used consistently when stating the symbol.
  2. The numerical section would benefit from explicit statements of the viscosity functions tested (e.g., discontinuous or highly oscillatory) and the precise mesh sizes employed.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the constructive comment and positive overall assessment. We agree that the distributional result for L^∞ coefficients is central and that the GLT arguments for the indefinite block structure merit additional explicit justification. The revised manuscript expands the relevant section with the requested details while preserving the original weak-regularity setting.

read point-by-point responses

  1. Referee: The distributional spectral result under merely measurable (L^∞) diffusion coefficients is load-bearing for the main claim, yet the argument that the symbol remains well-defined and that the essential spectrum localization carries over to the indefinite Stokes block structure is not fully detailed; standard GLT multiplication by a discontinuous coefficient requires additional approximation or density arguments that are not automatic here.

    Authors: We appreciate this observation. In the revised version we have added a new paragraph in Section 3.2 that first recalls the GLT symbol for the constant-viscosity Taylor-Hood pair and then justifies its extension to L^∞ viscosity via a density argument: continuous functions are dense in L^∞, the symbol map is continuous with respect to the L^∞ norm in the appropriate matrix-valued topology, and the essential spectrum of the approximating sequences converges to that of the limit operator. Because the saddle-point structure is preserved under this approximation (the pressure block remains unchanged and the velocity block remains symmetric positive definite), the localization of the essential spectrum carries over directly to the indefinite case. The same density argument handles the multiplication by a discontinuous coefficient without requiring extra regularity. We believe these additions fully address the concern while staying within the paper’s weak-assumption framework. revision: yes

Circularity Check

0 steps flagged

No circularity: spectral localization and distributional results derived via analytical GLT-style arguments under stated weak regularity assumptions.

full rationale

The paper claims localization and distributional spectral results for the Taylor-Hood P2-P1 stiffness matrix sequence of the variable-viscosity 2D Stokes problem. These are presented as consequences of the matrix sequence symbol under weak (L^infty or measurable) assumptions on the diffusion coefficient, accompanied by numerical tests. No step reduces a claimed prediction or uniqueness result to a fitted parameter, self-defined quantity, or load-bearing self-citation chain; the derivation remains self-contained against external matrix-sequence theory and does not rename known empirical patterns or smuggle ansatzes via prior work by the same authors. The central claims therefore do not collapse by construction to the inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review provides no explicit free parameters, axioms, or invented entities. The analysis rests on standard finite-element assumptions for Stokes problems and the stated weak regularity condition on the diffusion.

pith-pipeline@v0.9.0 · 5618 in / 978 out tokens · 32970 ms · 2026-05-18T03:35:52.231623+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

56 extracted references · 56 canonical work pages

  1. [1]

    Adriani, I

    A. Adriani, I. Furci, C. Garoni, S. Serra-Capizzano. Spectral and sin- gular value distribution of sequences of block matrices with rectangular Toeplitz blocks. Part I: Asymptotically rational block size ratios.Journal of Numerical Mathematics, 2026, in press

    work page 2026

  2. [2]

    Adriani, C

    A. Adriani, C. Garoni, S. Serra-Capizzano. Spectral and singular value dis- tribution of sequences of block matrices with rectangular Toeplitz blocks. Part I: Asymptotically irrational block size ratios.Journal of Numerical Mathematics, 2026, in press

    work page 2026

  3. [3]

    Adriani, A

    A. Adriani, A. J. A. Schiavoni-Piazza, S. Serra-Capizzano, C. Tablino- Possio. Revisiting the notion of approximating class of sequences for han- dling approximated PDEs on moving or unbounded domains.Electronic Transactions on Numerical Analysis, 63:424–451, 2025

    work page 2025

  4. [4]

    S. W. Armfield. Finite difference solutions of the Navier-Stokes equations on staggered and non-staggered grids.Computers & fluids, 20(1):1–17, 1991. 35

    work page 1991

  5. [5]

    Barakitis, M

    N. Barakitis, M. Donatelli, S. Ferri, V. Loi, S. Serra-Capizzano, R. Sor- mani. Blocking structures, approximation, and preconditioning.Numeri- cal Algorithms. 10.1007/s11075-025-02157-y. Online 8-7-2025

    work page doi:10.1007/s11075-025-02157-y 2025

  6. [6]

    Barbarino

    G. Barbarino. A systematic approach to reduced GLT.BIT Numerical Mathematics, 62(3):681–743, 2022

    work page 2022

  7. [7]

    Barbarino, D

    G. Barbarino, D. Bianchi, Carlo Garoni. Constructive approach to the monotone rearrangement of functions.Expositiones Mathematicae, 40(1):155–175, 2022

    work page 2022

  8. [8]

    Barbarino, C

    G. Barbarino, C. Garoni, M. Mazza, S. Serra-Capizzano. Rectangular GLT sequences.Electronic Transactions on Numerical Analysis.55: 585– 617, 2022

    work page 2022

  9. [9]

    Barbarino, C

    G. Barbarino, C. Garoni, S. Serra-Capizzano. Block generalized locally Toeplitz sequences: theory and applications in the unidimensional case. Electronic Transactions on Numerical Analysis, 53:28–112, 2020

    work page 2020

  10. [10]

    Barbarino, C

    G. Barbarino, C. Garoni, S. Serra-Capizzano. Block generalized locally Toeplitz sequences: theory and applications in the multidimensional case. Electronic Transactions on Numerical Analysis, 53:113–216, 2020

    work page 2020

  11. [11]

    Barbarino, S

    G. Barbarino, S. Serra-Capizzano. Non-Hermitian perturbations of Her- mitian matrix-sequences and applications to the spectral analysis of the numerical approximation of partial differential equations.Numerical Lin- ear Algebra with Applications.27(3), e2286, 31 pp., 2020

    work page 2020

  12. [12]

    Beckermann, A

    B. Beckermann, A. B. J. Kuijlaars. Superlinear convergence of conjugate gradients.SIAM Journal on Numerical Analysis, 39(1): 300—329, 2001

    work page 2001

  13. [13]

    Beckermann, S

    B. Beckermann, S. Serra-Capizzano. On the asymptotic spectrum of finite element matrix sequences.SIAM Journal on Numerical Analysis, 45(2): 746—769, 2007

    work page 2007

  14. [14]

    Benzi, G

    M. Benzi, G. H. Golub, J. Liesen. Numerical solution of saddle point problems.Acta numerica, 14:1–137, 2005

    work page 2005

  15. [15]

    Borzacchiello, E

    D. Borzacchiello, E. Leriche, B. Blotti` ere, J. Guillet. Box-relaxation based multigrid solvers for the variable viscosity Stokes problem.Computers & Fluids, 156:515–525, 2017

    work page 2017

  16. [16]

    Brezzi, R

    F. Brezzi, R. S. Falk. Stability of higher-order Hood-Taylor methods. SIAM Journal on Numerical Analysis, 28(3):581–590, 1991

    work page 1991

  17. [17]

    Burkhart, N

    A. Burkhart, N. Kohl, B. Wohlmuth, J. Zawallich. A robust matrix- free approach for large-scale non-isothermal high-contrast viscosity Stokes flow on blended domains with applications to geophysics.arXiv preprint arXiv:2506.04157, 2025. 36

    work page arXiv 2025

  18. [18]

    Dorostkar, M

    A. Dorostkar, M. Neytcheva, S. Serra-Capizzano. Spectral analysis of cou- pled PDEs and of their Schur complements via generalized locally Toeplitz sequences in 2D.Computer Methods in Applied Mechanics and Engineer- ing, 309:74–105, 2016

    work page 2016

  19. [19]

    Ekstr¨ om, S

    S.-E. Ekstr¨ om, S. Serra-Capizzano. Eigenvalues and eigenvectors of banded Toeplitz matrices and the related symbols.Numerical Linear Al- gebra with Applications, 25(5):e2137, 2018

    work page 2018

  20. [20]

    Elman, V

    H. Elman, V. E. Howle, J. Shadid, R. Shuttleworth, R. Tuminaro. Block preconditioners based on approximate commutators.SIAM Journal on Scientific Computing, 27(5):1651–1668, 2006

    work page 2006

  21. [21]

    Farrell, P

    P. Farrell, P. A. G. Orozco, E. S¨ uli. Finite element approximation and pre- conditioning for anisothermal flow of implicitly-constituted non-newtonian fluids.Mathematics of Computation, 91(334):659–697, 2022

    work page 2022

  22. [22]

    Furci, A

    I. Furci, A. Adriani, Stefano Serra-Capizzano. Block structured matrix- sequences and their spectral and singular value canonical distributions: a general theory.arXiv preprint arXiv:2409.06465, 2024

    work page arXiv 2024

  23. [23]

    Garoni, S

    C. Garoni, S. Serra-Capizzano.Generalized locally Toeplitz sequences: theory and applications, volume 1. Springer, 2017

    work page 2017

  24. [24]

    Garoni, S

    C. Garoni, S. Serra-Capizzano.Generalized locally Toeplitz sequences: theory and applications, volume 2. Springer, 2018

    work page 2018

  25. [25]

    Garoni, S

    C. Garoni, S. Serra-Capizzano. Generalized locally Toeplitz sequences: a spectral analysis tool for discretized differential equations.Splines and PDEs: from approximation theory to numerical linear algebra, Lecture Notes in Math., Fond. CIME/CIME Found. Subser., Springer, Cham, volume 2219:161–236, 2018

    work page 2018

  26. [26]

    Garoni, H

    C. Garoni, H. Speleers, S.-E. Ekstr¨ om, A. Reali, S. Serra-Capizzano, T. J. R. Hughes. Symbol-based analysis of finite element and isogeomet- ric B-spline discretizations of eigenvalue problems: exposition and re- view.Archives of Computational Methods in Engineering, 26(5):1639– 1690, 2019

    work page 2019

  27. [27]

    Gergelits, K

    T. Gergelits, K. Mardal, B. F. Nielsen, Z. Strakos. Laplacian precondi- tioning of elliptic PDEs: Localization of the eigenvalues of the discretized operator.SIAM Journal on Numerical Analysis, 57(3):1369–1394, 2019

    work page 2019

  28. [28]

    P. P. Grinevich, M. A. Olshanskii. An iterative method for the Stokes-type problem with variable viscosity.SIAM Journal on Scientific Computing, 31(5):3959–3978, 2009

    work page 2009

  29. [29]

    X. He, M. Neytcheva. Preconditioning the incompressible Navier-Stokes equations with variable viscosity.Journal of Computational Mathematics, 30(5): 461–482, 2012. 37

    work page 2012

  30. [30]

    X. He, M. Neytcheva, S. Serra-Capizzano. On an augmented Lagrangian- based preconditioning of Oseen type problems.BIT, 51(4): 865—888, 2011

    work page 2011

  31. [31]

    A. Hirn, M. Lanzend¨ orfer, J. Stebel. Finite element approximation of flow of fluids with shear-rate-and pressure-dependent viscosity.IMA Journal of Numerical Analysis, 32(4):1604–1634, 2012

    work page 2012

  32. [32]

    Kuchta, G

    M. Kuchta, G. Tobie, K. Miljkovi´ c, M. Bˇ ehounkov´ a, O. Souˇ cek, G. Choblet, O. ˇCadek. Despinning and shape evolution of Saturn’s moon Iapetus triggered by a giant impact.Icarus, 252:454–465, 2015

    work page 2015

  33. [33]

    A. B. J. Kuijlaars. Convergence analysis of Krylov subspace iterations with methods from potential theory.SIAM Review, 48(1): 3—40, 2006

    work page 2006

  34. [34]

    Mazza, A

    M. Mazza, A. Ratnani, S. Serra-Capizzano. Spectral analysis and spectral symbol for the 2D curl-curl stabilized) operator with applications to the related iterative solutions.Mathematics of Computation, 88(317):1155– 1188, 2019

    work page 2019

  35. [35]

    Mazza, M

    M. Mazza, M. Semplice, S. Serra-Capizzano, E. Travaglia. A matrix- theoretic spectral analysis of incompressible Navier-Stokes staggered DG approximations and a related spectrally based preconditioning approach. Numerische Mathematik, 149(4):933–971, 2021

    work page 2021

  36. [36]

    M. F. Murphy, G. H. Golub, A. J. Wathen. A note on precondition- ing for indefinite linear systems.SIAM Journal on Scientific Computing, 21(6):1969–1972, 2000

    work page 1969

  37. [37]

    Noutsos, S

    D. Noutsos, S. Serra-Capizzano, P. Vassalos. Matrix algebra precondi- tioners for multilevel Toeplitz systems do not insure optimal convergence rate.Theoretical Computer Science, 315(2-3):557-–579, 2004

    work page 2004

  38. [38]

    Pattyn, L

    F. Pattyn, L. Perichon, A. Aschwanden, B. Breuer, B. de Smedt, O. Gagliardini, G. H. Gudmundsson, R. C. A. Hindmarsh, A. Hubbard, J. V. Johnson, T. Kleiner, Y. Konovalov, C. Martin, A. J. Payne, D. Pollard, S. Price, M. R¨ uckamp, F. Saito, O. Souˇ cek, S. Sugiyama,. Benchmark experiments for higher-order and full Stokes ice sheet models. InAGU Fall Meeti...

    work page 2007

  39. [39]

    K. R. Rajagopal. On implicit constitutive theories for fluids.Journal of Fluid Mechanics, 550:243–249, 2006

    work page 2006

  40. [40]

    J. Rudi, A. Cristiano, I. Malossi, T. Isaac, G. Stadler, M. Gurnis, P. W. J. Staar, Y. Ineichen, C. Bekas, A. Curioni, O. Ghattas. An extreme-scale implicit solver for complex PDEs: Highly heterogeneous flow in earth’s mantle. InProceedings of the international conference for high perfor- mance computing, networking, storage and analysis, pages 1–12, 2015. 38

    work page 2015

  41. [41]

    J. Rudi, G. Stadler, O. Ghattas. Weighted BFBT preconditioner for Stokes flow problems with highly heterogeneous viscosity.SIAM Journal on Sci- entific Computing, 39(5):S272–S297, 2017

    work page 2017

  42. [42]

    N. S. Sarathkumar, S. Serra-Capizzano. GLT sequences and automatic computation of the symbol.Linear algebra and its applications, 697:468– 486, 2024

    work page 2024

  43. [43]

    Serra-Capizzano

    S. Serra-Capizzano. The rate of convergence of Toeplitz based PCG meth- ods for second order non-linear boundary value problems.Numerische Mathematik, 81(3):461-–495, 1999

    work page 1999

  44. [44]

    Serra-Capizzano

    S. Serra-Capizzano. Some theorems on linear positive operators and func- tionals and their applications.Computers and Mathematics with Applica- tion, 88(317):1155–1188, 2000. 39(7/8): 139-–167, (2000)

    work page 2000

  45. [45]

    Serra-Capizzano

    S. Serra-Capizzano. Distribution results on the algebra generated by Toeplitz sequences: a finite-dimensional approach.Linear algebra and its applications, 328:121–130, 2001

    work page 2001

  46. [46]

    Serra Capizzano

    S. Serra Capizzano. Matrix algebra preconditioners for multilevel Toeplitz matrices are not superlinear.Linear Algebra and its Applications, 343–344:303–319, 2002

    work page 2002

  47. [47]

    Serra Capizzano

    S. Serra Capizzano. Generalized locally Toeplitz sequences: spectral anal- ysis and applications to discretized partial differential equations.Linear Algebra and its Applications, 366:371–402, 2003

    work page 2003

  48. [48]

    Serra Capizzano

    S. Serra Capizzano. The GLT class as a generalized Fourier analysis and applications.Linear Algebra and its Applications, 419(1):180–233, 2006

    work page 2006

  49. [49]

    Y. Shih, G. Stadler, F. Wechsung. Robust multigrid techniques for aug- mented Lagrangian preconditioning of incompressible Stokes equations with extreme viscosity variations.SIAM Journal on Scientific Comput- ing, 45(3):S27–S53, 2022

    work page 2022

  50. [50]

    P. Tilli. A note on the spectral distribution of Toeplitz matrices.Linear and Multilinear Algebra, 45(2-3):147–159, 1998

    work page 1998

  51. [51]

    P. Tilli. Locally Toeplitz sequences: spectral properties and applications. Linear Algebra and its Applications, 278(1-3):91–120, 1998

    work page 1998

  52. [52]

    E. E. Tyrtyshnikov. A unifying approach to some old and new theorems on distribution and clustering.Linear algebra and its applications, 232:1–43, 1996

    work page 1996

  53. [53]

    E. E. Tyrtyshnikov, N. Zamarashkin. Spectra of multilevel Toeplitz ma- trices: advanced theory via simple matrix relationships.Linear algebra and its applications, 270:15–27, 1998. 39

    work page 1998

  54. [54]

    P. E. van Keken, C. Currie, S. D. King, M. D. Behn, A. Cagnioncle, J. He, R. F. Katz, S. Lin, E. M. Parmentier, M. Spiegelman, K. Wang. A community benchmark for subduction zone modeling.Physics of the Earth and Planetary Interiors, 171(1-4):187–197, 2008

    work page 2008

  55. [55]

    Wichrowski, P

    M. Wichrowski, P. Krzy˙ zanowski. A matrix-free multilevel preconditioner for the generalized Stokes problem with discontinuous viscosity.Journal of Computational Science, 63:101804, 2022

    work page 2022

  56. [56]

    Wichrowski, P

    M. Wichrowski, P. Krzy˙ zanowski, L. Heltai, S. Stupkiewicz. Exploit- ing high-contrast Stokes preconditioners to efficiently solve incompressible fluid-structure interaction problems.International Journal for Numerical Methods in Engineering, 124:5446–5470, 2023. 40

    work page 2023