Re: [deal.II] Reuse of the system_matrix in the Picard iteration for Inhomogenous Dirichlet boundary conditions

2019-04-27 Thread hanks0227
Oh, thank you very much!!! I was so stupid!! I thought that "system_matrix_tmp = system_matrix" is the same as "system_matrix_tmp.copy_from(system_matrix)" Thank you. Kyusik. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see

Re: [deal.II] Reuse of the system_matrix in the Picard iteration for Inhomogenous Dirichlet boundary conditions

2019-04-26 Thread Wolfgang Bangerth
On 4/26/19 8:04 PM, hanks0...@gmail.com wrote: > > > The several components in the matrix A_ij change to 0 or 1, and these changed > values are used in assembling the system_rhs. > (if it is wrong please let me know what is the reason) > > So, my question is reassembling the system_matrix