Actually I can't find the function AffineConstraints::compress(). I tried to use SparseMatrix::compress(::VectorOperation::values), on the other hand, but the error still remain. Maybe I don't use the right function??
Thanks Il giorno lun 21 ott 2019 alle ore 15:52 Wolfgang Bangerth < bange...@colostate.edu> ha scritto: > On 10/21/19 2:55 AM, Mauro Murer wrote: > > > > I tried to take constraints into account in the sparsity pattern, but I > don't > > know how to do it for a rectangular matrix such in my case. > > I have rows equal to the number of velocity dofs and columns equal to > pressure > > dofs. > > I didn't found any function for sparsity pattern in the documentation > that > > take in account different constraints (in my case velocity and pressure > > constraints). > > Then just build the sparsity pattern as you do and use > AffineConstraints::compress() after the fact -- I believe that there is > such a > function that can work on rectangular matrices. > > In general, you are discussing why it is simpler to just keep all > components > of the solution in the same DoFHandler object. In that case, the > rectangular > matrix is just one block in a square matrix, and all of the constraint > handling things work as expected. > > 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. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/146c2446-0d06-1cec-a908-0b474e4747dd%40colostate.edu > . > -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/CAK3LiWcu6Lya8rKQhfVv3o%3DMOitFH2hpD%2BRL_swYP7mr6rebCg%40mail.gmail.com.