On 11/19/21 6:26 PM, Marco Feder wrote:
indeed, as we have incompressibility, the last component of the rhs will
always be 0. I wrote a possible patch, I'll open the pr tomorrow so that I can
take another look at it!
Awesome, thanks!
W.
--
Thanks Wolfgang,
indeed, as we have incompressibility, the last component of the rhs will
always be 0. I wrote a possible patch, I'll open the pr tomorrow so that I
can take another look at it!
Best,
Marco
Il giorno venerdì 19 novembre 2021 alle 23:15:50 UTC+1 Wolfgang Bangerth ha
scritto:
On 11/19/21 12:57 AM, Ray Mclaren wrote:
I am solving a vector valued equation (mechanical equilibrium) using dealii's
matrix free infrastructure.
Although the library has several examples on how to apply neumann BC using
matrix free most are either scalar equations and using DG methods.
Marco,
I was going through step-22 and I saw that since you're using primitive
elements, then you compute the local_rhs contribution using
/fe.system_to_component_index(i).first/
However, as written in the program, we could as well multiply the dim+1 tensor
having the values of the i-th
Hello everyone!
This is deal.II newsletter #191.
It automatically reports recently merged features and discussions about the
deal.II finite element library.
## Below you find a list of recently proposed or merged features:
#12970: Replace cancel-builds-on-update by concurrency feature in the
