On 1/18/20 9:25 AM, David Eaton wrote: > Thank you for your explanations. Basically I formed a weak form of the PDE > for > one element and numerically integrate it at the Gaussian points based on the > interpolation from the local nodes. Subsequently, I assemble the weak forms > from all elements into a global system matrix based on a local-to-global > mapping of the nodes. After applying the boundary conditions, I solve this > linear system using a linear solver in Intel mkl.
Right -- that gives you the solution (u,p) of the problem. But then what do you do to visualize whatever it is that you find is/isn't discontinuous? Best W. -- ------------------------------------------------------------------------ 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/5c1ef574-52b6-8449-743c-aae9c355460d%40colostate.edu.