[deal.II] neumann boundary condition for a vector valued continuous FE equation using matrix free

2021-11-18 Thread Ray Mclaren
Hi, 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. Is there an example of continuous vect

[deal.II] Extractors for local_rhs in step-22

2021-11-18 Thread Marco Feder
Dear all, 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 sha