Skip to content

Mixed Poisson with efficient pre-conditioning#264

Merged
jorgensd merged 4 commits intomainfrom
dokken/mixed-poisson
Aug 12, 2025
Merged

Mixed Poisson with efficient pre-conditioning#264
jorgensd merged 4 commits intomainfrom
dokken/mixed-poisson

Conversation

@jorgensd
Copy link
Copy Markdown
Owner

@jorgensd jorgensd commented Aug 9, 2025

Resolve #145

@jorgensd
Copy link
Copy Markdown
Owner Author

Use the slightly altered form:

s = -(
    ufl.inner(ufl.grad(u), ufl.grad(v)) * ufl.dx
    - ufl.inner(ufl.avg(ufl.grad(v)), ufl.jump(u, n)) * ufl.dS
    - ufl.inner(ufl.jump(u, n), ufl.avg(ufl.grad(v))) * ufl.dS
    + (alpha * (degree + 1) ** 2 / ufl.avg(h))
    * ufl.inner(ufl.jump(u, n), ufl.jump(v, n))
    * ufl.dS
    - ufl.inner(ufl.grad(u), v * n) * dGammaD
    - ufl.inner(u * n, ufl.grad(v)) * dGammaD
    + (gamma * degree**2 / h) * u * v * dGammaD
)
@jorgensd
Copy link
Copy Markdown
Owner Author

@jorgensd jorgensd merged commit 7ba2fae into main Aug 12, 2025
1 check passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

1 participant