W dniu niedziela, 1 października 2017 21:20:06 UTC+2 użytkownik Daniel Arndt napisał: > > Michael, > > I implemented something similar to what you proposed. You might want to > have a look at https://github.com/dealii/dealii/pull/5175 > > Best, > Daniel > > Thanks, for now I'm using mine version (that seems to be working), I'll switch later to this one.
I've found another problem with off-diagonal blocks. The multiplication by B and BT is done by cell_loop and it triggers assertion in when used inside multigrid. The problem is that my divergence operator did had adjust_ghost_range_if_necessary method. I resolved the problem by copying from Base operator, but I think that it may be added to MatrixFree class so that it had not have to be copied every time. -- 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.
