Hi Daniel, I did indeed have this in.
I seem to have fixed it when I rewrote everything from scratch... seems sometimes it is quicker to start from the beginning. Thanks to all that helped - I have learned a great deal from this post alone and I hope I can contribute as much. On Saturday, September 8, 2018 at 5:36:28 PM UTC+1, Daniel Arndt wrote: > > Jane, > > >> For more info, the VectorTools flux conditions are in my setup_dofs code, >> within costraints.clear()... constraints.close(). After assembling the >> system, I don't use the distribute local to global function but I instead >> distribute the constraints onto the solution after I have solved the >> system. Is there something wrong with the order?? I am getting a similar >> issue with equations like step-22 but I will post separately as it is >> different. >> > > That sounds a bit strange to me. If you are using a ConstraintMatrix > object (or an AffineConstraints object), you should normally take the > constraints into account while assembling, i.e. use > > constraints.distrbiute_local_to_glocal(cell_matrix, cell_rhs, dof_indices, > system_matrix, system_rhs); > > and again after solving using > > constraints.distribute(solution); > > 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 dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
