Chucui, with a recent developer version your code fails for me in the 9th time step when the solver fails to converge within 1000 iterations. Hence, I can't reproduce the problem with a newer deal.II version.
Please provide a minimal example by removing all functionality that is not necessary to show the issue. It is simply not feasible to step through a code with almost 4000 lines of code. If I had to guess, I would say that you are calling clear() on one of your BlockSparseMatrix objects before the object is destroyed (without reinitializing the object in between) and this case is not handled correctly in version 8.5. Best, Daniel -- 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.
