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
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