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	565	December 26, 2021
Composite spaces, on() and tgv=-1 General Discussion	3	51	December 18, 2024
Application of Dirichlet BCs General Discussion	6	1718	June 17, 2020
Implementation of Dirichlet boundary condition when tgv=-1 Feature Request	4	3548	August 19, 2021
Exact Dirichlet BC General Discussion	3	1413	September 18, 2020

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

Related topics