Thank you, I'll take a look in detail at those examples. However, those examples don't use PETSc and the problem is that I can't find any method in PETScWrappers::SparseMatrix to copy a matrix. Which way would you suggest?
David. On Saturday, 20 August 2016 00:41:05 UTC+2, Wolfgang Bangerth wrote: > > On 08/19/2016 04:23 PM, David F wrote: > > > > P.S.: I tried to simply repeat the process to re-apply the BCs with a > new > > value using the same matrix, but I get always the first BCs I apply (I > think > > because once I apply it, the matrix knows it has been already condensed > and > > ignores the next calls). > > Yes. What you need to do is in every time step copy the matrix you have > previously assembled into a new matrix, call apply_boundary_values() to > this > matrix and your current right hand side, and then solve the linear system. > This way, you always start from the same matrix. > > I think one of steps 23, 24, 25, 26 does this. > > Best > Wolfgang > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: [email protected] > <javascript:> > 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 [email protected]. For more options, visit https://groups.google.com/d/optout.
