Dear Daniel,
Happy new year. I wish a wonderful new year for you !
This is only (directly) possible if you are not using inhomogeneous
> constraints. In that case, you can call
> ConstraintMatrix::distribute_local_to_global separately for matrix and
> right-hand side.
>
I am struggling to apply homogeneous periodic BC constraints in assembly
since ConstraintMatrix::condense is not usable in parallel and I can not
use ConstraintMatrix::distribute_local_to_global since I construct
system_matrix and system_rhs in global level as a linear combination of
pre-costructed mass_matrix(M), laplace_matrix(A), nonlinear
term(nl_term),nonlinear matrix(nl_matrix), solution (U) and
old_solution(U_old) vectors as follows:
system_rhs= [M+k*theta*A]*U-[M-k*(1-theta)*A]*U_old -k*nl_term(U,U_old)
in which theta is Crank_Nicolson parameter
system_matrix= M+k*theta*A-k*theta*nl_matrix(U, U_old)
However, to be able to use ConstraintMatrix::distribute_local_to_global, I
can assemble system_matrix in local level instead of global level. Having
that said, I can not do so for the system_rhs because solution and
old_solution vectors are global vectors.
Any hint regarding how to apply constraint on my system_rhs would be
appreciated.
Thanks,
Hamed
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
