On 10/16/2017 01:53 PM, Mark Ma wrote:
I think this problem lies in the time updating of solution using
old_solution, since the mass_matrix and laplace_matrix have already
eliminated the constraint node, /*mass_matrix_T.vmult
(system_rhs, old_solution_T_cal*//*);*/ is no longer valid for this
case.
Yes, this sounds reasonable from your description. Have you looked at
steps 24, 25, 26 to see how it is done there? I could imagine that you
need to think about which degrees of freedom you want eliminated/not
eliminated in the matrices you multiply with. To this end, write out the
weak form of the problem you need to solve in each time step, and think
specifically about the boundary values of the trial and test functions
and whether they should be part of the matrix you use on the right hand
side or not. It is possible that you need to use a different matrix for
the lhs and rhs operations.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
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 dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.