Mauro, apologies for not getting back to this earlier. Have you found a solution?
I think what you find is something that we had to learn over time and have learned since step-35 was written: That it's typically easiest to just put all solution variables into one large finite element system, and then build a square matrix that, if decomposed, may have rectangular matrices as building blocks. If we were to re-write step-35 today, that's the design we'd choose -- in fact, if anyone wanted to re-write the program in this style, I think that would make for an excellent project we would very gladly take! So that would be what I suggest you consider: Rewriting things in such a way that the various matrices become blocks in a large square matrix. If you want to go the cheap way, you can always just build a large square matrix for which you really only care about the one block in question and the rest is only there to make sure that you have a square matrix, rather than containing anything useful. Best W. On 10/21/19 9:22 AM, Mauro Murer wrote: > It seems it doesn't work for rectangular objects!!! > > An error occurred in line <1191> of file > </home/murer/Desktop/dealii/dealii-9.1.1/include/deal.II/lac/affine_constraints.templates.h> > > in function > void > dealii::AffineConstraints<number>::condense(dealii::DynamicSparsityPattern&) > const [with number = double] > The violated condition was: > sparsity.n_rows() == sparsity.n_cols() > Additional information: > This function only works for quadratic objects! > > Il giorno lun 21 ott 2019 alle ore 16:35 Wolfgang Bangerth > <bange...@colostate.edu <mailto:bange...@colostate.edu>> ha scritto: > > On 10/21/19 8:32 AM, Mauro Murer wrote: > > Actually I can't find the function AffineConstraints::compress(). > > Sorry, I meant ::condense(). > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: bange...@colostate.edu > <mailto: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 > <mailto:dealii%2bunsubscr...@googlegroups.com>. > To view this discussion on the web visit > > https://groups.google.com/d/msgid/dealii/4d17534b-c63e-5d49-d86e-0151b3ee0dfe%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 > <mailto:dealii+unsubscr...@googlegroups.com>. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/CAK3LiWecHAVGMq%2Bj3GXwwJ%3DJqat-1OY4PTFWZgH_X0FFA5EbbA%40mail.gmail.com > > <https://groups.google.com/d/msgid/dealii/CAK3LiWecHAVGMq%2Bj3GXwwJ%3DJqat-1OY4PTFWZgH_X0FFA5EbbA%40mail.gmail.com?utm_medium=email&utm_source=footer>. -- ------------------------------------------------------------------------ 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/16d80277-c38b-d33b-b07b-dabc895bb2a5%40colostate.edu.
