I created a ConstraintMatrix with both periodic and hanging node constraints,
and called close().
Then I pickeda constrained degree of freedom, lets say with global dof id =
“i” and printed the constraint equation pairs (j,a_ij) corresponding to “i” on
the processor for which “i” is
Dear dealii community,
I have a fluid-structure interaction type of problem to solve, where on one
domain I have the linearized compressible Navier-Stokes equations, and on
the other - equation of motion (elastic deformation). The PDEs that I have
are time-harmonic, so there is no time
Chucui,
Maybe a question more : how to make sure your period boundary coditions are
> added truly?
>
How would you do it for Dirichlet boundary conditions?
There are mutiple possibilities.
- You can just print the constraints (ConstraintMatrx.print(std::cout)),
print the locations of degrees
Dear Daniel,
Thank you very much! I use code like this:
constraints.distribute_local_to_global(local_matrix, local_dof_indices,
system_matrix) as you say,
and use codes like this
constraints.distribute(solution);
after every "solve()" step.
and now, my LittleCase0-3.cc can run.
Maybe a