Hello, I have a composite system with several overlapping FE spaces which I construct in blocks using independent variational forms for each space and then assembling them all into a large block matrix. I would like to enforce exact Dirichlet conditions without penalization (e.g. by using tgv = -1)…

What bookkeeping are you talking about? We could add a new tgv = -3 to handle this (I wanted to do this for a long time), but in the meantime, you can just impose on your off-diagonal blocks something like: varf vGamma(u, v) = on(gamma, u = 1); real[int] marker = vGamma(0, Vh, tgv = -1); for[i,j,v …

Preferred way to construct block matrices with exact Dirichlet BC (tgv = -1)?

General Discussion

arousta (Ali Roustaei) January 15, 2020, 11:14am 3

Also relevant, look at this post:
https://community.freefem.org/t/implementation-of-dirichlet-boundary-condition-when-tgv-1/113

1 Like

Topic		Replies	Views
Options to impose Dirichlet BC with "ON" General Discussion	4	480	December 26, 2021
Application of Dirichlet BCs General Discussion	6	1344	June 17, 2020
Implementation of Dirichlet boundary condition when tgv=-1 Feature Request	4	3297	August 19, 2021
Non-homogeneous Dirichlet BC at a point General Discussion	2	209	April 25, 2023
Constraining a specific element(s) in P0 space General Discussion	5	227	March 30, 2022

Preferred way to construct block matrices with exact Dirichlet BC (tgv = -1)?

Related Topics