Skip to content

GitLab

  • Menu
Projects Groups Snippets
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Sign in / Register
  • dumux dumux
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 91
    • Issues 91
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 82
    • Merge requests 82
  • Deployments
    • Deployments
    • Releases
  • Packages & Registries
    • Packages & Registries
    • Container Registry
  • Monitor
    • Monitor
    • Incidents
  • Analytics
    • Analytics
    • Value stream
    • Repository
  • Wiki
    • Wiki
  • Activity
  • Graph
  • Create a new issue
  • Commits
  • Issue Boards
Collapse sidebar
  • dumux-repositories
  • dumuxdumux
  • Issues
  • #804
Closed
Open
Created Jan 27, 2020 by Kilian Weishaupt@kweisOwner

[staggered] Matrix block arrangement does not comply with literature standard

Papers dealing with solvers for saddlepoint problems (like the incompressible Navier-Stokes equations, example) usually define the block matrix as

M = \begin{pmatrix}
  A & B^T\\
  B & C
\end{pmatrix}

where A is the block for the derivatives of the momentum balance equation residuals w.r.t velocity ("velocity block") and C is is the block for the derivatives of the mass balance equation residuals w.r.t pressure ("pressure block").

For incompressible fluids, C is zero.

However, the multitype blockmatrix in Dumux for staggered problems looks like this:

M = \begin{pmatrix}
  C & B^T\\
  B & A
\end{pmatrix}

which causes confusion and requires special care, e.g, for the implementation of the Uzawa solver !1827 (closed)

I suggest we adapt the matrix structure such that the velocity block is on the upper left. Are there any concerns / objections?

Edited Jan 27, 2020 by Kilian Weishaupt
Assignee
Assign to
Time tracking