I just use tecplot directly visualize the results. The vorticity contour from 
my simple code is continuous, and the results from deal.II is discontinuous 
(without L2 projection). Is it possible that the direct solver in Intel mkl did 
a similar projection step internally?

________________________________
From: Wolfgang Bangerth <bange...@colostate.edu>
Sent: Sunday, January 19, 2020, 12:49 AM
To: David Eaton; dealii@googlegroups.com
Subject: Re: [deal.II] discontinous contour over elements

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/SL2PR01MB2553A18C05BF7C941F6AFBC5C3330%40SL2PR01MB2553.apcprd01.prod.exchangelabs.com.

Reply via email to