Skip to content
GitLab
Projects Groups Topics Snippets
  • /
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Register
  • Sign in
  • dumux dumux
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributor statistics
    • Graph
    • Compare revisions
  • Issues 100
    • Issues 100
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 92
    • Merge requests 92
  • CI/CD
    • CI/CD
    • Pipelines
    • Jobs
    • Schedules
  • Deployments
    • Deployments
    • Environments
    • Releases
  • Packages and registries
    • Packages and registries
    • Container Registry
  • Monitor
    • Monitor
    • Incidents
  • Analytics
    • Analytics
    • Value stream
    • CI/CD
    • Repository
  • Wiki
    • Wiki
  • Activity
  • Graph
  • Create a new issue
  • Jobs
  • Commits
  • Issue Boards
Collapse sidebar
  • dumux-repositoriesdumux-repositories
  • dumuxdumux
  • Issues
  • #804
Closed
Open
Issue created Jan 27, 2020 by Kilian Weishaupt@kweisMaintainer

[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