Daniel, I've run the problem with a simple GridGenerator::hyper_rectangle and the results differ. I'm filing a bug report for this.
The norm of the matrix difference is 10^11, but at least the sparsity pattern remains the same, if this could be of any consolation! Should I look at something else? Thanks! Franco On Thursday, March 9, 2017 at 7:33:16 PM UTC+1, Daniel Arndt wrote: > > Franco, > > on the first glance your code looks good. Can you reduce the size of your > problem to a single cell? > Is there still a difference then? > > Best, > Daniel > > Am Mittwoch, 8. März 2017 13:04:01 UTC+1 schrieb Franco Milicchio: >> >> I have made available the snippet for this example here >> <https://bitbucket.org/snippets/fmilicchio/aebEd> with the possibility >> of using either the MeshWorker or not (just change the #define). >> >> As far as I see, the code is copied and pasted correctly. >> >> The right hand sides *are the same*, but *matrices differ*. >> >> Thanks for any pointers! >> Franco >> > -- 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. For more options, visit https://groups.google.com/d/optout.
